USGS Science⌕ Search

SEARCH · USGS Science

Results for “groundwater aquifer”

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

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

At least 1,441 records · Page 80Linked to original sources

Analyzing multi-year nitrate concentration evolution in Alabama aquatic systems using a machine learning model

Rising nitrate contamination in water systems poses significant risks to public health and ecosystem stability, necessitating advanced modeling to understand nitrate dynamics more accurately. This study applies the long short-term memory (LSTM) modeling to investigate the hydrologic and environmental factors influencing nitrate concentration dynamics in rivers and aquifers across the state of Alabama in the southeast of the United States. By integrating dynamic data such as streamflow and groundwater levels with static catchment attributes, the machine learning model identifies primary drivers of nitrate fluctuations, offering detailed insights into the complex interactions affecting multi-year nitrate concentrations in natural aquatic systems. In addition, a novel LSTM-based approach utilizes synthetic surface water nitrate data to predict groundwater nitrate levels, helping to address monitoring gaps in aquifers connected to these rivers. This method reveals potential correlations between surface water and groundwater nitrate dynamics, which is particularly meaningful given the lack of water quality observations in many aquifers. Field applications further show that, while the LSTM model effectively captures seasonal trends, limitations in representing extreme nitrate events suggest areas for further refinement. These findings contribute to data-driven water quality management, enhancing understanding of nitrate behavior in interconnected water systems.

Alabama↗

Simulation of regional groundwater flow and advective transport of per- and polyfluoroalkyl substances, Joint Base McGuire-Dix-Lakehurst and vicinity, New Jersey, 2018

A three-dimensional numerical model of groundwater flow was developed and calibrated for the unconsolidated New Jersey Coastal Plain aquifers underlying Joint Base McGuire-Dix-Lakehurst (JBMDL) and vicinity, New Jersey, to evaluate groundwater flow pathways of per- and polyfluoroalkyl substances (PFAS) contamination associated with use of aqueous film forming foam (AFFF) at the base. The regional subsurface flow model spans an area of approximately 518 square miles around JBMDL and is based on a previously developed hydrogeologic framework of the area. Steady-state flow in the unconsolidated aquifers was simulated using the MODFLOW 6 groundwater flow model, which is able to account for hydrostratigraphic pinchouts and discontinuities in the Coastal Plain aquifers underlying JBMDL. To account for local patterns of fluid flow driving advective subsurface migration of PFAS, the grid was refined using quadtree meshes spanning 21 areas where historical AFFF use was identified, five off-site reconnaissance areas identified by AFCEC as areas in which the occurrence of PFAS is most likely to pose a potential danger to local drinking water supplies, and along streams that behave as drains in the base-flow-dominated Coastal Plain. Following grid refinement, four physical processes known to govern subsurface flow were introduced to the model. These included effective precipitation recharge, discharge to streams and stream-connected wetlands, regional inflows and outflows along the model bottom, and withdrawals from wells, each of which were incorporated into the model as either external or internal boundary conditions. To account for effective precipitation recharge, a specified-flow boundary was assigned along the top of the model. Similarly, regional flows predicted using the modified U.S Geological Survey’s New Jersey Coastal Plain Regional Aquifer System Analysis model were treated as specified-flow boundary conditions along the bottom of the model. Base-flow losses were treated as drains along streams delineated using a 10-foot LiDAR dataset. Drains were also assigned to cells falling within stream-connected National Hydrologic Database wetlands. Finally, well-pumpage data mined from the New Jersey Water Transfer database were added to the model to account for extraction of groundwater through pumping from industrial-supply and drinking-water-supply wells. Along model edges established at groundwater divides, where the net flux of water across the boundary is equal to zero, natural no-flow boundary conditions were imposed. The refined flow model was calibrated using the parameter-estimation (PEST) program, which adjusts model parameters by performing a gradient search over the sum-of-squared-error objective function until the parameter set that produces simulated water levels and base flows most closely matches 544 water levels and 20 estimated base flows and closely adheres to initial parameter estimates. Based on the analysis of calibration residuals, the model did not appear to be affected by significant model structural error. The MODPATH particle-tracking algorithm was used to estimate advective transport paths of PFAS in the vicinity of JBMDL. Forward tracking was used to determine paths of PFAS away from AFFF source areas to streams, wetlands, pumping wells, and geographic areas that PFAS may contaminate. Additionally, reverse tracking was used to determine particle pathlines away from off-site PFAS reconnaissance areas, or areas within which all sources of PFAS might be advectively transported into subsurface drinking-water supplies, to locations at land surface that may indicate a source of PFAS. The coupled and calibrated groundwater flow and particle-tracking transport model provide valuable tools for predicting the relative extent of PFAS contamination from onsite legacy source areas. The calibrated model also provides measures of water-level and base-flow observation influence that can help guide future data-collection efforts related to groundwater and surface water sampling for PFAS.

