Geology topics

Groundwater & aquifers

Find USGS groundwater and aquifer publications. Search an aquifer name or place, compare source-supplied dates and open reports for hydrogeology, recharge and water-level studies.

Follow this topic with RSS

At least 883 records · Page 49Linked to original sources

USGS2002

Using chemical, hydrologic, and age dating analysis to delineate redox processes and flow paths in the riparian zone of a glacial outwash aquifer‐stream system

A combination of chemical and dissolved gas analyses, chlorofluorocarbon age dating, and hydrologic measurements were used to determine the degree to which biogeochemical processes in a riparian wetland were responsible for removing NO 3 − from groundwaters discharging to the Otter Tail River in west central Minnesota. An analysis of river chemistry and flow data revealed that NO 3 − concentrations in the river increased in the lower half of the 8.3 km study reach as the result of groundwater discharge to the river. Groundwater head measurements along a study transect through the riparian wetland revealed a zone of groundwater discharge extending out under the river. On the basis of combined chemical, dissolved gas, age date, and hydrologic results, it was determined that water chemistry under the riparian wetland was controlled largely by upgradient groundwaters that followed flow paths up to 16 m deep and discharged under the wetland, creating a pattern of progressively older, more chemically reduced, low NO 3 − water the farther one progressed from the edge of the wetland toward the river. These findings pose challenges for researchers investigating biogeochemical processes in riparian buffer zones because the progressively older groundwaters entered the aquifer in earlier years when less NO 3 − fertilizer was being used. NO 3 − concentrations originally present in the groundwater had also decreased in the upgradient aquifer as a result of denitrification and progressively stronger reducing conditions there. The resulting pattern of decreasing NO 3 − concentrations across the riparian zone may be incorrectly interpreted as evidence of denitrification losses there instead of in the upgradient aquifer. Consequently, it is important to understand the hydrogeologic setting and age structure of the groundwaters being sampled in order to avoid misinterpreting biogeochemical processes in riparian zones.

Water Resources Research↗

USGS2018

Relations between total phosphorus and orthophosphorus concentrations and rainfall, surface-water discharge, and groundwater levels in Big Cypress Seminole Indian Reservation, Florida, 2014–16

The Seminole Tribe of Florida (the Tribe) is partnering with the U.S. Environmental Protection Agency to develop a numeric phosphorus criterion for the 52,000-acre Big Cypress Seminole Indian Reservation (BCSIR), which is located downgradient of the Everglades Agricultural Area, and of other public and private lands, in southeastern Hendry County and northwestern Broward County in southern Florida. The U.S. Geological Survey (USGS), in cooperation with the Tribe, used water-quality data collected between October 2014 and September 2016 by the Tribe and the South Florida Water Management District (SFWMD), along with data from rainfall gages, surface-water stage and discharge gages, and groundwater monitoring wells, to (1) examine the relations between local hydrology and measured total phosphorus (TP) and orthophosphorus (OP) concentrations and (2) identify explanatory variables for TP concentrations. Of particular concern were conditions when TP exceeded 10 parts per billion (ppb) (0.01 milligram per liter [mg/L]) given that the State of Florida and the Miccosukee Tribe of Indians Alligator Alley Reservation (located downstream of the BCSIR) have adopted a 10-ppb maximum TP criterion for surface waters. From October 2014 to September 2016, the Tribe collected 47–52 samples at each of nine water-quality sites for analysis of TP and OP, except at one site where 28 samples were collected. For all sites sampled, concentrations of TP (as phosphorus [P]) ranged from less than 0.002 mg/L (2 ppb) to a maximum of nearly 0.50 mg/L (500 ppb), whereas concentrations of OP (as P), the reactive form of inorganic phosphorus readily absorbed by plants and (or) abiotically absorbed, ranged from less than 0.003 mg/L (3 ppb) to a maximum of 0.24 mg/L (240 ppb). The median and interquartile ranges of concentrations of TP and OP in the samples collected in 2014–16 by the Tribe were similar to the median and interquartile ranges of concentrations in samples collected by the SFWMD at nearby sites during the same period. Differences in concentrations can likely be explained by differences in sample collection methods, sampling locations, sample collection time, and the hydrology during sampling or by the number of samples collected. A major limitation of this study was the short duration of sample collection, which covers a limited range of hydrologic conditions within the BCSIR. The effect of surface-water and groundwater hydrologic conditions on TP and OP concentrations was assessed by using rainfall data and surface-water stage and discharge records. The highest TP and OP concentrations occurred during peak surface-water flows in the canals following long dry periods. Concentrations of TP and OP increased internal to the BCSIR in the western half of the BCSIR during wet periods, but increased concentrations tended to lag behind rainfall events, likely because control structures upstream of sampling sites do not release flows until the water levels in the canals reach predetermined levels. This pattern may indicate that bed sediments in the canals contain high concentrations of phosphorus that becomes resuspended during high flows or that phosphorus salts that had accumulated on dry land during dry periods are carried into the canals by runoff. The largest TP spikes usually occurred at the beginning of high-flow events, but then quickly tapered off even when flows remained high. Groundwater flows were assessed in the BCSIR by using groundwater level observations from two preexisting USGS monitoring well clusters, each characterized by a shallow well installed in the surficial aquifer system and a deeper well installed in the intermediate aquifer system. Groundwater levels were evaluated with respect to surface-water levels and discharge in the BCSIR during the period of surface-water sampling. During dry conditions water levels in canals were often higher than groundwater levels in the surficial aquifer, indicating the potential for surface water to recharge the surficial aquifer. During wetter conditions, this trend reversed, and there was potential for shallow groundwater discharge into the canals. From October 2014 to September 2016, concentrations of TP tended to decrease as surface-water inflows moved across the BCSIR from north to south. In both the western and eastern halves of the reservation, the mean concentration of TP was lower in the surface-water outflows from the BCSIR than in the inflows. The mean concentration of TP in the inflows to the western reservation was 0.04 mg/L (40 ppb), whereas the mean concentration of TP in the outflows was 0.03 mg/L (30 ppb). In the eastern reservation, the mean concentration of TP in the inflows was 0.07 mg/L (70 ppb), whereas the mean concentration of TP in the outflows was 0.04 mg/L (40 ppb). TP and OP concentrations were evaluated relative to other water-quality parameters, including turbidity, suspended solids, nitrate plus nitrite, dissolved oxygen, pH, and specific conductance, to determine if any relations existed between TP and other variables. Weak relations were indicated for turbidity and suspended solids at two sites, which indicates that there may be a relation of increased TP to mobilization of sediment.

