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 937 records · Page 52Linked to original sources
USGS2016

Simulation of groundwater storage changes in the eastern Pasco Basin, Washington

The Miocene Columbia River Basalt Group and younger sedimentary deposits of lacustrine, fluvial, eolian, and cataclysmic-flood origins compose the aquifer system of the Pasco Basin in eastern Washington. Irrigation return flow and canal leakage from the Columbia Basin Project have caused groundwater levels to rise substantially in some areas, contributing to landslides along the Columbia River. Water resource managers are considering extraction of additional stored groundwater to supply increasing demand and possibly mitigate problems caused by the increased water levels. To help address these concerns, the transient groundwater model of the Pasco Basin documented in this report was developed to quantify the changes in groundwater flow and storage. The MODFLOW model uses a 1-kilometer finite-difference grid and is constrained by logs and water levels from 846 wells in the study area. Eight model layers represent five sedimentary hydrogeologic units and underlying basalt formations. Head‑dependent flux boundaries represent the Columbia and Snake Rivers to the west and south, respectively, underflow to and (or) from adjacent areas to the northeast, and discharge to agricultural drains, springs, and groundwater withdrawal wells. Specified flux boundaries represent recharge from infiltrated precipitation and anthropogenic sources, including irrigation return flow and leakage from water-distribution canals. The model was calibrated with the parameter ‑ estimation code PEST++ to groundwater levels measured from 1907 through 2013 and measured discharge to springs and estimated discharge to agricultural drains. Increased recharge since pre-development resulted in a 6.8 million acre-feet increase in storage in the 508-14 administrative area of the Pasco Basin. Four groundwater-management scenarios simulate the 7-year drawdown resulting from withdrawals in different locations. Withdrawals of 2 million gallons per day (Mgal/d) from a hypothetical well field in the upper Ringold Formation along the Columbia River could generate 30–70 feet of drawdown, which may reduce landslide susceptibility along the White Bluffs. Drawdowns resulting from a 1 Mgal/d withdrawal from wells screened in either Pasco gravels, upper Ringold Formation, or both Ringold Formation and underlying basalt are simulated in the other three scenarios, and differ because of the contrasting hydraulic conductivities within the screened intervals.

Washington↗
USGS2016

Numerical simulation of the groundwater-flow system of the Kitsap Peninsula, west-central Washington

A groundwater-flow model was developed to improve understanding of water resources on the Kitsap Peninsula. The Kitsap Peninsula is in the Puget Sound lowland of west-central Washington, is bounded by Puget Sound on the east and by Hood Canal on the west, and covers an area of about 575 square miles. The peninsula encompasses all of Kitsap County, Mason County north of Hood Canal, and part of Pierce County west of Puget Sound. The peninsula is surrounded by saltwater, and the hydrologic setting is similar to that of an island. The study area is underlain by a thick sequence of unconsolidated glacial and interglacial deposits that overlie sedimentary and volcanic bedrock units that crop out in the central part of the study area. Twelve hydrogeologic units consisting of aquifers, confining units, and an underlying bedrock unit form the basis of the groundwater-flow model. Groundwater flow on the Kitsap Peninsula was simulated using the groundwater-flow model, MODFLOW‑NWT. The finite difference model grid comprises 536 rows, 362 columns, and 14 layers. Each model cell has a horizontal dimension of 500 by 500 feet, and the model contains a total of 1,227,772 active cells. Groundwater flow was simulated for transient conditions. Transient conditions were simulated for January 1985–December 2012 using annual stress periods for 1985–2004 and monthly stress periods for 2005–2012. During model calibration, variables were adjusted within probable ranges to minimize differences between measured and simulated groundwater levels and stream baseflows. As calibrated to transient conditions, the model has a standard deviation for heads and flows of 47.04 feet and 2.46 cubic feet per second, respectively. Simulated inflow to the model area for the 2005–2012 period from precipitation and secondary recharge was 585,323 acre-feet per year (acre-ft/yr) (93 percent of total simulated inflow ignoring changes in storage), and simulated inflow from stream and lake leakage was 43,905 acre-ft/yr (7 percent of total simulated inflow). Simulated outflow from the model primarily was through discharge to streams, lakes, springs, seeps, and Puget Sound (594,595 acre-ft/yr; 95 percent of total simulated outflow excluding changes in storage) and through withdrawals from wells (30,761 acre-ft/yr; 5 percent of total simulated outflow excluding changes in storage). Six scenarios were formulated with input from project stakeholders and were simulated using the calibrated model to provide representative examples of how the model could be used to evaluate the effects on water levels and stream baseflows of potential changes in groundwater withdrawals, in consumptive use, and in recharge. These included simulations of a steady-state system, no-pumping and return flows, 15-percent increase in current withdrawals in all wells, 80-percent decrease in outdoor water to simulate effects of conservation efforts, 15-percent decrease in recharge from precipitation to simulate a drought, and particle tracking to determine flow paths. Changes in water-level altitudes and baseflow amounts vary depending on the stress applied to the system in these various scenarios. Reducing recharge by 15 percent between 2005 and 2012 had the largest effect, with water-level altitudes declining throughout the model domain and baseflow amounts decreasing by as much as 18 percent compared to baseline conditions. Changes in pumping volumes had a smaller effect on the model. Removing all pumping and resulting return flows caused increased water-level altitudes in many areas and increased baseflow amounts of between 1 and 3 percent.

Washington↗
USGS2016

Hydrologic assessment of the shallow groundwater flow system beneath the Shinnecock Nation tribal lands, Suffolk County, New York

