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.

SourceResearch connectionFrom yearThrough year

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

At least 649 records · Page 36Linked to original sources

USGS2012

Shallow groundwater quality and geochemistry in the Fayetteville Shale gas-production area, north-central Arkansas, 2011

The Mississippian Fayetteville Shale serves as an unconventional gas reservoir across north-central Arkansas, ranging in thickness from approximately 50 to 550 feet and varying in depth from approximately 1,500 to 6,500 feet below the ground surface. Primary permeability in the Fayetteville Shale is severely limited, and successful extraction of the gas reservoir is the result of advances in horizontal drilling techniques and hydraulic fracturing to enhance and develop secondary fracture porosity and permeability. Drilling and production of gas wells began in 2004, with a steady increase in production thereafter. As of April 2012, approximately 4,000 producing wells had been completed in the Fayetteville Shale. In Van Buren and Faulkner Counties, 127 domestic water wells were sampled and analyzed for major ions and trace metals, with a subset of the samples analyzed for methane and carbon isotopes to describe general water quality and geochemistry and to investigate the potential effects of gas-production activities on shallow groundwater in the study area. Water-quality analyses from this study were compared to historical (pregas development) shallow groundwater quality collected in the gas-production area. An additional comparison was made using analyses from this study of groundwater quality in similar geologic and topographic areas for well sites less than and greater than 2 miles from active gas-production wells. Chloride concentrations for the 127 groundwater samples collected for this study ranged from approximately 1.0 milligram per liter (mg/L) to 70 mg/L, with a median concentration of 3.7 mg/L, as compared to maximum and median concentrations for the historical data of 378 mg/L and 20 mg/L, respectively. Statistical analysis of the data sets revealed statistically larger chloride concentrations (p-value <0.001) in the historical data compared to data collected for this study. Chloride serves as an important indicator parameter based on its conservative transport characteristics and relatively elevated concentrations in production waters associated with gas extraction activities. Major ions and trace metals additionally had lower concentrations in data gathered for this study than in the historical analyses. Additionally, no statistical difference existed between chloride concentrations from water-quality data collected for this study from 94 wells located less than 2 miles from a gas-production well and 33 wells located 2 miles or more from a gas-production well; a Wilcoxon rank-sum test showed a p-value of 0.71. Major ion chemistry was investigated to understand the effects of geochemical and reduction-oxidation (redox) processes on the shallow groundwater in the study area along a continuum of increased rock-water interaction represented by increases in dissolved solids concentration. Groundwater in sandstone formations is represented by a low dissolved solids concentration (less than 30 mg/L) and slightly acidic water type. Shallow shale aquifers were represented by dissolved solids concentrations ranging upward to 686 mg/L, and water types evolving from a dominantly mixed-bicarbonate and calcium-bicarbonate to a strongly sodium-bicarbonate water type. Methane concentration and carbon isotopic composition were analyzed in 51 of the 127 samples collected for this study. Methane occurred above a detection limit of 0.0002 mg/L in 32 of the 51 samples, with concentrations ranging upward to 28.5 mg/L. Seven samples had methane concentrations greater than or equal to 0.5 mg/L. The carbon isotopic composition of these higher concentration samples, including the highest concentration of 28.5 mg/L, shows the methane was likely biogenic in origin with carbon isotope ratio values ranging from -57.6 to -74.7 per mil. Methane concentrations increased with increases in dissolved solids concentrations, indicating more strongly reducing conditions with increasing rock-water interaction in the aquifer. As such, groundwater-quality data collected for this study indicate that groundwater chemistry in the shallow aquifer system in the study area is a result of natural processes, beginning with recharge of dilute atmospheric precipitation and evolution of observed groundwater chemistry through rock-water interaction and redox processes.

Arkansas↗

USGS2013

Groundwater withdrawals 1976, 1990, and 2000--10 and land-surface-elevation changes 2000--10 in Harris, Galveston, Fort Bend, Montgomery, and Brazoria Counties, Texas

The study area comprising Harris County and parts of Galveston, Fort Bend, Montgomery, and Brazoria Counties in southeastern Texas forms part of one of the largest areas of land-surface-elevation change in the United States. Land-surface-elevation change in the study area primarily is caused by the withdrawal of groundwater. Groundwater withdrawn from the Chicot and Evangeline aquifers has been the primary source of water for municipal supply, industrial and commercial use, and irrigation in the study area. Groundwater withdrawals cause compaction of clay and silt layers abundant in the aquifers, which has in turn resulted in the widespread, substantial land-surface-elevation changes in the region with increased flooding. To estimate land-surface-elevation changes, the U.S. Geological Survey (USGS), in cooperation with the Harris-Galveston Subsidence District (HGSD), documented land-surface-elevation changes in the study area that occurred during 2000–10 and 2005–10 based on elevation data measured by 11 USGS borehole-extensometer sites, a National Geodetic Survey Continuously Operating Reference Station, and Global Positioning System Port-A-Measure (PAM) sites operated by the HGSD and the Fort Bend Subsidence District. Groundwater withdrawals in the study area also were documented for 1976, 1990, and 2000–10. In 1976, about 428.9 million gallons per day (Mgal/d) were withdrawn from the aquifer system in Harris County, but by 2000, because of HGSD regulation, withdrawals had decreased to about 337.8 Mgal/d, or about a 21-percent reduction since 1976. By 2010, withdrawals had decreased to about 227.1 Mgal/d, or about a 47-percent reduction since 1976. Among the counties in the study area, the largest decrease in groundwater withdrawals has occurred in Galveston County since 1976. In 1976, about 27.4 Mgal/d were withdrawn from the aquifer system, and by 2000, withdrawals had decreased to about 4.12 Mgal/d, or about an 85-percent reduction since 1976. By 2010, withdrawals had decreased to about 0.626 Mgal/d, or about a 98-percent decrease since 1976. Since the mid-1970s, Fort Bend and Montgomery Counties have undergone extensive urban development and corresponding large increases in groundwater withdrawals. Total groundwater withdrawal for Fort Bend County in 1976 was about 16.0 Mgal/d, and by 2000, withdrawals had increased to about 86.5 Mgal/d, or about a 441-percent increase since 1976. By 2010, withdrawals in Fort Bend County had increased to about 99.8 Mgal/d, or about a 524-percent increase since 1976. Total groundwater withdrawal for Montgomery County in 1976 was about 7.84 Mgal/d, and by 2000, withdrawals had increased to about 43.6 Mgal/d, or about a 456-percent increase since 1976. By 2010, withdrawals in Montgomery County had increased to about 64.2 Mgal/d, or about a 719-percent increase since 1976. Total groundwater withdrawal in Brazoria County in 1976 was about 18.0 Mgal/d, and by 2000, withdrawals had increased to about 26.0 Mgal/d, or about a 44-percent increase. By 2010, withdrawals in Brazoria County had increased to about 24.7 Mgal/d, or about a 37-percent increase since 1976. Measured land-surface-elevation changes from December 31, 2000, to December 31, 2010, ranged from an elevation increase of 0.06 feet (ft), or an average increase in elevation of 0.006 ft per year, at the Seabrook borehole extensometer located near Seabrook, Tex., to an elevation decrease of 1.28 ft, or an average decrease in elevation of 0.128 ft per year, at a PAM station north of Jersey Village, Tex. (PAM 07). Measured land-surface-elevation changes from December 31, 2005, to December 31, 2010, ranged from an elevation increase of 0.07 ft, or an average increase in elevation of 0.014 ft per year, at PAM 09 in far northeastern Harris County to an elevation decrease of 0.51 ft, or an average decrease in elevation of 0.102 ft per year, at PAM 07.

