SEARCH · Geology Reports

Results for “Ground Water Modeling Investigation Report”

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.

1,134 records · Page 63Linked to original sources

USGS2014

Conceptual model of the uppermost principal aquifer systems in the Williston and Powder River structural basins, United States and Canada

The three uppermost principal aquifer systems of the Northern Great Plains—the glacial, lower Tertiary, and Upper Cretaceous aquifer systems—are described in this report and provide water for irrigation, mining, public and domestic supply, livestock, and industrial uses. These aquifer systems primarily are present in two nationally important fossil-fuelproducing areas: the Williston and Powder River structural basins in the United States and Canada. The glacial aquifer system is contained within glacial deposits that overlie the lower Tertiary and Upper Cretaceous aquifer systems in the northeastern part of the Williston structural basin. Productive sand and gravel aquifers exist within this aquifer system. The Upper Cretaceous aquifer system is contained within bedrock lithostratigraphic units as deep as 2,850 and 8,500 feet below land surface in the Williston and Powder River structural basins, respectively. Petroleum extraction from much deeper formations, such as the Bakken Formation, is rapidly increasing because of recently improved hydraulic fracturing methods that require large volumes of relatively freshwater from shallow aquifers or surface water. Extraction of coalbed natural gas from within the lower Tertiary aquifer system requires removal of large volumes of groundwater to allow degasification. Recognizing the importance of understanding water resources in these energy-rich basins, the U.S. Geological Survey (USGS) Groundwater Resources Program ( http://water.usgs.gov/ogw/gwrp/ ) began a groundwater study of the Williston and Powder River structural basins in 2011 to quantify this groundwater resource, the results of which are described in this report. The overall objective of this study was to characterize, quantify, and provide an improved conceptual understanding of the three uppermost and principal aquifer systems in energy-resource areas of the Northern Great Plains to assist in groundwater-resource management for multiple uses. The study area includes parts of Montana, North Dakota, South Dakota, and Wyoming in the United States and Manitoba and Saskatchewan in Canada. The glacial aquifer system is contained within glacial drift consisting primarily of till, with smaller amounts of glacial outwash sand and gravel deposits. The lower Tertiary and Upper Cretaceous aquifer systems are contained within several formations of the Tertiary and Cretaceous geologic systems, which are hydraulically separated from underlying aquifers by a basal confining unit. The lower Tertiary and Upper Cretaceous aquifer systems each were divided into three hydrogeologic units that correspond to one or more lithostratigraphic units. The period prior to 1960 is defined as the predevelopment period when little groundwater was extracted. From 1960 through 1990, numerous flowing wells were installed near the Yellowstone, Little Missouri and Knife Rivers, resulting in local groundwater declines. Recently developed technologies for the extraction of petroleum resources, which largely have been applied in the study area since about 2005, require millions of gallons of water for construction of each well, with additional water needed for long-term operation; therefore, the potential for an increase in groundwater extraction is high. In this study, groundwater recharge and discharge components were estimated for the period 1981–2005. Groundwater recharge primarily occurs from infiltration of rainfall and snowmelt (precipitation recharge) and infiltration of streams into the ground (stream infiltration). Total estimated recharge to the Williston and Powder River control volumes is 4,560 and 1,500 cubic feet per second, respectively. Estimated precipitation recharge is 26 and 15 percent of total recharge for the Williston and Powder River control volumes, respectively. Estimated stream infiltration is 71 and 80 percent of total recharge for the Williston and Powder River control volumes, respectively. Groundwater discharge primarily is to streams and springs and is estimated to be about 97 and 92 percent of total discharge for the Williston and Powder River control volumes, respectively. Most of the remaining discharge results from pumped and flowing wells. Groundwater flow in the Williston structural basin generally is from the west and southwest toward the east, where discharge to streams occurs. Locally, in the uppermost hydrogeologic units, groundwater generally is unconfined and flows from topographically high to low areas, where discharge to streams occurs. Groundwater flow in the Powder River structural basin generally is toward the north, with local variations, particularly in the upper Fort Union aquifer, where flow is toward streams.

Manitoba;Montana;North Dakota;Saskatchewan;South D↗

USGS1989

Geohydrology of the alluvial and terrace deposits of the North Canadian River from Oklahoma City to Eufaula Lake, central Oklahoma

This investigation was undertaken to describe the geohydrology of the alluvial and terrace deposits along the North Canadian River between Lake Overholser and Eufaula Lake, an area of about 1,835 square miles, and to determine the maximum annual yield of ground water. A 1982 water-level map of the alluvial and terrace aquifer was prepared using field data and published records. Data from test holes and other data from the files of the U.S. Geological Survey and the Oklahoma Water Resources Board were used to establish the approximate thickness of the alluvial and terrace deposits. The North Canadian River from Lake Overholser, near Oklahoma City, to Eufaula Lake is paralleled by a 2- to 3-mile wide band of alluvium. Scattered terrace deposits on either side of the alluvium reach an extreme width of 8 miles. Rocks of Permian age bound the alluvial and terrace deposits from the west to the midpoint of the study area; Pennsylvanian rocks bound the alluvial and terrace deposits from that point eastward. Three major aquifers are present in the study area: the alluvial and terrace aquifer, consisting of alluvium and terrace deposits of Quaternary age in a narrow band on either side of the North Canadian River; the Garber-Wellington aquifer of Permian age, consisting of an upper unconfined zone and a lower confined zone separated by relatively impermeable shales; and the Ada-Vamoosa aquifer of Pennsylvanian age. At locations were the alluvial and terrace aquifer overlies either of the other aquifers, there is hydraulic continuity between the alluvial and terrace aquifer and the other aquifers, and water levels are the same. Most large-scale municipal and industrial pumping from the Garber-Wellington aquifer is from the lower zone and has little discernible effect upon the alluvial and terrace aquifer. The total estimated base flow of the North Canadian River for the studied reach is 264 cubic feet per second. Evapotranspiration from the basin in August is about 60 cubic feet per second for the North Canadian River from Lake Overholser to a measuring station above Eufaula Lake. Estimated recharge rates to the alluvial and terrace aquifer in the basin range from 1.7 inches at the west edge of the study area to 7.0 inches at the east edge. Total permitted withdrawal from the aquifer, according to records of the Oklahoma Water Resources Board, ranged from 2,107 acre-feet per year in 1942 to about 21,415 acre-feet per year in 1982. Simulations of the alluvial and terrace aquifer from Lake Overholser to Eufaula Lake were made using a finite-difference model developed by McDonald and Harbaugh (1984). The area of the aquifers was subdivided into a finite-difference grid having 30 rows and 57 columns with cells measuring 1 mile in the north-south direction and 2 miles in the east-west direction. The model was calibrated in two steps: A steady-state calibration simulated head distribution prior to extensive pumping of the aquifer in 1942, and a transient calibration simulated head distribution after extensive pumpage. The final horizontal hydraulic conductivity used for the alluvial and terrace aquifer was 0.0036 feet per second (310 feet per day) at all locations. The recharge rate for the alluvial and terrace aquifer ranged from 1.7 inch per year in the west to 7.0 inches per year in the east, and averaged about 3.3 inches per year. A specific yield of 15 percent was used for the transient simulation. Permitted pumpage for 1942 through 1982 was used in the digital model to estimate the annual volume of water in storage in the alluvial and terrace aquifer for the years for this time period. The 1982 permitted pumpage rates were used for projections for 1983 to 2020. The estimated volume of water in storage was 1,940,000 acre-feet in 1982. Because the estimated recharge rate is equal to the allowed pumpage rate in 1982, the projected volume of water in storage in both 1993 and 2020 was 1,890,000 acre-feet.

Water-Resources Investigations Report↗

USGS2006

Diffuse-flow conceptualization and simulation of the Edwards aquifer, San Antonio region, Texas

A numerical ground-water-flow model (hereinafter, the conduit-flow Edwards aquifer model) of the karstic Edwards aquifer in south-central Texas was developed for a previous study on the basis of a conceptualization emphasizing conduit development and conduit flow, and included simulating conduits as one-cell-wide, continuously connected features. Uncertainties regarding the degree to which conduits pervade the Edwards aquifer and influence ground-water flow, as well as other uncertainties inherent in simulating conduits, raised the question of whether a model based on the conduit-flow conceptualization was the optimum model for the Edwards aquifer. Accordingly, a model with an alternative hydraulic conductivity distribution without conduits was developed in a study conducted during 2004-05 by the U.S. Geological Survey, in cooperation with the San Antonio Water System. The hydraulic conductivity distribution for the modified Edwards aquifer model (hereinafter, the diffuse-flow Edwards aquifer model), based primarily on a conceptualization in which flow in the aquifer predominantly is through a network of numerous small fractures and openings, includes 38 zones, with hydraulic conductivities ranging from 3 to 50,000 feet per day. Revision of model input data for the diffuse-flow Edwards aquifer model was limited to changes in the simulated hydraulic conductivity distribution. The root-mean-square error for 144 target wells for the calibrated steady-state simulation for the diffuse-flow Edwards aquifer model is 20.9 feet. This error represents about 3 percent of the total head difference across the model area. The simulated springflows for Comal and San Marcos Springs for the calibrated steady-state simulation were within 2.4 and 15 percent of the median springflows for the two springs, respectively. The transient calibration period for the diffuse-flow Edwards aquifer model was 1947-2000, with 648 monthly stress periods, the same as for the conduit-flow Edwards aquifer model. The root-mean-square error for a period of drought (May-November 1956) for the calibrated transient simulation for 171 target wells is 33.4 feet, which represents about 5 percent of the total head difference across the model area. The root-mean-square error for a period of above-normal rainfall (November 1974-July 1975) for the calibrated transient simulation for 169 target wells is 25.8 feet, which represents about 4 percent of the total head difference across the model area. The root-mean-square error ranged from 6.3 to 30.4 feet in 12 target wells with long-term water-level measurements for varying periods during 1947-2000 for the calibrated transient simulation for the diffuse-flow Edwards aquifer model, and these errors represent 5.0 to 31.3 percent of the range in water-level fluctuations of each of those wells. The root-mean-square errors for the five major springs in the San Antonio segment of the aquifer for the calibrated transient simulation, as a percentage of the range of discharge fluctuations measured at the springs, varied from 7.2 percent for San Marcos Springs and 8.1 percent for Comal Springs to 28.8 percent for Leona Springs. The root-mean-square errors for hydraulic heads for the conduit-flow Edwards aquifer model are 27, 76, and 30 percent greater than those for the diffuse-flow Edwards aquifer model for the steady-state, drought, and above-normal rainfall synoptic time periods, respectively. The goodness-of-fit between measured and simulated springflows is similar for Comal, San Marcos, and Leona Springs for the diffuse-flow Edwards aquifer model and the conduit-flow Edwards aquifer model. The root-mean-square errors for Comal and Leona Springs were 15.6 and 21.3 percent less, respectively, whereas the root-mean-square error for San Marcos Springs was 3.3 percent greater for the diffuse-flow Edwards aquifer model compared to the conduit-flow Edwards aquifer model. The root-mean-square errors for San Antonio and San Pedro Springs were appreciably greater, 80.2 and 51.0 percent, respectively, for the diffuse-flow Edwards aquifer model. The simulated water budgets for the diffuse-flow Edwards aquifer model are similar to those for the conduit-flow Edwards aquifer model. Differences in percentage of total sources or discharges for a budget component are 2.0 percent or less for all budget components for the steady-state and transient simulations. The largest difference in terms of the magnitude of water budget components for the transient simulation for 1956 was a decrease of about 10,730 acre-feet per year (about 2 per-cent) in springflow for the diffuse-flow Edwards aquifer model compared to the conduit-flow Edwards aquifer model. This decrease in springflow (a water budget discharge) was largely offset by the decreased net loss of water from storage (a water budget source) of about 10,500 acre-feet per year.

Texas↗

USGS1995

Geohydrology, water quality, and conceptual model of the hydrologic system Saco Landfill area, Saco, Maine

A geohydrologic study of the Saco Municipal Landfill in Saco, Maine, was done during 1993-94 to provide a preliminary interpretation of the geology and hydrology needed to guide additional studies at the landfill as part of the Superfund Program. The Saco Landfill, which was active from the early 1960's until 1986, includes three disposal areas on a 90-acre parcel. Sandy Brook, a small perennial stream, flows from north to south through the land-fill between the disposal areas. Discharge of leachate from the disposal areas to aquifers and streams has been documented since 1974. The landfill was declared a Superfund site in 1990 by the U.S. Environmental Protection Agency. Multiple lines of evidence are used in this study to indicate areas of ground-water contamination and sources of water flow in Sandy Brook. The geohydrologic system on the east side of Sandy Brook consists of an upper water-table aquifer and a lower aquifer, separated by a thick sequence of glaciomarine silt and clay. Depths to bedrock range from 60 to more than 200 ft (feet), on the basis of data from seismic-refraction studies and drilling. The upper aquifer, which is generally less than 15 ft thick, consists of fine-to medium-grained sand deposited in a shallow postglacial marine environment. The lower aquifer, which was deposited as a series of glaciomarine fans, contains two sediment types: Well-sorted sand and gravel and unsorted sediments called diamict sediments. East of Sandy Brook, the thickness of the lower aquifer ranges from 25 to 100 ft, based on drilling at the landfill. The glaciomarine silts and clays (known as the presumpscot Formation) range from 50 to more than 100 ft thick. West of Sandy Brook, the glaciomarine silt and clay is largely absent, and fractured bedrock is very close to land surface under one of the disposal areas in the northwestern part of the property. The lower aquifer is unconfined in the southwestern side of the study area; bedrock slopes towards the south, and the aquifer thickens to 100 ft at the southwestern end of the study area. Preliminary estimates of mean annual streamflow in Sandy Brook, based on a partial year of continuous record, indicate that runoff increases from approximately 2.1 ft3/s (cubic feet per second) upstream from the landfill to 2.7 ft3/s downstream from the landfill, although the drainage area down-stream is only 11 percent greater than the drainage area upstream. A water-budget estimate based on available streamflow and climatic data indicates that Sandy Brook below the landfill gains about 80 million gallons per year from sources outside the drainage-basin boundary. Possible sources include the lower aquifer north or west of the landfill area and the fractured bedrock northwest of Sandy Brook. Specific conductance of water in Sandy Brook increases downstream from the landfill. In September 1993, specific conductance was 184 liS/cm (microsiemens per centimeter at 25 degrees Celsius) upstream from the landfill and 496 uS/cm downstream from the landfill. Continuous monitoring of specific conductance in Sandy Brook shows that the downstream increase is less during periods of stormflow because of dilution. Electromagnetic terrain-conductivity surveys, results of ground-water chemical analyses, and changes in streamwater quality have been used to identify areas of likely ground-water contamination. The specific conductance of ground water exceeds 2,000 uS/cm in some areas near the landfills. This compares to specific conductances of less than 200 uS/cm in water from most shallow wells that are considered to represent background water quality. Ground water in the upper aquifer east of Sandy Brook and in the lower aquifer west of Sandy Brook has been affected by leachate flowing from the landfill areas. The extent of contamination in bedrock, if any, is unknown. Water levels measured in 16 wells were used to help determine the direction of ground-water flow. The electromagnetic terrain-conduct surveys and stream specific-conductance data sup port the interpretation that water in the upper aqui fer flows radially away from the two disposal areas east of Sandy Brook towards Sandy Brook and other small surface-water bodies in the area. West of Sandy Brook, ground water under the third dis posal area moves in the lower aquifer northeast and southeast towards Sandy Brook, where it discharges to the stream.

Maine↗

USGS2005

Summary of hydraulic properties of the Floridan Aquifer system in coastal Georgia and adjacent parts of South Carolina and Florida

Hydraulic-property data for the Floridan aquifer system and equivalent clastic sediments in a 67-county area of coastal Georgia and adjacent parts of South Carolina and Florida were evaluated to provide data necessary for development of ground-water flow and solute-transport models. Data include transmissivity at 324 wells, storage coefficient at 115 wells, and vertical hydraulic conductivity of 72 core samples from 27 sites. Hydraulic properties of the Upper Floridan aquifer vary greatly in the study area due to the heterogeneity (and locally to anisotropy) of the aquifer and to variations in the degree of confinement provided by confining units. Prominent structural features in the areathe Southeast Georgia Embayment, the Beaufort Arch, and the Gulf Troughinfluence the thickness and hydraulic properties of the sediments comprising the Floridan aquifer system. Transmissivity of the Upper Floridan aquifer and equivalent updip units was compiled for 239 wells and ranges from 530 feet squared per day (ft2/d) at Beaufort County, South Carolina, to 600,000 ft2/d in Coffee County, Georgia. In carbonate rock settings of the lower Coastal Plain, transmissivity of the Upper Floridan aquifer generally is greater than 20,000 ft2/d, with values exceeding 100,000 ft2/d in the southeastern and southwestern parts of the study area (generally coinciding with the area of greatest aquifer thickness). Transmissivity of the Upper Floridan aquifer generally is less than 10,000 ft2/d in and near the upper Coastal Plain, where the aquifer is thin and consists largely of clastic sediments, and in the vicinity of the Gulf Trough, where the aquifer consists of low permeability rocks and sediments. Large variability in the range of transmissivity in Camden and Glynn Counties, Georgia, and Nassau County, Florida, demonstrates the anisotropic distribution of hydraulic properties that may result from fractures or solution openings in the carbonate rocks. Storage coefficient of the Upper Floridan aquifer was compiled for 106 wells and ranges from about 0.00004 at Beaufort County, South Carolina, to 0.04 in Baker County, Florida. Transmissivity of the Lower Floridan aquifer and equivalent updip clastic units was compiled for 53 wells and ranges from about 170 ft2/d in Barnwell County, South Carolina, to about 43,000 ft2/d in Camden County, Georgia. Transmissivity of the Lower Floridan aquifer is greatest where the aquifer is thickest in southeastern Georgia and northeastern Floridawhere estimates are greater than 10,000 ft2/d; at one well in southeastern Georgia transmissivity was estimated to be as high as 200,000 ft2/d. Storage-coefficient data for the Lower Floridan aquifer are limited to three estimates in Barnwell and Allendale Counties, South Carolina, and to estimates determined from six multi-aquifer tests in Duval County, Florida. In the South Carolina tests, storage coefficient ranges from 0.0003 to 0.0004; this range is indicative of a confined aquifer. The storage coefficient for the combined Upper and Lower Floridan wells in Duval County, Florida, ranges from 0.00002 to 0.02. Vertical hydraulic conductivity was compiled from core samples collected at 27 sites. For the Upper Floridan confining unit, values from 39 core samples at 17 sites range from 0.0002 to 3 feet per day (ft/d). For the Lower Floridan confining unit, values from 10 core samples at 9 sites range from about 0.000004 to 0.16 ft/d. Vertical hydraulic conductivity of the Upper Floridan aquifer was compiled from 16 core samples at five sites, mostly in the Brunswick, Georgia, area and values range from 0.00134 to 160.4 ft/d. Vertical hydraulic conductivity for the semiconfining unit separating the upper and lower water-bearing zones of the Upper Floridan at Brunswick, Georgia, compiled from 6 core samples at three sites ranges from 0.000008 to 0.000134 ft/d. The vertical hydraulic conductivity of the Lower Floridan aquifer in a core sample from a well at Brunswick, G

Florida, Georgia, South Carolina↗

USGS1999

Hydrogeology, water use, and simulation of flow in the High Plains aquifer in northwestern Oklahoma, southeastern Colorado, southwestern Kansas, northeastern New Mexico, and northwestern Texas

The U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, began a three-year study of the High Plains aquifer in northwestern Oklahoma in 1996. The primary purpose of this study was to develop a ground-water flow model to provide the Water Board with the information it needs to manage the quantity of water withdrawn from the aquifer. The study area consists of about 7,100 square miles in Oklahoma and about 20,800 square miles in adjacent states to provide appropriate hydrologic boundaries for the flow model. The High Plains aquifer includes all sediments from the base of the Ogallala Formation to the potentiometric surface. The saturated thickness in Oklahoma ranges from more than 400 feet to less than 50 feet. Natural recharge to the aquifer from precipitation occurs throughout the area but is extremely variable. Dryland agricultural practices appear to enhance recharge from precipitation, and part of the water pumped for irrigation also recharges the aquifer. Natural discharge occurs as discharge to streams, evapotranspiration where the depth to water is shallow, and diffuse ground-water flow across the eastern boundary. Artificial discharge occurs as discharge to wells. Irrigation accounted for 96 percent of all use of water from the High Plains aquifer in the Oklahoma portion of the study area in 1992 and 93 percent in 1997. Total estimated water use in 1992 for the Oklahoma portion of the study area was 396,000 acre-feet and was about 3.2 million acre-feet for the entire study area. Since development of the aquifer, water levels have declined more than 100 feet in small areas of Texas County, Oklahoma, and more than 50 feet in areas of Cimarron County. Only a small area of Beaver County had declines of more than 10 feet, and Ellis County had rises of more than 10 feet. A flow model constructed using the MODFLOW computer code had 21,073 active cells in one layer and had a 6,000- foot grid in both the north-south and east-west directions. The model was used to simulate the period before major development of the aquifer and the period of development. The model was calibrated using observed conditions available as of 1998. The predevelopment-period model integrated data or estimates on the base of aquifer, hydraulic conductivity, streambed and drain conductances, and recharge from precipitation to calculate the predevelopment altitude of the water table, discharge to the rivers and streams, and other discharges. Hydraulic conductivity, recharge, and streambed conductance were varied during calibration so that the model produced a reasonable representation of the observed water table altitude and the estimated discharge to streams. Hydraulic conductivity was reduced in the area of salt dissolution in underlying Permianage rocks. Recharge from precipitation was estimated to be 4.0 percent of precipitation in greater recharge zones and 0.37 percent in lesser recharge zones. Within Oklahoma, the mean difference between water levels simulated by the model and measured water levels at 86 observation points is -2.8 feet, the mean absolute difference is 44.1 feet, and the root mean square difference is 52.0 feet. The simulated discharge is much larger than the estimated discharge for the Beaver River, is somewhat larger for Cimarron River and Wolf Creek, and is about the same for Crooked Creek. The development-period model added specific yield, pumpage, and recharge due to irrigation and dryland cultivation to simulate the period 1946 through 1997. During calibration, estimated specific yield was reduced by 15 percent in Oklahoma east of the Cimarron-Texas County line. Simulated recharge due to irrigation ranges from 24 percent for the 1940s and 1950s to 2 percent for the 1990s. Estimated recharge due to dryland cultivation is about 3.9 percent of precipitation. The mean difference between the simulated and observed waterlevel changes from predevelopment to 1998 at 162 observation points in Oklahoma is less than 0.01 foot, the mean absolute difference is 13.1 feet, and the root mean square difference is 17.9 feet. The model simulates 7.8 cubic feet per second discharge to the Beaver River above Optima Reservoir at the end of 1997, whereas the river was actually dry to this point by this time. The model simulates a decrease in discharge to the Cimarron River and Wolf Creek that appears to be reasonable. The sensitivity of the predevelopment-period model to changes in recharge and hydraulic conductivity was tested. Simulated water levels are sensitive to changes in both recharge and hydraulic conductivity. Simulated discharge to streams is sensitive to recharge but is insensitive to hydraulic conductivity. Recharge and hydraulic conductivity are closely related with respect to water levels but are not related with respect to discharge to streams. The sensitivity of the development-period model to changes in specific yield and recharge due to dryland cultivation was tested. The model appears more sensitive to specific yield, but the two inputs were not varied over the same range in percent change. The calibrated development-period model was used to simulate water-level changes from 1998 to 2020 using mean 1996-97 pumpage. The largest simulated water-level changes in Oklahoma occur in Texas County where water levels are simulated to decline 25 to 50 additional feet over a large area. Water levels also are simulated to decline 10 to 25 additional feet in two large areas of Cimarron County, two areas in Beaver County, and one area in Ellis County, all in Oklahoma. Water levels are simulated to decline more than 100 additional feet in several areas in Kansas and were simulated to decline 50 to 100 additional feet in several areas in Texas.

Colorado, Kansas, New Mexico, Oklahoma, Texas↗

USGS2003

Residence times and nitrate transport in ground water discharging to streams in the Chesapeake Bay Watershed

One of the major water-quality problems in the Chesapeake Bay is an overabundance of nutrients from the streams and rivers that discharge to the Bay. Some of these nutrients are from nonpoint sources such as atmospheric deposition, agricultural manure and fertilizer, and septic systems. The effects of efforts to control nonpoint sources, however, can be difficult to quantify because of the lag time between changes at the land surface and the response in the base-flow (ground water) component of streams. To help resource managers understand the lag time between implementation of management practices and subsequent response in the nutrient concentrations in the base-flow component of streamflow, a study of ground-water discharge, residence time, and nitrate transport in springs throughout the Chesapeake Bay Watershed and in four smaller watersheds in selected hydrogeomorphic regions (HGMRs) was conducted. The four watersheds were in the Coastal Plain Uplands, Piedmont crystalline, Valley and Ridge carbonate, and Valley and Ridge siliciclastic HGMRs. A study of springs to estimate an apparent age of the ground water was based on analyses for concentrations of chlorofluorocarbons in water samples collected from 48 springs in the Chesapeake Bay Watershed. Results of the analysis indicate that median age for all the samples was 10 years, with the 25th percentile having an age of 7 years and the 75th percentile having an age of 13 years. Although the number of samples collected in each HGMR was limited, there did not appear to be distinct differences in the ages between the HGMRs. The ranges were similar between the major HGMRs above the Fall Line (modern to about 50 years), with only two HGMRs of small geographic extent (Piedmont carbonate and Mesozoic Lowland) having ranges of modern to about 10 years. The median values of all the HGMRs ranged from 7 to 11 years. Not enough samples were collected in the Coastal Plain for comparison. Spring samples showed slightly younger water under wet conditions than under dry conditions. The apparent age of water from wells, springs, and other ground-water discharge points in the four targeted watersheds was modern to 60 years, which was similar to the apparent ages from the spring study. In the Pocomoke River Watershed in the Coastal Plain Uplands HGMR, the apparent age of ground-water samples ranged from 0 to 60 years; the ages in the vicinity of the streams ranged from 0 to 23 years. The apparent ages of ground water in the Polecat Creek Watershed in the Piedmont crystalline HGMR ranged from 2 to 30 years. The apparent ages of water from wells in the Muddy Creek Watershed in the Valley and Ridge carbonate HGMR ranged from 10 to 20 years (except for a single sample that was 45 years). The ages in the East Mahantango Creek Watershed in the Valley and Ridge siliciclastic HGMR ranged from 0 to 50 years. The distribution in apparent age of water from wells in the targeted watersheds, however, generally is older than that for water from the springs. The median age of water from wells in the Muddy Creek Watershed, for example, was 15 years, compared to 11 years for the water from the springs in that watershed, and less than 10 years for water from all springs in the spring study. The similarity in the ranges in apparent age of water from the wells and from the springs shows that the samples from the targeted watersheds and springs have bracketed the range of apparent ages that would be expected in the shallow ground-water-flow systems throughout the Chesapeake Bay Watershed. The apparent age of water from individual wells does not necessarily represent the entire distribution of ages of the discharging ground water, and it is this distribution of ages that affects the response of nutrient concentrations in stream base flow. Nutrient-reduction scenarios were modeled for two watersheds for which the distribution of apparent ground-water ages was available, the East Mahantango Creek Watershed in the Valley and Ridge siliciclastic HGMR and the Locust Grove Watershed in the Coastal Plain Uplands HGMR. A nutrient-reduction scenario was created for East Mahantango Creek, where the average residence time was determined to be approximately 10 years on the basis of the output of particle tracking from a ground-water-ﬂow model. This scenario showed decreases of nearly 50 percent in base-ﬂow concentrations of nitrate in streams within the ﬁrst year after the reduction in nitrogen input; smaller reductions in nitrate concentration occurred in each subsequent year. A second scenario for that same watershed, in which the same 10-year average residence time was assumed and an exponential model was used for analysis, showed that a 50-percent reduction in base-ﬂow concentrations of nitrate could take up to 5 years. For the Locust Grove Watershed, in which an average residence time of 32 years was assumed, simulation with the exponential model showed that it may take more than 20 years to achieve a 50-percent reduction in base-ﬂow concentra-tions of nitrate. Although it was not possible to construct such scenarios for all watersheds, these examples show the range of possible responses to changes in nutrient inputs in two very different types of watersheds. Findings from this study include information on factors that affect ground-water age, spatial distribution of ages, and nitrogen transport. In the East Mahantango Creek Watershed and the Polecat Creek Watershed, the residence time varied spatially depending on the position of the ﬂow path, and temporally depending on the recharge conditions. Generally, ground water in areas near the stream had short residence times and the water in upland areas had longer residence times. Water traveling through deep layers had longer residence times than water traveling through shallow layers, and residence times were faster under high recharge conditions than low recharge conditions. Ground water in the Pocomoke Watershed exhibits a similar pattern: younger water discharges to small order streams in headwater basins and older water discharges to larger streams near the basin outlet. Factors affecting nitrogen transport in ground water include spatial and temporal variation in input sources, ground-water age, and aquifer processes that lead to denitriﬁcation. Spatial and temporal variations in nitrogen sources affect all the watersheds. Tributaries with higher inputs of nitrogen have higher concentrations in stream base ﬂow. Areas where nitrogen application rates have increased over time show an age-nitrate relation in ground-water samples. The age-nitrate relation can be affected by denitriﬁcation, which occurs in Pocomoke and East Mahantango Creeks but is not evident in Polecat and Muddy Creeks. In East Mahantango Creek, the level of denitriﬁcation is signiﬁcant in water with residence times greater than 20 years, but because this is a small component of overall ground-water discharge to a stream, it may not remove a signiﬁcant quantity of nitrogen from the system. Denitriﬁcation in Pocomoke Creek is signiﬁcant and appears to affect mostly older water discharging to streams. Therefore, if most of the nitrogen entering these two streams is associated with the discharge of younger ground water, denitriﬁcation may not greatly affect the overall nitrogen delivery to these streams. Other ﬁndings of this study show that nitrate in ground water discharging along preferential ﬂow paths may not be affected by natural processes, such as denitriﬁcation or uptake by riparian vegetation. Seeps to swales and ditches beneath the north uplands at Polecat Creek indicate a shallow water table and discharge of young ground water whereas the absence of such seeps on the south side indicates a deep water table and a lack of young ground water. Similarly, discharge at the base of the slope and to the valley wetland south of the creek but not north of the creek indicates a different role for the riparian forest on the two sides of the creek. In many of the systems where water discharges at the base of slopes to wetlands, ditches have been dug to drain the valley. Such drainage circumvents possible removal of nitrate by riparian vegetation. Because ground-water residence times do not appear directly related to the HGMRs, the targeting of management practices will achieve the most rapid response in water quality if directed at 1) watersheds with large agricultural sources of nitrate, 2) areas with the shortest ground-water-ﬂow paths and 3) areas not affected by signiﬁcant denitriﬁcation. The fastest response in stream base-ﬂow concentrations of nitrogen to implementation of management practices would be to implement practices in those areas with the highest loads rather than attempt to target practices on the basis of HGMR stratiﬁcation. Overall ﬁndings of the study indicate that 1) ground-water contributions to nitrogen in streamﬂow are signiﬁcant, 2) some response to management practices should be evident in base-ﬂow concentrations of nitrogen and loads within 1 to 5 years in watersheds with the shortest average residence times, but response time may be closer to 20 years in watersheds with longer average ground-water residence times, 3) the majority of the response in ground-water discharge to any changes in management practices will be distributed over a 10-year time period even in the watersheds with the fastest response times, and 4) given that half the streamﬂow is from ground-water discharge and the other half is runoff or soil water, about 90 percent of total water being discharged to a stream will be less than about a decade old; therefore, full implementation of nutrient reductions may result in improved streamwater quality in about a decade. In the more-likely scenario of gradual source reduction, the reduction in concentrations of nitrate in streams and aquifers would take longer than the examples shown here.

