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.

SourceResearch connectionFrom yearThrough year

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

At least 1,405 records · Page 78Linked to original sources

USGS2008

Geochemical controls of elevated arsenic concentrations in groundwater, Ester Dome, Fairbanks district, Alaska

Ester Dome, an upland area near Fairbanks, Alaska, was chosen for a detailed hydrogeochemical study because of the previously reported elevated arsenic in groundwater, and the presence of a large set of wells amenable to detailed sampling. Ester Dome lies within the Fairbanks mining district, where gold-bearing quartz veins, typically containing 2–3 vol.% sulfide minerals (arsenopyrite, stibnite, and pyrite), have been mined both underground and in open cuts. Gold-bearing veins on Ester Dome occur in shear zones and the sulfide minerals in these veins have been crushed to fine-grained material by syn- or post-mineralization movement. Groundwater at Ester Dome is circumneutral, Ca–HCO 3 to Ca–SO 4 type, and ranges from dilute (specific conductance of 48 µS/cm) to more concentrated (specific conductance as high as 2070 µS/cm). In general, solute concentrations increase down hydrologic gradient. Redox species indicate that the groundwaters range from oxic to sub-oxic (low dissolved oxygen, Fe(III) reduction, no SO 4 reduction). Waters with the highest Fe concentrations, as high as 10.7 mg/L, are the most anoxic. Dissolved As concentrations range from < 1 to 1160 µg/L, with a median value of 146 µg/L. Arsenic concentrations are not correlated with specific conductance or Fe concentrations, suggesting that neither groundwater residence time, nor reductive dissolution of iron oxyhydroxides, control the arsenic chemistry. Furthermore, As concentrations do not covary with other constituents that form anions and oxyanions in solution (e.g., HCO 3 , Mo, F, or U) such that desorption of arsenic from clays or oxides also does not control arsenic mobility. Oxidation of arsenopyrite and dissolution of scorodite, in the near-surface environment appears to be the primary control of dissolved As in this upland area. More specifically, the elevated As concentrations are spatially associated with sulfidized shear zones and localities of gold-bearing quartz veins. Consistent with this interpretation, elevated dissolved Sb concentrations (as high as 59 µg/L), also correlated with occurrences of hypogene sulfide minerals, were measured in samples with high dissolved As concentrations.

Chemical Geology↗

USGS2013

Long-term groundwater contamination after source removal—The role of sorbed carbon and nitrogen on the rate of reoxygenation of a treated-wastewater plume on Cape Cod, MA, USA

The consequences of groundwater contamination can remain long after a contaminant source has been removed. Documentation of natural aquifer recoveries and empirical tools to predict recovery time frames and associated geochemical changes are generally lacking. This study characterized the long-term natural attenuation of a groundwater contaminant plume in a sand and gravel aquifer on Cape Cod, Massachusetts, after the removal of the treated-wastewater source. Although concentrations of dissolved organic carbon (DOC) and other soluble constituents have decreased substantially in the 15 years since the source was removed, the core of the plume remains anoxic and has sharp redox gradients and elevated concentrations of nitrate and ammonium. Aquifer sediment was collected from near the former disposal site at several points in time and space along a 0.5-km-long transect extending downgradient from the disposal site and analyses of the sediment was correlated with changes in plume composition. Total sediment carbon content was generally low (< 8 to 55.8 μmol (g dry wt)− 1) but was positively correlated with oxygen consumption rates in laboratory incubations, which ranged from 11.6 to 44.7 nmol (g dry wt)− 1 day− 1. Total water extractable organic carbon was < 10–50% of the total carbon content but was the most biodegradable portion of the carbon pool. Carbon/nitrogen (C/N) ratios in the extracts increased more than 10-fold with time, suggesting that organic carbon degradation and oxygen consumption could become N-limited as the sorbed C and dissolved inorganic nitrogen (DIN) pools produced by the degradation separate with time by differential transport. A 1-D model using total degradable organic carbon values was constructed to simulate oxygen consumption and transport and calibrated by using observed temporal changes in oxygen concentrations at selected wells. The simulated travel velocity of the oxygen gradient was 5–13% of the groundwater velocity. This suggests that the total sorbed carbon pool is large relative to the rate of oxygen entrainment and will be impacting groundwater geochemistry for many decades. This has implications for long-term oxidation of reduced constituents, such as ammonium, that are being transported downgradient away from the infiltration beds toward surface and coastal discharge zones.