New Jersey↗

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes significantly to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides significant volumes of agricultural products not just for California but the entire United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface water and groundwater resources in the valley. This has resulted in water management challenges in the Salinas Valley that are predominantly related to distribution of water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. To deal with the distribution issue, historically water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface water supplies are insufficient. Availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand on groundwater and surface water resources throughout the Salinas Valley has resulted in undesirable effects of unsustainable water use, such as surface water depletion, groundwater level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface water and groundwater resources. To meet water demands and reduce undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and a variable climate. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resource Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley Hydrogeologic system called the Salinas Valley System Model. The Salinas Valley Geologic Framework was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The Salinas Valley Watershed Model simulates the entire Salinas River watershed. Monthly surface water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

EarthArXiv↗

Ground-water resources investigation in the Amran Valley, Yeman Arab Republic

A program of hydrologic studies and exploratory drilling was conducted intermittently between 1974 and 1978 to evaluate the water-bearing properties of the unconsolidated alluvial sediments and associated rocks in the semi-arid Amran Valley basin, an 800-square-kilometer area in north-central Yemen Arab Republic. Inventory data from 395 wells were compiled, observation well and rain-gage networks were established and 16 standard complete chemical analyses were made for samples from selected wells. The water resources of the area were overexploited. The chemical quality of the water is generally good. Four aquifer tests were run to determine transmissivity and storage characteristics. The pumping tests show that groundwater occurs under semi-confined leaky-aquifer conditions in the valley fill. Wells drilled in the alluvial fill of the south-central part of the valley have the highest yields. Wells penetrating the limestone and volcanic rocks generally have little or no yield except in fracture zones. Basalt flows occur interbedded with the wadi alluvium at several depths. Cropping out rocks in the Amran Valley range in age from late Jurassic to Holocene. (USGS)

Open-File Report↗

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes substantially to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides substantial volumes of agricultural products not just for California but for the United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface-water and groundwater resources in the valley. This situation has resulted in water management challenges in the Salinas Valley that generally relate to the distribution of the water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. Historically, to deal with the distribution issue, water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface-water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface-water supplies are insufficient. The availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand for groundwater and surface-water resources throughout the Salinas Valley has resulted in undesirable effects from unsustainable water use, such as surface-water depletion, groundwater-level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface-water and groundwater resources. To meet water demands and reduce the undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and climate variability. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resources Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley hydrogeologic system called the Salinas Valley System Model. The geologic framework is known as the Salinas Valley Geologic Framework and was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The surface-water model is called the Salinas Valley Watershed Model and simulates the Salinas River watershed. Monthly surface-water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

California↗

Adaptation of the Carter-Tracy water influx calculation to groundwater flow simulation

