USGS Science⌕ Search

SEARCH · USGS Science

Results for “Groundwater”

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

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

At least 1,315 records · Page 73Linked to original sources
USGS2013

Hydrogeologic framework and estimates of groundwater storage for the Hualapai Valley, Detrital Valley, and Sacramento Valley basins, Mohave County, Arizona

We have investigated the hydrogeology of the Hualapai Valley, Detrital Valley, and Sacramento Valley basins of Mohave County in northwestern Arizona to develop a better understanding of groundwater storage within the basin fill aquifers. In our investigation we used geologic maps, well-log data, and geophysical surveys to delineate the sedimentary textures and lithology of the basin fill. We used gravity data to construct a basin geometry model that defines smaller subbasins within the larger basins, and airborne transient-electromagnetic modeled results along with well-log lithology data to infer the subsurface distribution of basin fill within the subbasins. Hydrogeologic units (HGUs) are delineated within the subbasins on the basis of the inferred lithology of saturated basin fill. We used the extent and size of HGUs to estimate groundwater storage to depths of 400 meters (m) below land surface (bls). The basin geometry model for the Hualapai Valley basin consists of three subbasins: the Kingman, Hualapai, and southern Gregg subbasins. In the Kingman subbasin, which is estimated to be 1,200 m deep, saturated basin fill consists of a mixture of fine- to coarse-grained sedimentary deposits. The Hualapai subbasin, which is the largest of the subbasins, contains a thick halite body from about 400 m to about 4,300 m bls. Saturated basin fill overlying the salt body consists predominately of fine-grained older playa deposits. In the southern Gregg subbasin, which is estimated to be 1,400 m deep, saturated basin fill is interpreted to consist primarily of fine- to coarse-grained sedimentary deposits. Groundwater storage to 400 m bls in the Hualapai Valley basin is estimated to be 14.1 cubic kilometers (km 3 ). The basin geometry model for the Detrital Valley basin consists of three subbasins: northern Detrital, central Detrital, and southern Detrital subbasins. The northern and central Detrital subbasins are characterized by a predominance of playa evaporite and fine-grained clastic deposits; evaporite deposits in the northern Detrital subbasin include halite. The northern Detrital subbasin is estimated to be 600 m deep and the middle Detrital subbasin is estimated to be 700 m deep. The southern Detrital subbasin, which is estimated to be 1,500 m deep, is characterized by a mixture of fine- to coarse-grained basin fill deposits. Groundwater storage to 400 m bls in the Detrital Valley basin is estimated to be 9.8 km 3 . The basin geometry model for the Sacramento Valley basin consists of three subbasins: the Chloride, Golden Valley, and Dutch Flat subbasins. The Chloride subbasin, which is estimated to be 900 m deep, is characterized by fine- to coarse-grained basin fill deposits. In the Golden Valley subbasin, which is elongated north-south, and is estimated to be 1,300 m deep, basin fill includes fine-grained sedimentary deposits overlain by coarse-grained sedimentary deposits in much of the subbasin. The Dutch Flat subbasin is estimated to be 2,600 m deep, and well-log lithologic data suggest that the basin fill consists of interlayers of gravel, sand, and clay. Groundwater storage to 400 m bls in the Sacramento Valley basin is estimated to be 35.1 km 3 .

Arizona↗
USGS2013

Groundwater depletion in the United States (1900−2008)

A natural consequence of groundwater withdrawals is the removal of water from subsurface storage, but the overall rates and magnitude of groundwater depletion in the United States are not well characterized. This study evaluates long-term cumulative depletion volumes in 40 separate aquifers or areas and one land use category in the United States, bringing together information from the literature and from new analyses. Depletion is directly calculated using calibrated groundwater models, analytical approaches, or volumetric budget analyses for multiple aquifer systems. Estimated groundwater depletion in the United States during 1900–2008 totals approximately 1,000 cubic kilometers (km 3 ). Furthermore, the rate of groundwater depletion has increased markedly since about 1950, with maximum rates occurring during the most recent period (2000–2008) when the depletion rate averaged almost 25 km 3 per year (compared to 9.2 km 3 per year averaged over the 1900–2008 timeframe).

Scientific Investigations Report↗
USGS2013

Relation between organic-wastewater compounds, groundwater geochemistry, and well characteristics for selected wells in Lansing, Michigan

In 2010, groundwater from 20 Lansing Board of Water and Light (BWL) production wells was tested for 69 organic-wastewater compounds (OWCs). The OWCs detected in one-half of the sampled wells are widely used in industrial and environmental applications and commonly occur in many wastes and stormwater. To identify factors that contribute to the occurrence of these constituents in BWL wells, the U.S. Geological Survey (USGS) interpreted the results of these analyses and related detections of OWCs to local characteristics and groundwater geochemistry. Analysis of groundwater-chemistry data collected by the BWL during routine monitoring from 1969 to 2011 indicates that the geochemistry of the BWL wells has changed over time, with the major difference being an increase in sodium and chloride. The concentrations of sodium and chloride were positively correlated to frequency of OWC detections. The BWL wells studied are all completed in the Saginaw aquifer, which consists of water-bearing sandstones of Pennsylvanian age. The Saginaw aquifer is underlain by the Parma-Bayport aquifer, and overlain by the Glacial aquifer. Two possible sources of sodium and chloride were evaluated: basin brines by way of the Parma-Bayport aquifer, and surficial sources by way of the Glacial aquifer. To determine if water from the underlying aquifer had influenced well-water geochemistry over time, the total dissolved solids concentration and changes in major ion concentrations were examined with respect to well depth, age, and pumping rate. To address a possible surficial source of sodium and chloride, 25 well, aquifer, or hydrologic characteristics, and 2 groundwater geochemistry variables that might influence whether, or the rate at which, water from the land surface could reach each well were compared to OWC detections and well chemistry. The statistical tests performed during this study, using available variables, indicated that reduced time of travel of water from the land surface to the well opening was significantly correlated with detections of OWCs. No specific well or aquifer characteristic was correlated with OWC detections; however, wells with detections tended to have less modeled confining material thickness (as simulated in the regional groundwater flow model), which is an estimate of the amount of clay or shale between the Glacial and Saginaw aquifers. Additional analyses and collection of other data would be required to more conclusively identify the source and to determine the potential vulnerability of other wells because each BWL well may have a somewhat unique set of characteristics that governs its response to pumping. Therefore, it is possible that a relevant explanatory variable was not included in this analysis. The current patterns of geochemistry, and the relation between these patterns and volume of pumpage for the BWL wells, indicates other wells may be susceptible to OWCs in the future.