Texas↗

USGS2013

Occurrence of methane in groundwater of south-central New York State, 2012-systematic evaluation of a glaciated region by hydrogeologic setting

A survey of methane in groundwater was undertaken to document methane occurrence on the basis hydrogeologic setting within a glaciated 1,810-square-mile area of south-central New York along the Pennsylvania border. Sixty-six wells were sampled during the summer of 2012. All wells were at least 1 mile from any known gas well (active, exploratory, or abandoned). Results indicate strong positive and negative associations between hydrogeologic settings and methane occurrence. The hydrogeologic setting classes are based on topographic position (valley and upland), confinement or non-confinement of groundwater by glacial deposits, well completion in fractured bedrock or sand and gravel, and hydrogeologic subcategories. Only domestic wells and similar purposed supply wells with well-construction and log information were selected for classification. Field water-quality characteristics (pH, specific conductance, dissolved oxygen, and temperature) were measured at each well, and samples were collected and analyzed for dissolved gases, including methane and short-chain hydrocarbons. Carbon and hydrogen isotopic ratios of methane were measured in 21 samples that had at least 0.3 milligram per liter (mg/L) of methane. Results of sampling indicate that occurrence of methane in groundwater of the region is common—greater than or equal to 0.001 mg/L in 78 percent of the groundwater samples. Concentrations of methane ranged over five orders of magnitude. Methane concentrations at which monitoring or mitigation are indicated (greater than or equal to 10 mg/L) were measured in 15 percent of the samples. Methane concentrations greater than 0.1 mg/L were associated with specific hydrogeologic settings. Wells completed in bedrock within valleys and under confined groundwater conditions were most closely associated with the highest methane concentrations. Fifty-seven percent of valley wells had greater than or equal to 0.1 mg/L of methane, whereas only 10 percent of upland wells equaled or exceeded that concentration. Isotopic signatures differed between these groups as well. Methane in valley wells was predominantly thermogenic in origin, likely as a result of close vertical proximity to underlying methane-bearing saline groundwater and brine and possibly as a result of enhanced bedrock fracture permeability beneath valleys that provides an avenue for upward gas migration. Isotopic signatures of methane from four upland well samples indicated a microbial origin (carbon-dioxide reduction) with one sample possibly altered by microbial methane oxidation. Water samples from wells in a valley setting that indicate a mix of thermogenic and microbial methane reflect the close proximity of regional groundwater flow and underlying saline water and brine in valley areas. The microbial methane is likely produced by bacteria that utilize carbon dioxide or formational organic matter in highly reducing environments within the subregional groundwater flow system. This characterization of groundwater methane shows the importance of subsurface information (hydrogeology, well construction) in understanding methane occurrence and provides an initial conceptual framework that can be utilized in investigation of stray gas in south-central New York.

New York↗

USGS2014

Hydrogeology and simulation of groundwater flow in the Central Oklahoma (Garber-Wellington) Aquifer, Oklahoma, 1987 to 2009, and simulation of available water in storage, 2010–2059

The Central Oklahoma (Garber-Wellington) aquifer underlies about 3,000 square miles of central Oklahoma. The study area for this investigation was the extent of the Central Oklahoma aquifer. Water from the Central Oklahoma aquifer is used for public, industrial, commercial, agricultural, and domestic supply. With the exception of Oklahoma City, all of the major communities in central Oklahoma rely either solely or partly on groundwater from this aquifer. The Oklahoma City metropolitan area, incorporating parts of Canadian, Cleveland, Grady, Lincoln, Logan, McClain, and Oklahoma Counties, has a population of approximately 1.2 million people. As areas are developed for groundwater supply, increased groundwater withdrawals may result in decreases in long-term aquifer storage. The U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, investigated the hydrogeology and simulated groundwater flow in the aquifer using a numerical groundwater-flow model. The purpose of this report is to describe an investigation of the Central Oklahoma aquifer that included analyses of the hydrogeology, hydrogeologic framework of the aquifer, and construction of a numerical groundwater-flow model. The groundwater-flow model was used to simulate groundwater levels and for water-budget analysis. A calibrated transient model was used to evaluate changes in groundwater storage associated with increased future water demands.

Oklahoma↗

USGS2014

Simulation of groundwater flow and saltwater movement in the Onslow County area, North Carolina: predevelopment-2010