The Carter-Tracy calculation for water influx is adapted to groundwater flow simulation with additional clarifying explanation not present in the original papers. The Van Everdingen and Hurst aquifer-influence functions for radial flow from an outer aquifer region are employed. This technique, based on convolution of unit-step response functions, offers a simple but approximate method for embedding an inner region of groundwater flow simulation within a much larger aquifer region where flow can be treated in an approximate fashion. The use of aquifer-influence functions in groundwater flow modeling reduces the size of the computational grid with a corresponding reduction in computer storage and execution time. The Carter-Tracy approximation to the convolution integral enables the aquifer influence function calculation to be made with an additional storage requirement of only two times the number of boundary nodes more than that required for the inner region simulation. It is a good approximation for constant flow rates but is poor for time-varying flow rates where the variation is large relative to the mean. A variety of outer aquifer region geometries, exterior boundary conditions, and flow rate versus potentiometric head relations can be used. The radial, transient-flow case presented is representative. An analytical approximation to the functions of Van Everdingen and Hurst for the dimensionless potentiometric head versus dimensionless time is given.

Water Resources Research↗

Groundwater modeling

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

Book chapter↗

Ground-water flow and solute transport at a municipal landfill site on Long Island, New York – Part 1: Hydrogeology and water quality

Hydrogeology and water quality in a 4 sq mi area surrounding the Brookhaven landfill site in the Town of Brookhaven, New York, were studied in 1981-83. The 60-acre sanitary landfill at the site was excavated in highly permeable glacial outwash that forms the upper glacial aquifer and is lined with a polyvinyl chloride membrane. Groundwater beneath the site is under water table conditions and flows southeast at approximately 1.1 ft/day. Samples from wells downgradient indicate that leachate has entered the aquifer despite the liner. A plume 3,700 ft long , 2,400 ft wide, and at least 90 ft thick was delineated based on specific conductance data. Water quality in the Magothy aquifer and in Beaverdam Creek, a groundwater fed stream 2,000 ft southeast of the landfill site, does not appear to be affected. (USGS)

New York↗

Estimation of regional-scale groundwater flow properties in the Bengal Basin of India and Bangladesh

Quantitative evaluation of management strategies for long-term supply of safe groundwater for drinking from the Bengal Basin aquifer (India and Bangladesh) requires estimation of the large-scale hydrogeologic properties that control flow. The Basin consists of a stratified, heterogeneous sequence of sediments with aquitards that may separate aquifers locally, but evidence does not support existence of regional confining units. Considered at a large scale, the Basin may be aptly described as a single aquifer with higher horizontal than vertical hydraulic conductivity. Though data are sparse, estimation of regional-scale aquifer properties is possible from three existing data types: hydraulic heads, 14C concentrations, and driller logs. Estimation is carried out with inverse groundwater modeling using measured heads, by model calibration using estimated water ages based on 14C, and by statistical analysis of driller logs. Similar estimates of hydraulic conductivities result from all three data types; a resulting typical value of vertical anisotropy (ratio of horizontal to vertical conductivity) is 104. The vertical anisotropy estimate is supported by simulation of flow through geostatistical fields consistent with driller log data. The high estimated value of vertical anisotropy in hydraulic conductivity indicates that even disconnected aquitards, if numerous, can strongly control the equivalent hydraulic parameters of an aquifer system. ?? US Government 2009.

Hydrogeology Journal↗

Variables that affect agricultural chemicals in groundwater in Nebraska