Michigan↗
USGS2015

Potential groundwater recharge for the State of Minnesota using the Soil-Water-Balance model, 1996-2010

Groundwater recharge is one of the most difficult components of a water budget to ascertain, yet is an important boundary condition necessary for the quantification of water resources. In Minnesota, improved estimates of recharge are necessary because approximately 75 percent of drinking water and 90 percent of agricultural irrigation water in Minnesota are supplied from groundwater. The water that is withdrawn must be supplied by some combination of (1) increased recharge, (2) decreased discharge to streams, lakes, and other surface-water bodies, and (3) removal of water that was stored in the system. Recent pressure on groundwater resources has highlighted the need to provide more accurate recharge estimates for various tools that can assess the sustainability of long-term water use. As part of this effort, the U.S. Geological Survey, in cooperation with the Minnesota Pollution Control Agency, used the Soil-Water-Balance model to calculate gridded estimates of potential groundwater recharge across Minnesota for 1996‒2010 at a 1-kilometer (0.621-mile) resolution. The potential groundwater recharge estimates calculated for Minnesota from the Soil-Water Balance model included gridded values (1-kilometer resolution) of annual mean estimates (that is, the means for individual years from 1996 through 2010) and mean annual estimates (that is, the mean for the 15-year period 1996−2010). The Soil-Water-Balance model uses a modified Thornthwaite-Mather soil-water-balance approach, with components of the soil-water balance calculated on a daily basis. A key advantage of this approach includes the use of commonly available geographic information system data layers that incorporate land cover, soil properties, and daily meteorological data to produce temporally and spatially variable gridded estimates of potential recharge. The Soil-Water-Balance model was calibrated by using a combination of parameter estimation techniques, making manual adjustments of model parameters, and using parameter values from previously published Soil-Water-Balance models. Each calibration simulation compared the potential recharge estimate from the model against base-flow estimates derived from three separate hydrograph separation techniques. A total of 35 Minnesota watersheds were selected for the model calibration. Meteorological data necessary for the model included daily precipitation, minimum daily temperature, and maximum daily temperature. All of the meteorological data were provided by the Daymet dataset, which included daily continuous surfaces of key climatological data. Land-cover data were provided by the 2001 and 2006 National Land Cover Database: the 2001 classification was used from 1994 through 2003, and the 2006 classification was used from 2004 through 2010. Soil data used in the model included hydrologic soils group and the available soil-water capacity. These soil data were obtained from the Natural Resources Conservation Service Soil Survey Geographic (SSURGO) database and the State Soil Geographic (STATSGO) database. The statewide mean annual potential recharge rate from 1996–2010 was 4.9 inches per year. Potential recharge estimates increased from west to east across Minnesota. The mean annual potential recharge estimates across Minnesota at a 1-km resolution for the overall simulation period (1996–2010) ranged from less than 0.1 to 17.8 inches per year. Some of the lowest potential recharge rates for the simulation period were in the Red River of the North Basin of northwestern Minnesota, and generally were between 1.0 and 1.5 inches per year. The highest potential recharge rates were in northeastern Minnesota and the Anoka Sand Plain in central Minnesota. Eighty-eight percent of the potential recharge rates (by grid cell) were between 2 and 8 inches per year from 1996–2010. Only about 3 percent of all the potential recharge estimates (by grid cell) were less than 2 inches per year, and 9 percent of estimates were greater than 8 inches per year. On an annual basis, however, potential recharge rates were as high as 27.2 inches per year. The highest annual mean recharge estimate across the State was for 2010, and the lowest mean recharge estimate was for 2003. Although precipitation variability partially explained the annual differences in potential recharge estimates, precipitation alone did not account for these differences, and other factors such as antecedent moisture conditions likely were important. Also, because precipitation gradients across the State can vary from year to year, the dominant land-cover class and hydrologic soil group combinations for a particular region had a large effect on the resulting potential recharge value. During 1996–2010, April had the greatest monthly mean potential recharge compared to all other months, accounting for a mean of 30 percent of annual potential recharge in this single month.

Minnesota↗
USGS2015

Groundwater flow in the Brunswick/Glynn County area, Georgia, 2000-04