Onslow County, North Carolina, is located within the designated Central Coastal Plain Capacity Use Area (CCPCUA). The CCPCUA was designated by law as a result of groundwater level declines of as much as 200 feet during the past four decades within aquifers in rocks of Cretaceous age in the central Coastal Plain of North Carolina and a depletion of water in storage from increased groundwater withdrawals in the area. The declines and depletion of water in storage within the Cretaceous aquifers increase the potential for saltwater migration—both lateral encroachment and upward leakage of brackish water. Within the CCPCUA, a reduction in groundwater withdrawals over a period of 16 years from 2003 to 2018 is mandated. Under the CCPCUA rules, withdrawals in excess of 100,000 gallons per day from any of the Cretaceous aquifer well systems are subject to water-use reductions of as much as 75 percent. To assess the effects of the CCPCUA rules and to assist with groundwater-management decisions, a numerical model was developed to simulate the groundwater flow and chloride concentrations in the surficial Castle Hayne, Beaufort, Peedee, and Black Creek aquifers in the Onslow County area. The model was used to (1) simulate groundwater flow from 1900 to 2010; (2) assess chloride movement throughout the aquifer system; and (3) create hypothetical scenarios of future groundwater development. After calibration of a groundwater flow model and conversion to a variable-density model, five scenarios were created to simulate future groundwater conditions in the Onslow County area: (1) full implementation of the CCPCUA rules with three phases of withdrawal reductions simulated through 2028; (2) implementation of only phase 1 withdrawal reductions of the CCPCUA rules and simulated through 2028; (3) implementation of only phases 1 and 2 withdrawal reductions of the CCPCUA rules and simulated through 2028; (4) full implementation of the CCPCUA rules with the addition of withdrawals from the Castle Hayne aquifer in Onslow County at the fully permitted amount in the final stress period and simulated through 2028; and (5) full implementation of the CCPCUA rules as in scenario 1 except simulated through 2100. Results from the scenarios give an indication of the water-level recovery in the Black Creek aquifer throughout each phase of the CCPCUA rules in Onslow County. Furthermore, as development of the Castle Hayne aquifers was increased in the scenarios, cones of depression were created around pumping centers. Additionally, the scenarios indicated little to no change in chloride concentrations for the time periods simulated.

North Carolina↗

USGS2014

Status of groundwater quality in the Borrego Valley, Central Desert, and Low-Use Basins of the Mojave and Sonoran Deserts study unit, 2008-2010: California GAMA Priority Basin Project

Groundwater quality in the approximately 963-square-mile Borrego Valley, Central Desert, and Low-Use Basins of the Mojave and Sonoran Deserts study unit was investigated as part of the Priority Basin Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is located in southern California in San Bernardino, Riverside, San Diego, and Imperial Counties. The GAMA Priority Basin Project is being conducted by the California State Water Resources Control Board in collaboration with the U.S. Geological Survey and the Lawrence Livermore National Laboratory. The GAMA Borrego Valley, Central Desert, and Low-Use Basins of the Mojave and Sonoran Deserts study was designed to provide a spatially unbiased assessment of the quality of untreated (raw) groundwater in the primary aquifer system. The assessment is based on water-quality and ancillary data collected by the U.S. Geological Survey from 52 wells (49 grid wells and 3 understanding wells) and on water-quality data from the California Department of Public Health database. The primary aquifer system was defined by the depth intervals of the wells listed in the California Department of Public Health database for the Borrego Valley, Central Desert, and Low-Use Basins of the Mojave and Sonoran Deserts study unit. The quality of groundwater in the primary aquifer system may be different from that in the shallower or deeper water-bearing zones; shallow groundwater may be more vulnerable to surficial contamination. This study assesses the status of the current quality of the groundwater resource by using data from samples analyzed for volatile organic compounds (VOCs), pesticides, and naturally occurring inorganic constituents, such as major ions and trace elements. This status assessment is intended to characterize the quality of groundwater resources in the primary aquifer system of the Borrego Valley, Central Desert, and Low-Use Basins of the Mojave and Sonoran Deserts study unit, not the treated drinking water delivered to consumers by water purveyors. Relative-concentrations (sample concentration divided by the health- or aesthetic-based benchmark concentration) were used for evaluating groundwater quality for those constituents that have Federal or California regulatory or non-regulatory benchmarks for drinking-water quality. A relative-concentration greater than 1.0 indicates a concentration greater than a benchmark, and a relative-concentration less than or equal to 1.0 indicates a concentration equal to or less than a benchmark. Relative-concentrations of organic constituents and special-interest constituents [perchlorate and N-nitrosodimethylamine (NDMA)] were classified as high (relative-concentration greater than 1.0), moderate (relative-concentration greater than 0.1 and less than or equal to 1.0), or low (relative-concentration less than or equal to 0.1). Relative-concentrations of inorganic constituents were classified as high (relative-concentration greater than 1.0), moderate (relative-concentration greater than 0.5 and less than or equal to 1.0), or low (relative-concentration less than or equal to 0.5). Aquifer-scale proportion was used as the primary metric in the status assessment for evaluating regional-scale groundwater quality. High aquifer-scale proportion is defined as the percentage of the area of the primary aquifer system with a high relative-concentration for a particular constituent or class of constituents; this percentage is based on an areal rather than a volumetric basis. Moderate and low aquifer-scale proportions were defined as the percentages of the primary aquifer system with moderate and low relative-concentrations, respectively, of a constituent or class of constituents. Two statistical approaches—grid-based and spatially weighted—were used to evaluate aquifer-scale proportions for individual constituents and classes of constituents. Grid-based and spatially weighted estimates were comparable to each other (within 90-percent confidence intervals) in the study unit. Inorganic constituents (one or more) with health-based benchmarks were detected at high relative-concentrations in 48 percent of the primary aquifer system and at moderate relative-concentrations in 26 percent of the primary aquifer system. The high aquifer-scale proportion of inorganic constituents primarily reflected high aquifer-scale proportions of fluoride (27 percent), arsenic (18 percent), molybdenum (16 percent), boron (10 percent), uranium (5.6 percent), gross alpha radioactivity (9.7 percent), and nitrate (2.7 percent). The inorganic constituents with secondary maximum contaminant levels (SMCLs) were detected at high relative-concentrations in 13 percent of the primary aquifer system and at moderate relative-concentrations in 39 percent. The high aquifer-scale proportion for SMCL constituents reflected high aquifer-scale proportions of total dissolved solids (TDS, 11 percent), manganese (2.8 percent), and chloride (2.8 percent). Organic constituents were not detected at high relative-concentrations in the primary aquifer system, and were present at moderate relative-concentrations in 5.0 percent, and at low relative-concentrations or were not detected in 95 percent of the primary aquifer system. Of the 148 organic constituents analyzed, 12 constituents were detected. Two organic constituents, chloroform and tetrachloroethene (PCE), were detected in more than 10 percent of samples, but were detected mostly at low relative-concentrations.

California↗

USGS2014

Groundwater and surface-water interaction and potential for underground water storage in the Buena Vista-Salida Basin, Chaffee County, Colorado, 2011