Chesapeake Bay watershed↗

USGS2015

Chloride concentrations, loads, and yields in four watersheds along Interstate 95, southeastern Connecticut, 2008-11: factors that affect peak chloride concentrations during winter storms

Chloride (Cl - ) concentrations and loads and other water chemistry characteristics were assessed to evaluate potential effects of road-deicer applications on streamwater quality in four watersheds along Interstate 95 (I–95) in southeastern Connecticut from November 1, 2008, through September 30, 2011. Streamflow and water quality were studied in the Four Mile River, Oil Mill Brook, Stony Brook, and Jordan Brook watersheds, where developed land ranged from 9 to 32 percent. Water-quality samples were collected and specific conductance was measured continuously at paired water-quality monitoring sites, upstream and downstream from I–95. Specific conductance values were related to Cl - concentrations to assist in determining the effects of road-deicing operations on the levels of Cl - in the streams. Streamflow and water-quality data were compared with weather data and with the timing, amount, and composition of deicers applied to State highways. Grab samples were collected during winter stormwater-runoff events, such as winter storms or periods of rain or warm temperatures in which melting takes place. Grab samples were also collected periodically during the spring and summer and during base-flow conditions. The estimated Cl - concentrations at the eight water-quality monitoring sites during winter storms peaked as high as 270 milligrams per liter (mg/L) and were well below the U.S. Environmental Protection Agency (EPA) recommended acute chloride toxicity criterion of 860 mg/L and the chronic 4-day average toxicity criterion of 230 mg/L. Cl - concentrations in streams, particularly at sites downstream from I–95, peaked during increased streamflow in the winter and early spring as a result of deicers applied to roads and washed off by stormwater or meltwater. Cl - concentrations during most of the nonwinter seasons decreased during increases in streamflow because storm runoff was more dilute than base flow. However, peaks in specific conductance and estimated chloride concentrations at streams in more urbanized watersheds corresponded to peaks in streamflow well after winter snow or ice events; these delayed peaks in Cl - concentration likely resulted from deicer residue that remained in melting snow piles and on roadsides and (or) that were flushed from soils and shallow groundwater, then discharged downstream. Estimated peak Cl - concentrations varied with the type of winter storm event and were highest during or after winter storms of frozen precipitation and rain, in which the rain or meltwater effectively washed off the deicers. Estimated peak Cl - concentrations correlated positively with the duration of deicer application but generally not with streamflow. Estimated peak Cl - concentrations during the winter season were highest during low streamflow at most sites. Chloride concentrations varied considerably in shallow groundwater as a result of land-use differences. Cl - concentrations were very high (as high as 800 mg/L) in shallow groundwater downstream from I–95 at the Four Mile River site. Chloride/bromide mass concentration ratios and the proximity of a former landfill and sewage lagoon upstream indicate a likely source of Cl - is landfill leachate and possibly sewage leachate. Cl - loads in streams generally were highest in the winter and early spring. The estimated daily Cl - yield for the four monitoring sites downstream from I–95 ranged from 0.0004 ton per day per square mile for one of the least developed watersheds to 0.052 ton per day per square mile for the watershed with the highest percentage of urban development and impervious area. The estimated median contribution of Cl - load from atmospheric deposition was small and ranged from 0.07 percent of Cl - load at the Jordan Brook watershed to 0.57 percent at the Oil Mill Brook watershed. The Cl - loads in streams (outputs) were compared with Cl - load inputs, which include atmospheric deposition and deicer applications; Cl - load inputs were slightly larger than the Cl - load outputs at most of the sites during most years but do not account for the Cl - load in groundwater leaving the watersheds. A multiple linear regression model was developed to describe the variability of the natural log of peak specific conductance, as well as estimated Cl - concentrations. Five significant variables best explained the variability in the natural log of the peak specific conductance after deicing events: (1) snow on ground before deicing event; (2) winter precipitation with rain; (3) specific conductance in base flow; (4) State-operated road lane miles divided by watershed area; and (5) amount of Cl - from deicers applied to State-operated roads per lane mile. In this report, winter precipitation is defined as any type of precipitation, including frozen precipitation and rain, that occurs during the active deicing season, typically November through March. Frozen precipitation is defined here as snow, sleet, freezing rain, or any winter precipitation except rain. The addition of a lane mile in both directions on I–95 would result in an estimate of approximately 2 to 11 percent increase in Cl - input from deicers applied to I–95 and other roads maintained by Connecticut Department of Transportation. The largest estimated increase in Cl - load was in the watersheds with the greatest number miles of I–95 corridor relative to the total lane miles maintained by Connecticut Department of Transportation. On the basis of these estimates and the estimated peak Cl - concentrations during the study period, it is unlikely that the increased use of deicers on the additional lanes would lead to Cl - concentrations that exceed the aquatic habitat criteria.

