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,117 records · Page 62Linked to original sources

Status of groundwater quality in the Santa Barbara Study Unit, 2011: California GAMA Priority Basin Project

Groundwater quality in the 48-square-mile Santa Barbara study unit was investigated in 2011 as part of the California State Water Resources Control Board’s Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project. The study unit is mostly in Santa Barbara County and is in the Transverse and Selected Peninsular Ranges hydrogeologic province. The GAMA Priority Basin Project is carried out by the U.S. Geological Survey in collaboration with the California State Water Resources Control Board and Lawrence Livermore National Laboratory. The GAMA Priority Basin Project was designed to provide a statistically unbiased, spatially distributed assessment of the quality of untreated groundwater in the primary aquifer system of California. The primary aquifer system is defined as that part of the aquifer corresponding to the perforation interval of wells listed in the California Department of Public Health database for the Santa Barbara study unit. This status assessment is intended to characterize the quality of groundwater resources in the primary aquifer system of the Santa Barbara study unit, not the treated drinking water delivered to consumers by water purveyors. The status assessment for the Santa Barbara study unit was based on water-quality and ancillary data collected in 2011 by the U.S. Geological Survey from 23 sites and on water-quality data from the California Department of Public Health database for January 24, 2008–January 23, 2011. The data used for the assessment included volatile organic compounds; pesticides; pharmaceutical compounds; two constituents of special interest, perchlorate and N -nitrosodimethylamine (NDMA); and naturally present inorganic constituents, such as major ions and trace elements. Relative-concentrations (sample concentration divided by the health- or aesthetic-based benchmark concentration) were used to evaluate groundwater quality for those constituents that have federal or California regulatory and non-regulatory benchmarks for drinking-water quality. For inorganic, organic, and special-interest constituents, a relative-concentration greater than 1.0 indicates a concentration greater than the benchmark and is classified as high. Inorganic constituents are classified as moderate if relative-concentrations are greater than 0.5 and less than or equal to 1.0 and are classified as low if relative-concentrations are less than or equal to 0.5. For organic and special-interest constituents, the boundary between moderate and low relative-concentrations was set at 0.1. Aquifer-scale proportion was used as the primary metric for evaluating regional-scale groundwater quality. High aquifer-scale proportion is defined as the areal percentage of the primary aquifer system with a relative-concentration greater than 1.0 for a particular constituent or class of constituents. Moderate and low aquifer-scale proportions were defined as the areal percentage of the primary aquifer system that had moderate and low relative-concentrations, respectively. Two statistical approaches—grid based and spatially weighted—were used to calculate aquifer-scale proportions for individual constituents and constituent classes. Grid-based and spatially weighted estimates were comparable in this the study (within 90-percent confidence intervals). Grid-based results were selected for use in the status assessment unless, as was observed in a few cases, a grid-based result was zero and the spatially weighted result was not zero, in which case, the spatially weighted result was used. Inorganic constituents that have human-health benchmarks were present at high relative-concentrations in 5.3 percent of the primary aquifer system and at moderate concentrations in 32 percent. High aquifer-scale proportions of inorganic constituents primarily were a result of high aquifer-scale proportions of boron (5.3 percent) and fluoride (5.3 percent). Inorganic constituents that have aesthetic-based benchmarks, referred to as secondary maximum contaminant levels, were present at high relative-concentrations in 58 percent of the primary aquifer system and at moderate concentrations in 37 percent. Iron, manganese, sulfate, and total dissolved solids were the inorganic constituents with secondary maximum contaminant levels present at high relative-concentrations. In contrast, organic and special-interest constituents that have health-based benchmarks were not detected at high relative-concentrations in the primary aquifer system. Of the 218 organic constituents analyzed, 10 were detected—9 that had human-health benchmarks. Organic constituents were present at moderate relative-concentrations in 11 percent of the primary aquifer system. The moderate aquifer-scale proportions were a result of moderate relative-concentrations of the volatile organic compounds methyl tert- butyl ether (MTBE, 11 percent) and 1,2-dichloroethane (5.6 percent). The volatile organic compounds 1,1,1-trichloroethane, 1,1-dichloroethane, bromodichloromethane, chloroform, MTBE, and perchloroethene (PCE); the pesticide simazine; and the special-interest constituent perchlorate were detected at more than 10 percent of the sites in the Santa Barbara study unit. Perchlorate was present at moderate relative-concentrations in 50 percent of the primary aquifer system. Pharmaceutical compounds and NDMA were not detected in the Santa Barbara study unit.

California↗

Delineation of areas contributing groundwater to selected receiving surface water bodies for long-term average hydrologic conditions from 1968 to 1983 for Long Island, New York

