USGS Science⌕ Search

SEARCH · USGS Science

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.

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

1,001 records · Page 56Linked to original sources
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-flow model. This scenario showed decreases of nearly 50 percent in base-flow concentrations of nitrate in streams within the first 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-flow 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-flow 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 flow 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 denitrification. Spatial and temporal variations in nitrogen sources affect all the watersheds. Tributaries with higher inputs of nitrogen have higher concentrations in stream base flow. 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 denitrification, which occurs in Pocomoke and East Mahantango Creeks but is not evident in Polecat and Muddy Creeks. In East Mahantango Creek, the level of denitrification is significant 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 significant quantity of nitrogen from the system. Denitrification in Pocomoke Creek is significant 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, denitrification may not greatly affect the overall nitrogen delivery to these streams. Other findings of this study show that nitrate in ground water discharging along preferential flow paths may not be affected by natural processes, such as denitrification 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-flow paths and 3) areas not affected by significant denitrification. The fastest response in stream base-flow 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 stratification. Overall findings of the study indicate that 1) ground-water contributions to nitrogen in streamflow are significant, 2) some response to management practices should be evident in base-flow 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 streamflow 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↗
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↗