Defining the distribution and flow of shallow groundwater beneath the Shinnecock Nation tribal lands in Suffolk County, New York, is a crucial first step in identifying sources of potential contamination to the surficial aquifer and coastal ecosystems. The surficial or water table aquifer beneath the tribal lands is the primary source of potable water supply for at least 6 percent of the households on the tribal lands. Oyster fisheries and other marine ecosystems are critical to the livelihood of many residents living on the tribal lands, but are susceptible to contamination from groundwater entering the embayment from the surficial aquifer. Contamination of the surficial aquifer from flooding during intense coastal storms, nutrient loading from fertilizers, and septic effluent have been identified as potential sources of human and ecological health concerns on tribal lands. The U.S. Geological Survey (USGS) facilitated the installation of 17 water table wells on and adjacent to the tribal lands during March 2014. These wells were combined with other existing wells to create a 32-well water table monitoring network that was used to assess local hydrologic conditions. Survey-grade, global-navigation-satellite systems provided centimeter-level accuracy for positioning wellhead surveys. Water levels were measured by the USGS during May (spring) and November (fall) 2014 to evaluate seasonal effects on the water table. Water level measurements were made at high and low tide during May 2014 to identify potential effects on the water table caused by changes in tidal stage (tidal flux) in Shinnecock Bay. Water level contour maps indicate that the surficial aquifer is recharged by precipitation and upgradient groundwater flow that moves from the recharge zone located generally beneath Sunrise Highway, to the discharge zone beneath the tribal lands, and eventually discharges into the embayment, tidal creeks, and estuaries that bound the tribal lands to the east, south, and west. Water levels in many of the wells in the network fluctuated in response to precipitation, upgradient groundwater flow, and tidal flux in Shinnecock Bay. Water level altitudes ranged from 6.66 to 0.47 feet (ft) above the North American Vertical Datum of 1988 during the spring measurement period, and from 5.25 to -0.24 ft (NAVD 88) during fall 2014. Historically, annual and seasonal precipitation seem to indicate long-term water level trends in an index well located in the town of Southampton, correlates with changes in storage in the upper glacial aquifer, but does not necessarily indicate water level extremes in the shallow groundwater system. To place the study period in perspective, calendar year 2014 was the 32d wettest year on record, with precipitation for the year totaling 48.1 inches, a 2.6-percent increase from the annual average (46.9 inches per year), based on 81 years of complete record at the National Oceanographic and Atmospheric Administration, National Weather Service cooperative meteorological station at Bridgehampton, New York. Estimated recharge to the water table beneath the tribal lands from precipitation for 2014 is 25.4 inches. Tidal flux caused water levels in wells to fluctuate from 0.30 to -0.24 ft during May 2014. Water levels in wells located north of Old Fort Pond and beneath the southernmost extent of the tribal lands were most influenced by tidal flux. During June 2014, hydrographs indicate that tidal flux influenced water levels by 0.48 ft in a well located near the southernmost extent of the tribal lands approximately 0.3 miles north of Shinnecock Bay, and was zero at a well located approximately 0.5 miles south of Montauk Highway, and 0.4 miles west of Heady Creek, near the geographic center of the tribal lands. Tidal-influence delay time (time interval between peak high-tide stage and corresponding peak high-water level) ranged from 1.75 hours at the well located near the southernmost extent of the tribal lands, to more than 4 hours at a well located north of Old Fort Pond, near the northwestern part of the tribal lands. Estimated hydraulic-conductivity values derived from the results of specific-capacity tests that were completed at nine observation wells during March 2015 were used to calculate average linear velocity. Average linear velocity along conceptualized flow-path segments of the upper glacial aquifer located beneath the tribal lands was estimated using an assumed effective porosity value, and hydraulic-conductivity and hydraulic-head values that were interpolated from measured values. Groundwater travel times were estimated by dividing the length of the flow-path segment by the average linear velocity along the flow-path segment. Total estimated groundwater travel time along a conceptualized flow path, beginning near Sunrise Highway and terminating at Shinnecock Bay, is approximately 45 years using a porosity value of 30 percent. A surficial-silty unit was identified from approximately 0 to 10 ft below land surface at multiple locations beneath the tribal lands. The lithology of the surficial unit was verified by interpreted gamma log results obtained from select wells, and auger-rig drill cuttings from an observation well located near the geographic center of the tribal lands. The altitude of the unit varies with topography and was delineated along a cross section line that trends north-south along the approximate centerline (spine) of the tribal lands. The altitude of the hydrogeologic contact between the upper glacial and the Magothy aquifers generally decreases from northwest to southeast, occurs at a depth ranging from about 150 to 200 ft beneath the tribal lands, and was identified at two locations north of the tribal lands, near Sunrise Highway and Sebonac Road. Results of electrical geophysical surveys indicate that the depth to the freshwater/saltwater interface decreases from north to south with decreasing water level altitude, and the Magothy and upper glacial aquifers contain saltwater at varying depths along the north-south trending section. Results of the surveys also indicate that the Magothy aquifer beneath the tribal lands contains brackish and salty water and is not considered a source of potable water supply. In general, depth to the interface increases with increasing geographic distance from the coastline. Low water table altitudes can result in increased saltwater encroachment into the surficial aquifer beneath the tribal lands. This upward movement and shallow depth of the freshwater/saltwater interface can jeopardize water quality in wells that supply water for domestic use.

New York↗
USGS2017

Groundwater discharge by evapotranspiration, flow of water in unsaturated soil, and stable isotope water sourcing in areas of sparse vegetation, Amargosa Desert, Nye County, Nevada

This report documents methodology and results of a study to evaluate groundwater discharge by evapotranspiration (GWET) in sparsely vegetated areas of Amargosa Desert and improve understanding of hydrologic-continuum processes controlling groundwater discharge. Evapotranspiration and GWET rates were computed and characterized at three sites over 2 years using a combination of micrometeorological, unsaturated zone, and stable-isotope measurements. One site (Amargosa Flat Shallow [AFS]) was in a sparse and isolated area of saltgrass ( Distichlis spicata ) where the depth to groundwater was 3.8 meters (m). The second site (Amargosa Flat Deep [AFD]) was in a sparse cover of predominantly shadscale ( Atriplex confertifolia ) where the depth to groundwater was 5.3 m. The third site (Amargosa Desert Research Site [ADRS]), selected as a control site where GWET is assumed to be zero, was located in sparse vegetation dominated by creosote bush ( Larrea tridentata ) where the depth to groundwater was 110 m. Results indicated that capillary rise brought groundwater to within 0.9 m (at AFS) and 3 m (at AFD) of land surface, and that GWET rates were largely controlled by the slow but relatively persistent upward flow of water through the unsaturated zone in response to atmospheric-evaporative demands. Greater GWET at AFS (50 ± 20 millimeters per year [mm/yr]) than at AFD (16 ± 15 mm/yr) corresponded with its shallower depth to the capillary fringe and constantly higher soil-water content. The stable-isotope dataset for hydrogen (δ 2 H) and oxygen (δ 18 O) illustrated a broad range of plant-water-uptake scenarios. The AFS saltgrass and AFD shadscale responded to changing environmental conditions and their opportunistic water use included the time- and depth-variable uptake of unsaturated-zone water derived from a combination of groundwater and precipitation. These results can be used to estimate GWET in other areas of Amargosa Desert where hydrologic conditions are similar.

Nevada↗
USGS2018

Hydrogeology of, simulation of groundwater flow in, and potential effects of sea-level rise on the Kirkwood-Cohansey aquifer system in the vicinity of Edwin B. Forsythe National Wildlife Refuge, New Jersey