To assist resource managers and planners in developing informed strategies to address nitrogen loading to coastal water bodies of Long Island, New York, the U.S. Geological Survey and the New York State Department of Environmental Conservation initiated a program to delineate a comprehensive dataset of groundwater recharge areas (or areas contributing groundwater), travel times, and outflows to streams and saline embayments on Long Island. A four-layer regional three-dimensional finite-difference groundwater-flow model of hydrologic conditions from 1968 to 1983 was used to provide delineations of 48 groundwater watersheds on Long Island. Sixteen particle starting points were evenly spaced within each of the 4,000- by 4,000-foot model cells that receive water-table recharge and tracked using forward particle-tracking analysis modeling software to outflow zones. For each particle, simulated travel times were grouped by age as follows: less than or equal to 10 years, greater than 10 years and less than or equal to 100 years, greater than 100 years and less than or equal to 1,000 years, and greater than 1,000 years; and simulated ending zones were grouped into 48 receiving water bodies, based on the New York State Department of Environmental Conservation Waterbody Inventory/Priority Waterbodies List. Areal delineation of travel time zones and groundwater contributing areas were generated and a table was prepared presenting the sum of groundwater outflow for each area.

New York↗

Groundwater resources of the Devils Postpile National Monument—Current conditions and future vulnerabilities

This study presents an extensive database on groundwater conditions in and around Devils Postpile National Monument. The database contains chemical analyses of springs and the monument water-supply well, including major-ion chemistry, trace element chemistry, and the first information on a list of organic compounds known as emerging contaminants. Diurnal, seasonal, and annual variations in groundwater discharge and chemistry are evaluated from data collected at five main monitoring sites, where streams carry the aggregate flow from entire groups of springs. These springs drain the Mammoth Mountain area and, during the fall months, contribute a significant fraction of the San Joaquin River flow within the monument. The period of this study, from fall 2012 to fall 2015, includes some of the driest years on record, though the seasonal variability observed in 2013 might have been near normal. The spring-fed streams generally flowed at rates well below those observed during a sequence of wet years in the late 1990s. However, persistence of flow and reasonably stable water chemistry through the recent dry years are indicative of a sizeable groundwater system that should provide a reliable resource during similar droughts in the future. Only a few emerging contaminants were detected at trace levels below 1 microgram per liter (μg/L), suggesting that local human visitation is not degrading groundwater quality. No indication of salt from the ski area on the north side of Mammoth Mountain could be found in any of the groundwaters. Chemical data instead show that natural mineral water, such as that discharged from local soda springs, is the main source of anomalous chloride in the monument supply well and in the San Joaquin River. The results of the study are used to develop a set of recommendations for future monitoring to enable detection of deleterious impacts to groundwater quality and quantity

California↗

Estimates of nitrate loads and yields from groundwater to streams in the Chesapeake Bay watershed based on land use and geology

The water quality of the Chesapeake Bay may be adversely affected by dissolved nitrate carried in groundwater discharge to streams. To estimate the concentrations, loads, and yields of nitrate from groundwater to streams for the Chesapeake Bay watershed, a regression model was developed based on measured nitrate concentrations from 156 small streams with watersheds less than 500 square miles (mi2 ) at baseflow. The regression model has three predictive variables: geologic unit, percent developed land, and percent agricultural land. Comparisons of estimated and actual values within geologic units were closely matched. The coefficient of determination (R2 ) for the model was 0.6906. The model was used to calculate baseflow nitrate concentrations at over 83,000 National Hydrography Dataset Plus Version 2 catchments and aggregated to 1,966 total 12-digit hydrologic units in the Chesapeake Bay watershed. The modeled output geospatial data layers provided estimated annual loads and yields of nitrate from groundwater into streams. The spatial distribution of annual nitrate yields from groundwater estimated by this method was compared to the total watershed yields of all sources estimated from a Chesapeake Bay SPAtially Referenced Regressions On Watershed attributes (SPARROW) water-quality model. The comparison showed similar spatial patterns. The regression model for groundwater contribution had similar but lower yields, suggesting that groundwater is an important source of nitrogen for streams in the Chesapeake Bay watershed.

Chesapeake Bay watershed↗

Hydrologic conditions and simulation of groundwater and surface water in the Great Dismal Swamp of Virginia and North Carolina

The U.S. Geological Survey (USGS), in cooperation with the U.S Fish and Wildlife Service, has investigated the hydrology of the Great Dismal Swamp (Swamp) National Wildlife Refuge (Refuge) in Virginia and North Carolina and developed a three-dimensional numerical model to simulate groundwater and surface-water hydrology. The model was developed with MODFLOW-NWT, a USGS numerical groundwater flow modeling program, in combination with the Surface-Water Routing Process, a software package that simulates dynamic surface-water flows, water control structure management, and groundwater/surface-water interactions. The steady-state model was calibrated to average spring conditions by using automated parameter estimation software (PEST) to reduce simulation errors and assess model parameter sensitivity. The model was then used to simulate wet and dry climatic conditions and a variety of hypothetical scenarios in which water levels in the Swamp were raised and lowered by simulated management of water control structures. Results of the model simulations indicate that, under average spring conditions, precipitation is the primary water input (92%); surface-water (5%) and groundwater (3%) inflows make up the remainder. The primary outflow (or loss) is evapotranspiration (55%), with surface outflows (about 41%) and groundwater outflow (about 4%) making up the remainder. Simulated adjustment of water control structure weir levels demonstrates that groundwater levels are affected by water levels in adjacent ditches and that surface-water and groundwater levels can be controlled through management of water control structures, allowing the Refuge to better manage fire risks and preserve forested-wetland ecosystems in the Refuge. The 13 water control structures proposed in the simulated scenario representing possible future conditions effectively raised simulated water levels in the northeastern corner of the study area, a goal of the Refuge management. Results of this study demonstrate use of MODFLOW with the Surface-Water Routing Process for simulating water management options in peat wetlands and will help Refuge managers to better understand existing hydrologic conditions, assess the hydrologic effects of planned changes to water control structures, and apply the new simulation tool to guide water management on the Refuge.