Florida↗

USGS2020

Water-table elevation maps for 2008 and 2016 and water-table elevation changes in the aquifer system underlying eastern Albuquerque, New Mexico

The addition of surface water from the San Juan-Chama Drinking Water Project to the Albuquerque water supply and the reduction in per capita water use has led to decreased groundwater withdrawals. This decrease in withdrawals has resulted in rising groundwater levels since 2008 in portions of the aquifer underlying Albuquerque. The wells used to assess the Kirtland Air Force Base Bulk Fuels Facility (KAFB BFF) ethylene dibromide (EDB) groundwater contamination were installed with well screens that crossed the water table in order to monitor and sample groundwater within the EDB plume. While replacement wells have been installed, an understanding of the water-table response to decreases in regional groundwater withdrawals is required to evaluate the monitoring well network. Water-table elevation maps of the aquifer underlying the Albuquerque metropolitan area east of the Rio Grande for 2008 and 2016 and a map of the change in elevations in this 8-year period provide an improved understanding of the water-table elevations and the changes that are occurring. The water-table elevation contours for both 2008 and 2016 show that groundwater generally flows from the Rio Grande and from the mountain-front recharge in the southeast toward the center of the study area, a major groundwater pumping center. The water-table elevation increased in most of the study area from 2008 to 2016. The area of greatest increase in the water-table elevation covers most of the northeastern part of the study area, where there has historically been pumping-related drawdown and subsequent groundwater-level rises in the production zone of the Santa Fe Group aquifer system.

New Mexico↗

USGS2023

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

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

New Jersey↗

USGS2009

Source, Transport, and Fate of Groundwater Contamination at Site 45, Marine Corps Recruit Depot, Parris Island, South Carolina

Groundwater contamination by tetrachloroethene and its dechlorination products is present in two partially intermingled plumes in the surficial aquifer near a former dry-cleaning facility at Site 45, Marine Corps Recruit Depot, Parris Island, South Carolina. The northern plume originates from the vicinity of former above-ground storage tanks. Free-phase tetrachloroethene from activities in this area entered the groundwater and the storm sewer. The southern plume originates at a nearby new dry-cleaning facility, but probably was the result of contamination released to the aquifer from a leaking sanitary sewer line from the former dry-cleaning facility. Discharge of dissolved groundwater contamination is primarily to leaking storm sewers below the water table. Extensive biodegradation of the contamination takes place in the surficial aquifer; however, the biodegradation is insufficient to reduce trichloroethene to less than milligram-per-liter concentrations prior to discharging into the storm sewers. The groundwater volatile organic compounds entering the storm sewers are substantially diluted by tidal flushing upon entry and are subject to volatilization as they are transported through the storm sewer to a discharge point in a tributary to Ballast Creek. TCE concentrations of about 2-6 micrograms per liter were present in storm-sewer water near the discharge point (sampled at manhole STS26). On three out of four sampling events at manhole STS14, the storm-sewer water contained no vinyl chloride. During a time of relatively high groundwater levels, however, 20 micrograms per liter of vinyl chloride was present in STS14 storm-sewer water. Because groundwater leaks into that storm sewer and because the storm sewer upgradient from manhole STS14 is adjacent to part of the aquifer where 2,290 micrograms per liter of vinyl chloride have been detected, there is a potential for substantially increased concentrations of vinyl chloride to discharge at the storm-sewer outfall under conditions of high groundwater levels and low tidal flushing. In addition, the observation that free-phase tetrachloroethene may have entered the storm-sewer system during the 1994 discharge means that dense nonaqueous phase liquid tetrachloroethene could have leaked from various parts of the storm sewer or discharged to surface water at the storm-sewer outfall.

South Carolina↗

USGS2011

Investigation of the potential source area, contamination pathway, and probable release history of chlorinated-solvent-contaminated groundwater at the Capital City Plume Site, Montgomery, Alabama, 2008-2010