The Edwin B. Forsythe National Wildlife Refuge encompasses more than 47,000 acres of New Jersey coastal habitats, including salt marshes, freshwater wetlands, tidal wetlands, barrier beaches, woodlands, and swamps. The refuge is along the Atlantic Flyway and provides breeding habitat for fish, migratory birds, and other wildlife species. The refuge area may be threatened by global climate change, including sea-level rise (SLR). The Kirkwood-Cohansey aquifer system underlies the Edwin B. Forsythe National Wildlife Refuge. Groundwater is an important source of freshwater flow into the refuge, but information about the interaction of surface water and groundwater in the refuge area and the potential effects of SLR on the underlying aquifer system is limited. The U.S. Geological Survey (USGS), in cooperation with the U.S. Fish and Wildlife Service (USFWS), conducted a hydrologic assessment of the refuge in New Jersey and developed a groundwater flow model to improve understanding of the geohydrology of the refuge area and to serve as a tool to evaluate changes in groundwater-level altitudes that may result from a rise in sea level. Groundwater flow simulations completed for this study include a calibrated baseline simulation that represents 2005–15 hydraulic conditions and three SLR scenarios―20, 40, and 60 centimeters (cm) (0.656, 1.312, and 1.968 feet, respectively). Results of the three SLR simulations indicate that the water table in the unconfined Kirkwood-Cohansey aquifer system in the refuge area will rise, resulting in increased discharge of fresh groundwater to freshwater wetlands and streams. As sea level rises, simulated groundwater discharge to the salt marsh, bay, and ocean is projected to decrease. Flow from the salt marsh, bay, and ocean to the overlying surface water is projected to increase as sea level rises. The simulated movement of the freshwater-seawater interface as sea level rises depends on the hydraulic-head gradient. In the center of the Forsythe model area, topographic relief is 23 feet (ft) and the hydraulic-head gradient is 0.0033. In the center of the Forsythe model area, the simulated interface moved inland about 600 ft and downward about 15 ft from the baseline simulation to scenario 3 as a result of a SLR of 60 cm. In the southern part of the Forsythe model area, the topography is flatter (relief of 8 ft) and the hydraulic-head gradient is smaller (0.001). In the southern part of the Forsythe model study area, the simulated interface in this area is projected to move inland about 200 ft from the baseline simulation to scenario 3 and does not move downward.

New Jersey↗
USGS2018

Effects of surface-water and groundwater inflows and outflows on the hydrology of the Tsala Apopka Lake Basin in Citrus County, Florida

The U.S. Geological Survey, in cooperation with the Southwest Florida Water Management District, initiated a study to quantify the inflows and outflows in the Floral City, Inverness, and Hernando pools of the Tsala Apopka Lake Basin in Citrus County, Florida. This study assesses hydrologic changes in pool stages, groundwater levels, spring flows, and streamflows caused by the diversion of streamflow from the Withlacoochee River to the Tsala Apopka Lake Basin through water-control structures. A surface-water/groundwater flow model was developed using hydraulic parameters for lakes, streams, the unsaturated zone, and the underlying surficial and Upper Floridan aquifers estimated using an inverse modeling calibration technique. After calibration, the model was used to assess the relation between inflows and outflows in the Tsala Apopka Lake Basin and changes in pool stages. Simulation results using the calibrated surface-water/groundwater flow model showed that leakage rates from the pools to the Upper Floridan aquifer were largest at the deep lake cells and that these leakage rates to the Upper Floridan aquifer were the highest in the model area. Downward leakage to the Upper Floridan aquifer occurred beneath most of the extent of the Floral City, Inverness, and Hernando pools. These leakage rates depended on the lakebed leakance and the difference between lake stages and heads in the Upper Floridan aquifer. Leakage rates were higher for the Floral City pool than for the Inverness pool, and higher for the Inverness pool than for the Hernando pool. Lakebed leakance was higher for the Floral City pool than for the Hernando pool, and higher for the Hernando pool than for the Inverness pool. Simulation results showed that the average recharge rate to the surficial aquifer was 10.3 inches per year for the 2004 to 2012 simulation period. Areas that recharge the surficial aquifer covered about 86 percent of the model area. Simulations identified areas along segments of the Withlacoochee River and within land-surface depressions that receive water from the surficial aquifer. Recharge rates were largest in physiographic regions having a deep water table. Simulated heads in the Upper Floridan aquifer indicated the general flow directions in the active flow model area were from the northeast toward the southwest and then westward toward the coast, and from the southeast toward the northwest and then westward toward the coast, consistent with flow directions inferred from the estimated potentiometric surface map for May 2010. The largest inflow in the water budget of the Upper Floridan aquifer was downward leakage from the overlying hydrogeologic unit. The largest outflow in the water budget of the Upper Floridan aquifer was spring flow. The calibrated surface-water and groundwater flow model was used to simulate hydrologic scenarios that included changes in rainfall rates, projected increases in groundwater pumping rates for 2025 and 2035, no flow for the 2004–12 period through the eight water-control structures in the Tsala Apopka Lake Basin, and the removal of the Inglis Dam and the Inglis Bypass Spillway on Lake Rousseau. Scenario simulation results were compared to annual average calibrated water levels and flows from 2004 to 2012. Simulated declines in the Tsala Apopka Lake pool stages under the 10-percent lower rainfall scenario were about 0.8, 0.3, and 1.3 feet (ft) for the Floral City, Inverness, and Hernando pools, respectively. Simulated groundwater levels under the same scenario declined up to 5.4 ft in the surficial aquifer and up to 2.9 ft in the Upper Floridan aquifer. Under the projected increases in groundwater pumping rates for 2035 that represented an increase of 36 percent from average 2004 to 2012 pumping rates, the simulated declines in the Floral City, Inverness, and Hernando pool stages were, in downstream order, 0.02, 0.06, and 0.04 ft. The largest drawdown under the projected increases in groundwater pumping rates for 2035 was 2.1 ft in the surficial aquifer and about 1.8 ft in the Upper Floridan aquifer. A scenario of decreased rainfall by 10 percent caused greater declines in water levels and pool stages than projected increases in groundwater pumping rates. The simulation with no flow through the eight Tsala Apopka Lake water-control structures resulted in simulated declines in average pool stage of 1.8, 1.9, and 0.5 ft in the Floral City, Inverness, and Hernando pools, respectively. The simulated removal of the two water-control structures in Lake Rousseau caused flow to increase at Rainbow Springs by 28 cubic feet per second, an increase of 4.7 percent from the average calibrated flow for 2004 to 2012.

Florida↗
USGS2018

Groundwater chemistry and water-level elevations in bedrock aquifers of the Piceance and Yellow Creek watersheds, Rio Blanco County, Colorado, 2013–16