By 2030, the population of the Arkansas Headwaters Region, which includes all of Chaffee and Lake Counties and parts of Custer, Fremont, and Park Counties, Colorado, is forecast to increase about 73 percent. As the region’s population increases, it is anticipated that groundwater will be used to meet much of the increased demand. In September 2009, the U.S. Geological Survey, in cooperation with the Upper Arkansas Water Conservancy District and with support from the Colorado Water Conservation Board; Chaffee, Custer, and Fremont Counties; Buena Vista, Cañon City, Poncha Springs, and Salida; and Round Mountain Water and Sanitation District, began a 3-year study of groundwater and surface-water conditions in the Buena Vista-Salida Basin. This report presents results from the study of the Buena Vista-Salida Basin including synoptic gain-loss measurements and water budgets of Cottonwood, Chalk, and Browns Creeks, changes in groundwater storage, estimates of specific yield, transmissivity and hydraulic conductivity from aquifer tests and slug tests, an evaluation of areas with potential for underground water storage, and estimates of stream-accretion response-time factors for hypothetical recharge and selected streams in the basin. The four synoptic measurements of flow of Cottonwood, Chalk, and Browns Creeks, suggest quantifiable groundwater gains and losses in selected segments in all three perennial streams. The synoptic measurements of flow of Cottonwood and Browns Creeks suggest a seasonal variability, where positive later-irrigation season values in these creeks suggest groundwater discharge, possibly as infiltrated irrigation water. The overall sum of gains and losses on Chalk Creek does not indicate a seasonal variability but indicates a gaining stream in April and August/September. Gains and losses in the measured upper segments of Chalk Creek likely are affected by the Chalk Cliffs Rearing Unit (fish hatchery). Monthly water budgets were estimated for selected segments of five perennial streams (Cottonwood, North Cottonwood, Chalk, and Browns Creeks, and South Arkansas River) in the Buena Vista-Salida Basin for calendar year 2011. Differences between reported diversions and estimated crop irrigation requirements were used to estimate groundwater recharge in the areas irrigated by water supplied from the diversions. The amount of groundwater recharge in all the basins varied monthly; however, the greatest amount of recharge was during June and July for Cottonwood, North Cottonwood, and Chalk Creeks and South Arkansas River. The greatest amount of recharge in 2011 in Browns Creek occurred in July and August. The large seasonal fluctuations of groundwater near irrigated areas in the Buena Vista-Salida Basin indicate that the increased groundwater storage resulting from infiltration of surface-water diversions has dissipated by the following spring. Areas within the Buena Vista-Salida Basin with the potential for underground storage were identified using geographic information system data, including topographic, geologic, and hydrologic data, excluding the mountainous areas that border the Buena Vista-Salida Basin and igneous and metamorphic rock outcrop areas. The areas that met the selection criteria for underground water storage are located on terrace deposits near the Arkansas River and adjacent to its major tributaries. The selected areas also contain much of the irrigated land within the basin; consequently, irrigation ditches and canals could provide a means of conveying water to potential recharge sites.

Colorado↗

USGS2014

A reconnaissance spatial and temporal assessment of methane and inorganic constituents in groundwater in bedrock aquifers, Pike County, Pennsylvania, 2012-13

Pike County in northeastern Pennsylvania is underlain by the Devonian-age Marcellus Shale and other shales, formations that have potential for natural gas development. During 2012–13, the U.S. Geological Survey in cooperation with the Pike County Conservation District conducted a reconnaissance study to assess baseline shallow groundwater quality in bedrock aquifers prior to possible shale-gas development in the county. For the spatial component of the assessment, 20 wells were sampled in summer 2012 to provide data on the occurrence of methane and other aspects of existing groundwater quality throughout the county, including concentrations of inorganic constituents commonly present at low levels in shallow, fresh groundwater but elevated in brines. For the temporal component of the assessment, 4 of the 20 wells sampled in summer 2012 were sampled monthly from July 2012 through June 2013 to provide data on seasonal variability in groundwater quality. All water samples were analyzed for major ions, nutrients, selected inorganic trace constituents (including metals and other elements), stable isotopes of water, radon-222, gross alpha- and gross beta-particle activity, dissolved gases (methane, ethane, and ethene), and, if possible, isotopic composition of methane. Additional analyses for boron and strontium isotopes, age-dating of water, and radium-226 were done on water samples collected from six wells in June 2013. Results of the summer 2012 sampling show that water from 16 (80 percent) of 20 wells had detectable concentrations of methane, but concentrations were less than 0.1 milligram per liter (mg/L) in most well-water samples; only two well-water samples had concentrations greater than 1 mg/L. The groundwater with elevated methane also had a chemical composition that differed in some respects (pH, selected major ions, and inorganic trace constituents) from groundwater with low methane concentrations. The two well-water samples with the highest methane concentrations (about 3.7 and 5.8 mg/L) also had the highest pH values (8.7 and 8.3, respectively) and the highest concentrations of sodium, lithium, boron, fluoride, and bromide. Elevated concentrations of some other constituents, such as barium, strontium, and chloride, were not limited to well-water samples with elevated methane, although the two samples with elevated methane also had among the highest concentrations of these constituents. One sample with elevated methane concentrations also had elevated arsenic concentrations, with the arsenic concentration of 30 micrograms per liter (μg/L) exceeding the drinking-water standard of 10 µg/L for arsenic. No other sample from the 20 wells sampled in summer 2012 had concentrations of constituents that exceeded any established primary drinking-water standards. However, radon-222 activities ranging up to 4,500 picocuries per liter (pCi/L) exceeded the proposed drinking-water standard of 300 pCi/L in 85 percent of the 20 well-water samples. The isotopic composition methane in the two high-methane samples (δC CH4 values of -64.55 and -64.41 per mil and δD CH4 values of -216.9 and -201.8 per mil, respectively) indicates a predominantly microbial source for the methane formed by a carbon dioxide reduction process. The stable isotopic composition of water (δD H20 and δ 18 O H20 ) in samples from all 20 wells falls on the local meteoric line, indicating water in the wells was of relatively recent meteoric origin (modern precipitation), including samples with elevated methane concentrations. Analytical results for 4 of the 20 wells sampled monthly for 1 year ending June 2013 in order to assess temporal variability in groundwater quality show that concentrations of major ions generally varied by less than 20 percent, with most differences less than 4 mg/L. Concentrations of methane varied by less than 1 μg/L (0.001 mg/L) in samples from three wells with low methane and by as much as 1 mg/L (1,000 μg/L) in samples from one well with relatively high methane. The isotopic composition of methane in the one well with relatively high methane varied slightly in the monthly samples, ranging from about -64.5 to -64.8 per mil for δ 13 C CH4 and from about -217 to -228 per mil for δD CH4 . The δ 13 C values for dissolved inorganic carbon (DIC) in water from this well were consistent with microbial methane formation by carbon dioxide reduction (drift-type methane) and varied little in the temporal samples, ranging from -10.5 to -10.1 per mil. Additional analyses of samples collected in late June 2013 from six wells with a range of methane and trace constituent concentrations provided baseline data on strontium and boron isotopic compositions ( 87 Sr/ 86 Sr ratios and δ 11 B, respectively) that potentially may be used to differentiate among sources of these constituents. The strontium and boron isotopic composition determined in the six shallow Pike County groundwater samples had 87 Sr/ 86 Sr ratios of 0.71426 to 0.71531 and δ 11 B values of 11.7 to 27.0 per mil, which differ from those reported for brines in Devonian-age formations in Pennsylvania. The June 2013 samples were also analyzed for radium-226 and age-dating dissolved gases. Activities of radium-226 ranged from 0.041 to 0.29 pCi/L in water samples from the six wells and were less than the drinking-water standard of 5 pCi/L for combined radium-226 and radium-228. Age-dating of groundwater using a method based on the presence of anthropogenic gases (chlorofluorocarbons and sulfur hexafluoride) released into the atmosphere yielded estimated recharge dates for water from these six wells that ranged from the 1940s to early 2000s. The oldest water was in samples from wells that had the highest methane concentrations and the youngest water was in samples from a continuously pumped 300-foot deep production well.