Connecticut↗

USGS2005

Response curves for phosphorus plume lengths from reactive-solute-transport simulations of onland disposal of wastewater in noncarbonate sand and gravel aquifers

Surface-water resources in Massachusetts often are affected by eutrophication, excessive plant growth, which has resulted in impaired use for a majority of the freshwater ponds and lakes and a substantial number of river-miles in the State. Because supply of phosphorus usually is limiting to plant growth in freshwater systems, control of phosphorus input to surface waters is critical to solving the impairment problem. Wastewater is a substantial source of phosphorus for surface water, and removal of phosphorus before disposal may be necessary. Wastewater disposed onland by infiltration loses phosphorus from the dissolved phase during transport through the subsurface and may be an effective disposal method; quantification of the phosphorus loss can be simulated to determine disposal feasibility. In 2003, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, initiated a project to simulate distance of phosphorus transport in the subsurface for plausible conditions of onland wastewater disposal and subsurface properties. A coupled one-dimensional unsaturated-zone and three-dimensional saturated-zone reactive-solute-transport model (PHAST) was used to simulate lengths of phosphorus plumes. Knowledge of phosphorus plume length could facilitate estimates of setback distances for wastewater-infiltration sites from surface water that would be sufficient to protect the surface water from eutrophication caused by phosphorus transport through the subsurface and ultimate discharge to surface water. The reactive-solute-transport model PHAST was used to simulate ground-water flow, solute transport, equilibrium chemistry for dissolved and sorbed species, and kinetic regulation of organic carbon decomposition and phosphate mineral formation. The phosphorus plume length was defined for the simulations as the maximum extent of the contour for the 0.015 milligram-per-liter concentration of dissolved phosphorus downgradient from the infiltration bed after disposal cessation. Duration of disposal before cessation was assumed to be 50 years into an infiltration bed of 20,000 square feet at the rate of 3 gallons per square foot per day. Time for the maximum extent of the phosphorus plume to develop is on the order of 100 years after disposal cessation. Simulations indicated that phosphorus transport beyond the extent of the 0.015 milligram-per-liter concentration contour was never more than 0.18 kilogram per year, an amount that would likely not alter the ecology of most surface water. Simulations of phosphorus plume lengths were summarized in a series of response curves. Simulated plume lengths ranged from 200 feet for low phosphorus-concentration effluents (0.25 milligram per liter) and thick (50 feet) unsaturated zones to 3,400 feet for high phosphorus-concentration effluents (14 milligrams per liter) discharged directly into the aquifer (unsaturated-zone thickness of 0 feet). Plume length was nearly independent of unsaturated-zone thickness at phosphorus concentrations in the wastewater that were less than 2 milligrams per liter because little or no phosphorus mineral formed at low phosphorus concentrations. For effluents of high phosphorus concentration, plume length varied from 3,400 feet for unsaturated-zone thickness of 0 to 2,550 feet for unsaturated-zone thickness of 50 feet. Model treatments of flow and equilibrium-controlled chemistry likely were more accurate than rates of kinetically controlled reactions, notably precipitation of iron-phosphate minerals; the kinetics of such reactions are less well known and thus less well defined in the model. Sensitivity analysis indicated that many chemical and physical aquifer properties, such as hydraulic gradient and model width, did not affect the simulated plume length appreciably, but duration of discharge, size of infiltration bed, amount of dispersion, and number of sorption sites on the aquifer sediments did affect plume length appreciably. Because simulation of plume length in carbonate-mineral sediments indicated that the plume would be substantially longer than in noncarbonate-mineral sediments, the application of the response curves in locations with carbonate-mineral sediments would be inappropriate. The effect of carbonate minerals in sediments is to increase pH, which causes decreased sorption of phosphorus on aquifer sediments. Phosphorus removal from solution by precipitation onto aquifer sediments is more efficient at high concentrations of disposed phosphorus than at low concentrations. At very low phosphorus concentrations, the solubility product of phosphorus minerals is not exceeded and no phosphorus mineral forms. An important consequence is that removal of dissolved phosphorus from the plume by processes in the subsurface is decreased the more that removal efforts are applied in treatment before wastewater is disposed. Model simulations indicate that removal of phosphorus from wastewater disposed through septic systems would have the advantage of efficient phosphorus removal in the subsurface because phosphorus concentrations are high in septic-system effluent. Short plume lengths result from wastewater disposal through septic systems because of the efficient phosphorus removal and because of the low volume of wastewater involved. The simulation results for small-volume systems are not quantitative, however, because wastewater-infiltration rates are much lower than those of the higher-volume system that was used to calibrate the model and to create the plume-length response curves. The response curves for phosphorus plume lengths, as defined by the maximum extent of the 0.015 milligram-per-liter concentration contour, is clearly defined in the model simulations, although the relation between simulated plume length and protective setback distance is subject to interpretation. Phosphorus does move beyond the point at which the simulated 0.015 milligram-per-liter concentration contour has stopped, so that a determination of protective setback distance must include a consideration of whether that continued flux, or some other flux amount, is appropriate. Also, simulations indicate that phosphorus plumes do not reach their full extent until 50 to 200 years after disposal cessation, depending on concentration of phosphorus disposed. No phosphorus plume has been monitored for that long after cessation, so there is no way to verify the long-term simulation results.