The Piceance and Yellow Creek watersheds in Rio Blanco County, Colorado, are known to contain important energy resources (oil shale and natural gas) and mineral resources (nahcolite). The primary sources of fresh groundwater in the Piceance and Yellow Creek watersheds are bedrock aquifers in the Uinta and Green River Formations. The aquifers are divided into an upper and lower aquifer separated by a regionally extensive semiconfining layer. These aquifers provide water to streams and springs in the watersheds and are an important resource to people living and working in the area. Development of energy and mineral resources has the potential to affect the quality of groundwater in several ways. The Bureau of Land Management and the U.S. Geological Survey began groundwater monitoring in 2010 to characterize the groundwater quality and water-level elevations of shallow bedrock aquifers in the Piceance and Yellow Creek watersheds. The purpose of this report is to present ground-water chemistry and water-level elevations collected during 2013–16. Comparisons are made to data that were collected from the bedrock aquifers from 2010 to 2012 to identify the potential for changes in water quality and water-level elevations. Appreciable changes in water-level elevations and hydraulic gradient were observed in early April 2015 in two wells completed in the upper and lower aquifers. The hydraulic gradient between the two wells was consistently downward from the upper aquifer to the lower aquifer during 2010–15; however, in early April 2015, the gradient changed from downward to upward between the two aquifers. Overall, water-level elevations declined by about 14 and 11 feet in the upper and lower aquifers, respectively, from 2013 to 2016. Previously published data estimated groundwater ages at 1,200 years old in the upper aquifer and 9,600 years old in the lower aquifer. These groundwater ages indicate that ground-water was recharged over thousands of years. With such long periods of time for aquifer recharge, declines in water-level elevation over short time steps (a few months) have important implications for sustainable management of this resource. Solution mining activities or drilling for oil and natural gas in the area could be related to the changes observed in water-level elevations in these wells; however, further investigation would be needed to evaluate causation. Changes in major-ion chemistry were evaluated in the bedrock aquifer using time series plots of select major-ion data from 2010 to 2016. Major-ion chemistry was variable for a single well from 2010 to 2016 where alkalinity and sulfate were the most variable constituents. One possible explanation for the observed changes in major-ion chemistry may be that the sample depth for that well no longer represents the most appreciable flow in the borehole. On a larger scale, potential changes in flow within the borehole may indicate changes in the regional flow system. Methane and volatile organic compound concentrations were evaluated using a similar approach to that of major ions and had similar findings. Methane concentrations in wells sampled from 2010 to 2016 were generally constant. The only exception was observed at a single well where the range of methane concentrations was from 57.4 (2010) to 4.02 milligrams per liter (2013). This is the same well where changes in water-level elevation, hydraulic gradient, and major-ion chemistry were observed, providing multiple lines of evidence to indicate change in the bedrock aquifers. Sampling of a well located in an area with little energy development but where faults or fractures could provide a path for the migration of fluids indicate mixing of groundwater between the upper and lower aquifers.

Colorado↗
USGS2020

Updated study reporting levels (SRLs) for trace-element data collected for the California Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project, October 2009–October 2018

Groundwater samples have been collected in California as part of statewide investigations of groundwater quality conducted by the U.S. Geological Survey for the Groundwater Ambient Monitoring and Assessment (GAMA) Priority Basin Project (PBP) since 2004. The GAMA-PBP is being conducted in cooperation with the California State Water Resources Control Board to assess and monitor the quality of groundwater resources used for public and domestic drinking-water supply and to improve public knowledge of groundwater quality in California. Quality-control samples (including but not limited to field, equipment, and source-solution blanks) were collected to evaluate and quantify the quality of the groundwater sample results. The GAMA-PBP previously determined study reporting levels (SRLs) for trace-element results based primarily on field blanks collected in California from May 2004 through March 2013. SRLs are raised reporting levels used to reduce the likelihood of reporting false detections attributable to contamination bias. The purpose of this report is to identify any changes in the pattern or magnitude of concentrations or detections in field blanks since the last evaluation that would require changing or ending the use of SRLs implemented in October 2009. Constituents analyzed were aluminum, antimony, arsenic, barium, beryllium, boron, cadmium, chromium, hexavalent chromium, cobalt, copper, iron, lead, lithium, manganese, molybdenum, nickel, selenium, silver, strontium, thallium, uranium, vanadium, and zinc. For this review, data from 167 field blanks collected from October 2009 through October 2018 by the GAMA-PBP for trace elements were compiled. Based on a consistent pattern of decreasing cobalt and manganese concentrations in field blanks from 2009 to 2013, the GAMA-PBP decided to reevaluate all trace-element SRLs, effectively setting an end date for previously defined SRLs. Beginning October 2013, SRLs would be determined from field-blank data collected through October 2018. The detection frequency and upper limit of potential contamination bias (BD-90/90) were determined from field blanks for each trace element. The BD-90/90, that is, the upper 90-percent confidence limit of the 90th percentile concentration of potential extrinsic contamination, was calculated by assuming the binomial probability distribution. These results were compared to each constituent’s detection limit to determine whether an SRL was necessary to minimize the potential for detections in the groundwater samples, attributed principally to contamination bias. Results of the evaluation were used to set SRLs for trace-element data collected by the GAMA-PBP between October 2013 and October 2018. Trace elements prescribed an SRL based on this review were hexavalent chromium, cobalt, copper, lead, and zinc. This review also resulted in the removal of SRLs from iron, manganese, molybdenum, and nickel. Although an SRL for hexavalent chromium could not be evaluated in the earlier reviews because the data were not collected regularly until 2015, one was established herein as 0.34 micrograms per liter (µg/L). The SRL for cobalt, as previously implemented, had been to reject all results; it was changed to 0.16 µg/L following a reduction in cobalt field-blank detection frequency resulting from mitigation steps, starting in 2014, aimed at reducing contamination bias introduced by high-capacity capsule filters used during sample collection. The SRL for copper did not change, and the SRL for lead changed very little based on this review. Lastly, the SRL for zinc was lowered from 6.2 µg/L to 3.9 µg/L.

California↗
USGS2020

Hydrogeology and groundwater flow in alluvial deposits, north Summerset, South Dakota