Massachusetts↗

USGS2017

Interannual to multidecadal climate forcings on groundwater resources of the U.S. West Coast

Study region The U.S. West Coast, including the Pacific Northwest and California Coastal Basins aquifer systems. Study focus Groundwater response to interannual to multidecadal climate variability has important implications for security within the water–energy–food nexus. Here we use Singular Spectrum Analysis to quantify the teleconnections between AMO, PDO, ENSO, and PNA and precipitation and groundwater level fluctuations. The computer program DAMP was used to provide insight on the influence of soil texture, depth to water, and mean and period of a surface infiltration flux on the damping of climate signals in the vadose zone. New hydrological insights for the region We find that PDO, ENSO, and PNA have significant influence on precipitation and groundwater fluctuations across a north-south gradient of the West Coast, but the lower frequency climate modes (PDO) have a greater influence on hydrologic patterns than higher frequency climate modes (ENSO and PNA). Low frequency signals tend to be preserved better in groundwater fluctuations than high frequency signals, which is a function of the degree of damping of surface variable fluxes related to soil texture, depth to water, mean and period of the infiltration flux. The teleconnection patterns that exist in surface hydrologic processes are not necessarily the same as those preserved in subsurface processes, which are affected by damping of some climate variability signals within infiltrating water.

California, Oregon, Washington↗

USGS2018

Long-term changes in pond permanence, size, and salinity in Prairie Pothole Region wetlands: The role of groundwater-pond interaction

Study Region Cottonwood Lake area wetlands, North Dakota, U.S.A. Study Focus Fluctuations in pond permanence, size, and salinity are key features of prairie-pothole wetlands that provide a variety of wetland habitats for waterfowl in the northern prairie of North America. Observation of water-level and salinity fluctuations in a semi-permanent wetland pond over a 20-year period, included periods when the wetland occasionally was dry, as well as wetter years when the pond depth and surface extent doubled while volume increased 10 times. New hydrological insights for the study region Compared to all other measured budget components, groundwater flow into the pond often contributed the least water (8–28 percent) but the largest amount (>90 percent) of specific solutes to the water and solute budgets of the pond. In drier years flow from the pond into groundwater represented > 10 percent of water loss, and in 1992 was approximately equal to evapotranspiration loss. Also during the drier years, export of calcium, magnesium, sodium, potassium, chloride, and sulfate by flow from the pond to groundwater was substantial compared with previous or subsequent years, a process that would have been undetected if groundwater flux had been calculated as a net value. Independent quantification of water and solute gains and losses were essential to understand controls on water-level and salinity fluctuations in the pond in response to variable climate conditions.

Journal of Hydrology: Regional Studies↗

USGS2021

Distributed memory parallel groundwater modeling for the Netherlands Hydrological Instrument

Worldwide, billions of people rely on fresh groundwater reserves for their domestic, agricultural and industrial water use. Extreme droughts and excessive groundwater pumping put pressure on water authorities in maintaining sustainable water usage. High-resolution integrated models are valuable assets in supporting them. The Netherlands Hydrological Instrument (NHI) provides the Dutch water authorities with open source modeling software and data. However, NHI integrated groundwater models often require long run times and large memory usage, therefore strongly limiting their application. As a solution, we present a distributed memory parallelization, focusing on the National Hydrological Model. Depending on the level of integration, we show that significant speedups can be obtained up to two orders of magnitude. As far as we know, this is the first reported integrated groundwater parallelization of an operational hydrological model used for national-scale integrated water management and policy making. The parallel model code and data are freely available.

Environmental Software & Modelling↗

USGS2005

Source and movement of helium in the eastern Morongo groundwater Basin: The influence of regional tectonics on crustal and mantle helium fluxes