North Carolina↗

Simulation of potential groundwater recharge for the glacial aquifer system east of the Rocky Mountains, 1980–2011, using the Soil-Water-Balance Model

An understanding of the spatial and temporal extent of groundwater recharge is critical for many types of hydrologic assessments involving water quality, contaminant transport, ecosystem health, and sustainable use of groundwater. Annual potential groundwater recharge was simulated at a 1-kilometer resolution with the Soil-Water-Balance (SWB) model for the glacial aquifer system east of the Rocky Mountains, from central Montana east to Maine, for calendar years 1980–2011. The SWB model used high resolution meteorological, land cover, and soil hydrology datasets that are nationally consistent and publicly available. The SWB model computed daily potential groundwater recharge as precipitation in excess of interception, runoff, evapotranspiration, and soil-water storage capacity. Daily potential recharge values within each year of the simulation were summed to produce annual potential recharge rates. Potential recharge as described in this report is water that infiltrates vertically below the plant rooting zone and is assumed to reach the water table. The calibrated SWB model in this report is called the glacial SWB model. Model calibration assumed that the area contributing to groundwater discharge equaled the surface watershed. The model was calibrated to stream base flows from 39 watersheds throughout the model domain that had hydrologic conditions appropriate for hydrograph separation. Base flows were calculated from daily streamflow records with the HYSEP local minimum hydrograph separation method The glacial SWB model reproduced the mean annual base-flow calibration targets well; the Nash-Sutcliffe efficiency coefficient was 0.94, and the root mean squared error was 1.28 inches per year. The glacial SWB model provides insight into the spatial and temporal variability in potential annual recharge across the glacial aquifer system. About 20 percent of the active model area had an average potential recharge rate of less than 1 inch per year. Total precipitation, total recharge, and recharge as a percentage of precipitation increased from west to east. A substantial amount of the recharge water (39 percent) entering the glacial aquifer system travels through developed (urbanized) and agricultural landscapes, which are known to cause water-quality impairments. Regional climatic events, such as the 1988 to 1989 drought, are apparent in the potential recharge time series. Potential recharge generally increased across the glacial aquifer system between 2001 and 2011. A comparison of the potential recharge from the glacial SWB model to previous broad-scale recharge estimates reveals several important considerations for future SWB modeling applications. Shifts in the overall distribution of potential recharge between separate models can be explained by methods used to generate base-flow calibration target datasets. Spatial patterns in potential recharge simulated by SWB models are strongly dependent on the data and assumptions used to assign model cells to hydrologic soil groups. A review of several SWB models used to estimate groundwater recharge (and not surface runoff) revealed that model results are most sensitive to input climatic data, followed by surface runoff (curve number) and root-zone depth parameters.

Scientific Investigations Report↗

Preliminary investigation of groundwater quality near a Michigan cemetery, 2016–17

The potential effect of cemetery leachate on groundwater quality in the United States has rarely been studied. Nutrients and other constituents associated with decomposition and burial processes (such as embalming) have the potential to reach shallow groundwater and could affect nearby drinking-water sources. The objective of this preliminary investigation was to evaluate the potential effect of cemetery leachate on shallow groundwater quality near Mt. Hope Cemetery in Ingham County, Lansing, Michigan, which is within the Well-head Protection Area for the City of Lansing. The constituents measured in this study include nutrients, trace metals, formaldehyde, fecal indicator bacteria, bacterial pathogen genes, contaminants of emerging concern (including pharmaceuticals, personal care products, and wastewater indicator compounds), and age-dating compounds. Three monitoring wells were installed 7 to 12 feet below land surface downgradient from the cemetery and sampled quarterly for 1 year. A fourth well (Fenner) was sampled to determine groundwater conditions outside the potential effects of cemetery leachate; samples from this well were collected near the water table. Nitrogen and phosphorus compounds were present at higher concentrations in two of the three monitoring wells (wells C1 and C3) than in the Fenner well. Formaldehyde and pharmaceuticals were not detected in any of the wells; however, several trace metals, including arsenic, manganese, and aluminum, were present in high concentrations, with arsenic concentrations typically exceeding the U.S. Environmental Protection Agency (EPA) drinking-water standard. Several wastewater indicator compounds, including atrazine, phenol, p-cresol, camphor, and skatole, were detected in the monitoring wells. Microbial data indicate the presence of staphylococci, enterococci, and Escherichia coli (E. coli) , with the highest concentrations being measured in the same two monitoring wells that exhibited elevated concentrations of nutrients in the groundwater (wells C1 and C3). Several bacterial pathogen genes were detected, including several Enterococcus species (spp.)— vanB (vancomycin-resistant enterococci), shiga-toxin-producing E. coli genes (including eaeA [attachment virulence trait] and stx1 [moderate toxin]), and the E. coli 16s ribosomal RNA (rDNA) gene ( E. coli species marker). These results were similar to results of studies conducted in Canada, Australia, and the United Kingdom, in which concentrations of bacteria, metals, and nutrients were elevated in groundwater near cemeteries.