Massachusetts↗

USGS2005

Evaluation of ground-water flow and land-surface subsidence caused by hypothetical withdrawals in the northern part of the Gulf Coast Aquifer system, Texas

During 2003–04 the U.S. Geological Survey, in cooperation with the Texas Water Development Board (TWDB) and the Harris-Galveston Coastal Subsidence District (HGCSD), used the previously developed Northern Gulf Coast Ground-Water Availability Modeling (NGC GAM) model to evaluate the effects of hypothetical projected withdrawals on ground-water flow in the northern part of the Gulf Coast aquifer system and land-surface subsidence in the NGC GAM model area of Texas. The Gulf Coast aquifer system comprises, from the surface, the Chicot and Evangeline aquifers, the Burkeville confining unit, the Jasper aquifer, and the Catahoula confining unit. Two withdrawal scenarios were simulated. The first scenario comprises historical withdrawals from the aquifer system for 1891–2000 and hypothetical projected withdrawals for 2001–50 compiled by the TWDB (TWDB scenario). The projected withdrawals compiled by the TWDB are based on ground-water demands estimated by regional water planning groups. The second scenario is a “merge” of the TWDB scenario with an alternate set of projected withdrawals from the Chicot and Evangeline aquifers in the Houston metropolitan area for 1995–2030 provided by the HGCSD (HGCSD scenario). Under the TWDB scenario withdrawals from the entire system are projected to be about the same in 2050 as in 2000. The simulated potentiometric surfaces of the Chicot aquifer for 2010, 2020, 2030, 2040, and 2050 show relatively little change in configuration from the simulated 2000 potentiometric surface (maximum water-level depths in southern Harris County 150–200 feet below NGVD 29). The simulated decadal potentiometric surfaces of the Evangeline aquifer show the most change between 2000 and 2010. The area of water levels 250– 400 feet below NGVD 29 in western Harris County in 2000 shifts southeastward to southern Harris County, and water levels recover to 200–250 feet below NGVD 29 by 2010. Water levels in southern Harris County recover to 150–200 feet below NGVD 29 by 2020 and remain in that range through 2050. A relatively small cone of depression in southern Montgomery County that did not appear in the 2000 surface develops and enlarges during the projected period, with a maximum depth of 250–300 feet below NGVD 29 in 2030, 2040, and 2050. The simulated decadal potentiometric surfaces of the Jasper aquifer each have a major cone of depression centered in southern Montgomery County that was minimally developed in 2000 but reaches depths of 550–650 feet below NGVD 29 in the 2020, 2030, 2040, and 2050 surfaces. Under the TWDB scenario the percentage of withdrawals supplied by net recharge increases from 75 percent in 2000 to 87 percent in 2050, and the percentage of withdrawals supplied by storage decreases from 25 percent in 2000 to 13 percent in 2050. Under the HGCSD scenario, withdrawals from the Chicot and Evangeline aquifers increase about 74 percent during 1995–2030; Jasper aquifer withdrawals are unchanged from those of the TWDB scenario. For the 2010, 2020, and 2030 potentiometric surfaces of the Chicot and Evangeline aquifers, the substantially greater withdrawals of the HGCSD scenario relative to those of the TWDB scenario result in progressively deeper cones of depression than those in the potentiometric surfaces associated with the TWDB scenario—for the Chicot aquifer in southern Harris County, 400–450 feet below NGVD 29 in 2030; for the Evangeline aquifer in southern Montgomery County, 700–750 feet below NGVD 29 in 2030. Although Jasper aquifer withdrawals are the same for both scenarios, the major cone of depression centered in southern Montgomery County in the 2030 potentiometric surface is 50 feet deeper at its center (600–700 feet below NGVD 29) than the cone in the 2030 surface under the TWDB scenario. Under the HGCSD scenario, the percentage of withdrawals supplied by net recharge decreases from 72 percent in 1995 to 57 percent in 2030, and the percentage of withdrawals supplied by storage increases from 28 percent in 2000 to 43 percent in 2030. About 85 percent of the increase supplied by storage is from the compaction of clay. Land-surface subsidence in the major area of subsidence centered in Harris and Galveston Counties during 2000–50 that results from simulating the TWDB withdrawal scenario expands slightly to the west and increases in places. The maximum change occurs in the Conroe area where subsidence increases from about 4 to about 13 feet during the projected period. Land-surface subsidence in the major area of subsidence during 1995–2030 that results from simulating the HGCSD withdrawal scenario increases substantially. For example, in east-central Harris County maximum subsidence increases from about 10–11 feet in 1995 to 22 feet in 2030. The hypothetical projected withdrawal scenarios are estimates of future withdrawals and might not represent actual future withdrawals. The simplifying assumptions that the downdip limit of freshwater flow in each hydrogeologic unit is a stable, sharp interface across which no flow occurs and that the base of the system is a no-flow boundary become less realistic and thus increase the uncertainty in results as drawdowns increase. The presence of uncertainty dictates that the results of the predictive simulations described in this report be used with caution in any decision-making process.