Pennsylvania↗

USGS2015

Simulation of groundwater flow and chloride transport in the “1,200-foot” sand with scenarios to mitigate saltwater migration in the “2,000-foot” sand in the Baton Rouge area, Louisiana

Groundwater withdrawals have caused saltwater to encroach into freshwater-bearing aquifers beneath Baton Rouge, Louisiana. The 10 aquifers beneath the Baton Rouge area, which includes East and West Baton Rouge Parishes, Pointe Coupee Parish, and East and West Feliciana Parishes, provided about 184.3 million gallons per day (Mgal/d) for public supply and industrial use in 2012. Groundwater withdrawals from the “1,200-foot” sand in East Baton Rouge Parish have caused water-level drawdown as large as 177 feet (ft) north of the Baton Rouge Fault and limited saltwater encroachment from south of the fault. The recently developed groundwater model for simulating transport in the “2,000-foot” sand was rediscretized to also enable transport simulation within the “1,200-foot” sand and was updated with groundwater withdrawal data through 2012. The model was recalibrated to water-level observation data through 2012 with the parameter-estimation code PEST and calibrated to observed chloride concentrations at observation wells within the “1,200-foot” sand and “2,000-foot” sand. The model is designed to evaluate strategies to control saltwater migration, including changes in the distribution of groundwater withdrawals and installation of scavenger wells to intercept saltwater before it reaches existing production wells. Seven hypothetical scenarios predict the effects of different groundwater withdrawal options on groundwater levels and the transport of chloride within the “1,200-foot” sand and the “2,000-foot” sand during 2015–2112. The predicted water levels and concentrations for all scenarios are depicted in maps for the years 2047 and 2112. The first scenario is a base case for comparison to the six other scenarios and simulates continuation of 2012 reported groundwater withdrawals through 2112 (100 years). The second scenario that simulates increased withdrawals from industrial wells in the “1,200-foot” sand predicts that water levels will be 12–25 ft lower by 2047 and that there will be a negligible difference in chloride concentrations within the “1,200-foot” sand. The five other scenarios simulate the effects of various withdrawal schemes on water levels and chloride concentrations within the “2,000-foot” sand. Amongst these five other scenarios, three of the scenarios simulate only various withdrawal reductions, whereas the two others also incorporate withdrawals from a scavenger well that is designed to extract salty water from the base of the “2,000-foot” sand. Two alternative pumping rates (2.5 Mgal/d and 1.25 Mgal/d) are simulated in each of the scavenger-well scenarios. For the “2,000-foot” sand scenarios, comparison of the predicted effects of the scenarios is facilitated by graphs of predicted chloride concentrations through time at selected observation wells, plots of salt mass in the aquifer through time, and a summary of the predicted plume area and average concentration. In all scenarios, water levels essentially equilibrate by 2047, after 30 years of simulated constant withdrawal rates. Although predicted water-level recovery within the “2,000-foot” sand is greatest for the scenario with the greatest reduction in groundwater withdrawal from that aquifer, the scavenger-well scenarios are most effective in mitigating the future extent and concentration of the chloride plume. The simulated scavenger-well withdrawal rate has more influence on the plume area and concentration than do differences among the scenarios in industrial and public-supply withdrawal rates.

Louisiana, Mississippi↗

USGS2016

Delineation of the Pahute Mesa–Oasis Valley groundwater basin, Nevada

This report delineates the Pahute Mesa–Oasis Valley (PMOV) groundwater basin, where recharge occurs, moves downgradient, and discharges to Oasis Valley, Nevada. About 5,900 acre-feet of water discharges annually from Oasis Valley, an area of springs and seeps near the town of Beatty in southern Nevada. Radionuclides in groundwater beneath Pahute Mesa, an area of historical underground nuclear testing at the Nevada National Security Site, are believed to be migrating toward Oasis Valley. Delineating the boundary of the PMOV groundwater basin is necessary to adequately assess the potential for transport of radionuclides from Pahute Mesa to Oasis Valley. The PMOV contributing area is defined based on regional water-level contours, geologic controls, and knowledge of adjacent flow systems. The viability of this area as the contributing area to Oasis Valley and the absence of significant interbasin flow between the PMOV groundwater basin and adjacent basins are shown regionally and locally. Regional constraints on the location of the contributing area boundary and on the absence of interbasin groundwater flow are shown by balancing groundwater discharges in the PMOV groundwater basin and adjacent basins against available water from precipitation. Internal consistency for the delineated contributing area is shown by matching measured water levels, groundwater discharges, and transmissivities with simulated results from a single-layer, steady-state, groundwater-flow model. An alternative basin boundary extending farther north than the final boundary was rejected based on a poor chloride mass balance and a large imbalance in the northern area between preferred and simulated recharge.

California, Nevada↗

USGS2016

Groundwater-flow model for the Wood River Valley aquifer system, south-central Idaho