Michigan↗

Tritium as an indicator of modern, mixed, and premodern groundwater age

Categorical classification of groundwater age is often used for the assessment and understanding of groundwater resources. This report presents a tritium-based age classification system for the conterminous United States based on tritium ( 3 H) thresholds that vary in space and time: modern (recharged in 1953 or later), if the measured value is larger than an upper threshold; premodern (recharged prior to 1953) if the measured value is smaller than a lower threshold; or mixed if the measured value is between the two thresholds. Inclusion of spatially varying thresholds, rather than a single threshold, accounts for the observed systematic variation in 3 H deposition across the United States. Inclusion of time-varying thresholds, rather than a single threshold, accounts for the date of sampling given the radioactive decay of 3 H. The efficacy of the tritium-based age classification system was evaluated at national and regional scales. The system was evaluated at a national scale by classifying samples from 1,788 public-supply wells distributed across 19 principal aquifers and comparing those results with expectations based on hydrogeologic principles. The regional-scale data are from five paired networks of shallow and deep wells (287 wells). As expected, modern groundwater is more prevalent in shallow wells than in deeper wells, in fractured-rock and carbonate aquifers as compared to clastic aquifers, in unconfined areas as compared to confined areas, and in humid climates as compared to arid climates. The results from a tritium-based age classification system compared favorably with the results of 14 previous studies of groundwater ages that used different age tracers and analytical methods. The wells and samples from the Cambrian-Ordovician aquifer that had been analyzed using a more complex multi-tracer analysis were also analyzed using the tritium-based age classification system, and there was a close match between the two methods. The results from these various studies suggest that the tritium-based age classification system may be informative as a screening tool prior to selecting more expensive and complex age-dating tracers and methods, or to provide an explanatory variable for other water-quality data where more complex methods or tracers are not available. This work improves on previous groundwater age classification using 3 H by developing methods that (1) determine 3 H thresholds for groundwater recharged in 1953 or later that minimize the misclassification of modern samples as mixed; (2) determine a pre-1953 threshold to estimate premodern background concentrations; and (3) add a mixed category to classify samples that are clearly neither entirely modern nor entirely premodern. As with any tritium-based approach, it can fail when the 3 H record in precipitation does not accurately reflect the record of 3 H in recharge

Conterminous United States↗

Trace metal and nutrient loads from groundwater seepage into the South Fork Coeur d’Alene River near Smelterville, northern Idaho, 2017

The Coeur d’Alene mining district in northern Idaho historically was a globally important source of lead, zinc, and silver, but over 100 years of mining has left a legacy of metals contamination in the Coeur d’Alene River valley. Previous studies by the U.S. Geological Survey (USGS) and others have indicated that groundwater discharging into the South Fork Coeur d’Alene River between Kellogg and Smelterville, Idaho, is a substantial source of dissolved zinc, dissolved cadmium, and total phosphorus. As part of its ongoing cleanup efforts, the U.S. Environmental Protection Agency is constructing a groundwater collection and treatment system to intercept and treat this contaminated water before it reaches the river. To establish conditions prior to construction, the USGS conducted a seepage study in September 2017 to quantify the rate and quality of groundwater discharging into the South Fork Coeur d’Alene River between Kellogg and Smelterville. Repeated measurements of streamflow were taken at multiple locations in the river and tributaries, and water-quality samples were collected and analyzed for trace metals and nutrients. Results showed consistent increases in streamflow (5.8 ± 1.3 cubic feet per second); and in dissolved zinc (85 ± 9.3 kilograms per day [kg/d]), dissolved cadmium (0.58 ± 0.10 kg/d) and total phosphorus (6.3 ± 0.45 kg/d) loads in a discrete segment of the reach. These gains exceeded tributary inputs, thereby implicating groundwater discharge as the main source of loading. Zinc and cadmium loads from groundwater in 2017 were less than those measured in 1999 but comparable to those measured from 2003 to 2008. This suggests that remedial actions in the late 1990s and early 2000s decreased trace-metal loading from 1999 to 2003, but that conditions remained similar from 2003 to 2017. A second seepage study will be conducted after construction and treatment plant system optimization are complete; this second study will evaluate changes in groundwater discharge to and water quality in the South Fork Coeur d’Alene River compared to the pre-construction conditions.

Idaho↗

Optimization assessment of a groundwater-level observation network in the Middle Rio Grande Basin, New Mexico