Texas↗

USGS2024

Nutrient chemistry in the Elizabeth Lake subwatershed—Effects of onsite wastewater treatment systems on groundwater and lake water quality, Los Angeles County, California

Nutrient (nitrogen [N] and phosphorus [P] chemistry) downgradient from onsite wastewater treatment system (OWTS) was evaluated with a groundwater study in the area surrounding Elizabeth Lake, the largest of three sag lakes within the Santa Clara River watershed of Los Angeles County, California. Elizabeth Lake is listed on the “303 (d) Impaired Waters List” for excess nutrients and is downgradient from more than 600 OWTS. The primary objective of this study was to develop a conceptual hydrogeological model to determine if discharge from OWTS is transported into shallow groundwater within the Elizabeth Lake subwatershed and contributes nutrients to Elizabeth Lake in excess of the total maximum daily load limit. An analysis of historical data and data collected for this study provided estimates of aquifer properties, such as hydraulic gradients and other parameters necessary to estimate boundary conditions. Electrical resistivity tomography (ERT) surveys were done to determine the best monitoring well locations and to estimate depth to groundwater. During 4 separate sampling events, 11 wells, 2 imported water tanks, 1 spring (sampled on March 17, 2019), and Elizabeth Lake were sampled, which occurred during February–September 2020. ERT transects and borehole geophysical measurements indicated that there were low to high resistivity materials in the subsurface and potential perched fresh water. Most of the aquifer material was characterized as sandy silt, occasionally with mixed clays and medium gravels, and was estimated to have a hydraulic conductivity from 3.28x10 −3 to 16.4 feet per day, a porosity from 0.34 to 0.42, and a hydraulic gradient from 0.01 to 0.03. Although bedrock was not obvious in ERT transects, all well depths were terminated at depths of an impassible confining layer observed to be a highly consolidated blue-gray clay. Depths to granitic bedrock, based on road outcrops and lithologic driller logs, varied throughout the study area. Depth to the bedrock was estimated to be shallow on the north side of Elizabeth Lake at approximately 30 feet below land surface (ft bls). Depth to bedrock is at 50 ft bls toward the east of the Elizabeth Lake subwatershed, which is at topographic ground surface to the north and south of the residential development. Groundwater levels ranged from approximately 0 to 12 ft bls during this study. Historical water levels ranged from 8 to 16 ft bls in the lower elevation of the study area and increased to depths of as much as 80 ft bls at higher elevations on the north and south boundaries of the Elizabeth Lake subwatershed. Water-quality samples were analyzed for major ions, nutrients, dissolved organic carbon, stable isotopes, and age-dating tracers. A principal component analysis was completed to determine organic matter sources. The proportion of recharge from imported waters, used for domestic consumption, was calculated using stable water isotopes, deuterium (δD) and oxygen (δ 18 O). Recharge from imported waters accounted for approximately 15–71 percent of the total recharge to groundwater within the study area. Total nitrogen concentrations ranged from 0.17 to 30.9 milligrams per liter (mg/L) as N, and phosphorus, measured in the soluble form as orthophosphate, ranged from 0.03 to 0.35 mg/L as P. Nitrate concentrations in groundwater samples ranged from less than the detection limit (0.01 mg/L as N) to approximately 24 mg/L as N. Nitrate was not detected in 3 of the 12 sites sampled during the study (2 wells and Elizabeth Lake). Dissolved organic carbon concentrations ranged from 0.4 to 27 mg/L in groundwater and from 9.9 to 100 mg/L in Elizabeth Lake. Ammonium and orthophosphate concentrations generally were low in groundwater. However, elevated concentrations of ammonium in Elizabeth Lake were assumed to be due to avian waste products or biological nitrogen fixation. Groundwater ages were mostly modern (recharged since 1952), with a median recharge temperature of 13 degrees Celsius. Redox conditions in groundwater indicated the likely occurrence of nitrate attenuation by denitrification downgradient from the wells to the south of Elizabeth Lake before groundwater discharges to the lake. Undetectable nitrate in Elizabeth Lake at the time of sampling was likely due to algal uptake. Most wells contained stable isotopes of nitrogen and oxygen in nitrate (δ 15 N-NO 3 and δ 18 O-NO 3 ) molecules with values consistent with denitrification. However, one monitoring well on the north of Elizabeth Lake (ELLA-8) had no evidence of denitrification, based on elevated concentrations of nitrate and a sufficient amount of dissolved oxygen such that the water was oxic and not favorable for the denitrification reaction. Consequently, this nitrate could be delivered to Elizabeth Lake through groundwater discharge if nitrate is not removed from the system by denitrifying bacteria downgradient from the well before the groundwater discharges into Elizabeth Lake. The principal component analysis demonstrated that dissolved organic matter optical properties track different sources of dissolved organic matter from decayed plants, animals, and animal-derived wastes. Two wells contained strong indicators of OWTS water presence, although geochemical evidence indicated other wells may also be affected by OWTS discharge.

California↗

USGS2000

Geohydrology and numerical simulation of groundwater flow in the central Virgin River Basin of Iron and Washington Counties, Utah

Because rapid growth of communities in Washington and Iron Counties, Utah, is expected to cause an increase in the future demand for water resources, a hydrologic investigation was done to better understand ground-water resources within the central Virgin River basin. This study focused on two of the principal ground-water reservoirs within the basin: the upper Ash Creek basin ground-water system and the Navajo and Kayenta aquifer system. The ground-water system of the upper Ash Creek drainage basin consists of three aquifers: the uppermost Quaternary basin-fill aquifer, the Tertiary alluvial-fan aquifer, and the Tertiary Pine Valley monzonite aquifer. These aquifers are naturally bounded by the Hurricane Fault and by drainage divides. On the basis of measurements, estimates, and numerical simulations of reasonable values for all inflow and outflow components, total water moving through the upper Ash Creek drainage basin ground-water system is estimated to be about 14,000 acre-feet per year. Recharge to the upper Ash Creek drainage basin ground-water system is mostly from infiltration of precipitation and seepage from ephemeral and perennial streams. The primary source of discharge is assumed to be evapotranspiration; however, subsurface discharge near Ash Creek Reservoir also may be important. The character of two of the hydrologic boundaries of the upper Ash Creek drainage basin ground-water system is speculative. The eastern boundary provided by the Hurricane Fault is assumed to be a no-flow boundary, and a substantial part of the ground-water discharge from the system is assumed to be subsurface outflow beneath Ash Creek Reservoir along the southern boundary. However, these assumptions might be incorrect because alternative numerical simulations that used different boundary conditions also proved to be feasible. The hydrogeologic character of the aquifers is uncertain because of limited data. Differences in well yield indicate that there is considerable variability in the transmissivity of the basin-fill aquifer. Field data also indicate that the basin-fill aquifer is more transmissive than the underlying alluvial-fan aquifer. Data from the Pine Valley monzonite aquifer indicate that its transmissivity may be highly variable and that it is strongly influenced by the connection of fractures. The Navajo and Kayenta aquifers provide most of the potable water to the municipalities of Washington County. Because of large outcrop exposures, uniform grain size, and large stratigraphic thickness, these formations are able to receive and store large amounts of water. In addition, structural forces have resulted in extensive fracture zones that enhance ground-water recharge and movement within these aquifers. Aquifer testing of the Navajo aquifer indicates that horizontal hydraulic-conductivity values range from 0.2 to 32 feet per day at different locations and may be primarily dependent on the extent of fracturing. Limited data indicate that the Kayenta aquifer generally is less transmissive than the Navajo aquifer. The aquifers are bounded to the south and west by the erosional extent of the formations and to the east by the Hurricane Fault, which completely offsets these formations and is assumed to be a lateral no-flow boundary. Like the Hurricane Fault, the Gunlock Fault is assumed to be a lateral no-flow boundary that divides the Navajo and Kayenta aquifers within the study area into two parts: the main part, between the Hurricane and Gunlock Faults; and the Gunlock part, west of the Gunlock Fault. Generally, the water in the Navajo and Kayenta aquifers contains few dissolved minerals. However, two distinct areas contain water with dissolved-solids concentrations greater than 500 milligrams per liter: a larger area north of the city of St. George and a smaller area a few miles west of the town of Hurricane. Mass-balance calculations indicate that in the higher-dissolved-solids area north of St. George, as much as 2.7 cubic feet per second may be entering the aquifer from underlying formations. For the area west of Hurricane, as much as 1.5 cubic feet per second may be entering the aquifer from underlying formations. On the basis of measurements, estimates, and numerical simulations, total water moving through the Navajo and Kayenta aquifers is estimated to be about 25,000 acre-feet per year for the main part and 5,000 acre-feet per year for the Gunlock part. The primary source of recharge is assumed to be infiltration of precipitation in the main part and seepage from the Santa Clara River in the Gunlock part. The primary source of discharge is assumed to be well discharge for both the main and Gunlock parts of the aquifers. Numerical simulations indicate that faults with major offset, such as the Washington Hollow Fault and an unnamed fault near Anderson Junction, may impede horizontal ground-water flow. Also, increased horizontal hydraulic conductivity along the orientation of predominant surface fracturing may be an important factor in regional ground-water flow. Simulations with increased north-south hydraulic conductivity substantially improved the match to measured water levels in the central area of the model between Snow Canyon and Mill Creek. Numerical simulation of the Gunlock part, using aquifer properties determined for the city of St. George municipal well field, resulted in a reasonable representation of regional water levels and estimated seepage from and to the Santa Clara River. To further quantify the Gunlock part of the Navajo and Kayenta aquifers, a better understanding of ground-water flow at the Gunlock Fault is needed.

Utah↗

USGS2008

Hydrogeology and water quality of the Leetown area, West Virginia