Detection of the organic solvent perchloroethylene (PCE) in a shallow public-supply well in 1991 and exposure of workers in 1993 to solvent vapors during excavation activities to depths near the water table provided evidence that the shallow aquifer beneath the capital city of Montgomery, Alabama, was contaminated. Investigations conducted from 1993 to 1999 by State and Federal agencies confirmed the detection of PCE in the shallow aquifer, as well as the detection of the organic solvent trichloroethylene (TCE) and various inorganic compounds, but the source of the groundwater contamination was not determined. In May 2000 the U.S. Environmental Protection Agency proposed that the site, called the Capital City Plume (CCP) Site, be a candidate for the National Priorities List. Between 2000 and 2007, numerous site-investigation activities also did not determine the source of the groundwater contamination. In 2008, additional assessments were conducted at the CCP Site to investigate the potential source area, contamination pathway, and the probable release history of the chlorinated-solvent-contaminated groundwater. The assessments included the collection of (1) pore water in 2008 from the hyporheic zone of a creek using passive-diffusion bag samplers; (2) tissue samples in 2008 and 2009 from trees growing in areas of downtown Montgomery characterized by groundwater contamination and from trees growing in riparian zones along the Alabama River and Cypress Creek; and (3) groundwater samples in 2009 and 2010. The data collected were used to investigate the potential source area of contaminants detected in groundwater, the pathway of groundwater contamination, and constraints on the probable contaminant-release history. The data collected between 2008 and 2010 indicate that the PCE and TCE contamination of the shallow aquifer beneath the CCP Site most likely resulted from the past use and disposal of industrial wastewater from printing operations containing chlorinated solvents into the sanitary sewer and (or) stormwater systems of Montgomery. Moreover, chlorinated-solvent use and disposal occurred at least between the 1940s and 1970s at several locations occupied by printing operations. The data also indicate that PCE and TCE contamination continues to occur in the shallow subsurface near potential release areas and that PCE and TCE have been transported to the intermediate part of the shallow aquifer.

Alabama↗

USGS2013

Groundwater quality and the relation between pH values and occurrence of trace elements and radionuclides in water samples collected from private wells in part of the Kickapoo Tribe of Oklahoma Jurisdictional Area, central Oklahoma, 2011

From 1999 to 2007, the Indian Health Service reported that gross alpha-particle activities and concentrations of uranium exceeded the Maximum Contaminant Levels for public drinking-water supplies in water samples from six private wells and two test wells in a rural residential neighborhood in the Kickapoo Tribe of Oklahoma Jurisdictional Area, in central Oklahoma. Residents in this rural area use groundwater from Quaternary-aged terrace deposits and the Permian-aged Garber-Wellington aquifer for domestic purposes. Uranium and other trace elements, specifically arsenic, chromium, and selenium, occur naturally in rocks composing the Garber-Wellington aquifer and in low concentrations in groundwater throughout its extent. Previous studies have shown that pH values above 8.0 from cation-exchange processes in the aquifer cause selected metals such as arsenic, chromium, selenium, and uranium to desorb (if present) from mineral surfaces and become mobile in water. On the basis of this information, the U.S. Geological Survey, in cooperation with the Kickapoo Tribe of Oklahoma, conducted a study in 2011 to describe the occurrence of selected trace elements and radionuclides in groundwater and to determine if pH could be used as a surrogate for laboratory analysis to quickly and inexpensively identify wells that might contain high concentrations of uranium and other trace elements. The pH and specific conductance of groundwater from 59 private wells were measured in the field in an area of about 18 square miles in Lincoln and Pottawatomie Counties. Twenty of the 59 wells also were sampled for dissolved concentrations of major ions, trace elements, gross alpha-particle and gross beta-particle activities, uranium, radium-226, radium-228, and radon-222 gas. Arsenic concentrations exceeded the Maximum Contaminant Level of 10 micrograms per liter in one sample having a concentration of 24.7 micrograms per liter. Selenium concentrations exceeded the Maximum Contaminant Level of 50 micrograms per liter in one sample having a concentration of 147 micrograms per liter. Both samples had alkaline pH values, 8.0 and 8.4, respectively. Uranium concentrations ranged from 0.02 to 383 micrograms per liter with 5 of 20 samples exceeding the Maximum Contaminant Level of 30 micrograms per liter; the five wells with uranium concentrations exceeding 30 micrograms per liter had pH values ranging from 8.0 to 8.5. Concentrations of uranium and radon-222 and gross alpha-particle activity showed a positive relation to pH, with the highest concentrations and activity in samples having pH values of 8.0 or above. The groundwater samples contained dissolved oxygen and high concentrations of bicarbonate; these characteristics are also factors in increasing uranium solubility. Concentrations of radium-226 and radium-228 (combined) ranged from 0.03 to 1.7 picocuries per liter, with a median concentration of 0.45 picocuries per liter for all samples. Radon-222 concentrations ranged from 95 to 3,600 picocuries per liter with a median concentration of 261 picocuries per liter. Eight samples having pH values ranging from 8.0 to 8.7 exceeded the proposed Maximum Contaminant Level of 300 picocuries per liter for radon-222. Eight samples exceeded the 15 picocuries per liter Maximum Contaminant Level for gross alpha-particle activity at 72 hours (after sample collection) and at 30 days (after the initial count); those samples had pH values ranging from 8.0 to 8.5. Gross beta-particle activity increased in 15 of 21 samples during the interval from 72 hours to 30 days. The increase in gross beta-particle activity over time probably was caused by the ingrowth and decay of uranium daughter products that emit beta particles. Water-quality data collected for this study indicate that pH values above 8.0 are associated with potentially high concentrations of uranium and radon-222 and high gross alpha-particle activity in the study area. High pH values also are associated with potentially high concentrations of arsenic, chromium, and selenium in groundwater when these elements occur in the aquifer matrix along groundwater-flow paths.