Agricultural chemicals from nonpoint sources in groundwater are present in the major provinces of the High Plains aquifer in Nebraska. Nitrate and triazine-herbicide concentrations in groundwater were assessed to establish preliminary relations between these constituents and selected hydrogeologic, climatic, and land-use variables. Also, macropore flow paths were measured in an attempt to delineate their contribution to non-point source pollution from the study areas. Water from 82 wells in six study areas was analyzed for nitrate; water from 57 of the 82 wells was analyzed for triazine herbicides. Twenty-one independent variables were identified that could potentially affect chemical concentrations in groundwater. Data for 9 of 21 independent variables suspected of affecting concentrations of nitrate and triazine herbicides in groundwater were collected from the well sites. The nine variables and their measured ranges were hydraulic gradient, 0.0006–0.0053; hydraulic conductivity, 1.5–45.4 m (5–149 ft) per day; specific discharge, 0.004–0.091 m (0.0128–0.2998 ft) per day; depth to water, 0.91–76 m (3–250 ft); well depth, 12–168 m (40–550 ft); annual precipitation, 30–100 cm (12.0–39.3 in.); soil permeability, 1.9–23 cm (0.76–9.0 in.); irrigation-well density, 0–8 irrigation wells per 2.59 km 2 (1 square mile); and annual nitrogen fertilizer use, 0–118 kg (0–260 lb) of nitrogen per acre. Macropore flow is listed in percent, average per study area based on determinations from dye studies. In this instance, macropore flow is used to also entail preferential flow paths. Nitrate concentrations ranged from 0.1 to 45 mgL −1 . Triazine-herbicide concentrations were detected in samples from five of the six study areas in concentrations ranging from 0.1 to 2.3 μL −1 . Analysis indicated that there were significant differences in nitrate concentrations (averages-at 95 % confidence level using Kendall Test) among the six study areas; no significant differences in triazineherbicide concentrations were found. Concentrations of nitrate and triazine herbicide were determined (using contingency-table analysis), to be significantly larger in more intensively irrigated areas compared to less intensively irrigated areas. Preliminary correlations with the independent variables and nitrate concentrations indicated significant relations at the 95%confidence level with variables hydraulic conductivity, well depth, and irrigation well density. Correlations with triazine-herbicide concentrations indicated significant relations with hydraulic conductivity, specific discharge, well depth, annual precipitation, and irrigation well density, as well as nitrate concentrations. Simple multiple-regression technique indicated that well depth and density and fertilizer use explained about 51 % of the variation in nitrate concentrations. Specific discharge and well depth explained about 60 % of the variation in triazine-herbicide concentrations. Macropore flow paths and specific discharge explained 84 % of the total variation in triazineherbicide concentrations. The use of trade names in this report is for identification purposes only and does not constitute endorsement by the U.S. Geological Survey.

Nebraska↗

Aquifers in Cretaceous rocks of the central Coastal Plain of North Carolina

Aquifers in rocks of Cretaceous age are the major source of groundwater for public supplies in the central Coastal Plain. These aquifers consist of sand, gravel, and limestone beds of the Peedee, Black Creek, and the upper and lower Cape Fear aquifers, each separated by a confining unit composed of clay and silt beds. The aquifers and confining units (1) rest upon crystalline basement rocks; (2) dip and thicken to the east-southeast; and (3) are overlain by younger aquifers and confining units in deposits of Quaternary and Tertiary age. The top of the uppermost aquifer, the Peedee, ranges from 122 ft above to 595 ft below sea level. The maximum thickness of all aquifers and confining units in Cretaceous rocks is more than 1,600 ft. Aquifers and confining units were defined and correlated for this report using 125 geophysical logs and accompanying drillers' logs, water level data, and water quality information regarding chloride concentrations in water. This analysis allowed the construction of seven hydrogeologic cross sections that depict the continuity of all the aquifers and confining units. These cross sections also show water levels and chloride concentrations in water from various test intervals and describe where chloride concentration in water exceeds 250 mg/L within each aquifer. Detailed maps of each Cretaceous aquifer show altitude of its top, thickness, areas of selected sand percentages, transition zones from fresh to saltwater, and the thickness of the confining unit. H ydrogeologic data for all aquifers and confining units are presented in tabular form.

North Carolina↗

Sources of groundwater based on Helium analyses in and near the freshwater/saline-water transition zone of the San Antonio segment of the Edwards Aquifer, South-Central Texas, 2002-03