We assess the role of fracturing and seismicity on fluid-driven mass transport of helium using groundwaters from the eastern Morongo Basin (EMB), California, USA. The EMB, located ???200 km east of Los Angeles, lies within a tectonically active region known as the Eastern California Shear Zone that exhibits both strike-slip and extensional deformation. Helium concentrations from 27 groundwaters range from 0.97 to 253.7 ?? 10-7 cm3 STP g-1 H2O, with corresponding 3He/4He ratios falling between 1.0 and 0.26 RA (where RA is the 3He/4He ratio of air). All groundwaters had helium isotope ratios significantly higher than the crustal production value of ???0.02 RA. Dissolved helium concentrations were resolved into components associated with solubility equilibration, air entrainment, in situ production within the aquifer, and extraneous fluxes (both crustal and mantle derived). All samples contained a mantle helium-3 (3Hem) flux in the range of 4.5 to 1351 ?? 10-14 cm3 STP 3He cm-2 yr-1 and a crustal flux (J0) between 0.03 and 300 ?? 10-7 cm3 STP 4He cm-2 yr-1. Groundwaters from the eastern part of the basin contained significantly higher 3Hem and deep crustal helium-4 (4Hedc) concentrations than other areas, suggesting a localized source for these components. 4Hedc and 3Hem are strongly correlated, and are associated with faults in the basin. A shallow thermal anomaly in a >3,000 m deep graben in the eastern basin suggests upflow of fluids through active faults associated with extensional tectonics. Regional tectonics appears to drive large scale crustal fluid transport, whereas episodic hydrofracturing provides an effective mechanism for mantle-crust volatile transport identified by variability in the magnitude of degassing fluxes (3Hem and J0) across the basin. Copyright ?? 2005 Elsevier Ltd.

Geochimica et Cosmochimica Acta↗

USGS2008

Evidence for crustal degassing of CF4 and SF6 in Mojave Desert groundwaters

Dissolved tetrafluoromethane (CF4) and sulfur hexafluoride (SF6) concentrations were measured in groundwater samples from the Eastern Morongo Basin (EMB) and Mojave River Basin (MRB) located in the southern Mojave Desert, California. Both CF4 and SF6 are supersaturated with respect to equilibrium with the preindustrial atmosphere at the recharge temperatures and elevations of the Mojave Desert. These observations provide the first in situ evidence for a flux of CF4 from the lithosphere. A gradual basin-wide enhancement in dissolved CF4 and SF6 concentrations with groundwater age is consistent with release of these gases during weathering of the surrounding granitic alluvium. Dissolved CF4 and SF6 concentrations in these groundwaters also contain a deeper crustal component associated with a lithospheric flux entering the EMB and MRB through the underlying basement. The crustal flux of CF4, but not of SF6, is enhanced in the vicinity of local active fault systems due to release of crustal fluids during episodic fracture events driven by local tectonic activity. When fluxes of CF4 and SF6 into Mojave Desert groundwaters are extrapolated to the global scale they are consistent, within large uncertainties, with the fluxes required to sustain the preindustrial atmospheric abundances of CF4 and SF6. ?? 2007 Elsevier Ltd. All rights reserved.

Geochimica et Cosmochimica Acta↗

USGS2012

Geochemical modeling of changes in shallow groundwater chemistry observed during the MSU-ZERT CO2 injection experiment