Oklahoma↗

USGS2013

Arsenic concentrations, related environmental factors, and the predicted probability of elevated arsenic in groundwater in Pennsylvania

Analytical results for arsenic in water samples from 5,023 wells obtained during 1969–2007 across Pennsylvania were compiled and related to other associated groundwater-quality and environmental factors and used to predict the probability of elevated arsenic concentrations, defined as greater than or equal to 4.0 micrograms per liter (µg/L), in groundwater. Arsenic concentrations of 4.0 µg/L or greater (elevated concentrations) were detected in 18 percent of samples across Pennsylvania; 8 percent of samples had concentrations that equaled or exceeded the U.S. Environmental Protection Agency’s drinking-water maximum contaminant level of 10.0 µg/L. The highest arsenic concentration was 490.0 µg/L. Comparison of arsenic concentrations in Pennsylvania groundwater by physiographic province indicates that the Central Lowland physiographic province had the highest median arsenic concentration (4.5 µg/L) and the highest percentage of sample records with arsenic concentrations greater than or equal to 4.0 µg/L (59 percent) and greater than or equal to 10.0 µg/L (43 percent). Evaluation of four major aquifer types (carbonate, crystalline, siliciclastic, and surficial) in Pennsylvania showed that all types had median arsenic concentrations less than 4.0 µg/L, and the highest arsenic concentration (490.0 µg/L) was in a siliciclastic aquifer. The siliciclastic and surficial aquifers had the highest percentage of sample records with arsenic concentrations greater than or equal to 4.0 µg/L and 10.0 µg/L. Elevated arsenic concentrations were associated with low pH (less than or equal to 4.0), high pH (greater than or equal to 8.0), or reducing conditions. For waters classified as anoxic (405 samples), 20 percent of sampled wells contained water with elevated concentrations of arsenic; for waters classified as oxic (1,530 samples) only 10 percent of sampled wells contained water with elevated arsenic concentrations. Nevertheless, regardless of the reduction-oxidation classification, 54 percent of samples with low pH (13 of 24 samples) and 25 percent of samples with high pH (57 of 230 samples) had elevated arsenic concentrations. Arsenic concentrations in groundwater in Pennsylvania were correlated with concentrations of several chemical constituents or properties, including (1) constituents associated with redox processes, (2) constituents that may have a similar origin or be mobilized under similar chemical conditions as arsenic, and (3) anions or oxyanions that have similar sorption behavior or compete for sorption sites on iron oxides. Logistic regression models were created to predict and map the probability of elevated arsenic concentrations in groundwater statewide in Pennsylvania and in three intrastate regions to further improve predictions for those three regions (glacial aquifer system, Gettysburg Basin, Newark Basin). Although the Pennsylvania and regional predictive models retained some different variables, they have common characteristics that can be grouped by (1) geologic and soils variables describing arsenic sources and mobilizers, (2) geochemical variables describing the geochemical environment of the groundwater, and (3) locally specific variables that are unique to each of the three regions studied and not applicable to statewide analysis. Maps of Pennsylvania and the three intrastate regions were produced that illustrate that areas most at risk are those with geology and soils capable of functioning as an arsenic source or mobilizer and geochemical groundwater conditions able to facilitate redox reactions. The models have limitations because they may not characterize areas that have localized controls on arsenic mobility. The probability maps associated with this report are intended for regional-scale use and may not be accurate for use at the field scale or when considering individual wells.

Pennsylvania↗

USGS2015

Hydrologic budget and conditions of Permian, Pennsylvanian, and Mississippian aquifers in the Appalachian Plateaus physiographic province