The city of Summerset is a growing community in west South Dakota. The Sun Valley Estates subdivision in the north part of the city was developed on unconsolidated deposits surrounded by steep terrain. During years with greater than normal precipitation, particularly in 2019, groundwater levels increased in the unconsolidated deposits and caused damage to stormwater systems, sewer infrastructure, and houses with basements. The U.S. Geological Survey, in cooperation with the City of Summerset, completed a study of the hydrogeology and groundwater flow in the alluvial aquifer part of the unconsolidated deposits in north Summerset to understand the groundwater system in the area and to provide hydrogeologic information in support of future development planning. The study area included most of the Sun Valley Estates subdivision in the north part of the city of Summerset in the east Black Hills of west South Dakota. About 0.7 square mile of water-bearing alluvial deposits is included in the study area. Precipitation in the study area from 2017 to 2019 was compared to the monthly normal values at a nearby climate site. The largest departure from normal was in May 2019 with precipitation exceeding the monthly normal by about 5 inches (in.). All months in 2019, except March, exceeded the monthly normal precipitation. Cumulative departure from normal precipitation in 2019 increased from about 4 in. greater than normal in January to about 18 in. greater than normal in December. The geologic setting of the study area is characterized by the surrounding Black Hills. Unconsolidated Quaternary-age deposits overlie consolidated to partially consolidated Mesozoic-age and Paleozoic-age shales, sandstones, and limestones. Surficial deposits of alluvium and other unconsolidated deposits are the primary surficial geologic units in the study area and form the components of the alluvium hydrogeologic unit of the study area. Results from previous studies of alluvium along nearby Rapid Creek estimated hydraulic conductivity to range from 89 to 2,292 feet per day (ft/d), transmissivity to range from 1,001 to 32,083 feet squared per day, and storage coefficients to range from 0.0002 to 0.16. Hydraulic conductivity and transmissivity generally decreased downstream along Rapid Creek (west to east). Slug tests were completed August 16, 2019, at two observation wells completed in the alluvial aquifer in the Sun Valley Estates subdivision to determine hydraulic conductivity. Hydraulic conductivity estimated from AQTESOLV curve-fitting analysis using the Bouwer-Rice method for all slug-in and slug-out trials from two observation wells in the study ranged from 0.20 to 0.26 ft/d for well 441318103220001 (SunValley1 well) and from 0.54 to 14 ft/d for well 441319103215701 (SunValley2 well). The mean, median, and standard deviation of all trials at both wells were 4.3 ft/d, 0.8 ft/d, and 5.6 ft/d, respectively. The extent of the alluvial aquifer was determined by geologic maps and lithologic logs. Alluvial deposits in the study area extend to about 1 mile in the north–south direction and about 1.5 miles in the southeast–northwest direction. The direction of groundwater flow was estimated using water-level records and topographic maps. The resulting potentiometric map indicated that groundwater in the alluvial aquifer under the Sun Valley Estates subdivision originates from higher elevations of the west part of the area of interest and from streams in the southeast part. Recharge and evapotranspiration estimates were results from a Soil-Water Balance model that calculated a matrix of recharge for 2019 with values ranging from 0 to 11.4 in. and an annual mean value of 5.1 in. across the study area. Soil-Water Balance-estimated potential evapotranspiration for 2019 ranged from 28.90 to 28.75 in. and the estimated annual mean was 28.86 in. across the study area. Estimated groundwater budget components for the alluvial aquifer in the area of interest included inflows and outflows. Total estimated groundwater budget components for inflows for 2019 were about 66 percent from recharge, 33 percent from streamflow, and 1 percent from inflow from adjacent aquifers. Total estimated outflows were about 99-percent evapotranspiration and less than 1-percent outflow to adjacent aquifers.

South Dakota↗
USGS2021

Trends in groundwater levels in and near the Rosebud Indian Reservation, South Dakota, water years 1956–2017

The U.S. Geological Survey (USGS), in cooperation with the Rosebud Sioux Tribe, completed a study to characterize water-level fluctuations in observation wells to examine driving factors that affect water levels in and near the Rosebud Indian Reservation, which comprises all of Todd County. The study investigates concerns regarding potential effects of groundwater withdrawals and climate conditions on groundwater levels within an area that includes Todd County and a surrounding area that extends 10 miles north, east, and west of the county border. Characterization of water-level fluctuations in observation wells and relative driving factors was accomplished by statistical trend analysis. Two statistical methods were used for analysis of temporal trends for climatic and hydrologic data. To determine which trend analysis to use, applicable datasets were tested for statistically significant short-term persistence (STP). In the absence of significant STP, existence of statistical trends was determined using the standard Mann-Kendall test for probability values less than or equal to 0.10 (90-percent confidence level); however, a modified Mann-Kendall test was used for datasets where statistically significant STP was detected. Trend magnitudes were computed using the Sen’s slope estimator. Monthly data from the Parameter-elevation Regressions on Independent Slopes Model (PRISM) were aggregated to obtain annual and seasonal datasets for total precipitation, minimum air temperature ( T min ), and maximum air temperature ( T max ) for the study area and a surrounding buffer area. Trend tests for total precipitation, T min , and T max were completed for annual and seasonal time series for water years 1956–2017, which is about 2 years before the earliest available water-level measurements. A 2-year offset was arbitrarily selected because scrutiny of water-level and precipitation data indicated that responses of groundwater levels for many of the observation wells lagged major changes in precipitation patterns by about 2 years. Statistically significant upward trends were detected for annual precipitation and annual T min for almost all of the study area and the surrounding buffer area. Statistically significant downward trends in T max were detected for a very small part of the study area; however, the sparse spatial coverage reduces confidence that these are true trends. Spatial distributions of statistically significant trends in seasonal climate data were generally similar to the annual trends, but with substantial differences in the spatial density of the trends. Groundwater trends for 58 observation wells were analyzed for three separate water-level parameters (minimum, median, and maximum) because wells are measured sporadically and data are biased towards more frequent measurements during periods of heaviest irrigation demand. Trends in the time series of annual precipitation (from PRISM) starting 2 years earlier than for the associated water-level trend also were analyzed for the location of each individual observation well. Sen’s slope and Mann-Kendall probability values (p-values) were computed for the three water-level parameters and for the annual precipitation time series. Graphs showing results of trend analyses for each observation well also showed changes over time in the sum of licensed groundwater withdrawals within six specified radii (0.5, 1, 2, 3, 4, and 5 miles) of each well as a qualitative indicator of proximal groundwater demand. Of all 58 observation wells considered, 28 wells had significant upward trends for at least one of the three water-level parameters, 11 wells had significant downward trends for at least one water-level parameter, and 19 wells did not have any significant trends. Significant upward trends in annual precipitation were detected for 48 of the 58 wells. Results of trend analyses likely show the effects of groundwater withdrawals on water levels in the Ogallala aquifer in areas of substantial demand. Precipitation trends are significantly upward for 43 of the 48 wells completed in the Ogallala aquifer that were analyzed. Of the 48 Ogallala aquifer wells, 24 had significant upward trends for at least one water-level parameter (17 with all 3); however, 10 wells had statistically significant downward trends for at least one water-level parameter (8 with all 3 parameters). All but one of the wells with significant downward trends are located in the south-central part of the study area where licensed irrigation withdrawals are concentrated.

South Dakota↗
USGS2021

Effects of climate and land-use change on thermal springs recharge—A system-based coupled surface-water and groundwater-flow model for Hot Springs National Park, Arkansas