An existing regional steady-state model for coastal Georgia, and parts of South Carolina and Florida, was revised to evaluate the local effects of pumping on the migration of high chloride (saline) water in the Upper Floridan aquifer located in the Brunswick/Glynn County, Georgia (Ga.) area. Revisions were focused on enhancing the horizontal and vertical resolution of the regional model grid in the vicinity of saline water. Modifications to the regional model consisted of (1) limiting grid size to a maximum of 500 feet (ft) per side in the vicinity of chloride contamination; (2) representing the upper and lower Brunswick aquifers with distinct model layers; (3) similarly, representing upper and lower water-bearing zones of the Upper Floridan aquifer with distinct model layers in Glynn and Camden Counties, Ga.; and (4) establishing new hydraulic-property zones in the Upper Floridan aquifer. The revised model simulated steady-state conditions that were assumed to exist during 2000 and 2004. Calibration of the revised steady-state model using pumping rates from 2000 indicates a "good" match (±10 ft) based on 181 observations, with median residuals (simulated minus observed water levels) in each of the active model layers ranging from -8.62 to 4.67 ft, and root mean square error (RMSE) ranging from 10.9 to 11.4 ft. In the Brunswick/Glynn County area, groundwater-level residuals in the upper water-bearing zone of the Upper Floridan aquifer (layer 7) indicate an "excellent" match (±5 ft) based on 41 observations with a median residual of -0.35 ft and RMSE of 4.32 ft. Calibration of the revised steady-state model using 2004 pumping rates and adjusted specified-head input values in the Floridan aquifer system indicates a "good" match (-10 ft) based on 88 observations, with median residuals in each of the active model layers ranging from -6.31 to -2.05 ft, and RMSE ranging from -6.95 to 14.5 ft. In the Brunswick/Glynn County area, groundwater-level residuals in the upper water-bearing zone of the Upper Floridan aquifer (layer 7) indicate an "excellent" match (±5 ft) based on 32 observations with a median residual of -1.50 ft and RMSE of 5.34 ft. Simulated potentiometric surfaces for 2000 and 2004 indicate coastward groundwater flow in the Upper and Lower Floridan aquifers influenced by pumping centers at Savannah, Jesup, and Brunswick, Ga., and indicate steep potentiometric gradients to the west and north of the Gulf Trough. In the Brunswick/Glynn County area, simulated industrial production wells located north of downtown Brunswick intercept local groundwater flow in the upper and lower water-bearing zones of the Upper Floridan aquifer and have created a cone of depression that locally alters the regional coastward flow direction. Maps of simulated water-level change during the 2000-04 period show differences in groundwater levels in the Upper Floridan aquifer that range from -2.5 ft to more than 5 ft in areas of coastal Georgia, and more than 20 ft near the Georgia-Florida State Line. Positive values indicate higher simulated water levels during 2004 than during 2000, which were caused by reduced pumping in the Upper Floridan aquifer prompted by the shutdown of a paper mill near the southern model boundary in 2002 and increased recharge following a prolonged drought during 1998-2002. Simulated potentiometric profiles for 2000 and 2004 were used to evaluate the potentiometric gradients in the upper water-bearing zone of the Upper Floridan aquifer (layer 7) near the chloride plume in the downtown Brunswick area. Four potentiometric profiles were constructed for 2000 to compare the simulated and observed water levels in 13 wells and were oriented outward from a primary well field. The simulated potentiometric gradients from the four profiles for 2000 ranged from 3.6 to 5.2 feet per mile (ft/mi) compared to observed values ranging from 4.1 to 5.6 ft/mi. The five potentiometric profiles constructed for 2004 allowed for a similar comparison using simulated and observed water levels in 18 wells. The simulated potentiometric gradients from the five profiles for 2000 ranged from 3.6 to 11.1 ft/mi compared to observed values ranging from 3.8 to 10.2 ft/mi. Simulated potentiometric gradients were higher for 2004 than for 2000 because of the inclusion of a well located within the cone of depression near downtown Brunswick. Composite-scaled sensitivities of the model parameters indicate the revised model is most sensitive to pumping rates, followed by the horizontal hydraulic conductivity in the Upper Floridan aquifer for zones along coastal Georgia. The revised model is least sensitive to the horizontal hydraulic conductivity of the confining units and vertical hydraulic conductivity of the aquifers. For parameters defined by hydraulic-property zones in the upper and lower water-bearing zones of the Upper Floridan aquifer, such as horizontal hydraulic conductivity, model sensitivity was not as great in the Brunswick/Glynn County area as other areas along coastal Georgia. The model exhibited more sensitivity to these parameters however, than to parameters representing the majority of zones defining the vertical hydraulic conductivity of the confining units, which originally were assumed to govern upward migration of chloride contamination into this aquifer. Analysis of simulated water-budget components for 2000 and 2004 indicate that specified-head boundaries in the Floridan aquifer system to the south and southwest of the regional model area control about 70 percent of inflows and nearly 50 percent of outflows to the model region. Other water-budget components indicate an 80-million-gallon-per-day decrease in pumping from the Floridan aquifer system during this period.

Georgia↗
USGS2016

Variability of surface-water quantity and quality and shallow groundwater levels and quality within the Rio Grande Project Area, New Mexico and Texas, 2009–13