In response to challenges to groundwater availability posed by historic land-use practices, expanding development of hydrocarbon resources, and drought, the U.S. Geological Survey Groundwater Resources Program began a regional assessment of the Appalachian Plateaus aquifers in 2013 that incorporated a hydrologic landscape approach to estimate all components of the hydrologic system: surface runoff, base flow from groundwater, and interaction with atmospheric water (precipitation and evapotranspiration). This assessment was intended to complement other Federal and State investigations and provide foundational groundwater-related datasets in the Appalachian Plateaus. A regional Soil-Water-Balance model was constructed for a 160,000-square-mile study area that extended to the topographic divide of all streams originating outside but flowing into areas underlain by Appalachian Plateaus aquifers. The model incorporated soil, landscape, and climate variables to estimate an annual water budget for the 32-year period from 1980 to 2011 and was calibrated using base-flow data estimated by hydrograph separation techniques from 20 streamflow gaging stations across the study area. Over this period, an average of 47 inches per year (in/yr) of precipitation fell on Appalachian Plateaus aquifers. Simulations from the regional Soil-Water-Balance model indicate that only 19 percent of the precipitation or an average 9 in/yr recharged aquifers, and 19 percent resulted in surface runoff to streams. The remaining 62 percent, an average of 27 in/yr of water, was returned to the atmosphere via evapotranspiration. Because withdrawals from aquifers due to pumping equated to less than 1 percent of the water budget, differences in predevelopment and postdevelopment regional water budgets of the Appalachian Plateaus were minimal. Storage changes caused by filling of abandoned coal-mine aquifers and long-term differences in aquifer storage resulting from climate fluctuations constitute a small portion of the overall water budget. The percentage of precipitation that results in recharge, runoff, or evapotranspiration from the landscape varies annually by up to a factor of two depending on temporal changes in prevailing climate conditions and spatial changes in basin characteristics, precipitation patterns, and sources of atmospheric moisture over a large study area. A comparison of water-budget estimates from the regional Soil-Water-Balance model for a dry year (1988) and wet year (2004) showed that evapotranspiration accounts for most of the annual differences in precipitation. As a portion of annual precipitation, evapotranspiration ranged from 69 percent (dry year) to 52 percent (wet year), a range four times greater than the 15 percent (dry year) to 18 percent (wet year) range estimated for recharge. Evapotranspiration as a percentage of precipitation peaks during dry periods, whereas base flow and runoff tend to reach minimum values. During wet periods, this relationship is reversed and base flow and runoff as a percentage of precipitation generally peak while evapotranspiration percentages reach minimum values. Annual recharge in the Appalachian Plateaus reaches a maximum at near 20 percent of annual precipitation, regardless of the severity of wet conditions. Hydrograph separation data from 849 streamflow gaging stations in the study area were used to assess trends in streamflow, base flow, surface runoff, and base-flow index, or ratio of base flow to streamflow, in the Appalachian Plateaus for the period from 1930 to 2011. Annual data anomalies for each of the four variables were individually defined as the annual standard deviation from the mean at all 849 streamflow gaging stations. Annual data anomalies confirm the close relation of annual precipitation to both base flow and runoff components of streamflow, and both components increased during the period of analysis. Around 1970, conditions shifted streamflow from values generally below to above long-term means. At a regional scale, increases in base flow account for most of these observed increases in mean annual streamflow. The independence of the base-flow index to annual climate trends indicate that changes in the components of streamflow of the Appalachian Plateaus are probably in response to shifts in seasonal precipitation or widespread land-use practices. A subset of 77 index streamgages, defined as having 60 or more years of complete record between the years 1930 and 2011 with no more than 20 percent missing data, was selected to show spatial patterns of change in the water budget. Data from the index streamgages showed that the overall trends in base flow are dependent upon the period of evaluation. Long-term (1930–2011) increases in base flow were observed throughout the study area. For two shorter periods (1930–1969 and 1970–2011) trends in base flow were largely negative. In general, spatial patterns of change in streamflow, base flow, and runoff were mixed but generally consistent with prevailing climate patterns and land-use changes.

Alabama, Kentucky, Maryland, Ohio, Pennslyvania, V↗

USGS1988

Ground-water resources and simulation of flow in aquifers containing freshwater and seawater, Island County, Washington

Aquifers in Island County, Washington, that are intruded by seawater from Puget Sound contain chloride concentrations that exceed 100 mg/L. Chloride concentrations exceeded 100 mg/L in 24 % of the wells that were drilled below sea level and sampled in August 1981, but most of the chloride concentrations did not exceed 1,000 mg/L. Groundwater occurs in glacial deposits that have a maximum thickness of 3,000 ft; the deposits were divided into five aquifers and five confining units. Four overlapping digital models were calibrated, using time-averaged data, to simulate three-dimensional steady flow of fresh groundwater in multiple aquifers containing freshwater and seawater separated by a sharp interface. Model simulations indicate that most of the recharge is discharged from aquifers C and D as springs below sea levels, and only a small fraction of the recharge moves downward below aquifer C. Simulations also indicate that aquifers beneath the islands are not recharged by groundwater that moves from the mainland through aquifers beneath Puget Sound except in the area of northeast Camano Island. Between Whidbey and Camano Islands, the freshwater-seawater interface intersects the bottom of Puget Sound and prevents movement of fresh groundwater between the two islands. (USGS)

Water-Resources Investigations Report↗

USGS1988

Hydrogeology and water-supply potential of the water-table aquifer on Dauphin Island, Alabama

The water table aquifer on Dauphin Island, Alabama, consists of a thin veneer of Holocene sand and an underlying Pleistocene unit locally known as the Gulfport Formation. The aquifer is from 28 to 35 ft thick with a thick marine clay at its base. Water in the aquifer generally is low in chloride content except near the coast. Excessively high iron concentrations in groundwater were found locally. A two-dimensional finite-difference groundwater flow model of the water table aquifer on Dauphin Island was used in the steady-state mode to evaluate the flow system under steady-state conditions. Model input data were obtained primarily from 40 test wells, 2 aquifer tests, continuous recording of groundwater levels, and rainfall. The model was calibrated to the low water-table conditions of July 1985 and high water table conditions of April 1985. The model was also used to simulate pumpage from the aquifer under transient conditions with no rainfall. Patterns of computed head changes compared favorably to the natural recession of water levels for the periods of April to May 1985 and May to June 1985. Simulation of groundwater withdrawals in the transient model showed the feasibility of producing 0.6 million gallons/day from eight wells that tap the water table aquifer without inducing lateral seawater encroachment. (USGS)

Water-Resources Investigations Report↗

USGS1989

Analysis of the effect of pumping on ground-water flow in the Springfield Plateau and Ozark aquifers near Springfield, Missouri