A three-dimensional hydrogeologic framework of the Hot Springs anticlinorium beneath Hot Springs National Park, Arkansas, was constructed to represent the complex hydrogeology of the park and surrounding areas to depths exceeding 9,000 feet below ground surface. The framework, composed of 6 rock formations and 1 vertical fault emplaced beneath the thermal springs, was discretized into 19 layers, 429 rows, and 576 columns and incorporated into a 3-dimensional steady-state groundwater-flow model constructed in MODFLOW-2005. Historical daily mean thermal spring flows were simulated for one stress period of approximately 34 years (1980–2014), chosen to represent the period of record for historical climate data used in the quantification of the boundary conditions. The groundwater-flow model was manually calibrated to historical daily mean thermal spring flows of 88,000 cubic feet per day observed over a 12-year period of record (1990–1995 and 1998–2005) at the thermal springs collection system. Calibration was achieved by calculating starting heads and general head boundary conditions from the Bernoulli equation and then adjusting the horizontal and vertical hydraulic conductivities of the rock formations and vertical fault and the hydraulic conductance of head-dependent flux boundaries. The groundwater-flow model was coupled to a surface-water model developed in the Precipitation-Runoff Modeling System (PRMS) by using PRMS-simulated gravity drainage as a specified flux recharge boundary condition in the groundwater-flow model. Together, the coupled models were used to (1) locate the areas of groundwater recharge to the thermal springs in the discretized hydrogeologic framework by using forward and reverse particle-tracking capabilities of MODPATH, (2) simulate the effects of variable recharge rates on the spring flows at the thermal springs, and (3) assess possible effects of climate and land-use change on the long-term variability of spring flows at the thermal springs. Forward and backward particle-tracking maps indicated that the most prevalent areas of recharge in the discretized hydrogeologic framework used in this study were within about 0.6–0.9 mile of the thermal springs. Forward particle tracking indicated a recharge area southwest of the thermal springs that corresponded to a location where the predominant lithologies are the Arkansas Novaculite, Hot Springs Sandstone, and Bigfork Chert. Backward particle tracking indicated a second localized area of recharge to the northeast of the thermal springs that corresponded to a location where the dominant lithology is the Bigfork Chert. The groundwater-flow model indicated that the most probable recharge formations are the Arkansas Novaculite, Bigfork Chert, and Hot Springs Sandstone. The simulated effects of climate and land-use changes on the variability of the spring-flow rates at the thermal springs generally resulted in reductions of thermal spring flow attributed to urban development and more extreme climates characterized by elevated mean surface air temperatures. The groundwater-flow model predicted a linear relation between the thermal spring discharge and the cumulative recharge volume applied to the hydrogeologic framework, and the positive slope of the predicted relation between recharge and simulated thermal spring flow indicates that more extreme precipitation events that supply more recharge may in fact increase the thermal spring-flow rates.

Arkansas↗
USGS2021

Delineation of areas contributing groundwater and travel times to receiving waters in Kings, Queens, Nassau, and Suffolk Counties, 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 New York State Department of Environmental Conservation initiated a program to delineate areas contributing groundwater to coastal water bodies by assembling a comprehensive dataset of areas contributing groundwater, travel times, and groundwater discharges to streams, lakes, marine surface waters, and subsea discharge boundaries. Steady-state, 25-layer regional, three-dimensional finite-difference groundwater-flow models of average regional hydrologic conditions were used for particle-tracking analysis to delineate areas contributing groundwater to 843 water bodies. Two steady-state conditions were simulated: recent conditions from 2005 to 2015 and predevelopment conditions of about 1900. About 14 million particles were evenly distributed across the water table and tracked forward to discharge zones. Using a uniform porosity of 25 percent, simulated recent condition travel times ranged from less than 2 years to greater than 10,000 years and were visualized in 11 travel time intervals. About 85 percent of particle travel times from the water table to points of discharge are less than 100 years. Simulated particle-tracking ending zones represented 843 receiving water bodies, based on the New York State Department of Environmental Conservation water body inventory and priority water bodies list. Areal delineation of travel-time intervals and areas contributing groundwater to water bodies were generated and are summarized with total groundwater outflow for each water body.

New York↗
USGS2021

Groundwater chemistry, hydrogeologic properties, bioremediation potential, and three-dimensional numerical simulation of the sand and gravel aquifer at Naval Air Station Whiting Field, near Milton, Florida, 2015–20

The U.S. Geological Survey completed a study between 2015 and 2020 of groundwater contamination in the sand and gravel aquifer at a Superfund site in northwestern Florida. Groundwater-quality samples were collected from representative monitoring wells located along a groundwater-flow pathway and analyzed in the field and laboratory. In general, ambient groundwater in the sand and gravel aquifer is acidic, dilute, and oxic. Groundwater age-dating results indicate recharge to the contaminated parts of the aquifer occurred between the 1970s and 1980s. Natural gamma, electromagnetic induction, and borehole nuclear magnetic resonance logs indicated that aquifer hydraulic conductivities generally increased with depth as the aquifer formation material became coarser, characteristic of a prograding marginal-marine delta depositional environment. Aquifer formation material incubated with radiocarbon (carbon-14) cis -1,2-Dichloroethylene demonstrated biodegradation directly to carbon dioxide in contaminated and uncontaminated parts of the aquifer. A three-dimensional, numerical groundwater-flow MODFLOW model of the sand and gravel aquifer in the study area was constructed. The calibrated model reasonably reproduced measured groundwater heads and streamflows. Moreover, the model can be used to run simulations of outcomes of potential remedial strategies, such as monitored natural attenuation, as part of future feasibility studies in the area.

Florida↗
USGS2023

Geology, hydrology, and groundwater contamination in the vicinity of Central Chemical facility, Hagerstown, Maryland

The soil and groundwater at the Central Chemical facility, Hagerstown, Maryland, are contaminated due to the blending and production of pesticides and fertilizers during much of the 20th century. Remedial investigations focus on two operable units (OU) consisting of the surface soils and waste disposal lagoon (OU-1) and the groundwater (OU-2). The contaminants of concern (COC) for groundwater include 41 compounds categorized within the subgroups of volatile organic compounds (VOCs), semi-volatile organic compounds (SVOCs), pesticides, and metals. The purpose of this report is to provide a conceptual site model of the hydrogeology and groundwater contaminant transport at and near the Central Chemical facility. The conceptual model was developed through review, synthesis, and interpretation of the results of hydrogeologic, soil, and other environmental investigations conducted at and in the vicinity of the facility in recent decades and is intended to support plans for environmental remediation of the groundwater in OU-2. The extent and nature of the groundwater contaminant plume associated with the bedrock was characterized for OU-2 of the site. Lithologic and structural comparisons between shallow soil, weathered rock, and epikarst and deeper competent but bedded, dipping, fractured, and karstic limestones illustrate two connected flow systems—a surficial flow system consisting of the unconsolidated overburden and epikarst and a structurally dominant bedrock flow system below the epikarst. Uncertainties exist regarding the nature and transport of contaminants within the epikarst system particularly within voids and perched epikarst water tables. Karst dissolution features are observed within the site including sinkholes and dissolution voids within wells at the site. Of interest, one well in the northern part of the study area (MW-J-71) appears to have a dissolution void connected to an offsite well (OW-2-115) farther to the north. This connection is supported by groundwater level data and elevated concentrations of total suspended solids (TSS) and chlorobenzene in both wells. The high level of TSS supports the possibility of offsite transport of particle-bound contaminants within the conduit system. Episodically elevated concentrations of COC from different groups also were observed within select wells in the epikarst near the waste disposal lagoon (particularly MW-A-51). The variability observed between different COC within the same well may be the result of additional contaminated source materials unrelated to the disposal lagoon. Storage and episodic transport of contaminated material within the epikarst system has the potential to hinder remediation efforts if not considered in the remedial action.