Drought conditions during the study period of January 1, 2009, to September 30, 2013, caused a reduction in surface-water releases from water-supply storage infrastructure of the Rio Grande Project, which led to changes in surface-water and groundwater (conjunctive) use in downstream agricultural alluvial valleys. Surface water and groundwater in the agriculturally dominated alluvial Rincon and Mesilla Valleys were investigated in this study to measure the influence of drought and subsequent change in conjunctive water use on quantity and quality of these water resources. In 2013, the U.S. Geological Survey, in cooperation with the New Mexico Environment Department and the New Mexico Interstate Stream Commission, began a study to (1) calculate dissolved-solids loads over the study period at streamgages in the study area where data are available, (2) assess the temporal variability of dissolved-solids loads at and between each streamgage where data are available, and (3) relate the spatiotemporal variability of shallow groundwater data (groundwater levels and quality) within the alluvial valleys of the study area to spatiotemporal variability of surface-water data over the study period. This assessment included the calculation of surface-water dissolved-solids loads at streamgages as well as a mass-balance approach to measure the change in salt load between these streamgages. Bimodal surface-water discharge data led to a temporally-dynamic volumetric definition of release and nonrelease seasons. Continuous surface-water discharge and water-quality data from three streamgages on the Rio Grande were used to calculate daily dissolved-solids loads over the study period, and the results were aggregated annually and seasonally. Results show the majority of dissolved-solids loading occurs during release season; however, decreased duration of the release season over the 5-year study period has resulted in a decrease of the total annual loads at each streamgage. Calculation of the change of salt loads using a mass-balance approach was applied between streamgages. Results from these calculations suggest differing responses to releases in the Rincon and Mesilla Valleys over the period of study; there is a decreasing sink of salt in the Rincon Valley whereas there is an increasing sink of salt in the Mesilla Valley. Daily groundwater-level and water-quality data from shallow wells within the two alluvial valleys show spatial heterogeneity of water quality over the study period. Mass-balance salt-loading trends during the study period are similar to previous trends during the 1950s drought as well as a wet period in the 1980s. The similarity of salt-loading trends from the 1950s, 1980s, and 2000s independent of the climate indicates salt loading in this hydrologic setting may be driven by water-use practices rather than a single climatic variable.

New Mexico, Texas↗
USGS2016

Simulating groundwater flow in karst aquifers with distributed parameter models—Comparison of porous-equivalent media and hybrid flow approaches

Understanding karst aquifers, for purposes of their management and protection, poses unique challenges. Karst aquifers are characterized by groundwater flow through conduits (tertiary porosity), and (or) layers with interconnected pores (secondary porosity) and through intergranular porosity (primary or matrix porosity). Since the late 1960s, advances have been made in the development of numerical computer codes and the use of mathematical model applications towards the understanding of dual (primary [matrix] and secondary [fractures and conduits]) porosity groundwater flow processes, as well as characterization and management of karst aquifers. The Floridan aquifer system (FAS) in Florida and parts of Alabama, Georgia, and South Carolina is composed of a thick sequence of predominantly carbonate rocks. Karst features are present over much of its area, especially in Florida where more than 30 first-magnitude springs occur, numerous sinkholes and submerged conduits have been mapped, and numerous circular lakes within sinkhole depressions are present. Different types of mathematical models have been applied for simulation of the FAS. Most of these models are distributed parameter models based on the assumption that, like a sponge, water flows through connected pores within the aquifer system and can be simulated with the same mathematical methods applied to flow through sand and gravel aquifers; these models are usually referred to as porous-equivalent media models. The partial differential equation solved for groundwater flow is the potential flow equation of fluid mechanics, which is used when flow is dominated by potential energy and has been applied for many fluid problems in which kinetic energy terms are dropped from the differential equation solved. In many groundwater model codes (basic MODFLOW), it is assumed that the water has a constant temperature and density and that flow is laminar, such that kinetic energy has minimal impact on flow. Some models have been developed that incorporate the submerged conduits as a one-dimensional pipe network within the aquifer rather than as discrete, extremely transmissive features in a porous-equivalent medium; these submerged conduit models are usually referred to as hybrid models and may include the capability to simulate both laminar and turbulent flow in the one-dimensional pipe network. Comparisons of the application of a porous-equivalent media model with and without turbulence (MODFLOW-Conduit Flow Process mode 2 and basic MODFLOW, respectively) and a hybrid (MODFLOW-Conduit Flow Process mode 1) model to the Woodville Karst Plain near Tallahassee, Florida, indicated that for annual, monthly, or seasonal average hydrologic conditions, all methods met calibration criteria (matched observed groundwater levels and average flows). Thus, the increased effort required, such as the collection of data on conduit location, to develop a hybrid model and its increased computational burden, is not necessary for simulation of average hydrologic conditions (non-laminar flow effects on simulated head and spring discharge were minimal). However, simulation of a large storm event in the Woodville Karst Plain with daily stress periods indicated that turbulence is important for matching daily springflow hydrographs. Thus, if matching streamflow hydrographs over a storm event is required, the simulation of non-laminar flow and the location of conduits are required. The main challenge in application of the methods and approaches for developing hybrid models relates to the difficulty of mapping conduit networks or having high-quality datasets to calibrate these models. Additionally, hybrid models have long simulation times, which can preclude the use of parameter estimation for calibration. Simulation of contaminant transport that does not account for preferential flow through conduits or extremely permeable zones in any approach is ill-advised. Simulation results in other karst aquifers or other parts of the FAS may differ from the comparison demonstrated herein.

Florida↗
USGS2016

Status of groundwater levels and storage volume in the Equus Beds aquifer near Wichita, Kansas, January 2016