A field experiment involving the release of carbon dioxide (CO 2 ) into a shallow aquifer was conducted near Bozeman, Montana, during the summer of 2008, to investigate the potential groundwater quality impacts in the case of leakage of CO 2 from deep geological storage. As an essential part of the Montana State University Zero Emission Research and Technology (MSU-ZERT) field program, food-grade CO 2 was injected over a 30 day period into a horizontal perforated pipe a few feet below the water table of a shallow aquifer. The impact of elevated CO 2 concentrations on groundwater quality was investigated by analyzing water samples taken before, during, and following CO 2 injection, from observation wells located in the vicinity of the injection pipe, and from two distant monitoring wells. Field measurements and laboratory analyses showed rapid and systematic changes in pH, alkalinity, and conductance, as well as increases in the aqueous concentrations of naturally occurring major and trace element species. The geochemical data were evaluated using principal component analysis (PCA) to (1) understand potential correlations between aqueous species, and (2) to identify minerals controlling the chemical composition of the groundwater prior to CO 2 injection. These evaluations were used to assess possible geochemical processes responsible for the observed increases in the concentrations of dissolved constituents, and to simulate these processes using a multicomponent reaction path model. Reasonable agreement between observed and modeled data suggests that (1) calcite dissolution was the primary pH buffer, yielding increased Ca +2 concentrations in the groundwater, (2) increases in the concentrations of most major and trace metal cations except Fe could be a result of Ca +2 -driven exchange reactions, (3) the release of anions from adsorption sites due to competitive adsorption of carbonate could explain the observed trends of most anions, and (4) the dissolution of reactive Fe minerals (presumed ferrihydrite and fougerite, from thermodynamic analyses) could explain increases in total Fe concentration. Highlights ► Because the possibility of CO 2 leakage cannot be completely ruled out, the potential impact of CO 2 intrusion on the quality of fresh water aquifers overlying CO 2 storage sites needs to be investigated. ► Geochemical data from a field experiment involving the release of carbon dioxide (CO 2 ) into a shallow aquifer were evaluated. ► Geochemical model used to assess possible geochemical processes responsible for the observed increases in the concentrations of dissolved constituents. ► Reasonable agreement between observed and modeled data suggests that increases in the concentrations of most major and trace metal cations except Fe could be a result of Ca +2 -driven exchange reactions and the release of anions from adsorption sites due to competitive adsorption of carbonate could explain the observed trends of most anions.

Montana↗

USGS2015

A comparison of estimates of basin-scale soil-moisture evapotranspiration and estimates of riparian groundwater evapotranspiration with implications for water budgets in the Verde Valley, Central Arizona, USA

Population growth in the Verde Valley in Arizona has led to efforts to better understand water availability in the watershed. Evapotranspiration (ET) is a substantial component of the water budget and a critical factor in estimating groundwater recharge in the area. In this study, four estimates of ET are compared and discussed with applications to the Verde Valley. Higher potential ET (PET) rates from the soil-water balance (SWB) recharge model resulted in an average annual ET volume about 17% greater than for ET from the basin characteristics (BCM) recharge model. Annual BCM PET volume, however, was greater by about a factor of 2 or more than SWB actual ET (AET) estimates, which are used in the SWB model to estimate groundwater recharge. ET also was estimated using a method that combines MODIS-EVI remote sensing data and geospatial information and by the MODFLOW-EVT ET package as part of a regional groundwater-flow model that includes the study area. Annual ET volumes were about same for upper-bound MODIS-EVI ET for perennial streams as for the MODFLOW ET estimates, with the small differences between the two methods having minimal impact on annual or longer groundwater budgets for the study area.

Arizona↗

USGS2006

Groundwater-surface water interaction in the riparian zone of an incised channel, Walnut Creek, Iowa

Riparian zones of many incised channels in agricultural regions are cropped to the channel edge leaving them unvegetated for large portions of the year. In this study we evaluated surface and groundwater interaction in the riparian zone of an incised stream during a spring high flow period using detailed stream stage and hydraulic head data from six wells, and water quality sampling to determine whether the riparian zone can be a source of nitrate pollution to streams. Study results indicated that bank storage of stream water from Walnut Creek during a large storm water runoff event was limited to a narrow 1.6 m zone immediately adjacent to the channel. Nitrate concentrations in riparian groundwater were highest near the incised stream where the unsaturated zone was thickest. Nitrate and dissolved oxygen concentrations and nitrate-chloride ratios increased during a spring recharge period then decreased in the latter portion of the study. We used MODFLOW and MT3DMS to evaluate dilution and denitrification processes that would contribute to decreasing nitrate concentrations in riparian groundwater over time. MT3DMS model simulations were improved with a denitrification rate of 0.02 1/d assigned to the floodplain sediments implying that denitrification plays an important role in reducing nitrate concentrations in groundwater. We conclude that riparian zones of incised channels can potentially be a source of nitrate to streams during spring recharge periods when the near-stream riparian zone is largely unvegetated. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