The U.S. Geological Survey, in cooperation with the Albuquerque Bernalillo County Water Utility Authority (ABCWUA), measures groundwater levels continuously (hourly) and discretely (semiannually and annually) at a network of wells and piezometers (hereafter called the observation network) within the Middle Rio Grande Basin in central New Mexico. Groundwater levels that are measured in this observation network provide a long-term hydrologic dataset that is heavily relied upon to make water management decisions. The desire to upgrade and perform maintenance on this observation network initiated this study, which assesses the spatial and temporal importance of measurements towards optimization of monitoring the observation network to reduce or redirect monitoring costs. This report describes the methods and results of the optimization assessment of this observation network, which included separate spatial and temporal methodologies and an evaluation using principal component analysis (PCA). Results from the spatial optimization assessment can be used to help identify observation network sites that do not significantly affect the generation of winter groundwater-elevation contour maps of the production zone. Results from the temporal optimization assessment and PCA can also be consulted when deciding which sites to remove from the network, especially for sites that are monitored more frequently than annually. Results from the temporal optimization assessment can be used to inform the minimum monitoring frequency at the observation network required to capture the trends shown in higher frequency monitoring. The PCA results distinguish spatially distributed characteristic water-level trends that can inform the management decisions that are made when using the spatial and temporal optimization assessment results. Reducing the temporal frequency or spatial density of monitoring is ultimately a management decision that depends on the amount of data loss or degradation that is deemed acceptable while still meeting the network objectives of the ABCWUA. This study can also serve as a starting point to a data gap analysis of local aquifer characteristics and help guide enhanced observation network design as needs arise or in advance of future water management decisions. The results of the spatial optimization assessment indicate that as many as about 20 specified sites can be removed from the observation network with a relatively small loss in the ability to represent the kriged groundwater-elevation surfaces of the production zone that were generated by using median groundwater elevations for two periods: the 2001 time interval and 2015 time interval. This analysis also demonstrated the importance of wells at the margin of the study area and in areas where there are large hydrologic gradients. At many of the 47 hourly monitored sites analyzed in the temporal optimization assessment, temporal trends were well represented for at least one of the reduced monitoring frequencies tested, indicating that a reduced frequency may be sufficient to adequately characterize seasonal and long-term trends. PCA and k-means clustering analysis of the 15 hourly monitored sites that are screened within the production zone indicate that the sites can be categorized into four groups, or clusters, of differing groundwater-level hydrograph characteristics. Except for one cluster, all of the clusters have the potential to be well represented by fewer index monitoring sites.

New Mexico↗

Geophysical surveys, hydrogeologic characterization, and groundwater flow model for the Truxton basin and Hualapai Plateau, northwestern Arizona

This is the third and final report in a series that describe the groundwater resources of the Hualapai Indian Reservation. These reports document the findings of a comprehensive groundwater study conducted on the reservation and adjacent areas from 2015 through 2018 by the U.S. Geological Survey in cooperation with the Bureau of Reclamation. The first report described the hydrologic framework and characterization of the Truxton aquifer on the Hualapai Indian Reservation (Bills and Macy, 2016). The second report described the hydrogeologic characterization of the Hualapai Plateau part of the reservation (Mason, Macy, and others, 2020). This report includes five chapters. Chapter A (Mason, Knight, and others, 2020) is a summary of this multichapter volume and briefly describes the study area. Chapter B (Mason, Bills, and Macy, 2020) describes the geology and hydrology of the Truxton basin and Hualapai Plateau. Chapter C (Kennedy, 2020) describes the results of a gravity geophysical survey of the Truxton basin. Chapter D (Ball, 2020) describes the findings of an airborne electromagnetic survey of the Truxton aquifer and Hualapai Plateau. Chapter E (Knight, 2020) describes the results of a transient groundwater model created for the entire Truxton aquifer both on and off the reservation. The groundwater-flow model is used to estimate projected groundwater levels based on future groundwater withdrawal scenarios.

Arizona↗

Simulated effects of sea-level rise on the shallow, fresh groundwater system of Assateague Island, Maryland and Virginia

The U.S. Geological Survey, in cooperation with the National Park Service, developed a three-dimensional groundwater-flow model for Assateague Island in eastern Maryland and Virginia to assess the effects of sea-level rise on the groundwater system. Sea-level rise is expected to increase the altitude of the water table in barrier island aquifer systems, possibly leading to adverse effects to ecosystems on the barrier islands. The potential effects of sea-level rise were evaluated by simulating groundwater conditions under sea-level-rise scenarios of 20 centimeters (cm), 40 cm, and 60 cm. Results show that as sea level rises, low-lying areas of the island originally represented as receiving freshwater recharge in the baseline scenario are inundated by saltwater. This change from freshwater recharge to saltwater decreases the overall amount of freshwater recharging the system. As the water table rises in response to the higher sea levels, freshwater flow out of the system changes, with more freshwater leaving as submarine groundwater discharge and less freshwater leaving as seeps and evapotranspiration. At the current land-surface altitude, as much as 50 percent of the island may be inundated with a 60-cm rise in sea level, and the low-lying marshes may change from freshwater to saltwater. Groundwater levels at 32 wells were monitored for as long as 12 months between October 2014 and September 2015 on Assateague Island. Results from objective classification analysis of 14 shallow monitoring wells show two dominant processes affecting groundwater levels in two different settings on the island. On the western side of the island, between the primary dune and the inland bays, water levels clearly respond to precipitation events. This side of the island is more protected from ocean tides and typically is more vegetated than the eastern side. On the eastern side of the island, between the Atlantic Ocean and the primary dune, water levels clearly respond to tidal events. Specific conductance was measured at four wells, two on the western part of the island and two on the eastern part of the island. Specific conductance values in the two wells west of the primary dune show episodic decreases, coinciding with precipitation events. Specific conductance values in the two wells on the eastern side of the primary dune show episodic increases, coinciding with high-tide events. These high frequency monitoring data are intended to aid in designing a monitoring network that can document both short-term and long-term hydrologic processes on Assateague Island National Seashore. This study uses a modeling approach consistent with models developed for Gateway National Recreation Area, Sandy Hook Unit (New Jersey) and Fire Island National Seashore (New York). Combined, these models are meant to improve the regional capabilities for predicting climate-change effects on barrier islands and provide resource managers with a common set of tools for adaptation and mitigation of potentially adverse effects of sea-level rise.