The Equus Beds aquifer in south-central Kansas, which is part of the High Plains aquifer, serves as a source of water for municipal and agricultural users in the area. The city of Wichita has used the Equus Beds aquifer as one of its primary water sources since the 1940s. The aquifer in and around Wichita’s well field reached historically low water levels in 1993, prompting the city to adopt new water-use and conservation strategies to ensure future water supply needs were met. Part of the plan was to initiate a managed aquifer recharge program called the Equus Beds Aquifer Storage and Recovery project. The goal of the managed aquifer recharge program is to artificially recharge the Equus Beds aquifer with treated water from the Little Arkansas River. As part of the Equus Beds Aquifer Storage and Recovery project, the city of Wichita and the U.S. Geological Survey have partnered in a long-term cooperative study to monitor and describe the quantity and quality of the water in the Equus Beds aquifer and the Little Arkansas River. The city of Wichita, the Equus Beds Groundwater Management District No. 2, the Kansas Department of Agriculture–Division of Water Resources, and the U.S. Geological Survey collected groundwater levels in numerous wells screened in the Equus Beds aquifer in the area in and around Wichita’s well field in January 2016. The measurements were used to interpolate potentiometric surfaces for shallow and deep parts of the aquifer in the study area. These potentiometric surfaces were compared with potentiometric surfaces from previous years to estimate changes in water levels and storage volume in the study area. Groundwater levels were generally higher in January 2016 than they were in January 2015. On average, in January 2016, groundwater levels in the shallow part of the aquifer were about 3.4 feet higher and groundwater levels in the deep part of the aquifer were about 3.8 feet higher than in January 2015. The volume of water stored in the study area decreased by about 74,000 acre-feet between predevelopment (the time period before substantial pumpage began in the 1940s) and January 2016; increased by about 121,000 acre-feet between the historic low in 1993 and January 2016; and increased by about 61,000 acre-feet between January 2015 and January 2016. About 62 percent of the storage volume lost between predevelopment and 1993 has been recovered. The increase in storage volume from January 2015 to January 2016 can probably be attributed to less pumping by the city of Wichita and irrigators, more recharge due to higher-than-average precipitation, and higher volumes of artificial recharge in 2015.

Kansas↗
USGS2017

Simulation of groundwater flow in the glacial aquifer system of northeastern Wisconsin with variable model complexity

The U.S. Geological Survey, National Water-Quality Assessment seeks to map estimated intrinsic susceptibility of the glacial aquifer system of the conterminous United States. Improved understanding of the hydrogeologic characteristics that explain spatial patterns of intrinsic susceptibility, commonly inferred from estimates of groundwater age distributions, is sought so that methods used for the estimation process are properly equipped. An important step beyond identifying relevant hydrogeologic datasets, such as glacial geology maps, is to evaluate how incorporation of these resources into process-based models using differing levels of detail could affect resulting simulations of groundwater age distributions and, thus, estimates of intrinsic susceptibility. This report describes the construction and calibration of three groundwater-flow models of northeastern Wisconsin that were developed with differing levels of complexity to provide a framework for subsequent evaluations of the effects of process-based model complexity on estimations of groundwater age distributions for withdrawal wells and streams. Preliminary assessments, which focused on the effects of model complexity on simulated water levels and base flows in the glacial aquifer system, illustrate that simulation of vertical gradients using multiple model layers improves simulated heads more in low-permeability units than in high-permeability units. Moreover, simulation of heterogeneous hydraulic conductivity fields in coarse-grained and some fine-grained glacial materials produced a larger improvement in simulated water levels in the glacial aquifer system compared with simulation of uniform hydraulic conductivity within zones. The relation between base flows and model complexity was less clear; however, the relation generally seemed to follow a similar pattern as water levels. Although increased model complexity resulted in improved calibrations, future application of the models using simulated particle tracking is anticipated to evaluate if these model design considerations are similarly important for understanding the primary modeling objective - to simulate reasonable groundwater age distributions.

Wisconsin↗
USGS2017

An update of hydrologic conditions and distribution of selected constituents in water, eastern Snake River Plain aquifer and perched groundwater zones, Idaho National Laboratory, Idaho, emphasis 2012-15

Since 1952, wastewater discharged to in ltration ponds (also called percolation ponds) and disposal wells at the Idaho National Laboratory (INL) has affected water quality in the eastern Snake River Plain (ESRP) aquifer and perched groundwater zones underlying the INL. The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy, maintains groundwater-monitoring networks at the INL to determine hydrologic trends and to delineate the movement of radiochemical and chemical wastes in the aquifer and in perched groundwater zones. This report presents an analysis of water-level and water-quality data collected from the ESRP aquifer, multilevel monitoring system (MLMS) wells in the ESRP aquifer, and perched groundwater wells in the USGS groundwater monitoring networks during 2012-15.

Idaho↗
USGS2018

Effects of groundwater withdrawals from the Hurricane Fault zone on discharge of saline water from Pah Tempe Springs, Washington County, Utah

Pah Tempe Springs, located in Washington County, Utah, contribute about 95,000 tons of dissolved solids annually along a 1,500-foot gaining reach of the Virgin River. The river gains more than 10 cubic feet per second along the reach as thermal, saline springwater discharges from dozens of orifices located along the riverbed and above the river on both banks. The spring complex discharges from fractured Permian Toroweap Limestone where the river crosses the north-south trending Hurricane Fault. The Bureau of Reclamation Colorado River Basin Salinity Control Program is evaluating the feasibility of capturing and desalinizing the discharge of Pah Tempe Springs to improve downstream water quality in the Virgin River. The most viable plan, identified by the Bureau of Reclamation in early studies, is to capture spring discharge by pumping thermal groundwater from within the Hurricane Fault footwall damage zone and to treat this water prior to returning it to the river. Three multiple-day interference tests were conducted between November 2013 and November 2014, wherein thermal groundwater was pumped from fractured carbonate rock in the fault damage zone at rates of up to 7 cubic feet per second. Pumping periods for these tests lasted approximately 66, 74, and 67 hours, respectively, and the tests occurred with controlled streamflows of approximately 2.0, 3.5, and 24.5 cubic feet per second, respectively, in the Virgin River upstream from the springs reach. Specific conductance, water temperature, and discharge were monitored continuously in the river (upstream and downstream of the springs reach) at selected individual springs, and in the pumping discharge during each of the tests. Water levels were monitored in three observation wells screened in the thermal system. Periodic stream and groundwater samples were analyzed for dissolved-solids concentration and the stable isotopes of oxygen and hydrogen. Additional discrete measurements of field parameters (specific conductance, water temperature, pH, and discharge) were made at up to 26 sites along the springs reach. These data demonstrate the interaction between the saline, thermal groundwater system and the Virgin River, and provide estimates of reductions in dissolved-solids loads to the river. The interference tests show that pumping thermal groundwater from the shallow carbonate aquifer adjacent to the springs is effective at capturing high dissolved-solids loads discharging from Pah Tempe Springs before they enter the Virgin River. Discharge measurements made in the Virgin River downstream of the springs reach show that streamflow is reduced by approximately the amount pumped, indicating that complete capture of thermal discharge is possible. During the February 2014 test, the dissolved-solids load removed by pumping (190 tons per day) was approximately equal to the dissolved-solids load reduction observed in the river below the springs reach, indicating near 100-percent efficient capture of spring-sourced dissolved solids. However, an observed decrease in temperature and specific conductance of the pumping discharge during the high-flow test in November 2014 showed that capture of the cool, fresh river water can occur and is more likely at a higher stage in the Virgin River.