Maryland↗
USGS2022

Simulation of regional groundwater flow and groundwater/lake interactions in the Central Sands, Wisconsin

A multiscale, multiprocess modeling approach was applied to the Wisconsin Central Sands region in central Wisconsin to quantify the connections between the groundwater system, land use, and lake levels in three seepage lakes in Waushara County, Wisconsin: Long and Plainfield (The Plainfield Tunnel Channel Lakes), and Pleasant Lakes. A regional groundwater-flow model, the Newton Raphson formulation of the U.S. Geological Survey modular finite-difference flow model groundwater-modeling package (MODFLOW-NWT), centered on the lakes, was used to extend regional surface-water boundaries to provide boundary conditions for two focused inset models, in the hydrologic simulation modeling package (MODFLOW 6), at higher resolution around the lakes. Land use and groundwater use were simulated at a regional scale using the Soil Water Balance model, which provided recharge and water-use boundary conditions for the MODFLOW models. Agricultural irrigation is the primary groundwater use in the area. Land and groundwater-use scenarios representing no irrigation, current (2018) irrigation, and potential future irrigation were simulated with the groundwater-flow model and the lake levels over a 38-year representative climate period.

Wisconsin↗
USGS2023

Groundwater-flow model of the Treasure Valley, southwestern Idaho, 1986–2015