The U.S. Geological Survey’s Leetown Science Center and the co-located U.S. Department of Agriculture’s National Center for Cool and Cold Water Aquaculture both depend on large volumes of cold clean ground water to support research operations at their facilities. Currently, ground-water demands are provided by three springs and two standby production wells used to augment supplies during periods of low spring flow. Future expansion of research operations at the Leetown Science Center is dependent on assessing the availability and quality of water to the facilities and in locating prospective sites for additional wells to augment existing water supplies. The hydrogeology of the Leetown area, West Virginia, is a structurally complex karst aquifer. Although the aquifer is a karst system, it is not typical of most highly cavernous karst systems, but is dominated by broad areas of fractured rock drained by a relatively small number of solution conduits. Characterization of the aquifer by use of fluorometric tracer tests, a common approach in most karst terranes, therefore only partly defines the hydrogeologic setting of the area. In order to fully assess the hydrogeology and water quality in the vicinity of Leetown, a multi-disciplinary approach that included both fractured rock and karst research components was needed. The U.S. Geological Survey developed this multi-disciplinary research effort to include geologic, hydrologic, geophysical, geographic, water-quality, and microbiological investigations in order to fully characterize the hydrogeology and water quality of the Leetown area, West Virginia. Detailed geologic and karst mapping provided the framework on which hydrologic investigations were based. Fracture trace and lineament analysis helped locate potential water-bearing fractures and guided installation of monitoring wells. Monitoring wells were drilled for borehole geophysical surveys, water-quality sampling, water-level measurements, and aquifer tests to characterize the quality of water and the hydraulic properties of the aquifer. Surface geophysical surveys provided a 3-dimensional view of bedrock resistivity in order to assess geologic and lithologic controls on ground-water flow. Borehole geophysical surveys were conducted in monitoring wells to assess the storage and movement of water in subsurface fractures. Numerous single-well, multi-well, and straddle packer aquifer tests and step-drawdown tests were conducted to define the hydraulic properties of the aquifer and to assess the role of bedrock fractures and solution conduits in the flow of ground water. Water samples collected from wells and springs were analyzed to assess the current quality of ground water and provide a baseline for future assessment. Microbiological sampling of wells for indicator bacteria and human and animal DNA provided an analysis of agricultural and suburban development impacts on ground-water quality. Light detection and ranging (LiDAR) data were analyzed to develop digital elevation models (DEMs) for assessing sinkhole distribution, to provide elevation data for development of a ground-water flow model, and to assess the distribution of major fractures and faults in the Leetown area. The flow of ground water in the study area is controlled by lithology and geologic structure. Bedrock, especially low permeability units such as the shale Martinsburg Formation and the Conococheague Limestone, act as barriers to water flowing down gradient and across bedding. This retardation of cross-strike flow is especially pronounced in the Leetown area, where bedding typically dips at steep angles. Highly permeable fault and fracture zones that disrupt the rocks in cross-strike directions provide avenues through which ground water can flow laterally across or through strata of low primary permeability. Significant strike parallel thrust faults and cross-strike faults typically coincide with larger solution conduits and act as drains for the more pervasive network of interconnected diffuse fractures. Results of borehole geophysical surveys indicate that although numerous fractures may intersect a borehole, only one or two of the fractures typically transmit most of the water to a well. The diffuse-flow dominated network of fractures that provides the majority of storage occupies only a small proportion of the total aquifer volume but constitutes the majority of porosity within the aquifer. Solution conduits, while occupying a relatively small volume of the overall aquifer, are especially important because they serve as primary drains for the ground-water flow system. Surface resistivity maps and cross-sectionsshow anomalous areas of low resistivities coincident with the prevailing geologic strike at N. 20º E., with major cross-strike faults, and with major springs in the region. Transmissivity derived from straddle packer tests was highly variable, and ranged over three orders of magnitude (1.8 x 10 -6 to 5.9 x 10 -3 ft 2 /d) in diffuse-flow fractures. A similar large variability in transmissivity was documented by single- and multi-well aquifer tests conducted in conduit-flow dominated portions of the aquifer (2.0 x 10 3 to 1.4 x 10 4 ft 2 /d) in lowland areas immediately adjacent to the Leetown Science Center. A stream-gaging station installed on Hopewell Run near the point where the stream exits the Leetown watershed indicates average daily streamflow for the Hopewell Run of approximately 11.2 ft 3 /s, and ranged from a minimum of 1.80 ft 3 /s on September 28, 2005, to a maximum of 73.0 ft 3 /s on December 11, 2003. Base-flow (ground-water) discharge surveys identified numerous small seeps adjacent to streams in the area. Hydrographs of the stage of Balch Spring show rapid response to individual storms. Strong correlation of the flow of Hopewell Run and Balch Spring indicates the nearby losing stream reach is partly responsible for higher fluctuations in the stage of Balch Spring. A water budget for the study period (2003-2005), based on measured precipitation and hydrograph analyses, is expressed as Precipitation (38.60 in/yr) = Surface Runoff (1.36 in/yr) + Ground-Water Discharge (17.73 in/yr) + Evapotranspiration (24.23 in/yr) – Change in storage (4.72 in/yr). Flow of ground water through the epikarst, a shallow zone of intensely weathered rock and regolith, can be rapid (on the order of days or weeks) as flow is concentrated in solution conduits. Flow within the intermediate and deeper zones is typically much slower. Eight dye-tracer tests conducted in the Leetown area found ground-water flow patterns to be divergent, with velocities ranging from about 12.5 to 610 ft/day and a median velocity of 50 ft/day. Estimates of ground-water age in carbonate rocks in the region are on the order of 15 years in the shallower portions of the aquifer to 50 years or older for deeper portions of the aquifer. Shallow springs can have a significant component of fairly young water (< 5 years in age). Ground-water samples collected from 16 sites (12 wells and 4 springs) in the Leetown area were analyzed for more than 340 constituents. Only turbidity, indicator bacteria, and radon were typically present in concentrations exceeding U.S. Environmental Protection Agency (USEPA) drinking-water or aquatic life standards.

West Virginia↗

USGS2019

Quantifying the eroded and deposited mass of mercury-contaminated sediment by using terrestrial laser scanning at the confluence of Humbug Creek and the South Yuba River, Nevada County, California, 2011–13

High-resolution, terrestrial laser scanning, also known as ground-based lidar (light detection and ranging), was used to quantify the volume of mercury-contaminated sediment eroded from an outcrop of historical placer-mining debris at the confluence of Humbug Creek and the South Yuba River in the Sierra Nevada foothills, about 17 kilometers northeast of Grass Valley, California, and delivered to a zone below an observed flood stage of the South Yuba River. Substantial quantities of mercury were used and lost to the environment from historical placer gold mining activities on the western slope of the Sierra Nevada, California, and recent studies have documented continued persistence of mercury and methylmercury concentrations in water, sediment, fish, and predatory invertebrates in the Yuba River drainage basin in relation to suspected mercury sources. To identify areas that have high levels of mercury contamination as possible remediation targets in the Yuba River drainage basin and other areas in the Sierra Nevada, the U.S. Geological Survey worked in cooperation with the Bureau of Land Management on this and other detailed studies. Malakoff Diggings, one of the largest hydraulic gold mines in the Sierra Nevada, is 3.5 kilometers north of the study site in the Humbug Creek subbasin. Terrestrial laser scanning was used to produce centimeter-scale, three-dimensional maps of the complex outcrop surface, which was composed of an upper erosional area (cliff and over-steepened slope) and a lower depositional area (colluvial slope). The outcrop could not be mapped non-destructively or in sufficient detail by traditional surveying techniques. The study site, which was approximately 70 meters long, 30 meters wide and 20 meters high, was surveyed four times in 2 years (December 15, 2011; October 25, 2012; January 4, 2013; and November 22, 2013) to determine volumetric differences in the upper erosional and lower depositional areas between surveys. Measured changes in volume for the upper erosional area and lower depositional area were multiplied by the corresponding sediment density so that a mass-balance relationship, between the eroded and deposited sediment during each period, could be used to estimate the amount of mercury-contaminated sediment that was transported to below the base of the colluvial slope, where it could be mobilized by the South Yuba River during a flood having a 5-to-10-year recurrence interval. On December 2, 2012, a flood of this estimated magnitude reached the base of the colluvial slope. Between the first and second surveys (December 15, 2011–October 25, 2012), an estimated mass of 18±9.2 kilograms of sediment was transported from steeper slopes to the gently sloping river bank below the base of the colluvial slope. Between the second and third surveys (October 25, 2012–January 4, 2013), an atmospheric river caused heavy precipitation at the study site during late November and early December 2012. This short-duration, high-intensity rain resulted in a large amount of erosion and deposition at the study site and also caused high streamflow (flood stage) in the South Yuba River. From October 2012 to January 2013, 51±31 kilograms of sediment was transported to below the base of the colluvial slope, that is, below the high-water mark of December 2, 2012. Between the third and fourth surveys (January 4, 2013–November 22, 2013), an additional 10±26 kilograms of sediment was transported to below the base of the colluvial slope. During the 24 months of the study, the total mass of sediment delivered below the base of the colluvial slope and the high-water mark of December 2, 2012, was 79±66 kilograms. In any given year there is a 10–20-percent chance (5-to-10-year recurrence interval) of a flood equal to or greater than that of the December 2, 2012, flood, which could transport mercury-contaminated sediment at the study site into the South Yuba River. Hydraulically modeled estimates of the South Yuba River stage during floods having a 50- and 100-year recurrence interval (2- and 1-percent annual exceedance probability, respectively) indicated that resulting river stages could be 2.2–3.0 meters above the base of the colluvial slope, or 2.2–3.0 meters above the high-water mark of December 2, 2012. Such high river stages would be likely to inundate the lower half of the colluvial slope and mobilize a substantial volume of mercury-contaminated sediment to downstream areas.

California↗

USGS1998

Ground-water levels, predevelopment ground-water flow, and stream-aquifer relations in the vicinity of the Savannah River Site, Georgia and South Carolina