Utah↗
USGS2019

Simulation of the regional groundwater-flow system in the St. Louis River basin, Minnesota

The St. Louis River Basin (SLRB) covers 3,600 square miles in northeastern Minnesota, with headwaters in the Mesabi Range and extensive wetlands and lakes throughout the basin. To better understand the regional groundwater system in the SLRB, a two-dimensional, steady-state groundwater-flow model of the SLRB was developed by the U.S. Geological Survey, in cooperation with the Minnesota Ojibwe Bands, using the analytic-element computer code GFLOW. The parameter-estimation software suite PEST was used to obtain a best fit of the modeled to measured groundwater levels and streamflows. The calibrated regional model was locally refined to create a smaller version of the model, the central SLRB model, that was used to evaluate hydrologic effects from extensive ditching in wetlands of the central SLRB. The refinements included adding ditches that were not represented in the regional model and modifying the aquifer base elevation to be more representative of the localized area. The central SLRB model was recalibrated to better match the distribution of mapped wetlands. Two scenarios were run of the central SLRB model: one with ditches and one without ditches. The model results were compared between the two scenarios to assess the effect of ditching on the groundwater system and potential changes to hydrologic conditions that support wetlands. Calibration of the regional SRLB model resulted in average horizontal hydraulic conductivity values of 6–39 feet per day for the glacial deposits and 3–4 feet per day for the uppermost fractured bedrock in the Biwabik Iron-Formation on the Mesabi Range. Average recharge across the calibrated model was 5.9 inches per year. Linesink resistance for the routed stream network was calibrated by using resistance categories based on the mapped soil hydrologic groups. The modeled regional groundwater-flow direction was generally to the south near the Mesabi Range topographic high and south or southwest across the rest of the basin. The updated calibration of the central SLRB model resulted in average horizontal hydraulic conductivity values of 5–36 feet per day for the glacial deposits and 3 feet per day for the uppermost fractured Biwabik Iron-Formation of the Mesabi Range. Average recharge across the ditch scenarios was 4.1 inches per year. Comparison of the preditch and postditch model scenarios showed that ditching reduced the area where the modeled water table was within 1 foot of the land surface (a wetland hydrology indicator) in as much as 40,000 acres, or 37 percent, of mapped permanent wetlands in the SLRB. An increase in the depth to the water table in wetland areas has the potential to degrade wetland persistence or function.

Minnesota↗
USGS2019

Hydrologic site assessment for passive treatment of groundwater nitrogen with permeable reactive barriers, Cape Cod, Massachusetts

Wastewater disposal associated with rapid population growth and development on Cape Cod, Massachusetts, during the past several decades has resulted in widespread contamination of groundwater with nitrogen. As a result, water quality in many of the streams, lakes, and coastal embayments on Cape Cod is impaired by excess nitrogen. To reduce nitrogen loads to these impaired water bodies, watershed-based planning is currently [2019] underway following a regional strategy, the section 208 areawide water-quality management plan update for Cape Cod. In the updated plan, traditional (sewering) and alternative wastewater management options are under consideration for restoring water quality in impaired surface-water bodies. Permeable reactive barriers, which are reactive zones emplaced below the water table for passive treatment of groundwater contaminants, are one of the alternatives being considered by Cape Cod towns as a potentially cost-effective technology for the removal of nitrogen from groundwater. However, the effectiveness of permeable reactive barriers depends on local conditions, and site-specific hydrologic and water-quality data are needed to inform the decision to install a permeable reactive barrier in a given location. These data are not available in most locations on Cape Cod; consequently, site assessments are needed before selecting this treatment option. To address this need, the U.S. Environmental Protection Agency, U.S. Geological Survey, and Cape Cod Commission formed a technical team in 2015 to develop and evaluate a hydrologic site-assessment approach for permeable reactive barrier installation. The approach developed by the technical team includes a preliminary regional assessment followed by a phased onsite investigation. The approach was intended to provide the hydrologic data needed to make informed decisions on site suitability and to support installation and monitoring should the site be deemed appropriate for a permeable reactive barrier. The factors that were evaluated to characterize local hydrologic conditions and inform site selection included groundwater flow directions and rates, depth to the water table, hydraulic conductivity and degree of heterogeneity of the aquifer, spatial distribution and concentration of nitrate and oxidation-reduction-sensitive constituents, thickness and depth of the treatment zone, distance to downgradient water bodies, and access for drilling and permeable reactive barrier installation. The approach was demonstrated on Cape Cod by conducting a preliminary assessment of 27 sites, from which 5 sites were selected for onsite investigations. Results indicated that the site-assessment approach was successful for screening sites and characterizing the geologic, hydrologic, and water-quality conditions at the sites selected for onsite investigations. Overall, the phased assessment evaluated in this study provided an efficient means of obtaining the hydrologic information needed to determine if a site was suitable for permeable reactive barrier installation on Cape Cod for the passive treatment of nitrogen in groundwater.