USGS2017

Prediction and visualization of redox conditions in the groundwater of Central Valley, California

Regional-scale, three-dimensional continuous probability models, were constructed for aspects of redox conditions in the groundwater system of the Central Valley, California. These models yield grids depicting the probability that groundwater in a particular location will have dissolved oxygen (DO) concentrations less than selected threshold values representing anoxic groundwater conditions, or will have dissolved manganese (Mn) concentrations greater than selected threshold values representing secondary drinking water-quality contaminant levels (SMCL) and health-based screening levels (HBSL). The probability models were constrained by the alluvial boundary of the Central Valley to a depth of approximately 300 m. Probability distribution grids can be extracted from the 3-D models at any desired depth, and are of interest to water-resource managers, water-quality researchers, and groundwater modelers concerned with the occurrence of natural and anthropogenic contaminants related to anoxic conditions. Models were constructed using a Boosted Regression Trees (BRT) machine learning technique that produces many trees as part of an additive model and has the ability to handle many variables, automatically incorporate interactions, and is resistant to collinearity. Machine learning methods for statistical prediction are becoming increasing popular in that they do not require assumptions associated with traditional hypothesis testing. Models were constructed using measured dissolved oxygen and manganese concentrations sampled from 2767 wells within the alluvial boundary of the Central Valley, and over 60 explanatory variables representing regional-scale soil properties, soil chemistry, land use, aquifer textures, and aquifer hydrologic properties. Models were trained on a USGS dataset of 932 wells, and evaluated on an independent hold-out dataset of 1835 wells from the California Division of Drinking Water. We used cross-validation to assess the predictive performance of models of varying complexity, as a basis for selecting final models. Trained models were applied to cross-validation testing data and a separate hold-out dataset to evaluate model predictive performance by emphasizing three model metrics of fit: Kappa; accuracy; and the area under the receiver operator characteristic curve (ROC). The final trained models were used for mapping predictions at discrete depths to a depth of 304.8 m. Trained DO and Mn models had accuracies of 86–100%, Kappa values of 0.69–0.99, and ROC values of 0.92–1.0. Model accuracies for cross-validation testing datasets were 82–95% and ROC values were 0.87–0.91, indicating good predictive performance. Kappas for the cross-validation testing dataset were 0.30–0.69, indicating fair to substantial agreement between testing observations and model predictions. Hold-out data were available for the manganese model only and indicated accuracies of 89–97%, ROC values of 0.73–0.75, and Kappa values of 0.06–0.30. The predictive performance of both the DO and Mn models was reasonable, considering all three of these fit metrics and the low percentages of low-DO and high-Mn events in the data.

California↗

USGS2024

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine↗

USGS2007

Investigation of the groundwater system at Masaya Caldera, Nicaragua, using transient electromagnetics and numerical simulation

The distribution of groundwater beneath Masaya Volcano, in Nicaragua, and its surrounding caldera was characterized using the transient electromagnetic method (TEM). Multiple soundings were conducted at 30 sites. Models of the TEM data consistently indicate a resistive layer that is underlain by one or more conductive layers. These two layers represent the unsaturated and saturated zones, respectively, with the boundary between them indicating the water-table elevation. A map of the TEM data shows that the water table in the caldera is a subdued replica of the topography, with higher elevations beneath the edifice in the south-central caldera and lower elevations in the eastern caldera, coinciding with the elevation of Laguna de Masaya. These TEM data, combined with regional hydrologic data, indicate that the caldera in hydrologically isolated from the surrounding region, with as much as 60??m of difference in elevation of the groundwater table across caldera-bounding faults. The water-table information and estimates of fluxes of water through the system were used to constrain a numerical simulation of groundwater flow. The simulation results indicate that basalt flows in the outer parts of the caldera have a relatively high transmissivity, whereas the central edifice has a substantially lower transmissivity. A layer of relatively high transmissivity must be present at depth within the edifice in order to deliver the observed flux of water and steam to the active vent. This hydrologic information about the caldera provides a baseline for assessing the response of this isolated groundwater system to future changes in magmatic activity. ?? 2007.