Pumpage of water from the Ozark aquifer for public supply and industry use by the city of Springfield and surrounding communities in southwestern Missouri has significantly altered the potentiometric surface of the aquifer. Springfield is located on a regional groundwater divide that trends east and west across southern Missouri. Groundwater that once flowed north and south from the divide now moves toward Springfield. Drawdown in the Ozark aquifer beneath Springfield has increased about 50 ft near the center of the city since 1974. The area of well influence also has increased, most notably to the south and southwest, because of increased pumpage by Springfield and new groundwater withdrawals in rapidly increasing communities, such as Republic and Nixa. Changes in the potentiometric surface of the Ozark aquifer, and to a lesser extent the Springfield Plateau aquifer, resulting from stresses applied by pumpage of water supply wells has altered the hydrologic budget of the Springfield area. Downward leakage of groundwater through the Ozark confining unit has increased from about 10 cu ft/sec to about 18 cu ft/sec because drawdown in the Ozark aquifer has resulted in an increased vertical hydraulic gradient across the confining unit. Minimal quantities of water are supplied by increased upward leakage through the St. Francois confining unit. Model simulations indicate substantial quantities of water are still (1987) being removed from storage within the Ozark and Springfield Plateau aquifers and the hydrologic system is not in equilibrium at this time. (USGS)

Water-Resources Investigations Report↗

USGS2013

Pathways of coupled arsenic and iron cycling in high arsenic groundwater of the Hetao basin, Inner Mongolia, China: an iron isotope approach

High As groundwater is widely distributed all over the world, which has posed a significant health impact on millions of people. Iron isotopes have recently been used to characterize Fe cycling in aqueous environments, but there is no information on Fe isotope characteristics in the groundwater. Since groundwater As behavior is closely associated with Fe cycling in the aquifers, Fe isotope signatures may help to characterize geochemical processes controlling As concentrations of shallow groundwaters. This study provides the first observation of Fe isotope fractionation in high As groundwater and evaluation of Fe cycling and As behaviors in shallow aquifers in terms of Fe isotope signatures. Thirty groundwater samples were taken for chemical and isotopic analysis in the Hetao basin, Inner Mongolia. Thirty-two sediments were sampled as well from shallow aquifers for Fe isotope analysis. Results showed that groundwater was normally enriched in isotopically light Fe with δ 56 Fe values between −3.40&permil; and 0.58&permil; and median of −1.14&permil;, while heavier δ 56 Fe values were observed in the sediments (between −1.10&permil; and 0.75&permil;, median +0.36&permil;). In reducing conditions, groundwaters generally had higher δ 56 Fe values, in comparison with oxic conditions. High As groundwaters, generally occurring in reducing conditions, had high δ 56 Fe values, while low As groundwaters normally had low δ 56 Fe values. Although sediment δ 56 Fe values were generally independent of lithological conditions, a large variation in sediment δ 56 Fe values was observed in the oxidation–reduction transition zone. Three pathways were identified for Fe cycling in shallow groundwater, including dissimilatory reduction of Fe(III) oxides, re-adsorption of Fe(II), and precipitation of pyrite and siderite. Dissimilatory reduction of Fe(III) oxides resulted in light δ 56 Fe values (around −1.0&permil;) and high As concentration (>50 μg/L) in groundwater in anoxic conditions. Re-adsorption of isotopically heavy Fe(II) produced by microbially mediated reduction of Fe(III) oxides led to further enrichment of isotopically light Fe in groundwater (up to −3.4&permil; of δ 56 Fe) in anoxic–suboxic conditions. Arsenic re-adsorption was expected to occur along with Fe(II) re-adsorption, decreasing groundwater As concentrations. In strongly reducing conditions, precipitation of isotopically light Fe-pyrite and/or siderite increased groundwater δ 56 Fe values, reaching +0.58&permil; δ 56 Fe, with a subsequent decrease in As concentrations via co-precipitation. The mixed effect of those pathways would regulate As and Fe cycling in most groundwaters.

Inner Mongolia↗

USGS2023

Evaluation of stream capture related to groundwater pumping, Lower Humboldt River Basin, Nevada

The Humboldt River Basin is the only river basin that is contained entirely within the State of Nevada. The effect of groundwater pumping on the Humboldt River is not well understood. Tools are needed to determine stream capture and manage groundwater pumping in the Humboldt River Basin. The objective of this study is to estimate capture and storage change caused by groundwater withdrawals in the lower Humboldt River Basin that can provide the Nevada State Engineer with data and information needed to manage groundwater and surface-water resources. A numerical groundwater flow model was developed for the purpose of estimating stream capture from pre-2016 and future pumping as well as for any location of potential future pumping within the lower Humboldt River Basin. This model was developed using MODFLOW-NWT to represent the lower Humboldt River Basin hydrologic system, including Humboldt River; Rye Patch Reservoir; groundwater evapotranspiration; pumping from municipal, agricultural, mining, and domestic wells; and agricultural drains. Aquifer properties were calibrated using results from numerous single- and multi-well aquifer tests (Nadler, 2020) and through the process of model calibration. Historical capture was estimated for 1960–2016 and predictive capture for the system was projected 100 years into the future (2017–2116) based on historical pumping patterns. Stream capture and drain capture are relatively low for the historical and predictive periods. During the historical period, increased pumping during dry years caused increased connections with capture sources and less water sourced to wells from aquifer storage. Storage and groundwater levels generally recovered during subsequent wet years. Overall, storage change has been the main source of water to wells in the lower Humboldt River Basin, followed by groundwater evapotranspiration capture. During the predictive period, pumping is projected to remain constant and capture 9 percent of stream water after 100 years. Capture and storage change maps were created to visualize spatial variability in potential capture and storage change through time and to provide a database of results that can be used to manage groundwater and surface-water resources. These maps show that potential stream capture would be a minor source of water to wells located across most of the simulated area, except for locations close to the Humboldt River and Rye Patch Reservoir. Drains also would be a minor potential source of water to wells except for those directly adjacent to the drains. In general, the potential supply of water to wells is storage-dominated and over time groundwater evapotranspiration-dominated in the agricultural area. Capture difference maps were generated to visualize where potential capture results might have greater limitations associated with nonlinear flow processes, such as head-dependent boundary conditions. Higher capture differences indicate larger capture map bias and therefore greater capture map uncertainty due to the inability of capture maps to account for nonlinear flow processes. Stream capture differences are highest directly adjacent to the river but are otherwise minimal. Drain capture differences are highest in the region of the agricultural drain network but are otherwise minimal. The Humboldt River, Rye Patch Reservoir, and drains introduce very little nonlinearity to the model, and their associated capture map bias is minimal. Potential groundwater evapotranspiration capture introduces a fair amount of nonlinearity to the model and has the potential to result in significant, localized groundwater evapotranspiration capture map bias over time. Groundwater evapotranspiration capture differences are the result of higher pumping rates lowering the water table below the root zone faster than lower pumping rates and essentially removing groundwater evapotranspiration as a potential source of capture faster than lower pumping rates. Wells that can no longer source their supply through groundwater evapotranspiration capture then generally source more of their water from storage. Thus, storage change bias increases over time as well. Capture prediction uncertainty due to parameter estimation was evaluated using a covariance matrix adaptation-evolution strategy. One hundred Monte Carlo realizations of model parameters were applied to the model to assess capture uncertainty at 13 grid cell locations within the model domain. In general, results indicated that greater capture uncertainty for a given source (river, drains, or evapotranspiration) is associated with proximity of a pumping well to that source. The magnitude of maximum capture fraction uncertainties after 100 years of pumping for stream capture, drain capture, groundwater evapotranspiration capture, and storage change were plus or minus (±) 0.17, ±0.10, ±0.20, and ±0.22, respectively.