This report evaluates dissolved noble gas data, specifically helium-3 and helium-4, collected by the U.S. Geological Survey, in cooperation with the San Antonio Water System, during 2002-03. Helium analyses are used to provide insight into the sources of groundwater in the freshwater/saline-water transition zone of the San Antonio segment of the Edwards aquifer. Sixty-nine dissolved gas samples were collected from 19 monitoring wells (categorized as fresh, transitional, or saline on the basis of dissolved solids concentration in samples from the wells or from fluid-profile logging of the boreholes) arranged in five transects, with one exception, across the freshwater/saline-water interface (the 1,000-milligrams-per-liter dissolved solids concentration threshold) of the Edwards aquifer. The concentration of helium-4 (the dominant isotope in atmospheric and terrigenic helium) in samples ranged from 63 microcubic centimeters per kilogram at standard temperature (20 degrees Celsius) and pressure (1 atmosphere) in a well in the East Uvalde transect to 160,587 microcubic centimeters per kilogram at standard temperature and pressure in a well in the Kyle transect. Helium-4 concentrations in the 10 saline wells generally increase from the western transects to the eastern transects. Increasing helium-4 concentrations from southwest to northeast in the transition zone, indicating increasing residence time of groundwater from southwest to northeast, is consistent with the longstanding conceptualization of the Edwards aquifer in which water recharges in the southwest, flows generally northeasterly (including in the transition zone, although more slowly than in the fresh-water zone), and discharges at major springs in the northeast. Excess helium-4 was greater than 1,000 percent for 60 of the 69 samples, indicating that terrigenic helium is largely present and that most of the excess helium-4 comes from sources other than the atmosphere. The helium data of this report cannot be used to identify sources of groundwater in and near the transition zone of the Edwards aquifer in terms of specific geologic (stratigraphic) units or hydrogeologic units (aquifers or confining units). However, the data indicate that the source or sources of the helium, and thus the water in which the helium is dissolved, in the transition zone are mostly terrigenic in origin rather than atmospheric. Whether most helium in and near the transition zone of the Edwards aquifer originated either in rocks outside the transition zone and at depth or in the adjacent Trinity aquifer is uncertain; but most of the helium in the transition zone had to enter the transition zone from the Trinity aquifer because the Trinity aquifer is the hydrogeologic unit immediately beneath and laterally adjacent to the transition zone of the Edwards aquifer. Thus the helium data support a hypothesis of sufficient hydraulic connection between the Trinity and Edwards aquifers to allow movement of water from the Trinity aquifer to the transition zone of the Edwards aquifer.

Texas↗

Chemical evolution of groundwater near a sinkhole lake, northern Florida: 2. Chemical patterns, mass-transfer modeling, and rates of chemical reactions

Chemical patterns along evolutionary groundwater flow paths in silicate and carbonate aquifers were interpreted using solute tracers, carbon and sulfur isotopes, and mass balance reaction modeling for a complex hydrologic system involving groundwater inflow to and outflow from a sinkhole lake in northern Florida. Rates of dominant reactions along defined flow paths were estimated from modeled mass transfer and ages obtained from CFC-modeled recharge dates. Groundwater upgradient from Lake Barco remains oxic as it moves downward, reacting with silicate minerals in a system open to carbon dioxide (CO 2 ), producing only small increases in dissolved species. Beneath and downgradient of Lake Barco the oxic groundwater mixes with lake water leakage in a highly reducing, silicate-carbonate mineral environment. A mixing model, developed for anoxic groundwater downgradient from the lake, accounted for the observed chemical and isotopic composition by combining different proportions of lake water leakage and infiltrating meteoric water. The evolution of major ion chemistry and the 13 C isotopic composition of dissolved carbon species in groundwater downgradient from the lake can be explained by the aerobic oxidation of organic matter in the lake, anaerobic microbial oxidation of organic carbon, and incongruent dissolution of smectite minerals to kaolinite. The dominant process for the generation of methane was by the CO 2 reduction pathway based on the isotopic composition of hydrogen (δ 2 H(CH 4 ) = −186 to −234‰) and carbon (δ 13 C(CH 4 ) = −65.7 to −72.3‰). Rates of microbial metabolism of organic matter, estimated from the mass transfer reaction models, ranged from 0.0047 to 0.039 mmol L −1 yr −1 for groundwater downgradient from the lake.

Water Resources Research↗

Boron contents and isotopic compositions of hog manure, selected fertilizers, and water in Minnesota