Ground-water levels, predevelopment ground-water flow, and stream-aquifer relations in the vicinity of the U.S. Department of Energy Savannah River Site, Georgia and South Carolina, were evaluated as part of a cooperative study between the U.S. Geological Survey, U.S. Department of Energy, and Georgia Department of Natural Resources. As part of this evaluation: (1) ground-water-level fluctuations and trends in three aquifer systems in sediment of Cretaceous and Tertiary age were described and related to patterns of ground-water use and precipitations; (2) a conceptual model ofthe stream-aquifer flow system was developed; (3) the predevelopment ground-water flow system, configuration of potentiometric surfaces, trans-river flow, and recharge-discharge relations were described; and (4) stream-aquifer relations and the influence of river incision on ground-water flow and stream-aquifer relations were described. The 5,147-square mile study area is located in the northern part of the Coastal Plain physiographic province of Georgia and South Carolina. Coastal Plain sediments comprise three aquifer systems consisting of seven aquifers that are separated hydraulically by confining units. The aquifer systems are, in descending order: (1) the Floridan aquifer system—consisting of the Upper Three Runs and Gordon aquifers in sediments of Eocene age; (2) the Dublin aquifer system—consisting of the Millers Pond, upper Dublin, and lower Dublin aquifers in sediments of Paleocene-Late Cretaceous age; and (3) the Midville aquifer system—consisting of the upper Midville and lower Midville aquifers in sediments of Late Cretaceous age. The Upper Three Runs aquifer is the shallowest aquifer and is unconfined to semi-confined throughout most of the study area. Ground-water levels in the Upper Three Runs aquifer respond to a local flow system and are affected mostly by topography and climate. Ground-water flow in the deeper, Gordon aquifer and Dublin and Midville aquifer systems is characterized by local flow near outcrop areas to the north, changing to intermediate flow and then regional flow downdip (southeastward) as the aquifers become more deeply buried. Water levels in these deeper aquifers show a pronounced response to topography and climate in the vicinity of outcrops, and diminish southeastward where the aquifer is more deeply buried. Stream stage and pumpage affect ground-water levels in these deeper aquifers to varying degrees throughout the study area. The geologic characteristics of the Savannah River alluvial valley substantially control the configuration of potentiometric surfaces, ground-water-flow directions, and stream-aquifer relations. Data from 18 shallow borings indicate incision into each aquifer by the paleo Savannah River channel and subsequent infill of permeable alluvium, allowing for direct hydraulic connection between aquifers and the Savannah River along parts of its reach. This hydraulic connection may be the cause of large ground-water discharge to the river near Jackson, S.C., where the Gordon aquifer is in contact with Savannah River alluvium, and also the cause of lows or depressions formed in the potentiometric surfaces of confined aquifers that are in contact with the alluvium. Ground water in these aquifers flows toward the depressions. The influence of the river is diminished downstream where the aquifers are deeply buried, and upstream and downstream ground-water flow is possibly separated by a water divide or "saddle". Water-level data indicate that saddle features probably exist in the Gordon aquifer and Dublin aquifer system, and also might be present in the Midville aquifer system. Ground-water levels respond seasonally or in long term to changes in precipitation, evapotranspiration, pumpage, and river stage. Continuous water-level data and water-levels measured in a network of 271 wells during the Spring (May) and Fall (October) in 1992, indicate that seasonal water-level changes generally are less than 4 feet, and that larger changes occur near supply wells. Long-term data (more than 10 years of record) from 283 wells indicate that water levels in most aquifers declined during some period prior to 1993. Although most declines were less than 15 feet and were limited to scattered areas influenced by supply wells, widespread declines of as much as 59 feet occurred in the lower Midville aquifer in Richmond County, Georgia. Water-level data at well-cluster sites indicate that in the vicinity of major ground-water divides, head decreases with depth; and in the vicinity of a regional drain, such as the Savannah River, head increases with depth. Although this vertical head distribution holds true over much of the study area, water-level data indicate that the Gordon, Millers Pond, and lower Dublin aquifers are apparent hydrologic "sinks" in parts of the study area, indicating a potential for vertical leakage from both above and below. Reasons for anomalously low heads in these three aquifers are unclear, but probably are related to: (1) subsurface pinchout of the aquifer, that influences flow patterns in the ground-water flow system, (2) hydraulic connection of the aquifer to river alluvium and associated large ground-water discharge, or (3) water-level declines as a result of pumpage. Estimated average ground-water discharge to the Savannah River based on data from 1941, 1942, and 1949, was 1,220 cubic feet per second; of which 46 percent was from the local flow system, 41 percent was from the intermediate flow system, and 13 percent was from the regional flow system. Estimated ground-water discharge during the 1954 and 1986 droughts indicates that the contribution from the intermediate flow system decreased in a downstream direction, and that most ground-water discharge was from the local flow system. The decreased contribution from the intermediate flow system in downstream reaches may be related to a downdip (downstream) thickening of confining units underlying the stream, that reduces upward leakage of water. Flow lines based on contours of the estimated predevelopment potentiometric surface of the confined Dublin and Midville aquifer systems, suggest ground-water flow beneath the floodplain of the Savannah River from one side of the river to the other side. This phenomenon, termed trans-river flow, is assumed to occur for a short distance into Georgia prior to discharge into the Savannah River alluvial valley. Trans-river flow is influenced by changes in hydraulic gradient (ground-water levels) near the river. Trans-river flow cannot be inferred from contours of the potentiometric surface of the Gordon or Upper Three Runs aquifers.

Georgia, South Carolina↗

USGS2019

Water-budget analysis of the Upper Big Sandy Designated Ground-water Basin alluvial aquifer, Elbert, El Paso, and Lincoln Counties, Colorado, 2016

The U.S. Geological Survey in cooperation with the Colorado Water Conservation Board and the Upper Big Sandy Groundwater Management District carried out a study in 2016 to evaluate potential groundwater storage changes within the Upper Big Sandy Designated Groundwater Basin (UBSDGB) alluvial aquifer, including groundwater flow between the UBSDGB alluvial aquifer and the Denver Basin bedrock aquifers. The UBSDGB alluvial aquifer is located along the ephemeral Big Sandy Creek on the east-central edge of the Denver Basin aquifer system and covers an area of about 66,560 acres within the UBSDGB. The UBSDGB alluvial aquifer consists of unconsolidated Quaternary sand and gravel deposits that contain an unconfined (water table) groundwater system. The western three-fourths of the UBSDGB alluvial aquifer overlies the Tertiary and Cretaceous bedrock formations that compose the Denver Basin aquifer system. The updated water budget for the UBSDGB alluvial aquifer, including annual change in groundwater storage in 2016, was determined by combining water-budget information from an existing Denver Basin model for about three-fourths of the study area with best estimates for the major water-budget components for the area outside the Denver Basin aquifer system. The western part of the UBSDGB was included in the Denver Basin model (modeled area), whereas the eastern part of the UBSDGB was not included in the Denver Basin model (unmodeled area). The water-budget components were first estimated for the modeled area using outputs from the Denver Basin model, which uses the modular finite-difference groundwater flow computer model MODFLOW-2000 with 1-mile grid cells. For this study, the Denver Basin model was updated with additional data from 2004 through 2016 to generate current (2016) estimates of water consumption in the UBSDGB alluvial aquifer. A basin-specific water budget for the UBSDGB alluvial aquifer from the Denver Basin model was computed using a modeling tool called ZONEBUDGET. The modeled area groundwater budget, along with previous studies, was used to estimate a groundwater budget for the unmodeled area, and results for the modeled and unmodeled areas were combined for an overall water-budget estimate for the entire UBSDGB alluvial aquifer. The net groundwater flow into the basin from adjacent alluvial aquifers was positive with flow entering the UBSDGB alluvial aquifer. Combining the total inflow from adjacent alluvial and the total outflow to adjacent alluvial aquifers resulted in a net flow from adjacent alluvial aquifers to UBSDGB alluvial aquifer of 5,125 acre-feet (ac-ft) in 2016. The net flow between the underlying bedrock aquifers and the UBSDGB alluvial aquifer was positive with flow entering the UBSDGB alluvial aquifer from the bedrock aquifers. The net flow from the bedrock aquifers to the UBSDGB alluvial aquifer was 347 ac-ft in 2016. Net recharge (precipitation and irrigation return flows minus evaporation) into the UBSDGB alluvial aquifer was negative with groundwater being removed from the UBSDGB alluvial aquifer over the total area of the basin. Combining the total inflow from recharge to the UBSDGB alluvial aquifer of 11,153 ac-ft in 2016 and the total evapo-transpiration of −11,656 ac-ft from the UBSDGB alluvial aquifer in 2016 resulted in a net recharge from UBSDGB alluvial aquifer of −503 ac-ft in 2016. Combining the modeled and unmodeled well pumping resulted in a total well pumping volume of −3,735 ac-ft in 2016 from the UBSDGB alluvial aquifer. The net groundwater flow to the stream network in the basin was negative with flow discharging from the UBSDGB alluvial aquifer into streams. Combining the total inflow from streams and the total outflow to streams for the UBSDGB alluvial aquifer resulted in −1,032 ac-ft in 2016 that was lost to the stream network in the UBSDGB. The net groundwater flow out of the UBSDGB was negative with flow leaving the UBSDGB alluvial aquifer. Combining the total area inflow to the basin from upgradient areas and the total area outflow from the basin for the UBSDGB alluvial aquifer resulted in a net flow out of the basin of −2,300 ac-ft. In the annual groundwater budget for 2016, groundwater storage in the UBSDGB alluvial aquifer system was removed because annual groundwater outflows from storage exceeded groundwater inflows to storage; in other words, water was removed from storage to balance the annual water budget. Combining the net flow from storage for the modeled area of 73 ac-ft and the inflow from storage for the unmodeled area of 2,025 ac-ft resulted in a net positive flow from storage of the UBSDGB alluvial aquifer of 2,098 ac-ft. Increased pumping since 1958 in the Denver and upper Arapahoe aquifers, not necessarily in the UBSDGB, has caused a change in flow from bedrock units, which were minor or non-contributors of inflow to the UBSDGB alluvial aquifer, to receiving outflow from the UBSDGB alluvial aquifer. Since 2000, aquifer storage has been an inflow component of the water budget, which means that outflow from the modeled area exceeded inflow for the UBSDGB alluvial aquifer. Increased recharge from wetter than average years could replenish the UBSDGB alluvial aquifer. From 2003 through 2016, 13 of the 25 observation wells completed in the UBSDGB alluvial aquifer had a decline in the groundwater-level elevation with an average decline of −2.21 feet, and 12 of the 25 observation wells had an increase in the groundwater-level elevation with an average increase of 1.54 feet. In general, wells at the eastern and western edges of the UBSDGB showed an increase in groundwater-level elevation that appears related to areas of groundwater discharge from the lower Dawson and Laramie-Fox Hills bedrock aquifers to the UBSDGB alluvial aquifer. The remaining wells exhibited water-level declines. Future work could include the development of a basin-specific model to serve as a basin management tool for modeling changes in groundwater levels and storage under various future groundwater recharge and withdrawal scenarios.

Colorado↗

USGS1993

Ground-water-quality assessment of the Central Oklahoma Aquifer, Oklahoma: Geochemical and geohydrologic investigations