Massachusetts↗
USGS2020

Variable-density groundwater flow and contaminant transport, Operable Unit 1, Naval Base Kitsap, Keyport, Washington

Chlorinated volatile organic compounds (CVOCs) have migrated to groundwater beneath a former 9-acre landfill at Operable Unit 1 (OU-1) on Naval Base Kitsap, which was active from the 1930s through 1973 on the Keyport Peninsula, in Kitsap County, Washington. Biodegradation of CVOCs at OU-1 limits the mass of dissolved-phase CVOCs in groundwater that discharges to surface water, but contaminant concentrations up to 630 milligrams per liter persist in localized areas, likely from the dissolution of residual, non-aqueous phase liquids. Variable-density groundwater-flow and contaminant-transport models were developed using the SEAWAT-Version 4 computer program to simulate the direction and rate of groundwater flow in a 5.9 square-mile (mi 2 ) - area surrounding the Keyport Peninsula, to estimate the CVOC mass in groundwater and the rate of mass loading, and to assess possible remedial activities at OU-1. The study area is underlain by Quaternary deposits consisting of alternating glacial and interglacial sediments ranging from 500 to 1,500 feet (ft) thick. A hydrogeologic model delineated a sequence of 10 units including a relatively thin package (less than 100 ft) of recent sediments (Vashon Stade and younger) beneath the Keyport Peninsula that are underlain by the much thicker (more than 300 ft) Clover Park Aquitard, which overlies a confined, sea-level aquifer.

Washington↗
USGS2022

Geology and hydrogeology of the Yucaipa groundwater subbasin, San Bernardino and Riverside Counties, California

The Yucaipa groundwater subbasin (referred to in this report as the Yucaipa subbasin) is located about 75 miles (mi) east of of Los Angeles and about 12 mi southeast of the City of San Bernardino. In the Yucaipa subbasin, as in much of southern California, limited annual rainfall and large water demands can strain existing water supplies; therefore, understanding local surface water and groundwater conditions is essential for managing these resources. To better understand the hydrogeology and water resources in the Yucaipa subbasin, especially groundwater, the San Bernardino Valley Municipal Water District and the U.S. Geological Survey initiated a cooperative study to evaluate the hydrogeologic system of the Yucaipa subbasin and the encompassing Yucaipa Valley watershed. Previous studies of the area provided information on general geologic and hydrologic conditions, but this study provides the first comprehensive definition of the hydrogeology of the subsurface throughout the entire subbasin. The Yucaipa subbasin is located between the northwest trending San Andreas fault zone and San Jacinto fault. Several northeast-trending dip-slip faults dissect the Yucaipa subbasin, providing the mechanism for structural relief within the sediment-filled subbasin and between the subbasin and surrounding mountains and highlands. Several of these dip-slip faults have been previously identified as potential barriers to groundwater flow. This report provides a synthesis of previous studies and a discussion of the geologic interpretations that were used as the foundation for hydrogeologic classification of the Yucaipa subbasin. Notably, this report (1) adopts the recently named and classified sedimentary deposits of Live Oak Canyon geologic formation and extends the mapped distribution of the formation into the Yucaipa subbasin, and (2) adopts the interpretation that activity along the Banning fault predates the deposition of most basin-fill sedimentary materials in the Yucaipa subbasin. Four hydrogeologic units were classified in the Yucaipa subbasin: (1) crystalline basement, (2) consolidated sedimentary materials, (3) unconsolidated sediment, and (4) surficial materials. The crystalline basement unit forms the bottom boundary of the aquifer system, and the three other units comprise the basin-fill aquifer system. The four hydrogeologic units vary in extent, thickness, and structural relief across the subbasin, with the unconsolidated sediment unit serving as the primary aquifer unit. A three-dimensional hydrogeologic framework model was developed for the Yucaipa subbasin and surrounding area to characterize the thickness, extent, and hydrogeologic variability of the aquifer system. Geologic maps, borehole geophysical logs, drillers’ lithology logs, and depth-to-basement gravity data were used to map and interpolate the subsurface extent and structure of the hydrogeologic units within the subbasin. Faults and structures of geologic and (or) hydrogeologic importance were included in the model for future evaluation of their potential effects on groundwater flow. The resulting hydrogeologic framework is consistent with existing geologic concepts and the tectonic and structural history of the Yucaipa subbasin and surrounding area. The framework is also suitable for use in basin-scale hydrogeologic investigations.

California↗
USGS2022

Characterization of and temporal changes in groundwater quality of the Upper Black Squirrel Creek Basin, El Paso County, Colorado, 2018–20