A three-dimensional numerical model of groundwater flow was developed for the Wood River Valley (WRV) aquifer system, Idaho, to evaluate groundwater and surface-water availability at the regional scale. This mountain valley is located in Blaine County and has a drainage area of about 2,300 square kilometers (888 square miles). The model described in this report can serve as a tool for water-rights administration and water-resource management and planning. The model was completed with support from the Idaho Department of Water Resources, and is part of an ongoing U.S. Geological Survey effort to characterize the groundwater resources of the WRV. A highly reproducible approach was taken for constructing the WRV groundwater-flow model. The collection of datasets, source code, and processing instructions used to construct and analyze the model was distributed as an R statistical-computing and graphics package. Flow in the WRV aquifer was simulated using the MODFLOW-USG groundwater flow model. The transient flow model simulates groundwater flow between 1995 and 2010. The model uses a 100-meter (328-feet) uniform grid spacing with 54,922 active model cells distributed over three model layers. A confining unit in the south-central part of the Bellevue fan necessitated the use of a multi-layer model. Specified-flow boundaries were used to simulate the groundwater inflows from each of the major tributary basins (also known as tributary basin underflow) and the areal recharge of precipitation and applied irrigation. Head‑dependent flow boundaries were used to simulate the stream-aquifer flow exchange in river reaches and the groundwater discharge at the outlet boundaries of Stanton Crossing and Silver Creek. The model was calibrated by adjusting aquifer hydraulic properties to match simulated and measured water levels and stream-aquifer flow exchange, using the parameter-estimation program PEST. The model reasonably simulated the measured water-table elevation, orientation, and gradients. Stream-aquifer flow exchange along river reaches also was reasonably simulated by the model. Inflow into the WRV aquifer system originates from three sources (from largest to smallest): Streamflow loss from the Big Wood River and Silver Creek, Areal recharge of precipitation and applied irrigation, and Tributary basin underflow. Outflow from the WRV aquifer system originates from five sources (from largest to smallest): Aquifer discharge into the Big Wood River and Silver Creek, Production-well pumping, Discharge across the water table into the vadose zone, Subsurface outflow beneath Silver Creek near Picabo, and Subsurface outflow beneath the Big Wood River near Stanton Crossing. Temporal changes in aquifer storage are most affected by areal recharge and groundwater pumping, and also contribute to changes in streamflow gains.

Idaho↗

USGS2017

Groundwater flux and nutrient loading in the northeast section of Bear Lake, Muskegon County, Michigan, 2015

Bear Lake in North Muskegon, Michigan, is listed as part of the Muskegon Lake area of concern as designated by the U.S. Environmental Protection Agency. This area of concern was designated as a result of eutrophication and beneficial use impairments. On the northeast end of Bear Lake, two man-made retention ponds (Willbrandt Pond East and Willbrandt Pond West), formerly used for celery farming, may contribute nutrients to Bear Lake. Willbrandt Ponds (East and West) were previously muck fields that were actively used for celery farming from the early 1900s until 2002. The restoration and reconnection of the Willbrandt Ponds into Bear Lake prompted concerns of groundwater nutrient loading into Bear Lake. Studies done by the State of Michigan and Grand Valley State University revised initial internal phosphorus load estimates and indicated an imbalance in the phosphorus budget in Bear Lake. From June through November 2015, the U.S. Geological Survey (USGS) did an investigative study to quantify the load of nutrients from shallow groundwater around the Willbrandt Ponds in an effort to update the phosphorus budget to Bear Lake. Seven sampling locations were established, including five shallow groundwater wells and two surface-water sites, in the Willbrandt pond study area and Bear Lake. A total of 12 nutrient samples and discrete water-level measurements were collected from each site from June through November 2015. Continuous water-level data were recorded for both surface-water monitoring locations for the entire sampling period. Water-level data indicated that Willbrandt Pond West had the highest average water-level elevation of all sites monitored, which indicated the general direction of flux is from Willbrandt Pond West to Bear Lake. Nutrient and chloride loading from Willbrandt Pond West to Bear Lake was calculated using two distinct methods: Dupuit and direct seepage methods. Shallow groundwater loading calculations were determined by using groundwater levels to first determine a flux of shallow groundwater, then nutrient concentrations to determine a load. It was determined that Willbrandt Pond East and Willbrandt Pond West contributed between 2 to 4 percent of the total annual phosphorus load to Bear Lake by way of shallow groundwater flow. Annual loads calculated for other constituents include orthophosphate (40–100 pounds per year [lb P/yr]), total nitrogen (200–830 lb/yr), chloride (12,700–32,100 lb/yr), and ammonia (130–670 lb N/yr). Study results indicated that mean groundwater and surface-water nutrient concentrations calculated in this study were higher than reported Michigan statewide values. The data collected in this study allow understanding of groundwater nutrient loading into Bear Lake in an effort to help inform future restoration and management decisions.

Michigan↗

USGS2018

Groundwater-storage change and land-surface elevation change in Tucson Basin and Avra Valley, south-central Arizona--2003-2016

The U.S. Geological Survey monitors groundwater-storage change and land-surface elevation change caused by groundwater withdrawal in Tucson Basin and Avra Valley—the two most populated alluvial basins within the Tucson Active Management Area. The Tucson Active Management Area is one of five active management areas in Arizona established by the 1980 Groundwater Management Act and governed by the Arizona Department of Water Resources. Gravity and land-surface elevation change were monitored every 1 to 3 years at wells and benchmarks in Tucson Basin and Avra Valley from 2003 to 2016. Monitoring resulted in estimates of land-surface elevation change and groundwater-storage change. Interferometric synthetic aperture radar (InSAR) interferograms showing land-surface elevation change were constructed for the Tucson metropolitan area from (1) May 2003 to July 2006, (2) July 2006 to June 2008, (3) June 2008 to April 2011, (4) April 2011 to November 2014, and (5) November 2014 to March 2016. For the Tucson metropolitan area, maximum subsidence of about 2 inches occurred during May 2003 to July 2006. From July 2006 to June 2008, maximum subsidence of approximately 0.8 inches occurred in two regions in the Tucson metropolitan area. From June 2008 to April 2011, about 0.8 inches of subsidence also occurred in two regions. Additionally, for the period April 2011 to November 2014, a maximum of about 0.9 inches of subsidence occurred in the same two regions of Tucson Basin. For the entire monitoring period from May 2003 to March 2016, maximum subsidence of as much as 5.3 inches occurred in the Tucson metropolitan area south of Irvington Road between south 12th Avenue and south Park Avenue, and as much as 4 inches in central Tucson south of Broadway between Country Club Road and Craycroft Road. The InSAR data indicated that there was no significant land-surface deformation from 2003 to 2016 in Avra Valley, and no change in either basin from 2014 to 2016. The volume of stored groundwater in the monitored part of Tucson Basin showed net zero change from spring 2003 to summer 2006. From summer 2006 to summer 2008 the volume of stored groundwater in the monitored part of Tucson Basin increased approximately 50,000 acre-feet; however, overdraft conditions resumed from summer 2008 to spring 2011, resulting in decreased storage of approximately 178,000 acre-feet. From spring 2011 to fall 2014, the volume of stored groundwater in Tucson Basin decreased about 200,000 acre-feet, following a period of lower than average rainfall in 2012 and 2013. The volume of stored groundwater in the monitored part of Tucson Basin increased approximately 167,000 acre-feet from fall 2014 to spring 2016. Groundwater storage in Avra Valley increased during the entire monitoring period from spring 2003 to spring 2016, largely as a result of managed recharge of Central Arizona Project water in the monitored region. From 2003 to 2016, artificial recharge in Avra Valley totaled approximately 1,788,000 acre-feet, and in Tucson Basin artificial recharge for the entire period was about 636,790 acre-feet. Artificial recharge exceeded pumping in Avra Valley for each time interval. Pumping in Tucson Basin exceeded artificial recharge for every period except 2014 to 2016. Overall, long-term water-level declines have stabilized or reversed since 2000 at most areas in Tucson Basin and Avra Valley.