Journal of Volcanology and Geothermal Research↗

USGS2008

Seasonal changes in submarine groundwater discharge to coastal salt ponds estimated using 226Ra and 228Ra as tracers

Submarine groundwater discharge (SGD) to coastal southern Rhode Island was estimated from measurements of the naturally-occurring radioisotopes 226Ra (t1/2 = 1600??y) and 228Ra (t1/2 = 5.75??y). Surface water and porewater samples were collected quarterly in Winnapaug, Quonochontaug, Ninigret, Green Hill, and Pt. Judith-Potter Ponds, as well as nearly monthly in the surface water of Rhode Island Sound, from January 2002 to August 2003; additional porewater samples were collected in August 2005. Surface water activities ranged from 12-83??dpm 100??L- 1 (60??dpm = 1??Bq) and 21-256??dpm 100??L- 1 for 226Ra and 228Ra, respectively. Porewater 226Ra activities ranged from 16-736??dpm 100??L- 1 (2002-2003) and 95-815??dpm 100??L- 1 (2005), while porewater 228Ra activities ranged from 23-1265??dpm 100??L- 1. Combining these data with a simple box model provided average 226Ra-based submarine groundwater fluxes ranging from 11-159??L m- 2 d- 1 and average 228Ra-derived fluxes of 15-259??L m- 2 d- 1. Seasonal changes in Ra-derived SGD were apparent in all ponds as well as between ponds, with SGD values of 30-472??L m- 2 d- 1 (Winnapaug Pond), 6-20??L m- 2 d- 1 (Quonochontaug Pond), 36-273??L m- 2 d- 1 (Ninigret Pond), 29-76??L m- 2 d- 1 (Green Hill Pond), and 19-83??L m- 2 d- 1 (Pt. Judith-Potter Pond). These Ra-derived fluxes are up to two orders of magnitude higher than results predicted by a numerical model of groundwater flow, estimates of aquifer recharge for the study period, and values published in previous Ra-based SGD studies in Rhode Island. This disparity may result from differences in the type of flow (recirculated seawater versus fresh groundwater) determined using each technique, as well as variability in porewater Ra activity. ?? 2007 Elsevier B.V. All rights reserved.

Marine Chemistry↗

USGS2008

Eddy correlation measurements of submarine groundwater discharge

This paper presents a new, non-invasive means of quantifying groundwater discharge into marine waters using an eddy correlation approach. The method takes advantage of the fact that, in virtually all aquatic environments, the dominant mode of vertical transport near the sediment–water interface is turbulent mixing. The technique thus relies on measuring simultaneously the fluctuating vertical velocity using an acoustic Doppler velocimeter and the fluctuating salinity and/or temperature using rapid-response conductivity and/or temperature sensors. The measurements are typically done at a height of 5–15 cm above the sediment surface, at a frequency of 16 to 64 Hz, and for a period of 15 to 60 min. If the groundwater salinity and/or temperature differ from that of the water column, the groundwater specific discharge (cm d− 1) can be quantified from either a heat or salt balance. Groundwater discharge was estimated with this new approach in Salt Pond, a small estuary on Cape Cod (MA, USA). Estimates agreed well with previous estimates of discharge measured using seepage meters and 222Rn as a tracer. The eddy correlation technique has several desirable characteristics: 1) discharge is quantified under in-situ hydrodynamic conditions; 2) salinity and temperature can serve as two semi-independent tracers of discharge; 3) discharge can be quantified at high temporal resolution, and 4) long-term records of discharge may be possible, due to the low power requirements of the instrumentation.

Massachussetts↗

USGS2023

Physicochemical coastal groundwater dynamics between Kauhakō Crater lake and Kalaupapa settlement, Moloka‘i, Hawai‘i