Maryland, Virginia↗

Surface-water and groundwater interactions in the Big Lost River, south-central Idaho

The Big Lost River of south-central Idaho interacts with the underlying aquifer by gaining and losing streamflow throughout various areas in the Big Lost River Valley. Surface-water and groundwater resources are used throughout the valley to sustain domestic, agricultural, and livestock needs. The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources, evaluated streamflow gains and losses by differential streamgaging in the lower Big Lost River, Idaho, during four measurement events: March 27–28, 2019; October 16–17, 2019; October 6–7, 2020; and March 30, 2021. This report presents and analyzes streamflow measurement and uncertainty data from each measurement event to describe surface-water/groundwater interactions. This report is the second chapter of a multi-chapter volume that characterizes water resources in the Big Lost River Basin. During the four measurement events, 100 streamflow measurements were made at 46 unique sites on the Big Lost River, James Creek, and diversions or tributaries between Mackay Reservoir near Mackay and Arco, Idaho. Aquifer lithology and dimensions affected spatial patterns of streamflow gains and losses between the upper, middle, and lower reaches; changes in water supply, groundwater levels, and surface-water management affected seasonal differences within reaches. In the upper reach of the Big Lost River, streamflow losses and gains were greater during the wetter 2019 events and lesser during the drier 2020 and 2021 events. The middle reach includes the largest losses from the Big Lost River to groundwater; these losses occurred in the Darlington Sinks where 42 percent or more of streamflow was lost as the aquifer widens and groundwater deepens. These results suggest that changing surface-water supply, irrigation use, and recharge affect interannual groundwater levels and, in turn, affect patterns of streamflow gains and losses in the middle reach. Finally, surface-water management is the primary control on surface-water/groundwater interactions in the lower reach. Overall patterns of streamflow gains and losses in this study generally were consistent with previous reports. However, paired with the related hydrogeologic framework and water budget, this investigation provides new insights into how hydrogeologic conditions and interannual variability in water supply, groundwater levels, and surface-water management affect surface-water/groundwater interactions in the Big Lost River Valley.

Idaho↗

Hydrogeologic characterization of the Yucaipa groundwater subbasin

Introduction Water management in the Santa Ana River watershed in San Bernardino and Riverside Counties in southern California (fig. A1) is complex with various water purveyors navigating geographic, geologic, hydrologic, and political challenges to provide a reliable water supply to stakeholders. As the population has increased throughout southern California, so has the demand for water. The Yucaipa groundwater subbasin (hereafter referred to as “Yucaipa subbasin”), one of nine groundwater subbasins in what the California Department of Water Resources (DWR) refers to as the Upper Santa Ana Valley groundwater basin (California Department of Water Resources, 2016; fig. A1; the DWR naming convention is used within this report), is no exception; steady population growth since the 1940s and changes in water use have forced local water purveyors to regularly adapt their water infrastructure. Water demands within the Yucaipa subbasin have historically been supplied by groundwater, but water imported via the California State Water Project has augmented the total water supply through direct use and through anthropogenic recharge at the Wilson Creek and Oak Glen Creek spreading basins since 2002. Overall demand for groundwater continues to rise, and local water managers are concerned that despite the influx of imported water, groundwater levels may decline to a point where producing water will be uneconomical, severely limiting the ability of local agencies to meet water-supply demand. To better understand the hydrogeology and water resources in the Yucaipa subbasin, the U.S. Geological Survey (USGS) initiated a study in cooperation with the San Bernardino Valley Municipal Water District (SBVMWD) to characterize and model the hydrologic system of the Yucaipa subbasin and the surrounding Yucaipa Valley watershed (YVW; fig. A2). To gain this comprehensive understanding, a three-dimensional (3D) hydrogeologic framework model (HFM; Cromwell and Matti, 2022) was constructed to quantify the structure and extent of hydrogeologic units in the YVW; the hydrologic system was conceptualized and quantified (described in chapter A); and the Yucaipa Integrated Hydrological Model (YIHM; described in chapter B ) was developed to simulate the integrated surface-water and aquifer systems, including natural and anthropogenic recharge and discharge throughout the study area during 1947–2014.

Scientific Investigations Report↗

Trends in groundwater levels, and orthophosphate and nitrate concentrations in the Middle Snake River Region, south-central Idaho