Arizona↗

USGS2020

Hydrogeologic characterization, groundwater chemistry, and vulnerability assessment, Ute Mountain Ute Reservation, Colorado and Utah

The U.S. Geological Survey, in cooperation with the Ute Mountain Ute Tribe (UMUT), initiated a study in 2016 to increase understanding of the hydrogeology and chemistry of groundwater within select areas of the Ute Mountain Ute Reservation (UMUR) in Colorado and Utah, identify vulnerabilities to the system and other natural resources, and outline information needs to aid in the understanding and protection of groundwater resources. The results presented for this study can be used to support the UMUT’s goal of protecting their vital groundwater resources on the UMUR. Hydrogeologic conditions were characterized for the surficial aquifer contained in Quaternary-age unconsolidated surficial deposits and the Dakota aquifer contained in the Cretaceous-age Dakota Sandstone. In the surficial aquifer, median depth to water ranges from about 5.4 to 17.2 feet below land surface in the Farm and Ranch Enterprise area and 11 to 34 feet below land surface in the Towaoc area, and the water table slopes generally southwest or south. A map of depth to the top of the Dakota Sandstone was constructed from existing well data. Depths range from zero in outcrop areas to more than 3,000 feet below land surface on mesas in the southeastern part of the UMUR. Groundwater-chemistry data were collected by the UMUT from 13 springs and 31 wells from 1996 through 2017. Specific conductance was much lower for samples from springs than from wells; median values were 512 and 6,024 microsiemens per centimeter at 25 degrees Celsius, respectively. Spring samples were well oxygenated. A few well samples were anoxic (dissolved oxygen concentrations less than 0.5 milligrams per liter [mg/L]), indicating reducing conditions in the aquifer. About 75 percent of spring samples had fresh water (total dissolved solids concentrations less than 1,000 mg/L), and about 85 percent of well samples had brackish or highly saline water (total dissolved solids concentrations greater than 1,000 mg/L). Water type for springs on the Ute Mountains was calcium bicarbonate. Lower-altitude springs had a calcium-sulfate water type. Most well samples had sodium as the dominant cation, and sulfate, bicarbonate, and chloride as the dominant anions. Fluoride concentrations in about 45 percent of well samples were greater than an agricultural-use standard of 2 mg/L. Nitrate plus nitrite concentrations in most spring and well samples were less than about 1.6 mg/L per liter. Concentrations in samples from wells in the irrigated agricultural area were elevated; the maximum concentration was 78.5 mg/L. About one-half of the trace-element samples had concentrations that were less than laboratory reporting limits. Only aluminum, arsenic, and selenium in spring samples, and boron and selenium in well samples, were detected at concentrations greater than surface-water standards or water-quality standards for agricultural use of groundwater. Only three organic compounds, the pesticides alachlor and atrazine and the volatile organic compound di(2-ethylhexyl) phthalate, were detected in well samples. The Escherichia coli bacteria was detected in 47 and 23 percent of samples from wells and springs, respectively. The E. coli detections included samples from three culturally significant springs, which did not meet the UMUT cultural-use standard of total absence of E. coli. Tritium and carbon-14 were the primary environmental tracers used for interpreting groundwater ages for Lopez 2 Spring and five wells (AP–1, 5000 Block, Cottonwood Spring, Goodknight, and SE Toe). Water from the AP–1 well contained a mixture of pre- and post-1950s recharge. Tritium and carbon-14 recharge ages for Lopez 2 Spring (post-1950s in age), Goodknight and SE Toe wells (pre-1950s in age), and Cottonwood Spring well (primarily pre-1950s in age) are supported by helium-4 data. The helium-4 data for the 5000 Block well are inconsistent with the tritium and carbon-14 age of pre-1950s recharge because of interference caused by high methane concentrations in the water. Springs and surficial deposits are more vulnerable to contamination from anthropogenic chemicals than deeper bedrock wells. Bedrock aquifers are vulnerable in areas where the geologic formations containing the aquifers are exposed at the land surface. Groundwater in deep bedrock aquifers is likely thousands of years old and is not currently affected by present-day land uses. Both shallow and deep groundwater are vulnerable to naturally occurring salts and minerals, such as of total dissolved solids, major ions, nitrate, and trace elements. Effects of a changing climate on water resources and other ecological characteristics of the UMUR could include changes in evapotranspiration, a decrease in snowpack, decreased aquifer recharge and flow of springs, a decrease in soil moisture, and increased occurrence of wildfires and forest mortality. Of particular interest for the UMUT are possible effects of a changing climate on medicinal and culturally important plants and springs Several information needs were identified during this study that would aid in the understanding and protection of groundwater resources on the UMUR. These include well-completion information for bedrock wells, the collection of environmental tracer data at additional wells, the addition of methane and hydrocarbon analysis to well sampling plans, and the resampling of springs and wells that were last sampled in 2002 or earlier.

Colorado↗

USGS2020

Simulation of groundwater flow in the regional aquifer system on Long Island, New York, for pumping and recharge conditions in 2005–15