The National Water-Quality Assessment pilot project for the Central Oklahoma aquifer examined the chemical and isotopic composition of ground water, the abundances and textures of minerals in core samples, and water levels and hydraulic properties in the flow system to identify geochemical reactions occurring in the aquifer and rates and directions of ground-water flow. The aquifer underlies 3,000 square miles of central Oklahoma and consists of Permian red beds, including parts of the Permian Garber Sandstone, Wellington Formation, and Chase, Council Grove, and Admire Groups, and Quaternary alluvium and terrace deposits. In the part of the Garber Sandstone and Wellington Formation that is not confined by the Permian Hennessey Group, calcium, magnesium, and bicarbonate are the dominant ions in ground water; in the confined part of the Garber Sandstone and Wellington Formation and in the Chase, Council Grove, and Admire Groups, sodium and bicarbonate are the dominant ions in ground water. Nearly all of the Central Oklahoma aquifer has an oxic or post-oxic environment as indicated by the large dissolved concentrations of oxygen, nitrate, arsenic(V), chromium(VI), selenium(VI), vanadium, and uranium. Sulfidic and methanic environments are virtually absent. Petrographic textures indicate dolomite, calcite, sodic plagioclase, potassium feldspars, chlorite, rock fragments, and micas are dissolving, and iron oxides, manganese oxides, kaolinite, and quartz are precipitating. Variations in the quantity of exchangeable sodium in clays indicate that cation exchange is occurring within the aquifer. Gypsum may dissolve locally within the aquifer, as indicated by ground water with large concentra-tions of sulfate, but gypsum was not observed in core samples. Rainwater is not a major source for most elements in ground water, but evapotranspiration could cause rainwater to be a significant source of potassium, sulfate, phosphate and nitrogen species. Brines derived from seawater are the most likely source of bromide and chloride in the aquifer. The dominant reaction in recharge is the uptake of carbon dioxide gas from the unsaturated zone (about 2.0 to 4.0 millimoles per liter) and the dissolution of dolomite (about 0.3 to 1.0 millimoles per liter). This reaction generates calcium, magnesium, and bicarbonate water composition. If dolomite does not dissolve to equilibrium, pH values range from 6.0 to 7.3; if dolomite dissolves to equilibrium, pH values are about 7.5. By the time recharge enters the deeper flow system, all ground water is saturated or supersaturated with dolomite and calcite. After carbonate-mineral equilibration has occurred, cation exchange of calcium and magnesium for sodium is the dominant geochemical reaction, which occurs to a substantial extent only in parts of the aquifer. Mass transfers of cation exchange greater than 2.0 millimoles per liter occur in the confined part of the Garber Sandstone and Wellington Formation and in parts of the Chase, Council Grove, and Admire Groups. Associated with cation exchange is dissolution of small quantities of dolomite, calcite, biotite, chlorite, plagioclase, or potassium feldspar, which produces pH values that range from 8.6 to 9.1. Large tritium concentrations indicate ground-water ages of less than about 40 years for most samples of recharge. Carbon-14 ages for samples from the unconfined aquifer generally are less than 10,000 years. Carbon-14 ages of ground water in the confined part of the aquifer range from about 10,000 to 30,000 years or older. These ages produce a time trend in deuterium values that qualitatively is consistent with the timing of the transition from the last glacial maximum to the present interglacial period. The most transmissive geologic units in the Central Oklahoma aquifer are the Garber Sandstone and Wellington Formation and the alluvium and terrace deposits; the Chase, Council Grove, and Admire Groups are less transmissive on the basis of available specific-capacity data. The transmissivities of the Permian geologic units depend largely on the percentage of sandstone; the percentage is greatest in the central part of the aquifer and decreases in all directions from this central part. Because of large mudstone and siltstone contents, the Hennessey Group and the Vanoss Formation are assumed to be confining units above and below the aquifer. The Cimarron and Canadian Rivers are defined to be the northern and southern extent of the aquifer because of decreases in transmissivity beyond the rivers and because there is no indication of ground-water underflow at these rivers. The eastern boundary of the aquifer is the limit of the outcrop of the Chase, Council Grove, and Admire Groups. The presence of brines in the western part of the study unit and below the aquifer indicate the extent of the freshwater flow system in these directions. Regional ground-water flow is west to east; the Deep Fork is a major discharge area for the regional flow system. Local flow systems are present within the unconfined part of the study unit. Most streams are gaining streams, and very few losing streams are evident. Median values of aquifer properties were estimated as follows: recharge to the saturated zone, 1.6 inches per year; evapotranspiration of water that never reaches the saturated zone, 25 to 30 inches per year; porosity, 0.22; storage coefficient, 0.0002; transmissivity, 260 to 450 feet squared per day; horizontal hydraulic conductivity, 4.5 feet per day; and the ratio of horizontal to vertical hydraulic conductivity, 10,000. Reported ground-water withdrawals peaked in 1985 at 13,900 million gallons but had decreased to 7,860 million gallons by 1989. Unreported domestic withdrawals were estimated to be 1,685 million gallons in 1980. The flow system in the aquifer can be considered to have three major components: (1) A shallow, local flow system in the unconfined part of the aquifer, (2) a deep, regional flow system in the unconfined part of the aquifer, and (3) a deep, regional flow system in the confined part of the aquifer. In the shallow, local flow system, water flows relatively quickly along short flowlines from the point of recharge to the point of discharge at the nearest stream. Many water samples from shallow wells contain large concentrations of tritium, which indicate ground-water ages of less than 40 years. In the deep, regional flow system in the unconfined part of the aquifer, water takes more time to flow along longer flowlines than in the shallow, local flow system. Much of the water in this flow system is recharged along ridges that correspond to ground-water divides between drainage basins. Transit times for water recharging the aquifer along ridges is greater than 5,000 years, computed using a numerical flow model in conjunction with a particle-tracking model. The deep, regional flow system in the confined part of the Garber Sandstone and Wellington Formation is recharged from a small part of the outcrop area of the Garber Sandstone. From the recharge area, water flows west under the confining unit to discharge to streams as far away as the Cimarron River. Flowpaths are relatively long, as much as 50 miles. The transit times in this flow system range from thousands to tens of thousands of years. The long-term hydrogeochemical process occurring in the Central Oklahoma aquifer is removal of unstable minerals, including dolomite, calcite, biotite, chlorite, and feldspars, and the replacement of exchangeable sodium on clays with calcium and magnesium. Over geologic time, the flux of water through the rapidly moving, local flow system has been sufficient to remove most of the dolomite, calcite, and exchangeable sodium. In places, chlorite and feldspars have been removed. In the deep, regional flow system of the unconfined part of the Garber Sandstone and Wellington Formation, the flux of water has been sufficient to remove most of the exchangeable sodium, but carbonate minerals remain sufficiently abundant to maintain dolomite and calcite equilibrium. In the confined part of the Garber Sandstone and Wellington Formation and in the less transmissive parts of the unconfined aquifer, including the Chase, Council Grove, and Admire Groups, ground-water flow is slowest, and the flux of water and extent of reaction have been insufficient to remove either the carbonate minerals or the exchangeable sodium on clays.

Oklahoma↗

USGS2007

Recharge area, base-flow and quick-flow discharge rates and ages, and general water quality of Big Spring in Carter County, Missouri, 2000-04

Exploration for lead deposits has occurred in a mature karst area of southeast Missouri that is highly valued for its scenic beauty and recreational opportunities. The area contains the two largest springs in Missouri (Big Spring and Greer Spring), both of which flow into federally designated scenic rivers. Concerns about potential mining effects on the area ground water and aquatic biota prompted an investigation of Big Spring. Water-level measurements made during 2000 helped define the recharge area of Big Spring, Greer Spring, Mammoth Spring, and Boze Mill Spring. The data infer two distinct potentiometric surfaces. The shallow potentiometric surface, where the depth-to-water is less than about 250 feet, tends to mimic topographic features and is strongly controlled by streams. The deep potentiometric surface, where the depth-to-water is greater than about 250 feet represents ground-water hydraulic heads within the more mature karst areas. A highly permeable zone extends about 20 mile west of Big Spring toward the upper Hurricane Creek Basin. Deeper flowing water in the Big Spring recharge area is directed toward this permeable zone. The estimated sizes of the spring recharge areas are 426 square miles for Big Spring, 352 square miles for Greer Spring, 290 square miles for Mammoth Spring, and 54 square miles for Boze Mill Spring. A discharge accumulation curve using Big Spring daily mean discharge data shows no substantial change in the discharge pattern of Big Spring during the period of record (water years 1922 through 2004). The extended periods when the spring flow deviated from the trend line can be attributed to prolonged departures from normal precipitation. The maximum possible instantaneous flow from Big Spring has not been adequately defined because of backwater effects from the Current River during high-flow conditions. Physical constraints within the spring conduit system may restrict its maximum flow. The largest discharge measured at Big Spring during the period of record (water years 1922 through 2004) was 1,170 cubic feet per second on December 7, 1982. The daily mean water temperature of Big Spring was monitored during water years 2001 through 2004 and showed little variability, ranging from 13 to 15° C (degree Celsius). Water temperatures generally vary less than 1° C throughout the year. The warmest temperatures occur during October and November and decrease until April, indicating Big Spring water temperature does show a slight seasonal variation. The use of the traditional hydrograph separation program HYSEP to determine the base flow and quick flow or runoff components at Big Spring failed to yield base-flow and quick-flow discharge curves that matched observations of spring characteristics. Big Spring discharge data were used in combination with specific conductance data to develop an improved hydrograph separation method for the spring. The estimated annual mean quick flow ranged from 15 to 48 cubic feet per second for the HYSEP analysis and ranged from 26 to 154 cubic feet per second for the discharge and specific conductance method for water years 2001 to 2004. Using the discharge and specific conductance method, the estimated base-flow component rises abruptly as the spring hydrograph rises, attains a peak value on the same day as the discharge peak, and then declines abruptly from its peak value. Several days later, base flow begins to increase again at an approximately linear trend, coinciding with the time at which the percentage of quick flow has reached a maximum after each recharge-induced discharge peak. The interval between the discharge peak and the peak in percentage quick flow ranges from 8 to 11 days for seven hydrograph peaks, consistent with quick-flow traveltime estimates by dye-trace tests from the mature karst Hurricane Creek Basin in the central part of the recharge area. Concentrations of environmental tracers chlorofluorocarbons (CFCs: CFC-11, CFC-12, CFC-113), and sulfur hexafluoride in discharge from Big Spring vary approximately linearly with percent quick flow from about 5 to 45 percent of discharge. Linear extrapolation to 100 percent quick flow implies CFC and SF 6 concentrations nearly identical to those in the 2002 atmosphere and indicates a modern age for the quick-flow component. Tracer concentrations for less than about 5 percent quick flow are increasingly lower than those expected from linear extrapolation to zero percent quick flow, indicating that the reservoir of older water in the Big Spring watershed may be a series of water mixtures with piston-flow ages greater than those obtained by extrapolation to zero percent quick flow. Each sample point with a low percentage of quick flow (less than 5 percent) may be a unique mixture. Environmental tracer data from Big Spring plot intermediate to the simple binary mixing of modern and old, pre-tracer water and results from the exponential mixture model. The mean ages of waters in the base-flow component approximately range from 30 to 200 years. The mean age of the base-flow component is youngest (30 to 40 years) in samples containing the highest quick-flow component (45 percent quick flow) and increases to 200 years or more as the fraction of quick flow decreases to less than 5 percent. Tritium data are consistent with a model of dilution of a modern component with an old, pre-tracer component and indicates that the old fraction is mostly pre-1960s in age with mean residence time of more than several hundred years. All of the samples from Big Spring and Greer Spring have water temperatures warmer than their nitrogen-argon recharge temperature, which range from approximately 10.5 to 14° C, suggesting recharge to the Big Spring watershed occurs primarily in late winter to early spring. The water temperatures at Big Spring are consistent with relatively shallow circulation (less than about 600 feet), and the water does not appear to be warmed by deep circulation along a geothermal gradient. Specific conductance values and concentrations of most inorganic constituents in water samples from Big Spring generally decrease with increasing discharge because of dilution with quick-flow water of lower ionic strength. Concentrations of some constituents such as chloride and nitrite plus nitrate, and fecal coliform densities, however, did not decrease with increasing discharge, indicating that quick flow probably is a more important source of these constituents compared to base flow. Water samples from Big Spring plot along the line of dolomite dissolution by carbonic acid, are at equilibrium with dolomite and calcite, and have a molar ratio of Ca:Mg of near 1, indicating dissolution of the mineral dolomite as the primary control on concentrations of calcium, magnesium, and bicarbonate. The flux of calcium and magnesium from Big Spring represents the dissolution of about 1,950 cubic feet of dolomite per day. The suspended sediment load of Big Spring was estimated to range from about 1 to about 70 tons per day, and the sediment load during base-flow periods ranged from about 1 to about 7 tons per day.

Missouri↗