The U.S. Geological Survey (USGS) evaluated nitrate and orthophosphate concentrations in groundwater for temporal trends (monotonic and step trends) for the middle Snake River region (Cassia, Gooding, Jerome, Lincoln, Minidoka, and Twin Falls Counties) in south-central Idaho using the Regional Kendall test (monotonic trends) and the Wilcoxon signed rank test (step trends). The study evaluated two trend periods: 2000–09 and 2010–19/20. The study area was divided into six hydrogeologic zones (HZs) that had similar geologic and hydrologic characteristics and that correlated with county boundaries where possible. Two well networks sampled by the USGS National Water Quality Program within the HZs were also evaluated. The northern Gooding County HZ had statistically significant increasing nitrate concentration trends for both the monotonic and step trends in the early trend period, while the Cassia and Jerome/Southern Gooding County HZs only had one of the statistical tests with statistically significant increasing nitrate concentrations. The Minidoka County HZ had conflicting results between the two statistical tests for the early time period with a statistically significant increasing monotonic trend in nitrate concentration and a statistically significant decreasing step trend. The differing results between these two statistical tests indicates the significance of concentration data during the middle of the time period. Both the Lincoln and Twin Falls County HZs did not have statistically significant trends for either test during either time period as well as the Northern Gooding County HZ for the latter time period. The Minidoka County HZ had statistically significant nitrate trends for both tests in the latter time period along with one of the trend tests for the Cassia and Jerome/Southern Gooding County HZ. Most of the nitrate concentration trend rates are low from 0.01 to 0.12 milligram per liter per year (mg/L/year) with the northern Gooding County HZ having the highest trend rate during the early time period of 0.28 mg/L/year for the step trend and 0.55 mg/L/year for the monotonic trend. All the HZs and both well networks had statistically significant increasing orthophosphate-concentrations trends in groundwater for the early time period except for the Lincoln County HZ and the step-trend for the Minidoka County HZ. Orthophosphate concentration trend rates for the early period were low, ranging from 0.001 to 0.015 mg/L/year. Only two HZs and the well networks had enough orthophosphate concentration data available in the latter time period to do statistical analysis. The two HZs (Minidoka and Southern Gooding/Jerome County) both have decreasing orthophosphate concentration trends, with only the monotonic trend for the Southern Gooding/Jerome County HZ being statistically significant at 90 percent with a rate of −0.001 mg/L/year. Groundwater levels in two well networks in the eastern Snake River Plain aquifer were also evaluated for trends (monotonic and step), with both networks having statistically significant declining groundwater levels for the 1993–2009 trend period. The latter trend period (2010–20) had statistically significant declining groundwater levels for the A&B well network and statistically significant increasing groundwater levels for the Jerome/Gooding well network, which is downgradient from an aquifer recharge area.

Idaho↗

Examination of dissolved uranium concentrations in regional shallow groundwater relative to Operable Unit 8 of the Denver Radium Superfund Site

A radium industry existed between about 1914 and 1920 in Denver, Colorado, with operations located along the South Platte River. Sites associated with that industry were contaminated with radium and uranium processing residues and were incorporated into clean-up efforts as Operating Units (OUs) of the Denver Radium Superfund Site. Concentrations of uranium exceeding the U.S. Environmental Protection Agency maximum contaminant level of 0.03 milligrams per liter for drinking water are present in shallow groundwater at OU8. However, previous studies have shown concentrations of dissolved uranium can be naturally high in shallow groundwater of the South Platte River valley compared to other rivers of the world. This report compares dissolved uranium concentrations measured by the U.S. Geological Survey across the South Platte River valley to data collected at the OU8 of the Denver Radium Superfund Site. The U.S. Geological Survey data represent 5 distinct urban or agricultural geographic areas and included 230 sampling events at 114 wells during 1993 to 2013. The OU8 data represent 13 wells and groundwater discharge locations sampled during the years 2017 and 2018. Dissolved uranium concentrations were statistically significantly greater for both years of the OU8 data compared to three datasets from shallow groundwater beneath urban areas in the Denver metropolitan area. However, compared to OU8, concentrations were significantly greater in shallow groundwater from an agricultural area of the South Platte River valley distant from Denver. Additionally, each of the urban area datasets contained some individual dissolved uranium concentrations greater than the greatest concentrations from the two OU8 datasets. Thus, naturally occurring concentrations of dissolved uranium in shallow groundwater that are greater than those observed at OU8 are not uncommon in the South Platte River valley.

Colorado↗

Evaluation of alternative groundwater-withdrawal scenarios on water levels in Kingsbury Pond, upper Charles River Basin, eastern Massachusetts