Most of the population of the Treasure Valley and the surrounding area of southwestern Idaho and easternmost Oregon depends on groundwater for domestic supply, either from domestic or municipal-supply wells. Current and projected rapid population growth in the area has caused concern about the long-term sustainability of the groundwater resource. In 2016, the U.S. Geological Survey, in cooperation with the Idaho Water Resource Board and the Idaho Department of Water Resources, began a project to construct a numerical groundwater-flow model of the westernmost portion of the western Snake River Plain aquifer system, called the Treasure Valley. The development of the model was guided by several objectives, including: to improve the understanding of groundwater and surface water interactions; to facilitate conjunctive water management; to provide a tool for water resources planning; and to provide a tool for water allocation. The model was constructed with a spatial scale and level of detail that aimed to meet these objectives while balancing the sometimes-competing goals of fast runtimes, numerical stability, usability, and parsimony. The Treasure Valley Groundwater Flow Model (TVGWFM) is a three-dimensional finite-difference numerical model constructed using MODFLOW 6 (Langevin and others, 2017, Documentation for the MODFLOW 6 Groundwater Flow Model: U.S. Geological Survey Techniques and Methods, book 6, chap. A55, 197 p., https://doi.org/10.3133/tm6A55 ). The model covers the westernmost portion of the western Snake River Plain and is discretized into a regular grid of 64 by 65 cells with a side length of 1 mile and 6 layers of varying depth and active area. A historical model period was developed consisting of 360 month-long stress periods for 1986–2015. The model builds upon previous modeling efforts by adding a transient period, incorporating new head and discharge observations to constrain parameters, incorporating information from the hydrogeologic framework model (HFM) of Bartolino (2019, Hydrogeologic framework of the Treasure Valley and surrounding area, Idaho and Oregon: U.S. Geological Survey Scientific Investigations Report 2019–5138, https://doi.org/10.3133/sir20195138 ) and incorporating refined estimates of evapotranspiration and irrigation classification of lands in the study area. The TVGWFM includes all significant components of recharge to and discharge from the aquifer. Inflows include canal seepage, irrigation and precipitation recharge, mountain-front recharge, rivers and stream seepage, and seepage from Lake Lowell. Outflows include discharge to agricultural drainage ditches, discharge to rivers and streams, pumping, and discharge to Lake Lowell. Each recharge or discharge component is represented separately using individual MODFLOW 6 packages. Parameter values were derived with a combination of trial-and-error steps and automated parameter estimation using PEST software (Doherty, J.E., 2005, PEST, model-independent parameter estimation–User manual: Watermark Numerical Computing, https://pesthomepage.org/documentation ). Parameter estimates were constrained with several types of observation data, including water levels, water level changes, vertical water level differences, drain discharges, change in drain discharges, river seepage, seepage from Lake Lowell, and change in seepage from Lake Lowell. Material properties from the hydrogeologic framework were also used to assign the minimum and maximum values of some parameters. A final parameter realization was reached that minimized residuals between the observed and modelled values for the various observation groups. Mean residuals for the observation groups were 15.4 feet (ft) for water levels, 0.2 ft for water level changes, 19.4 ft for vertical water level differences, −3.9 cubic feet per second (ft 3 /s) for drain discharges, 0.0 ft 3 /s for changes in drain discharge, 45.0 ft 3 /s for river seepage, −40.1 ft 3 /s for Lake Lowell seepage, and 126.3 ft 3 /s for changes in Lake Lowell seepage. The quality of the model’s fit to observations varied spatially, with notable areas of under- or over-simulation of water levels present to the northwest and southwest of Lake Lowell, in the foothills along the eastern model boundary, and near the City of Eagle. Trends were observed in the residuals of many of the observation groups, indicating that the model is missing or not fully reproducing some phenomena that are observed in the system. The TVGWFM can be used as a tool for water resource planning, for understanding the interactions of groundwater and surface water at a basin scale, and for facilitating conjunctive management, but may lack the precision needed for water rights administration at a local scale. Additional sources of uncertainty or limitations of the model are noted. The quantity and spatial distribution of canal seepage and infiltration of irrigation water recharge, the largest sources of recharge to the system, are unknown and approximated indirectly. There is poor understanding of how canal seepage and incidental recharge change as land is converted from agricultural (irrigated) to suburban (semi-irrigated). These uncertainties will affect any scenarios that investigate changes to land use or irrigation practices. Finally, the model has relatively high water-level residuals around and to the southwest of Lake Lowell and should not be used to estimate water level effects in that region. The model was built with multiple, broadly expressed objectives and did not optimize performance for specific uses. However, the model and the tools included in an associated data release provide ample flexibility to improve the model for future uses. Adjustments and improvements could be made by refining the model in an area of interest, collecting additional calibration data, applying more rigorous boundary conditions, or re-estimating model parameters to optimize model performance for a specific model forecast.

Idaho↗
USGS2024

Simulation of groundwater flow in the Long Island, New York regional aquifer system for pumping and recharge conditions from 1900 to 2019

The U.S. Geological Survey has developed a transient, groundwater-flow model that simulates hydrologic conditions in the Long Island aquifer system as part of an ongoing (since 2016) multiyear, cooperative investigation with the New York State Department of Environmental Conservation. The goals of this investigation are to assist stakeholders and resource managers to evaluate the response of the hydrologic system to changes in future hydraulic stresses. Responses in the hydrologic system include changes in water levels in the hydrogeologic units; discharge to streams, coastal waters, and subsurface infrastructure; and the extent of saline groundwater in the aquifers. Hydraulic stresses include future water-supply management and changes in land use and infrastructure. The numerical model synthesizes a diverse set of physiographic, geologic, climatic, land-use, and historical population, water use, and infrastructure data to physically represent the Long Island aquifer system from land surface to bedrock and to simulate annual hydrologic conditions between 1900 and 2019. A three-dimensional hydrogeologic framework was developed from existing and recently collected borehole geologic and geophysical data collected as part of a companion drilling program. Water-transmitting properties of the principal aquifer sediments were defined in three dimensions from new and existing lithologic logs. The distribution of recharge from precipitation was estimated from landscape characteristics and climate data. Anthropogenic recharge from wastewater, leaky infrastructure, and storm runoff were estimated from population, infrastructure, and pumping data. Water-use data, including well locations, depths, and pumping rates, were obtained from historical sources and records and used to estimate pumping stresses continuously in time and space, at an annual average time scale. The data were incorporated into a three-dimensional numerical model using the U.S. Geological Survey finite difference modeling code MODFLOW 6; the model encompassed all of Long Island and surrounding surface waters and simulated historical hydrologic conditions from 1900 to 2019. The calibration process involved trial and error adjustments using prior knowledge to improve general fit to observations followed by an inverse calibration to update and optimize input parameters, using an iterative ensemble smoother algorithm implemented in PEST++ version 5.0. This resulted in a model that generally was in good agreement with observed, dynamically varying hydrologic conditions from 1900 to 2019. The calibrated model was used to develop two base-case models for scenario testing of future, hypothetical conditions where one represented average-annual conditions, and one represented average-seasonal conditions from 2010 to 2019. The model representing average-annual conditions was modified further to represent an alternate sea-level position of 6 feet above the North American Vertical Datum of 1988, and the model representing average-seasonal conditions was modified to represent the average seasonal effects of a 5-year drought imposed upon current hydrologic conditions. Recharge is the sole source of water to the aquifer system; groundwater discharges to coastal water and streams and is withdrawn by pumped wells. Model-estimated annual recharge ranged from about 11 inches in 1965 to 41 inches in 1983. On average, from 2010 to 2019, about 23 percent of water was pumped from wells, and about 47 and 27 percent discharged to coastal waters and streams, respectively; the remaining 4 percent was water that moved into storage in the aquifer matrix. Water levels on Long Island vary naturally during time in response to changes in recharge; the amount of variation is largest in the interior of the island, in areas with highest water table altitudes near groundwater divides and lowest near streams and the coastal waters. The total range of water table altitudes on Long Island between 1900 and 2019 ranged from near 0 to more than 70 feet in western parts of Long Island. The largest range in altitudes is in New York City and is associated with areas of large historical withdrawals between the 1920s and the late 1980s. Water table altitudes generally varied by less than 10 feet in eastern Suffolk County, where the aquifer is under more natural conditions. Saltwater intrusion is of great concern on Long Island, particularly in western Long Island where both the unconfined and confined parts of the aquifer system have been intruded in response to large-scale groundwater withdrawals; however, the volume of freshwater in the islandwide aquifer system only has changed by about 5 percent between 1900 and 2019. The decadal change in the freshwater volume was largest during the early and mid-20th century, corresponding to the largest historical pumping, but that volume change did not exceed 1 percent. The negligible change in freshwater volume suggests that saltwater intrusion as of 2019 was limited at an islandwide scale but continues to occur in local areas of Queens and Nassau Counties, adversely affecting current water supplies and limiting future water supplies for affected communities. The regional groundwater model developed for this investigation is a tool that can be used to help determine the viability of current and future water supplies at a regional scale and can be used to support development of additional models at finer scale to support more focused assessments of groundwater sustainability.

New York↗
USGS2025

Groundwater response to managed aquifer recharge at the Southeast Houghton Artificial Recharge Project in Tucson, Arizona

Managed aquifer recharge is a widespread practice for storing water in the subsurface as groundwater. At a managed aquifer recharge facility in southern Arizona, groundwater-level and repeat microgravity data were collected to monitor aquifer response. These data were used to inform parameter identification for an unsaturated-zone flow model used to simulate the recharge process. The facility, the Southeast Houghton Artificial Recharge Project (SHARP), consists of 3 surface basins (about 27,600 square meters [6.8 acres] total surface area) where recycled water is distributed in recharge cycles lasting several months, with dry periods in between. During the study period, December 2020–December 2022, Tucson Water (the City of Tucson’s water utility) reported 6.56×10 6 cubic meters of water (5,320 acre-feet) recharged. Monitoring included groundwater-level observations at 3 monitoring wells and repeat microgravity measurements at as many as 22 locations (some stations were destroyed between surveys). Six gravity surveys were carried out using absolute- and relative-gravity meters. Large gravity increases, more than 250 microgals, were observed during the first repeat survey, 3.5 months after the start of recharge, but only in the immediate vicinity of the recharge basins. Data show that water moved downward to the water table, and storage changes in the unsaturated zone away from the facility were likely minimal. Gravity decreased at stations more than 1 kilometer from the facility, consistent with regional groundwater-level changes. Groundwater-level increases in wells adjacent to the recharge basins began 2 months after the second repeat gravity survey, and 5.5 months after recharge began. Unsaturated-zone flow modeling was carried out using software that simulates water movement and parameter estimation. Model calibration was carried out by minimizing an objective function calculated from the differences between simulated and observed groundwater levels, and between simulated and observed repeat microgravity data. Including repeat microgravity data in the objective function reduced the uncertainty in estimated parameter values for saturated hydraulic conductivity and saturated water content. Modeling indicated that the unsaturated zone between the recharge basins and the water table does not become saturated even after 685 days of simulated infiltration. This gradual wetting may account for increasing infiltration rates over time, as hydraulic conductivity increases with increasing water content. Unsaturated-zone water content decreased rapidly between recharge cycles. Model-simulated groundwater mounding extended about 1 kilometer from the center of SHARP after the 685-day period following the onset of recharge.

Arizona↗