Land-based sources of groundwater pollution can be a critical threat to coral reefs, and a better understanding of “ridge-to-reef” water movement is required to advance management and coral survival in the Anthropocene. In this study a more complete understanding of the geological, atmospheric, and oceanic drivers behind coastal groundwater exchange on the Kalaupapa peninsula, on Moloka‘i, Hawai‘i, is obtained by analyzing high resolution geochemical and geophysical time-series data. In concert with multiyear water level analyses, a tidally and precipitation-driven groundwater connection between Kauhakō Crater lake and submarine groundwater discharge (SGD) fluxes are demonstrated. Results include an average discharge rate of 190 cm d −1 and the detection of water-flow pathways past cesspools that likely contribute to higher nutrient loading near the SGD sites. This underlines the importance of managing anthropogenic nutrients that enter the shallow freshwater lens such as through cesspools and are consequently discharged via SGD onto coral reef habitats.

Hawai'i↗

USGS2014

Factors affecting temporal variability of arsenic in groundwater used for drinking water supply in the United States

The occurrence of arsenic in groundwater is a recognized environmental hazard with worldwide importance and much effort has been focused on surveying and predicting where arsenic occurs. Temporal variability is one aspect of this environmental hazard that has until recently received less attention than other aspects. For this study, we analyzed 1245 wells with two samples per well. We suggest that temporal variability, often reported as affecting very few wells, is perhaps a larger issue than it appears and has been masked by datasets with large numbers of non-detect data. Although there was only a slight difference in arsenic concentration variability among samples from public and private wells (p = 0.0452), the range of variability was larger for public than for private wells. Further, we relate the variability we see to geochemical factors—primarily variability in redox—but also variability in pH and major-ion chemistry. We also show that in New England there is a weak but statistically significant indication that seasonality may have an effect on concentrations, whereby concentrations in the first two quarters of the year (January–June) are significantly lower than in the second two quarters (July–December) (p < 0.0001). In the Central Valley of California, though not statistically significant (p = 0.4169), arsenic concentration is lower in the first quarter of the year but increases in subsequent quarters. In both regions, these changes appear to follow groundwater levels. It is possible that this difference in arsenic concentrations is related to groundwater level changes, pumping stresses, evapotranspiration effects, or perhaps mixing of more oxidizing, lower pH recharge water in wetter months. Focusing on the understanding the geochemical conditions in aquifers where arsenic concentrations are concerns and causes of geochemical changes in the groundwater environment may lead to a better understanding of where and by how much arsenic will vary over time.

Science of the Total Environment↗

USGS2019

Predicting attenuation of salinized surface- and groundwater-resources from legacy energy development in the Prairie Pothole Region

Oil and gas (energy) development in the Williston Basin, which partly underlies the Prairie Pothole Region in central North America, has helped meet U.S. energy demand for decades. Historical handling and disposal practices of saline wastewater co-produced during energy development resulted in salinization of surface and groundwater at numerous legacy energy sites. Thirty years of monitoring (1988–2018) at Goose Lake, which has been producing since the 1960s, documents long-term spatial and temporal changes in water quality from legacy energy development. Surface water quality was highly variable and decoupled from changes in groundwater quality, likely due to annual and regional climatic fluctuations. Therefore, changes in surface water-quality were not considered a reliable indicator of subsurface chloride migration. However, chloride concentrations in monitoring wells near wastewater sources exhibited systematic temporal reductions allowing for estimates of the time required for natural attenuation of groundwater to U.S. Environmental Protection Agency acute and chronic chloride toxicity benchmarks and a local background level. Point attenuation rates differed based on sediment type (outwash vs till) and yielded a range of predicted years when water-quality targets will be reached: acute – 2045 to 2113; chronic – 2069 to 2160; background – 2126 to 2275. Bulk attenuation rates from four separate years of data were used to calculate the distances chloride could migrate downgradient from the largest wastewater source. Potential distances of downgradient migration before dilution to water-quality targets decreased from 1989 to 2018: acute – 949 to 673 m; chronic – 1220 to 922 m; background – 1878 to 1525 m. Several downgradient wetlands are within these distances and will continue to receive saline contaminated groundwater for years. While these results demonstrate chloride attenuation at a legacy energy site, they also highlight the persistence of saline wastewater contamination and the need to mitigate future spills to prevent long-term salinization from energy development.

Montana↗