Kingsbury Pond is a glacial kettle pond in the town of Norfolk, Massachusetts, in the Mill River Basin, which is part of the Upper Charles River Basin in eastern Massachusetts. The pond is hydraulically connected to the surrounding groundwater-flow system, and water levels in the pond fluctuate in response to recharge to the aquifer from precipitation and wastewater return flows through septic systems, to withdrawals from the aquifer at nearby wells, and to precipitation directly on the pond surface. Concerns about the effects of withdrawals on water levels in the pond prompted an investigation to better understand the hydrology of Kingsbury Pond and its response to groundwater withdrawals and to determine if withdrawals from wells in Franklin, Mass., can be modified to simultaneously reduce the effect on water levels in the pond and yet meet the water-supply demands of the Town of Franklin. An existing, transient groundwater-flow model of the Upper Charles River Basin was modified for this study in the area near Kingsbury Pond to improve representation of the hydrologic system near the pond. The mean annual water-level altitude simulated for the pond for nonpumping conditions using the modified model is 136 feet (ft), which falls within the range of likely annual pond-altitude fluctuations of 135 to 140 ft estimated for average hydrologic conditions before the beginning of withdrawals at two nearby wells operated by the Town of Franklin (wells FR–04 and FR–05). The mean annual water-level altitude at the pond decreased by 3.8 ft to 132.2 ft for simulated mean monthly withdrawal rates at all wells within the Upper Charles River Basin from 2010 to 2019 (referred to as the baseline withdrawal condition). A groundwater management model that links the groundwater-flow model with a mathematical optimization method was developed to evaluate the effects of three alternative groundwater-withdrawal scenarios for the Franklin public-water system on water levels in Kingsbury Pond. In the first scenario, monthly withdrawal rates at wells FR–04 and FR–05 were increased from the baseline withdrawal rates to their maximum authorized rates for all months of the year; all other Franklin wells were specified at their baseline withdrawal rates. This scenario resulted in a mean annual water-level altitude at the pond of 129.3 ft, or a mean annual decline of 6.7 ft compared with nonpumping conditions and a decline of 2.9 ft compared with baseline conditions. The results of the second scenario showed that water levels in the pond can be increased relative to 2010–19 conditions while meeting Franklin’s 2010–19 monthly water-supply demands if withdrawals at wells FR–04 and FR–05 were shifted to other Franklin wells. In this scenario, monthly withdrawal rates at wells FR–04 and FR–05 were decreased from their baseline rates to one-third their maximum practical rates for all months of the year; increased withdrawal rates at other Franklin wells were determined by the management model. The decrease in withdrawal rates at wells FR–04 and FR–05 resulted in a mean water-level altitude at the pond of 134.1 ft, which was equivalent to a 51 percent increase (improvement) in the mean annual water level of the pond relative to the baseline condition. A third scenario was done to determine if Franklin’s existing water-supply system has the capacity to meet the mean annual maximum permitted withdrawal rate of the system of 3.45 million gallons per day while maintaining monthly withdrawal rates at wells FR–04 and FR–05 at their 2010–19 rates and water levels in Kingsbury Pond at baseline conditions. The analysis indicated that the capacity of the system cannot meet the increased demand during some months of the year with withdrawal rates at the two wells fixed at their monthly 2010–19 rates; however, the existing system is capable of meeting about 90 percent of the maximum permitted rate (3.10 million gallons per day) by increasing withdrawal rates at other Franklin wells above their 2010–19 rates.

Massachusetts↗

Hydrogeologic conceptual model of groundwater occurrence and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River is a major concern in the southwestern United States where the river provides water to about 40 million people for municipal and industrial use and is used to irrigate about 5.5 million acres of land. Much of the salinity in the Colorado River Basin is derived from natural interactions of surface water and groundwater with various geologic materials (rocks, soils, and alluvial deposits). The Dolores River in southwest Colorado is a major tributary of the Colorado River that historically accounts for about 6 percent of the salinity load to the Upper Colorado River Basin with the Paradox Valley being the primary source of salinity to the Dolores River. The Paradox Valley, one of several salt-anticline valleys in the region, is a fault-bounded topographic basin aligned with and exposing an underlying salt-anticline core. Salt deposits in the Pennsylvanian Paradox Formation of the Hermosa Group form an elongated salt diapir oriented northwest to southeast that is up to 12,000 feet (ft) thick beneath the present valley floor. Surface erosion, groundwater circulation, and weathering during Tertiary and Quaternary valley formation contributed to development of a cap rock, collapse features, breccia, and brine at the top of the exposed salt diapir. Today (2023), brine occurring in the brecciated cap rock and underlying salt deposits is in hydraulic connection with an overlying freshwater alluvial aquifer, and depending on seasonal river stage and hydrologic conditions, the brine discharges to the Dolores River causing the observed increase in salinity as the river crosses the Paradox Valley. To reduce salinity concentrations in the Dolores River, the Bureau of Reclamation (Reclamation) operates the Paradox Valley Unit (PVU). The PVU project consists of nine shallow brine pumping wells near the Dolores River and one deep disposal well where the brine is injected for disposal. When operational, the PVU pumping wells extract brine from the base of the alluvial aquifer that is piped and injected into a deep disposal well about 3 miles southwest of the PVU. The PVU became fully operational July 1, 1996, and by 2015, operation of the PVU had reduced salinity concentrations in the Dolores River by as much as 70 percent compared to pre-PVU conditions. In response to a 4.5 magnitude earthquake, injection operations, and thus PVU pumping, were ceased from March 2019 to June 2022. A trial period of PVU operation began in June 2022 with a reduced injection rate, and thus PVU pumping rate, of about two-thirds capacity to gather additional information and guide future operational decisions. In cooperation with Reclamation, the U.S. Geological Survey (USGS) developed this report to present the current (2023) understanding of groundwater and brine occurrence and discharge to the Dolores River in the Paradox Valley. Results from the compilation of spatial datasets, groundwater sampling and age dating, and aquifer tests are presented to provide improved understanding of the Paradox Valley hydrogeology, to supply datasets for a numerical groundwater-flow and brine-transport model, and to support future operations of the PVU. The hydrogeologic data provided herein, along with the most recent loading analysis for the Dolores River in the Paradox Valley, and a previous conceptual model for brine discharge to the river are used to present a conceptual understanding of groundwater occurrence in the Paradox Valley.

Colorado↗