Boron-isotope (δ 11 B) values may be useful as surrogate tracers of contaminants and indicators of water mixing in agricultural settings. This paper characterizes the B contents and isotopic compositions of hog manure and selected fertilizers, and presents δ 11 B data for ground and surface water from two agricultural areas. Boron concentrations in dry hog manure averaged 61 mg/kg and in commercial fertilizers ranged from below detection limits in some brands of ammonium nitrate and urea to 382 mg/kg in magnesium sulfate. Values of δ 11 B of untreated hog manure ranged from 7.2 to 11.2o/oo and of N fertilizers were −2.0 to 0.7o/oo. In 22 groundwater samples from a sand-plain aquifer in east-central Minnesota, B concentrations averaged 0.04 mg/L and δ 11 B values ranged from 2.3 to 41.5o/oo. Groundwater beneath a hog feedlot and a cultivated field where hog manure was applied had B-isotope compositions consistent with the water containing hog-manure leachate. In a 775-km 2 watershed with silty-loam soils in southcentral Minnesota: 18 samples of subsurface drainage from corn ( Zea mays L.) and soybean ( Glycine max L. Merr.) fields had average B concentrations of 0.06 mg/L and δ 11 B values of 5.3 to 15.1o/oo; 27 stream samples had average B concentrations of 0.05 mg/L and δ 11 B values of 1.0 to 19.0o/oo; and eight groundwater samples had average B concentrations of 0.09 mg/L and δ 11 B values of −0.3 to 23.0o/oo. Values of δ 11 B and B concentrations, when plotted against one another, define a curved mixing trend that suggests subsurface drainage and stream water contain mixtures of B from shallow and deep groundwater.

Minnesota↗

Application of a soil-water-balance model to estimate annual groundwater recharge for Long Island, New York, 1900–2019

A soil-water-balance (SWB) model was developed for Long Island, New York, to estimate the potential amount of annual groundwater recharge to the Long Island aquifer system from 1900 to 2019. The SWB model program is a computer code based on a modified Thornthwaite-Mather SWB approach and uses spatially and temporally distributed meteorological, land-cover, and soil properties as input to compute potential daily groundwater recharge. Simulated outputs indicate that island-wide potential groundwater recharge trends, as a percentage of precipitation, have increased approximately 3 percent during the 120-year period. The simulated results account for both climatic and land-cover changes that have occurred during the period. A change from undeveloped (forested land cover) to low- and medium-density residential land cover or land use increased potential groundwater recharge because of a decrease in evapotranspiration. During the 30-year period from 1900 to 1930, the simulated potential average groundwater recharge rate on Long Island was estimated to be 18.50 inches per year (in/yr), or a total of 1,243 million gallons per day, during the 30-year period from 1985 to 2015, the simulated potential average groundwater recharge rate estimate increased to 20.73 in/yr (a total of around 1,393 million gallons per day). During the 1900–2019 simulation period, the potential average annual groundwater recharge rate was about 19.24 in/yr. The data for that period included values for a 3-year meteorological drought from 1963 to 1965, where the mean precipitation was about 26.5 percent lower than the long-term average of 46.7 in/yr, and the potential groundwater recharge rate was about 12.3 in/yr. During a 3-year wet period from 1982 to 1984, where mean precipitation was about 19.6 percent higher than the long-term average, the estimated potential groundwater recharge rate was about 26.8 in/yr.

New York↗

Influence of water chemistry and travel distance on bacteriophage PRD-1 transport in a sandy aquifer

Experiments were conducted to evaluate the impact of groundwater chemistry and travel distance on the transport and fate behavior of PRD-1, a bacteriophage employed as a surrogate tracer for pathogenic enteric viruses. The experiments were conducted in the unconfined aquifer at the United States Geological Survey Cape Cod Toxic-Substances Hydrology Research Site in Falmouth, Massachusetts. The transport behavior of bromide (Br − ) and PRD-1 were evaluated in a sewage-effluent contaminated zone and a shallower uncontaminated zone at this site. Several multilevel sampling devices located along a 13-m transect were used to collect vertically discrete samples to examine longitudinal and vertical variability of PRD-1 retardation and attenuation. The concentration of viable bacteriophage in the aqueous phase decreased greatly during the first few meters of transport. This decrease is attributed to a combination of colloid filtration (attachment) and inactivation. The removal was greater (10 −12 relative recovery) and occurred within the first meter for the uncontaminated zone, whereas it was lesser (10 −9 relative recovery) and occurred over 4 m in the contaminated zone. The lesser removal observed for the contaminated zone is attributed to the influence of sorbed and dissolved organic matter, phosphate, and other anions, which are present in higher concentrations in the contaminated zone, on PRD-1 attachment. After the initial decrease, the aqueous PRD-1 concentrations remained essentially constant in both zones for the remainder of the tests (total travel distances of 13 m), irrespective of variations in geochemical properties within and between the two zones. The viable, mobile PRD-1 particles traveled at nearly the rate of bromide, which was used as a non-reactive tracer. The results of this study indicate that a small fraction of viable virus particles may persist in the aqueous phase and travel significant distances in the subsurface environment.