Nevada↗

USGS2020

Farmer behavior under groundwater management scenarios: Implications for groundwater conservation in the Mississippi Alluvial Plain

Concern about sustained availability of fresh groundwater for agricultural use in the Mississippi Alluvial Plain (MAP) mounts as groundwater levels decline. We evaluate elasticities of demand for groundwater and other agricultural inputs, as well as overall and output specific economies of scale for four major irrigated commodities (rice, corn, soybeans, and cotton) in the MAP region. Additionally, we investigate impacts of two groundwater management policy scenarios, including increasing pumping cost and groundwater use restrictions, on irrigation behavior. The results show price elasticity of demand for groundwater to be -0.13, indicating that it is inelastic, and an increasing cost of pumping will not significantly decrease relative demand for groundwater in the region. Even with policy scenarios that either increase the costs of pumping significantly or restrict groundwater use in the region, groundwater demand still appears to be inelastic. We also document significant overall economies of scale in the region. Our findings have implications for potential policy options aimed at reducing groundwater use. Efficient management practices are important to increase aquifer recharge, and at the same time, incorporation of human behavior via economic analysis will improve projections of groundwater availability in the MAP region.

Missouri, Tennessee, Arkansas, Louisiana, Mississi↗

USGS2019

Comparison of groundwater-model construction methods, representations of glacial geology, model designs, and groundwater-model flow simulations within Elkhart County, Indiana

Automated data-processing methods allow hydrologists to efficiently incorporate digital well-record datasets into the construction of hydrostratigraphic frameworks for groundwater-flow models. The method selected to construct the hydrostratigraphic framework can affect the extent of geologic heterogeneity that can be included in the model. The detail generated from a hydrostratigraphic framework can affect groundwater simulation results. The effects of detail on model accuracy, groundwater-flow simulations, and particle-tracking simulations are described in this study. This report compares differences in hydrostratigraphic frameworks and results of groundwater models using (1) a method that incorporates more hydrologic judgment at the expense of using limited lithologic data and (2) a method that is more automated and uses all available lithologic data. The study additionally evaluates the effect of model discretization and inclusion of more (or less) geologic detail on simulation results. Two methods were used to create hydrostratigraphic frameworks of glacial deposits in the St. Joseph River Basin. One method, referred to as the subjective method, manually identifies stratigraphic boundaries using a sample of well logs from State databases and uses two-dimensional kriging to create three model layers of the study area. Indicator kriging is used to define aquifer extent in each layer. The second method, referred to as the objective method, uses three-dimensional kriging to automatically create a detailed heterogeneous model of the study area using all wells logs from the State database. The objective method increases detail in the vertical by greatly increasing the number of computer groundwater model layers from 3 to 30. In Elkhart County, Indiana, a previously published model represents the product of the subjective method, and a newly calibrated model of the same area represents the product of the objective method. An automated calibration procedure was used with the objective model (derived from the objective method) for Elkhart County. The two most-sensitive parameters for the Elkhart County objective model are horizontal hydraulic conductivity of the sand and the combined sand and gravel/gravel deposits. Vertical hydraulic conductivity of the fine-grained and intermediate-sized deposits could not be estimated, possibly indicating major flow paths are along a continuously connected series of sand and gravel deposits and not through a confining layer. The statistics measuring model calibration accuracy for the objective model were slightly better than statistics for the subjective model (model derived from the subjective method) of Elkhart County, but the hydraulic conductivities and flow rates for the two models were different. The mean absolute errors between simulated and measured groundwater levels are 2.04 and 2.16 feet for the objective and subjective models, respectively. Simulated seepage losses from and groundwater discharges to measured stream reaches in the objective model were evenly balanced in terms of over and under simulations of measured values; the subjective model tended to overpredict measured groundwater discharge to streams. The overprediction may be related to the 58 percent greater total inflow and outflow through the subjective model. The greater flow rate through the subjective model results from higher horizontal hydraulic conductivities in the subjective model than in the objective model. Horizontal hydraulic conductivity ranged from 23.9 to 111 feet per day in the objective model and generally ranged from 170 to 370 feet per day in the subjective model. The improvement in calibration statistics for the objective model relative to the subjective model may be from increased detail in how the objective model represents the distribution of fine- and coarse-grained deposits. The improvement also could be associated with the difference in methods used to represent the continuity of the confining unit. The effect of differences in horizontal hydraulic conductivity distributions between the two models for Elkhart County is evident in the groundwater-flow paths simulated by the objective and subjective models. At a withdrawal well location, the flow lines produced by the objective model indicate a wider contributing area than that for the subjective model. The discontinuous confining unit represented in the objective model provided the opportunity for groundwater flow to split into an upper and lower path. The split in flow simulated by the objective model at one location was independently supported by bromide concentrations in groundwater; the subjective model did not duplicate the split in flow.