A three-dimensional groundwater-flow model was developed for the aquifer system of Long Island, New York, to evaluate (1) responses of the hydrologic system to changes in natural and anthropogenic hydraulic stresses, (2) the subsurface distribution of groundwater age, and (3) the regional-scale distribution of groundwater travel times and the source of water to fresh surface waters and coastal receiving waters. The model also provides the groundwater flow components used to define model boundaries for possible inset models used for local-scale analyses. The three-dimensional, groundwater flow model developed for this investigation uses the numerical code MODFLOW–NWT to represent steady-state conditions for average groundwater pumping and aquifer recharge for 2005–15. The particle-tracking algorithm MODPATH, which simulates advective transport in the aquifer, was used to estimate groundwater age, delineate the areas at the water table that contribute recharge to coastal and freshwater bodies, and estimate total travel times of water from the water table to discharge locations. A three-dimensional, 1-meter (3.3-foot) topobathymetric model was used to determine land-surface altitudes for the island and seabed altitudes for the surrounding coastal waters. The mapped extents and surface altitudes of major geologic units were compiled and used to develop a three-dimensional hydrogeologic framework of the aquifer system, including aquifers and confining units. Lithologic data from deep boreholes and previous aquifer-test results were used to estimate the three-dimensional distribution of hydraulic conductivity in principal aquifers. Natural recharge from precipitation was estimated for 2005–15 using a modified Thornthwaite-Mather methodology as implemented in a soil-water balance model. Components of anthropogenic recharge—wastewater return flow, storm water inflow, and inflow from leaky infrastructure—also were estimated for 2005–15. Groundwater withdrawals for various sources, including public water supply, industrial, remediation, and agricultural, were compiled or estimated for the same period. These data were incorporated into the model development to represent the aquifer system geometry, boundaries, and initial hydraulic properties of the regional aquifers and confining units within the Long Island aquifer system. Average hydraulic conditions—water levels and streamflows—for 2005–15 were estimated using existing data from the U.S. Geological Survey National Water Information System database. Model inputs were adjusted to best match average hydrologic conditions using inverse methods as implemented in the parameter-estimating software PEST. The calibrated model was used to simulate average hydrologic conditions in the aquifer system for 2005–15. About 656 cubic feet per second of water was withdrawn on average annually for 2005–15 for water supply and an average of about 349 cubic feet per second of water recharged the aquifer annually from return flow and leaky infrastructure. Parts of New York City have drawdowns exceeding 25 feet, mostly because of urbanization and associated large decreases in recharge rates. Large areas in the western part of the island have drawdowns exceeding 10 feet, mostly from large groundwater withdrawals and sewering, which largely eliminates wastewater return flow. Water-table altitudes in eastern parts of the island increased by more than 2 feet in some areas as a result of wastewater return flow in unsewered areas and changes in land use. Changes in streamflows show a similar pattern as water-table altitudes. Streamflows generally decrease in western parts of the island where there are large drawdowns and increase in eastern parts of the island where water-table altitudes increase.

New York↗

USGS2021

Volcanic Aquifers of Hawai‘i—Construction and calibration of numerical models for assessing groundwater availability on Kaua‘i, O‘ahu, and Maui

Steady-state numerical groundwater-flow models were constructed for the islands of Kaua‘i, O‘ahu, and Maui to enable quantification of the hydrologic consequences of withdrawals and other stresses that can place limits on groundwater availability. The volcanic aquifers of Hawai‘i supply nearly all drinking water for the islands’ residents, freshwater for diverse industries, and natural discharge to springs, streams, and nearshore areas that support ecosystems, cultural practices, aesthetics, and recreation. Increases in groundwater withdrawal and changes in climate can cause water-table depression, saltwater rise, and reduction of natural groundwater discharge—all of which can limit fresh groundwater availability. The numerical models described in this report are designed to quantify these consequences. Separate models were created for each island using MODFLOW-2005 with the Seawater Intrusion package, which allows simulation of freshwater and saltwater in ocean-island aquifers. Calibration resulted in models that generally replicate observed water-level, stream base-flow, and spring-flow data, and simulate groundwater-flow directions and fresh groundwater thicknesses that are consistent with conceptual models. The calibrated models use hydraulic properties that are consistent with the ranges reported in previous studies. The models show that the relative distribution of fresh groundwater discharge to the ocean, streams, and springs and withdrawals for human use differ substantially among the three islands studied here. These differences indicate that consequences that limit the availability of fresh groundwater for human use are likely to differ among the three islands.

Hawaii↗

USGS2023

Groundwater residence times in glacial aquifers—A new general simulation-model approach compared to conventional inset models

Groundwater is important as a drinking-water source and for maintaining base flow in rivers, streams, and lakes. Groundwater quality can be predicted, in part, by its residence time in the subsurface, but the residence-time distribution cannot be measured directly and must be inferred from models. This report compares residence-time distributions from four areas where groundwater flow and travel time were simulated with conventional simulation-inset models (IMs) and with a new automated model-construction method called general simulation models (GSMs). The comparison provides an opportunity to explore controls on travel time and improve the methods used in the creation of GSMs. These models can be useful for three main-use cases: (1) rapid testing of relationships that govern groundwater flow and age, (2) generation of consistent examples for training a machine-learning metamodel, and (3) serving as a starting point for more detailed models. Comparison of the GSMs to IMs indicated a qualified pattern of agreement for residence-time distributions as indicated by the Nash-Sutcliffe efficiency and Spearman’s correlation coefficient. The agreement was best for the median values of the simulated residence times in young fractions of groundwater (defined as the fractions of groundwater in samples less than 65 years old) at the scale of the eight-digit hydrologic-unit code. Generally, the median values of the young fractions in the IMs were correlated with the median values from the GSMs. The relative trends across the four areas also were similar for the other residence-time metrics. The medians of residence-time metrics at finer scales show a fair degree of scatter. The GSM results compared most poorly for median travel times in the older fraction of groundwater (older than 65 years). The GSM approach is intended as a flexible framework for developing models that can be useful individually as screening tools or collectively to support projects in statistical learning. Although one set of GSM algorithms was presented here, the approach can accommodate many types of data and also different categories of prior information. Comparison of GSMs and IMs suggests ways in which the GSMs, while remaining easy to construct and calibrate, can be improved for estimating groundwater travel times. IMs do not yield exact travel times, and matching GSMs to IMs does not guarantee an improvement; however, IMs provide a convenient benchmark against which to explore relations between physical characteristics of watersheds and the distribution of travel times within them. This effort was undertaken as part of the National Water Quality Program of the U.S. Geological Survey to assist in determining the susceptibility of groundwater in glacial aquifers to a variety of natural and anthropogenic contaminants.

Illinois, Indiana, Michigan, Wisconsin↗

USGS2022

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

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

New York↗