Water Research↗

Characterization of potential transport pathways and implications for groundwater management near an anticline in the Central Basin area, Los Angeles County, California

The Central Groundwater Basin (Central Basin) of southern Los Angeles County includes ~280 mi 2 of the Los Angeles Coastal Plain and serves as the primary source of water for more than two million residents. In the Santa Fe Springs–Whittier–Norwalk area, located in the northeastern part of the basin, several sources of volatile organic compounds have been identified. The volatile organic compunds are thought to have contributed to a large, commingled contaminant plume in groundwater that extends south-southwest downgradient from the Omega Chemical Corporation Superfund Site across folded geologic strata, known as the Santa Fe Springs Anticline. A multifaceted study—that incorporated a three-dimensional sequence-stratigraphic geologic model, two-dimensional groundwater particle-tracking simulations, and new groundwater chemistry data—was conducted to gain insight into the geologic and hydrologic controls on contaminant migration in the study area and to assess the potential for this shallow groundwater contamination to migrate into producing aquifer zones. Conceptual flow models were developed along a flow-parallel cross section based on the modeled stratigraphic architecture, observed geochemistry, and numerical model simulations that generally agree with observed water levels and contaminant distributions. These models predict that contaminants introduced into groundwater at shallow depths near the Omega Chemical Corporation Superfund Site and along the study cross section will likely migrate downgradient to depths intercepted by public supply wells. These conclusions, however, are subject to limitations and simplifications inherent in the modeling approaches used, as well as a significant scarcity of available geologic and hydrogeochemical information at depth and in the downgradient parts of the study area.

California↗

A finite-element simulation model for saturated-unsaturated, fluid-density-dependent ground-water flow with energy transport or chemically- reactive single-species solute transport

SUTRA (Saturated-Unsaturated Transport) is a computer program which can be used to simulate the movement of fluid and the transport of either energy or dissolved substances in a subsurface environment. The model employs a two-dimensional hybrid finite-element and integrated-finite-difference method to approximate the governing equations that describe the two interdependent processes that are simulated by SUTRA: (1) fluid density-dependent saturated or unsaturated groundwater flow, and either (2a) transport of a solute in the groundwater, in which the solute may be subject to: equilibrium adsorption on the porous matrix, and both first-order and zero-order production or decay, or, (2b) transport of thermal energy in the groundwater and solid matrix of the aquifer. SUTRA provides, as the primary calculated results, fluid pressures and either solute concentrations or temperatures, as they vary with time, everywhere in the simulated subsurface system. SUTRA may also be used to simulate simpler subsets of the above process. SUTRA may be employed for areal and cross-sectional models of saturated groundwater flow systems, and for cross-sectional models of unsaturated zone flow. Solute transport simulation using SUTRA may be used to simulate natural or man-induced chemical transport, solute sorption, production and decay. SUTRA may be used for simulation of variable density leachate movement, and for cross-sectional simulation of salt-water intrusion in aquifers at near-well or regional scales, with either dispersed or relatively sharp transition zones between fresh water and salt water. SUTRA energy transport simulation may be employed to model thermal regimes in aquifers, subsurface heat conduction, aquifer thermal energy storage systems, geothermal reservoirs, thermal pollution of aquifers, and natural hydrogeologic convection systems. (USGS)

Water-Resources Investigations Report↗