Indiana↗

USGS2024

Estimating groundwater level records using MOVE.1 and computing monthly percentiles from estimated groundwater records in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, performed record extensions on groundwater levels at select wells using the Maintenance of Variance Extension type 1 (MOVE.1) method. The groundwater levels estimated from these record extensions were used to compute monthly percentiles to improve future determinations of a groundwater index. In Massachusetts, 27 of 29 short-record study wells with continuous groundwater levels between 0.8 and 8.1 years were suitable for record extensions; 37 long-record index wells were used to extend the groundwater level records at the study wells. The index well selected to pair with a study well was chosen based on Pearson correlation coefficient values; cross-correlation between the two wells; geologic and topographic similarity; and smallest distance spanning the wells. Each study well and its corresponding index well have 1 or more years of concurrent, overlapping data; a Pearson correlation coefficient that exceeded a threshold value of 0.8; and a similar aquifer type and hydrologic characteristics. Of the 29 study wells, 2 showed poor correlations with all index wells and were not considered for record extensions. Performance metrics used to assess the accuracy of the MOVE.1 models indicated that most models provided reasonable estimates of groundwater levels. Root mean square error values ranged from 0.097 to 2.292 feet, with a median of 0.536 foot. Nash-Sutcliffe efficiency coefficient values ranged from 0.623 to 0.996, with a median value of 0.759. Generally, study wells in close geographical proximity to their index well resulted in stronger model performance. The average length of groundwater level records was extended by 14.1 years to a new average of 18.1 years. The estimated groundwater level records from the MOVE.1 models resulted in an increase in the range of highest and lowest groundwater levels at 23 of 27 wells. The increase in range of groundwater levels was between 0.08 to 7.95 feet. Monthly percentiles for State drought indices were computed from the estimated MOVE.1 records and observed records through December 31, 2021. Percentiles computed from estimated records show an average groundwater level about 1.0 foot lower than observed data at the 2d percentile and 0.1 foot lower at the 30th percentile.

Massachusetts↗

USGS2012

Effects of groundwater withdrawals associated with combined-cycle combustion turbine plants in west Tennessee and northern Mississippi

The Mississippi Embayment Regional Aquifer Study groundwater-flow model was used to simulate the potential effects on future groundwater withdrawals at five powerplant sites-Gleason, Weakley County, Tennessee; Tenaska, Haywood County, Tennessee; Jackson, Madison County, Tennessee; Southaven, DeSoto County, Mississippi; and Magnolia, Benton County, Mississippi. The scenario used in the simulation consisted of a 30-year average water-use period followed by a 30-day peak water-demand period. Effects of the powerplants on the aquifer system were evaluated by comparing the difference in simulated water levels in the aquifers at the end of the scenario (30 years plus 30 days) with and without the combined-cycle-plant withdrawals. Simulated potentiometric surface declines in source aquifers at potential combined-cycle-plant sites ranged from 56 feet in the upper Wilcox aquifer at the Magnolia site to 20 feet in the Memphis aquifer at the Tenaska site. The affected areas in the source aquifers at the sites delineated by the 4-foot potentiometric surface-decline contour ranged from 11,362 acres at Jackson to 535,143 acres at Southaven. The extent of areas affected by potentiometric surface declines was similar at the Gleason and Magnolia sites. The affected area at the Tenaska site was smaller than the affected areas at the other sites, most likely as a result of lower withdrawal rates and greater aquifer thickness. The extent of effect was smallest at the Jackson site, where the nearby Middle Fork Forked Deer River may act as a recharge boundary. Additionally, the Jackson site lies in the Memphis aquifer outcrop area where model-simulated recharge rates are higher than in areas where the Memphis aquifer underlies less permeable deposits. The potentiometric surface decline in aquifers overlying or underlying a source aquifer was generally 2 feet or less at all the sites except Gleason. At the Gleason site, withdrawals from the Memphis aquifer resulted in declines of as much as 9 feet in the underlying Fort Pillow aquifer. The simulated potentiometric surface change occurring in the Fort Pillow aquifer appears to be the result of leakage through the Flour Island Formation separating the Memphis and Fort Pillow aquifers where this confining unit is thin, sandy, or absent.

Arkansas;Kentucky;Mississippi;Missouri;Tennessee↗