In 2018–20, the U.S. Geological Survey, in cooperation with Upper Black Squirrel Creek Ground Water Management District, sampled 48 wells for Phase III of a multiphase plan investigating groundwater quality in the alluvial aquifer of the Upper Black Squirrel Creek Basin (UBSB), El Paso County, Colorado. Results for samples collected from October to December each year were used to assess spatial and temporal changes in groundwater quality and to differentiate sources of nitrate. Groundwater was predominantly classified as mixed-cation and mixed-anion water type in the aquifer, with variable chemistry along the periphery. Concentrations of constituents in groundwater were generally less than regulatory standards, except for nitrate in four wells. Isotopes of nitrogen and oxygen in nitrate identified four different potential sources or processes affecting nitrate in the alluvial aquifer: naturally occurring nitrate from soils, nitrate from animal and (or) human waste, and an unknown source, along with evidence of denitrification. Pharmaceutical compounds and personal-care products were detected in seven wells, with three wells having multiple detections. Stable isotopes of water indicated variability in seasonality of recharge throughout the UBSB alluvial aquifer. Nitrate concentrations from the 1984 study and the 1996 study were compared to the more recent concentrations in the 2013 study and the 2018–20 study. The northern one-third of the UBSB alluvial aquifer had a statistically significant increase in nitrate concentration from the 2013 study to the 2018–20 study, but no change was shown from the 1984 study to the 1996 study. The opposite was found true for the southern two-thirds of the UBSB alluvial aquifer with no statistically significant difference in nitrate concentration from the 2013 study to the 2018–20 study. Analysis of temporal changes indicated an increase in median and maximum nitrate concentrations from the 2013 study to the 2018–20 study throughout the UBSB alluvial aquifer. Continued sampling of wells in the UBSB would be beneficial to better determine temporal changes in groundwater quality, characterize human effects on water quality, and understand characteristics of the alluvial aquifer pertaining to sustainability of the resource.

Colorado↗
USGS2023

Preliminary machine learning models of manganese and 1,4-dioxane in groundwater on Long Island, New York

Manganese and 1,4-dioxane in groundwater underlying Long Island, New York, were modeled with machine learning methods to demonstrate the use of these methods for mapping contaminants in groundwater in the Long Island aquifer system. XGBoost, a gradient boosted, ensemble tree method, was applied to data from 910 wells for manganese and 553 wells for 1,4-dioxane. Explanatory variables included soil properties, groundwater flow, land use, and other features that describe the hydrogeology and geochemistry of the aquifer system. Four models were developed to predict the probability of manganese concentrations greater than a detection level of 10 micrograms per liter (μg/L) and greater than three threshold concentrations (50, 150, and 300 μg/L) relevant to drinking-water quality. One model was developed to predict the probability of 1,4-dioxane concentrations greater than a detection level of 0.07 μg/L. The 1,4-dioxane model was limited geographically to Suffolk County because of data availability. Predictions were made for two layers in the upper glacial aquifer and three layers in the Magothy aquifer, which are the upper two of the three major aquifers of the Long Island aquifer system. The objective of the study described in this report was to demonstrate the application of the methods rather than to develop precise estimates of manganese or 1,4-dioxane concentrations at any given location. The predictive models developed in the study are considered preliminary in the sense that they are an initial effort at developing these kinds of models specifically for Long Island. The models could be improved by the inclusion of additional data, by the use of methods to improve the modeling of infrequent high concentrations of manganese and 1,4-dioxane (above threshold concentrations), and by including more explanatory variables that specifically describe conditions and contaminant sources on Long Island. Nonetheless, the distribution of model predictions and the influence of explanatory variables in the models were consistent with the expected relations between contaminant concentrations and groundwater-flow-system characteristics and the distribution of manmade sources. Mapped predictions indicated that manganese detections were more probable in the upper glacial aquifer and along the southern shore of Long Island, consistent with the distribution of anoxic conditions in groundwater in the Long Island aquifer system. Manganese was infrequently predicted at concentrations greater than thresholds of concern for drinking-water quality in any of the aquifer layers. Detections of 1,4-dioxane were predicted in the western, more highly developed parts of Suffolk County, in the upper glacial aquifer and the top and middle layers of the Magothy aquifer, and in northwestern Suffolk County in the bottom layer of the Magothy aquifer. Although preliminary in nature and based on limited data, these mapped predictions can be used to generally identify areas where manganese and 1,4-dioxane may be present at concentrations of concern to prioritize areas for future monitoring and to guide future modeling and mapping efforts.

New York↗
USGS2023

Assessing the effects of chloride deicer applications on groundwater near the Siskiyou Pass, southwestern Oregon, July 2018–February 2021

The U.S. Geological Survey, in cooperation with the Oregon Department of Transportation (ODOT), evaluated the effects of cold-weather chloride deicers (road deicing chemicals) on groundwater quality, with a focus on chloride, near the Siskiyou Pass in southwestern Oregon. The study covered the period during July 2018 through February 2021. Between the years 2016 and 2020 ODOT applied up to 16,000 gallons per mile of chloride deicer and 143,000 pounds per mile of road salt along an 11-mile stretch of Interstate 5 (I-5) through the Siskiyou Pass. Despite the benefit of safer driving conditions, there are potentially negative environmental effects associated with the use of chloride-based deicers (such as magnesium chloride and sodium chloride). The results from this study are intended to help ODOT assess the water-quality effects from the application of chloride deicers at the Siskiyou Pass and inform decisions on how those chemicals are used. Dissolved chloride concentrations tended to be greater in groundwater downgradient from I-5 compared to groundwater upgradient from the interstate. Specific conductance was a good predictor of dissolved chloride concentration (R 2 = 0.905). Continuous monitoring showed that specific conductance measurements were greater at four downgradient spring-fed sites at the end of the study period compared with measurements at the beginning of the study. The study results indicate that chloride levels in shallow groundwater downgradient from I-5 are increasing, but dissolved chloride concentrations in domestic wells are not above the U.S. Environmental Protection Agency drinking water recommendations. The approach and methods used in this study, with modifications as site conditions warrant, can be applied in other areas of chloride deicer application to determine if groundwater is affected.

Oregon↗