Geology topics

Groundwater & aquifers

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

Follow this topic with RSS

At least 685 records · Page 38Linked to original sources

USGS2019

Simulation of groundwater flow and chloride transport in the “1,500-foot” sand, “2,400-foot” sand, and “2,800-foot” sand of the Baton Rouge area, Louisiana

Groundwater withdrawals since the 1940s have lowered water levels, altered groundwater-flow directions, and caused saltwater to intrude within some freshwater-containing sands of the fluvial-deltaic Southern Hills regional aquifer system beneath Baton Rouge, Louisiana. New interpretations of stratigraphic correlations amongst geophysical well logs were utilized to revise a hydrogeologic framework that delineates the depth and thickness variations of aquifers and confining units in the Southern Hills regional aquifer system. A groundwater-flow and chloride-transport model incorporating the revised framework was constructed to assess the effects of groundwater withdrawals on the rate and pathways of saltwater migration in the “1,500-foot” sand, “2,400-foot” sand, and the “2,800-foot” sand. Groundwater withdrawals reported since 1940 were compiled to specify annual average withdrawal rates through 2016 for 722 wells. Regional groundwater flow throughout the Southern Hills regional aquifer system was first simulated with MODFLOW, and flow-model parameters were calibrated to 8,810 water levels observed through 2016 by using the parameter-estimation code PEST++. Saltwater transport was subsequently simulated for the “1,500-foot” sand, “2,400-foot” sand, and the “2,800-foot” sand by using the variable-density code, SEAWAT. Chloride-concentration measurements were used as a proxy for saltwater to formulate the concentration initial conditions and calibrate the transport-model parameters. Three groundwater-management scenarios were simulated to evaluate the effects of different groundwater withdrawals on future groundwater levels and saltwater concentrations in the “1,500-foot” sand, “2,400-foot” sand, and “2,800-foot” sand. All three scenarios simulated the period from 2017 through 2112 (96 years), and the water levels and concentrations simulated for 2047 and 2112 were compared among the scenarios. The first scenario simulated a continuation of groundwater withdrawals at 2016 rates and represents the “status quo” of groundwater withdrawals. The second scenario simulated the effects of discontinuing 10,620 gallons per minute (gal/min) of withdrawals from the “2,800-foot” sand, and the third scenario simulated reallocating 2,000 gal/min of withdrawals from the “1,500-foot” sand to the “2,800-foot” sand. Continuation of the “status quo” withdrawals results in lower water levels by 2047 around groundwater-withdrawal locations in the “1,500-foot” sand, “2,400-foot” sand, and “2,800-foot” sand. By 2112, water levels recover to higher levels as flow in the aquifer approaches equilibrium. Saltwater within the “1,500-foot” sand would continue migrating toward public-supply wells located 2.4 miles (mi) north of the Baton Rouge Fault, but a “scavenger well” that removes relatively concentrated water from the base of the “1,500-foot” sand attenuates chloride concentrations at the public-supply wells. Saltwater within the “2,400-foot” sand would continue to encroach on a well with large withdrawals and farther east within an area about 1 mi north of the Baton Rouge Fault. Saltwater within the “2,800-foot” sand would migrate northward toward withdrawal wells located about 3 mi north of the industrial district. Cessation of 10,620 gal/min of industrial withdrawals from the “2,800-foot” sand about 12 mi northwest of the industrial district (scenario 2) would cause a substantial water-level recovery in the “2,800-foot” sand in the area of discontinued withdrawals. Groundwater levels 3 mi north of the industrial district would be 25–30 feet higher in 2047 than predicted for the “status quo” withdrawals. Saltwater encroachment toward wells north of the industrial district would be slowed because of the decreased hydraulic gradient. Reallocating 2,000 gal/min of withdrawals from the “1,500-foot” sand to the “2,800-foot” sand 12 mi northwest of the industrial district (scenario 3) would have a negligible effect on water levels and chloride concentrations in the “1,500-foot” sand 15 mi to the south-southeast where saltwater is encroaching toward wells in the “1,500-foot” sand. Within the “2,800-foot” sand, the area of saltwater encroachment is only 3 mi from increased withdrawals in the “2,800-foot” sand, and water levels would be about 5 feet lower in 2047 than for the “status quo” scenario. A larger hydraulic gradient would cause slightly faster saltwater transport and higher chloride concentrations within this area of the “2,800-foot” sand.

Louisiana↗

USGS2015

Geochemical conditions and the occurrence of selected trace elements in groundwater basins used for public drinking-water supply, Desert and Basin and Range hydrogeologic provinces, 2006-11: California GAMA Priority Basin Project

The geochemical conditions, occurrence of selected trace elements, and processes controlling the occurrence of selected trace elements in groundwater were investigated in groundwater basins of the Desert and Basin and Range (DBR) hydrogeologic provinces in southeastern California as part of the Priority Basin Project (PBP) of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The GAMA PBP is designed to provide an assessment of the quality of untreated (raw) groundwater in the aquifer systems that are used for public drinking-water supply. The GAMA PBP is being conducted by the California State Water Resources Control Board in collaboration with the U.S. Geological Survey and the Lawrence Livermore National Laboratory. The DBR hydrogeologic provinces consist of 141 defined groundwater basins separated by mountain ranges, faults, and other features. This report presents analyses of data collected from nine study areas within the DBR hydrogeologic provinces: Antelope Valley, Borrego Valley, the Central Desert area, Coachella Valley, Colorado River, Indian Wells Valley, Low-Use Basins of the Mojave and Sonoran Deserts, the Mojave, and Owens Valley. Collectively, these nine study areas are referred to as the DBR study unit. The study unit covers approximately 7,000 square miles and includes the 63 groundwater basins in the DBR hydrogeologic provinces in which groundwater is used for public drinking-water supply. The vast majority of the 223 wells sampled for this study were long-screened production wells used primarily for public supply. Uncorrected carbon-14 ( 14 C) groundwater ages for samples collected in the DBR study unit ranged from less than (<) 100 to 33,700 years before present (BP). Sixty-six percent of sample ages were greater than (>) 100 years BP, and 40 percent were >3,800 years BP. Samples collected from wells located adjacent to mountain-front recharge areas or major surface-water features generally had younger groundwater ages than did samples collected from wells located away from mountain fronts or towards the distal ends of basin groundwater flow paths. Most groundwater sampled in the DBR study unit had alkaline pH: 89 percent of sample pH values ranged from 7.1 to 9.8, with 37 percent greater than or equal to (≥) 7.9. Groundwater age was significantly correlated (positively) with pH, likely because silicate weathering is a primary control on groundwater pH and is a slow process. The oxidation-reduction (redox) condition of the groundwater sampled in the DBR study unit was predominantly oxic (71 percent), except in the Colorado River study area where organic-rich fluvial aquifers provide the electron donors necessary to support iron-reducing (anoxic-Fe) redox processes. The cation type of 78 percent of the samples was either sodium- or mixed-type, and the anion type of 83 percent of the samples was either bicarbonate- or mixed-type. Sodium-type groundwaters generally were older and more alkaline than calcium-type groundwaters, consistent with the change in water chemistry expected from cation exchange between groundwater and aquifer sediments over long periods of time. Because of the correlation with young groundwater, calcium-type groundwater was predominantly from wells located adjacent to mountain-front recharge areas. Arsenic (As), boron (B), fluoride (F), molybdenum (Mo), strontium (Sr), uranium (U), and vanadium (V) were selected for assessment in this study because they occurred at concentrations greater than California Department of Public Health or U.S. Environmental Protection Agency regulatory or non-regulatory drinking-water-quality benchmarks in more than 2 percent of the 223 samples collected in the DBR study unit. As and F were detected most commonly (18 and 13 percent, respectively) at concentrations above associated water-quality benchmarks and Sr and V least frequently (both at 3 percent). Given that 14 C groundwater ages are predominantly >100 years BP, land use in the study unit is primarily undeveloped, and chemicals derived from anthropogenic sources, such as volatile organic compounds, were infrequently detected, high concentrations of these trace elements in groundwater were most likely the result of natural factors and not anthropogenic factors. As, F, Mo, and V concentrations showed significant positive correlations to groundwater age and to pH. This relation is partly due to the sources of trace elements likely being the weathering of primary minerals, such as silicate minerals, which is a slow process that takes place over hundreds to thousands of years. This relation also reflects the positive correlation between groundwater age and pH. Geochemical modeling predicted that the dominant species of As, Mo, and V in solution were oxyanions (HAsO 4 2– , MoO 4 2– , and H 2 VO 4– ), which are likely to be mobile in alkaline groundwater because mineral surfaces composing aquifer matrices have a predominantly negative surface charge under alkaline conditions. F also exists predominantly as a negatively charged ion (F – ). At pH values >7.5, saturation indices generated by the geochemical modeling program PHREEQC indicated that F solubility may be somewhat limited by the precipitation of the mineral fluorapatite [Ca 5 (PO 4 ) 3 F]. Speciation modeling of As in anoxic-Fe groundwater (iron-reducing conditions) showed that samples were supersaturated with orpiment (As 2 S 3 ), indicating that mineral precipitation may be responsible for low As concentrations observed in reducing groundwater. In contrast, U concentrations showed significant negative correlations to groundwater age and to pH. Higher U concentrations generally occurred in samples for which geochemical modeling indicated that the uncharged ternary complex Ca 2 UO 2 (CO 3 ) 3 was the dominant aqueous U species. This uncharged complex is not attracted to the charged surfaces of minerals and thus increases U solubility. Formation of Ca 2 UO 2 (CO 3 ) 3 was greater in younger groundwaters because calcium and uranium concentrations generally were lower in older groundwaters, likely due to cation-exchange processes and precipitation of the mineral calcite as groundwater pH increased. Co-precipitation of U with the calcite (CaCO 3 ) may remove U from the aqueous phase. Saturation indices indicated that the anoxic-Fe groundwaters from the Colorado River study area were supersaturated with the mineral uraninite (UO 2 ), suggesting that UO 2 precipitation may be responsible for the low concentrations of U observed in these samples. Concentrations of strontium, which exists primarily in a cationic form (Sr 2+ ), were not significantly correlated with either groundwater age or pH. Strontium concentrations showed a strong positive correlation with total dissolved solids (TDS). Dissolved constituents, such as Sr, that interact with mineral surfaces through outer-sphere complexation become increasingly soluble with increasing TDS concentrations of groundwater. Boron concentrations also showed a significant positive correlation with TDS, indicating the B may interact to a large degree with mineral surfaces through outer-sphere complexation.

California↗

USGS2023

Determining three-dimensional hydrologic processes in the eastern Snake River Plain aquifer using geochemical mass-balance modeling, Idaho National Laboratory, eastern Idaho, with contributions by Treinen, K.C.

Waste constituents discharged to the eastern Snake River Plain aquifer at the U.S. Department of Energy (DOE) Idaho National Laboratory (INL) pose risks to the water quality of the aquifer. To understand these risks, the U.S. Geological Survey, in cooperation with the DOE, used geochemical mass-balance modeling to identify three-dimensional hydrologic processes in that portion of the aquifer underlying the southwestern part of the INL that affect the movement of groundwater and waste constituents. Modeling was performed using water chemistry of 74 water samples collected from 30 wells. Fifty-four of the water samples were collected from 11 wells equipped with multilevel monitoring systems with vertically discrete sampling zones that encompass the upper 750 feet of the aquifer. Water samples from these multilevel wells were collected during 2007‒13, a period when conditions in the aquifer were approximately steady-state because there was little or no recharge from the Big Lost River. The primary source of water in groundwater at the multilevel wells during 2007‒13 was the Big Lost River. Other sources of water include groundwater from the Little Lost River valley, precipitation, and wastewater. Horizontal groundwater-flow directions appear to be similar in both the shallow and deep parts of the aquifer, and surface-water sources of water in most deep groundwater shows that groundwater moves downward. Surface-water sources of water in deep groundwater noticeably decrease within and below the Matuyama flow and associated sedimentary interbeds, which indicates that these units are semi-impermeable and retard the downward movement of groundwater.

Idaho↗

USGS2019

Potentiometric surface of the Mississippi River Valley alluvial aquifer, spring 2016

A potentiometric surface map for spring 2016 was created for the Mississippi River Valley alluvial (MRVA) aquifer using selected available groundwater-altitude data from wells and surface-water-altitude data from streamgages. Most of the wells were measured annually or one time after installation, but some wells were measured more than one time or continually; streamgages are typically operated continuously. Personnel from the Arkansas Natural Resources Commission, Arkansas Department of Health, Arkansas Geological Survey, Illinois Department of Agriculture, Illinois State Water Survey, Louisiana Department of Natural Resources, Louisiana Department of Transportation and Development, Mississippi Department of Environmental Quality, Yazoo Mississippi Delta Joint Water Management District, U.S. Department of Agriculture–Natural Resources Conservation Service, and the U.S. Geological Survey (USGS) routinely collect groundwater data from wells screened in the MRVA aquifer. The USGS and the U.S. Army Corps of Engineers routinely collect data on river stage and discharge for the rivers overlying the MRVA aquifer. The potentiometric surface map for 2016 was created using existing data as part of the USGS Water Availability and Use Science Program to support investigations that characterize the MRVA aquifer. Sufficient groundwater-altitude data were available to create a potentiometric-surface map for spring 2016 for about 81 percent of the aquifer area. The potentiometric contours ranged from 10 to 340 feet. The regional direction of groundwater flow in the MRVA aquifer was generally towards the south-southwest, except in areas of groundwater-altitude depressions, where groundwater flows into the depressions, and near rivers, where groundwater flow generally parallels the flow in the rivers. There are large depressions in the potentiometric surface of the MRVA aquifer in the lower half of the Cache region and in most of the Grand Prairie and Delta regions.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗

USGS2009

Simulations of Groundwater Flow and Particle Tracking Analysis in the Area Contributing Recharge to a Public-Supply Well near Tampa, Florida, 2002-05

Shallow ground water in the north-central Tampa Bay region, Florida, is affected by elevated nitrate concentrations, the presence of volatile organic compounds, and pesticides as a result of groundwater development and intensive urban land use. The region relies primarily on groundwater for drinking-water supplies. Sustainability of groundwater quality for public supply requires monitoring and understanding of the mechanisms controlling the vulnerability of public-supply wells to contamination. A single public-supply well was selected for intensive study based on the need to evaluate the dominant processes affecting the vulnerability of public-supply wells in the Upper Floridan aquifer in the City of Temple Terrace near Tampa, Florida, and the presence of a variety of chemical constituents in water from the well. A network of 29 monitoring wells was installed, and water and sediment samples were collected within the area contributing recharge to the selected public-supply well to support a detailed analysis of physical and chemical conditions and processes affecting the water chemistry in the well. A three-dimensional, steady-state groundwater flow model was developed to evaluate the age of groundwater reaching the well and to test hypotheses on the vulnerability of the well to nonpoint source input of nitrate. Particle tracking data were used to calculate environmental tracer concentrations of tritium and sulfur hexafluoride and to calibrate traveltimes and compute flow paths and advective travel times in the model area. The traveltime of particles reaching the selected public-supply well ranged from less than 1 day to 127.0 years, with a median of 13.1 years; nearly 45 percent of the simulated particle ages were less than about 10 years. Nitrate concentrations, derived primarily from residential/commercial fertilizer use and atmospheric deposition, were highest (2.4 and 6.11 milligrams per liter as nitrogen, median and maximum, respectively) in shallow groundwater from the surficial aquifer system and lowest (less than the detection level of 0.06 milligram per liter) in the deeper Upper Floridan aquifer. Denitrification occurred near the interface of the surficial aquifer system and the underlying intermediate confining unit, within the intermediate confining unit, and within the Upper Floridan aquifer because of reducing conditions in this part of the flow system. However, simulations indicate that the rapid movement of water from the surficial aquifer system to the selected public-supply well through karst features (sinkholes) and conduit layers that bypass the denitrifying zones (short-circuits), coupled with high pumping rates, allow nitrate to reach the selected public-supply well in concentrations that resemble those of the overlying surficial aquifer system. Water from the surficial aquifer system with elevated concentrations of nitrate and low concentrations of some volatile organic compounds and pesticides is expected to continue moving into the selected public-supply well, because calculated flux-weighted concentrations indicate the proportion of young affected water contributing to the well is likely to remain relatively stable over time. The calculated nitrate concentration in the selected public-supply well indicates a lag of 1 to 10 years between peak concentrations of nonpoint source contaminants in recharge and appearance in the well.

Scientific Investigations Report↗

USGS2020

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

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

South Dakota↗

USGS2007

Bedrock aquifers and population growth in the Denver Basin, Colorado, USA

The Denver Basin bedrock aquifer system consists of Tertiary and Cretaceous age sedimentary rocks known as the Dawson, Denver, Arapahoe and Laramie-Fox Hills aquifers. The number of bedrock wells has increased from 12,000 in 1985 to over 33,700 in 2001 and the withdrawal of groundwater has caused water level declines in excess of 75 meters. Water level declines now range from 3 to 12 meters per year in the critical Arapahoe Aquifer. The groundwater supplies were once thought to be sufficient for 100 years but there is concern that they may be depleted in 10 to 15 years in areas on the west side of the basin. Groundwater is being mined from the aquifer system because the withdrawal through wells exceeds the rate of recharge. Increased groundwater withdrawal will cause further water level declines, increased costs to pump groundwater, and reduced yield from existing wells. In the Denver Basin, hydrologists have some capability to monitor declines in water levels for the Arapaho Aquifer, but generally have a limited ability to monitor water use. More complete and accurate water use data are needed to predict groundwater longevity for the Arapahoe Aquifer. The life of the Arapahoe Aquifer can be extended with artificial recharge using imported surface water, water reuse, restrictions on lawn watering, well permit restrictions and other conservation measures. Availability of groundwater may limit growth in the Denver Basin over the next 20 years unless residents are willing to pay for additional new sources of supply.

Episodes↗

USGS2016

Geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers within northern Bexar and Comal Counties, Texas

During 2014–16, the U.S. Geological Survey, in cooperation with the Edwards Aquifer Authority, documented the geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers within northern Bexar and Comal Counties, Texas. The Edwards and Trinity aquifers are major sources of water for agriculture, industry, and urban and rural communities in south-central Texas. Both the Edwards and Trinity are classified as major aquifers by the State of Texas. The purpose of this report is to present the geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers within northern Bexar and Comal Counties, Tex. The report includes a detailed 1:24,000-scale hydrostratigraphic map, names, and descriptions of the geology and hydrostratigraphic units (HSUs) in the study area. The scope of the report is focused on geologic framework and hydrostratigraphy of the outcrops and hydrostratigraphy of the Edwards and Trinity aquifers within northern Bexar and Comal Counties, Tex. In addition, parts of the adjacent upper confining unit to the Edwards aquifer are included. The study area, approximately 866 square miles, is within the outcrops of the Edwards and Trinity aquifers and overlying confining units (Washita, Eagle Ford, Austin, and Taylor Groups) in northern Bexar and Comal Counties, Tex. The rocks within the study area are sedimentary and range in age from Early to Late Cretaceous. The Miocene-age Balcones fault zone is the primary structural feature within the study area. The fault zone is an extensional system of faults that generally trends southwest to northeast in south-central Texas. The faults have normal throw, are en echelon, and are mostly downthrown to the southeast. The Early Cretaceous Edwards Group rocks were deposited in an open marine to supratidal flats environment during two marine transgressions. The Edwards Group is composed of the Kainer and Person Formations. Following tectonic uplift, subaerial exposure, and erosion near the end of Early Cretaceous time, the area of present-day south-central Texas was again submerged during the Late Cretaceous by a marine transgression resulting in deposition of the Georgetown Formation of the Washita Group. The Early Cretaceous Edwards Group, which overlies the Trinity Group, is composed of mudstone to boundstone, dolomitic limestone, argillaceous limestone, evaporite, shale, and chert. The Kainer Formation is subdivided into (bottom to top) the basal nodular, dolomitic, Kirschberg Evaporite, and grainstone members. The Person Formation is subdivided into (bottom to top) the regional dense, leached and collapsed (undivided), and cyclic and marine (undivided) members. Hydrostratigraphically the rocks exposed in the study area represent a section of the upper confining unit to the Edwards aquifer, the Edwards aquifer, the upper zone of the Trinity aquifer, and the middle zone of the Trinity aquifer. The Pecan Gap Formation (Taylor Group), Austin Group, Eagle Ford Group, Buda Limestone, and Del Rio Clay are generally considered to be the upper confining unit to the Edwards aquifer. The Edwards aquifer was subdivided into HSUs I to VIII. The Georgetown Formation of the Washita Group contains HSU I. The Person Formation of the Edwards Group contains HSUs II (cyclic and marine members [Kpcm], undivided), III (leached and collapsed members [Kplc,] undivided), and IV (regional dense member [Kprd]), and the Kainer Formation of the Edwards Group contains HSUs V (grainstone member [Kkg]), VI (Kirschberg Evaporite Member [Kkke]), VII (dolomitic member [Kkd]), and VIII (basal nodular member [Kkbn]). The Trinity aquifer is separated into upper, middle, and lower aquifer units (hereinafter referred to as “zones”). The upper zone of the Trinity aquifer is in the upper member of the Glen Rose Limestone. The middle zone of the Trinity aquifer is formed in the lower member of the Glen Rose Limestone, Hensell Sand, and Cow Creek Limestone. The regionally extensive Hammett Shale forms a confining unit between the middle and lower zones of the Trinity aquifer. The lower zone of the Trinity aquifer consists of the Sligo and Hosston Formations, which do not crop out in the study area. The upper zone of the Trinity aquifer is subdivided into five informal HSUs (top to bottom): cavernous, Camp Bullis, upper evaporite, fossiliferous, and lower evaporite. The middle zone of the Trinity aquifer is composed of the (top to bottom) Bulverde, Little Blanco, Twin Sisters, Doeppenschmidt, Rust, Honey Creek, Hensell, and Cow Creek HSUs. The underlying Hammett HSU is a regional confining unit between the middle and lower zones of the Trinity aquifer. The lower zone of the Trinity aquifer is not exposed in the study area. Groundwater recharge and flow paths in the study area are influenced not only by the hydrostratigraphic characteristics of the individual HSUs but also by faults and fractures and geologic structure. Faulting associated with the Balcones fault zone (1) might affect groundwater flow paths by forming a barrier to flow that results in water moving parallel to the fault plane, (2) might affect groundwater flow paths by increasing flow across the fault because of fracturing and juxtaposing porous and permeable units, or (3) might have no effect on the groundwater flow paths. The hydrologic connection between the Edwards and Trinity aquifers and the various HSUs is complex. The complexity of the aquifer system is a combination of the original depositional history, bioturbation, primary and secondary porosity, diagenesis, and fracturing of the area from faulting. All of these factors have resulted in development of modified porosity, permeability, and transmissivity within and between the aquifers. Faulting produced highly fractured areas that have allowed for rapid infiltration of water and subsequently formed solutionally enhanced fractures, bedding planes, channels, and caves that are highly permeable and transmissive. The juxtaposition resulting from faulting has resulted in areas of interconnectedness between the Edwards and Trinity aquifers and the various HSUs that form the aquifers.

Texas↗

USGS2021

Water-level conditions in the confined aquifers of the New Jersey Coastal Plain, 2013

The Coastal Plain aquifers of New Jersey provide an important source of water for more than 3.5 million people. In 2013, groundwater withdrawals from 10 confined aquifers of the New Jersey Coastal Plain totaled about 190 million gallons per day. Steadily increasing withdrawals from the late 1800s to the early 1990s resulted in declining water levels and the formation of regional cones of depression in many confined Coastal Plain aquifers. Starting in 1978, the U.S. Geological Survey (USGS) began mapping the potentiometric surfaces of the major confined Coastal Plain aquifers every 5 years to provide a regional assessment of groundwater conditions. In a study conducted by the USGS, in cooperation with the New Jersey Department of Environmental Protection, water levels in 10 confined aquifers of the New Jersey Coastal Plain were measured and evaluated to provide a regional overview of groundwater conditions during fall 2013. Water levels were measured in 987 wells in New Jersey, and parts of Pennsylvania and Delaware. Potentiometric-surface maps were prepared for, in ascending order of age, the confined Cohansey aquifer of Cape May County, Rio Grande water-bearing zone, Atlantic City 800-foot sand, Piney Point aquifer, Vincentown aquifer, Wenonah-Mount Laurel aquifer, Englishtown aquifer system, and the Upper, Middle, and Lower aquifers of the Potomac-Raritan-Magothy (PRM) aquifer system. Persistent, regionally extensive cones of depression were present in the potentiometric surfaces of the Englishtown aquifer system and Wenonah-Mount Laurel aquifer in Ocean and Monmouth Counties; Wenonah-Mount Laurel and Upper, Middle, and Lower PRM aquifers in Camden County; and Atlantic City 800-foot sand in Atlantic County. Changes in water levels from 2008 to 2013 were measured in many Coastal Plain aquifers in New Jersey. In some areas, water levels continued to decline as a result of pumping, but in other areas water levels continued to recover as a result of regulated decreases in groundwater withdrawals. Since 2008, in the confined Cohansey aquifer in Cape May County, water levels generally did not change; however, cones of depression in the potentiometric surface of the Piney Point aquifer in some areas of Cumberland County deepened by more than 20 feet (ft). In Critical Area 1, an area of restricted withdrawals, measured water levels in the Wenonah-Mount Laurel aquifer declined in parts of southern Monmouth County by more than 10 ft; however, rises in water levels of more than 10 ft were measured in parts of northern Ocean and Monmouth Counties. Since 2008, in Critical Area 2, also an area of restricted withdrawals, measured water levels in the Wenonah-Mount Laurel aquifer rose more than 20 ft in parts of western Burlington County and more than 20 ft in parts of western Camden County. Since 2008, in Critical Area 1, measured water levels in the Englishtown aquifer system declined in parts of eastern Ocean County by more than 10 ft and in southeastern Monmouth County by more than 20 ft; however, rises in water levels of more than 10 ft were measured in other parts of Ocean and Monmouth Counties. In general, since 2008 in Critical Area 2, in the Upper PRM aquifer, measured water levels continued to rise by 10 ft or more in central and western Burlington and central Camden Counties. In the Middle PRM aquifer in Critical Area 2, measured water levels rose in parts of central Camden County by 10 ft or more. However, measured water levels in the Lower PRM aquifer in Critical Area 2 were more than 10 ft lower in the center of the cone of depression in central Camden County, but measured water levels continued to rise updip from this area in Critical Area 2. Seasonal water-level fluctuations are presented in time-series hydrographs for 77 wells during 1978–2013. Analyses of long-term water-level changes for the period 2008–13 indicate downward water-level trends at 14 wells (18 percent), upward trends at 34 wells (44 percent), and no substantial change at 29 wells (38 percent). Downward trends were most often observed for wells screened in the Piney Point aquifer and the Atlantic City 800-foot sand. Upward water-level trends were most often measured for wells screened in the PRM aquifer system. Upward water-level trends also were measured for wells in the Englishtown aquifer system and the Wenonah-Mount Laurel aquifer in Critical Area 1 in some areas; however, downward trends and no substantial changes were measured in other areas.

New Jersey↗

USGS2012

Tabulated Transmissivity and Storage Properties of the Floridan Aquifer System in Florida and Parts of Georgia, South Carolina, and Alabama

A goal of the U.S. Geological Survey Groundwater Resources Program is to assess the availability of fresh water within each of the principal aquifers in the United States with the greatest groundwater withdrawals. The Floridan aquifer system (FAS), which covers an area of approximately 100,000 square miles in Florida and parts of Georgia, Alabama, Mississippi, and South Carolina, is one such principal aquifer, having the fifth largest groundwater withdrawals in the Nation, totaling 3.64 billion gallons per day in 2000. Compilation of FAS hydraulic properties is critical to the development and calibration of groundwater flow models that can be used to develop water budgets spatially and temporally, as well as to evaluate resource changes over time. Wells with aquifer test data were identified as Upper Floridan aquifer (UFA), Lower Floridan aquifer (LFA), Floridan aquifer system (FAS, Upper Floridan with some middle and/or Lower Floridan), or middle Floridan confining unit (MCU), based on the identification from the original database or report description, or comparison of the open interval of the well with previously published maps. This report consolidates aquifer hydraulic property data obtained from multiple databases and reports of the U.S. Geological Survey, various State agencies, and the Water Management Districts of Florida, that are compiled into tables to provide a single information source for transmissivity and storage properties of the FAS as of October 2011. Transmissivity calculated from aquifer pumping tests and specific-capacity data are included. Values for transmissivity and storage coefficients are intended for use in regional or sub regional groundwater flow models; thus, any tests (aquifer pumping tests and specific capacity data) that were conducted with packers or for open intervals less than 30 feet in length are excluded from the summary statistics and tables of this report, but are included in the database. The transmissivity distribution from the aquifer pumping tests is highly variable. The transmissivity based on aquifer pumping tests (from 1,045 values for the UFA and FAS) ranges from 8 to about 9,300,000 square feet per day (ft 2 /d) and values of storage coefficient (646 reported) range from 3x10 -9 to 0.41. The 64 transmissivity values for the LFA range from about 130 to 4,500,000 ft 2 /d, and the 17 storage coefficient values range from 7x10 -8 to 0.03. The 14 transmissivity values for the MCU range from 1 to about 600,000 ft 2 /d and the 10 storage coefficient values range from 8x10 -8 to 0.03. Transmissivity estimates for the UFA and FAS for 442 specific capacity tests range from approximately 200 to 1,000,000 ft 2 /d.

Alabama, Florida, Georgia, South Carolina↗

USGS2014

Survival of bacterial indicators and the functional diversity of native microbial communities in the Floridan aquifer system, south Florida

The Upper Floridan aquifer in the southern region of Florida is a multi-use, regional scale aquifer that is used as a potable water source and as a repository for passively recharged untreated surface waters, and injected treated surface water and wastewater, industrial wastes, including those which contain greenhouse gases (for example, carbon dioxide). The presence of confined zones within the Floridan aquifer that range in salinity from fresh to brackish allow regulatory agencies to permit the injection of these different types of product waters into specific zones without detrimental effects to humans and terrestrial and aquatic ecosystems. The type of recharge that has received the most regulatory attention in south Florida is aquifer storage and recovery (ASR). The treated water, prior to injection and during recovery, must meet primary and secondary drinking water standards. The primary microbiology drinking water standard is total coliforms, which have been shown to be difficult to inactivate below the regulatory standard during the treatment process at some ASR facilities. The inefficient inactivation of this group of indicator bacteria permits their direct injection into the storage zones of the Floridan aquifer. Prior to this study, the inactivation rates for any member of the total coliform group during exposure to native geochemical conditions in groundwater from any zone of the Floridan aquifer had not been derived. Aboveground flow through mesocosms and diffusion chambers were used to quantify the inactivation rates of two bacterial indicators, Escherichia coli and Pseudomonas aeruginosa, during exposure to groundwater from six wells. These wells collect water from two ASR storage zones: the Upper Floridan aquifer (UFA) and Avon Park Permeable Zone (APPZ). Both bacterial strains followed a biphasic inactivation model. The E. coli populations had slower inactivation rates in the UFA (range: 0.217–0.628 per hour (h -1 )) during the first phase of the model than when exposed to groundwater from the APPZ (range: 0.540–0.684 h -1 ). The inactivation rates for the first phase of the models for P. aeruginosa were not significantly different between the UFA (range: 0.144–0.770 h -1 ) and APPZ (range: 0.159–0.772 h -1 ) aquifer zones. The inactivation rates for the second phase of the model for this P. aeruginosa were also similar between UFA (range: 0.003–0.008 h -1 ) and APPZ (0.004–0.005 h -1 ) zones, although significantly slower than the model’s first phase rates for this bacterial species. Geochemical data were used to determine which dissimilatory biogeochemical reactions were most likely to occur under the native conditions in the UFA and APPZ zones using thermodynamics principles to calculate free energy yields and other cell-related energetics data. The biogeochemical processes of acetotrophic and hydrogenotrophic sulfate reduction, methanogenesis and anaerobic oxidation of methane dominated in all six groundwater sites. A high throughput DNA microarray sequencing technology was used to characterize the diversity in the native aquifer bacterial communities (bacteria and archaea) and assign putative physiological capabilities to the members of those communities. The bacterial communities in both zones of the aquifer were shown to possess the capabilities for primary and secondary fermentation, acetogenesis, methanogenesis, anaerobic methane oxidation, syntrophy with methanogens, ammonification, and sulfate reduction. The data from this study provide the first determination of bacterial indicator survival during exposure to native geochemical conditions of the Floridan aquifer in south Florida. Additionally, the energetics and functional bacterial diversity characterizations are the first descriptions of native bacterial communities in this region of the Floridan aquifer and reveal how these communities persist under such extreme conditions. Collectively, these types of data can be used to develop and refine groundwater models.

Florida↗

USGS2011

Simulations of groundwater flow and particle-tracking analysis in the zone of contribution to a public-supply well in San Antonio, Texas

In 2006, a public-supply well in San Antonio, Texas, was selected for intensive study to assess the vulnerability of public-supply wells in the Edwards aquifer to contamination by a variety of compounds. A local-scale, steady-state, three-dimensional numerical groundwater-flow model was developed and used in this study to evaluate the movement of water and solutes from recharge areas to the selected public-supply well. Particle tracking was used to compute flow paths and advective traveltimes throughout the model area and to delineate the areas contributing recharge and zone of contribution for the selected public-supply well. The local-scale model grid has a finer vertical discretization than do previous regional Edwards aquifer models and incorporates refined parameter zones corresponding with multiple (10) hydrogeologic units representing the Edwards aquifer. In the Edwards aquifer, high matrix porosity and permeability likely are overshadowed by high permeability developed in structurally influenced karstic conduit systems that transmit water into, through, and out of the aquifer system. The complexity of the aquifer system in the local-scale study area is further increased by numerous faults with varying vertical displacements. The extensive faulting results in the juxtaposition of hydrogeologic units with differing hydraulic properties and has appreciable effects on groundwater flow in the Edwards aquifer. The local-scale model simulations use the MODFLOW Hydrogeologic-Unit Flow Package and include two hydrogeologic units with high hydraulic conductivities (one or more orders of magnitude higher than for the other simulated hydrogeologic units) that are intended to simulate fast flow paths attributable to karst features. The two “conduit” hydrogeologic units of the Edwards aquifer represent the lower 8 meters of the leached and collapsed members and the Kirschberg evaporite member of the Edwards Group. The MODFLOW Horizontal-Flow Barrier Package was used to simulate faults in the local-scale model. The assumption was made that the degree to which a fault acts as a barrier to groundwater flow is proportional to the fault displacement. The final calibrated hydraulic-conductance values ranged from 0.01 to 0.2 per day for fault displacements ranging from 0 to more than 100 percent of the total aquifer thickness. The calibrated steady-state simulation generally reproduces the spatial distribution of measured water-level altitudes. Simulated water-level altitudes were within 9.0 meters of measured water-level altitudes at 74 of the 84 wells used as targets for the local-scale model for the calibrated steady-state simulation. The overall mean absolute difference between simulated and measured water-level altitudes is 4.2 meters, and the mean algebraic difference is 1.9 meters. The simulated springflow for San Antonio Springs was 7.7 percent greater and for San Pedro Springs was 4.2 percent less than the median measured springflow. Simulated tritium concentrations were within 0.14 tritium units of measured tritium concentrations for 11 of the 13 local-scale study tritium observations from the 10 local-scale study wells used to calibrate the steady-state local-scale model, with a mean absolute difference between simulated and measured tritium concentrations of 0.11 tritium units and a mean algebraic difference of -0.04 tritium units. Simulated tritium concentrations in the selected public-supply well during November 2007 were within 0.09 tritium units of the measured concentrations, with the exception of the shallowest observation from the well. The steady-state simulation water budget indicates that recharge occurring in the local-scale study area accounts for 31.8 percent of the sources of water to the Edwards aquifer in the local-scale model area and that inflow through the model boundaries contributes 68.2 percent. Most of the flow into the local-scale model area through the model boundaries occurs through the western and southern boundaries, 58.2 and 39.6 percent, respectively. The largest discharges from the Edwards aquifer in the local-scale model area are boundary outflow (71.4 percent) and withdrawals by wells (24.9 percent). Most of the flow out of the local-scale model area through the model boundaries occurs through the southern and eastern boundaries, 54.2 and 39.6 percent, respectively. The simulated zones of contribution for the selected public-supply well, Timberhill well nest, and Zarzamora well nest extend to the north, northeast, and northwest from each site in the confined zone of the aquifer into the recharge zone, where all recharge to the aquifer occurs. The area contributing recharge for the selected public-supply well has the greatest extent. The area contributing recharge for the Timberhill well nest encompasses approximately the western one-half of the area contributing recharge for the selected public-supply well, and that for the Zarzamora well nest encompasses approximately the eastern two-thirds of the area contributing recharge for the selected public-supply well. Simulated particle ages ranged from less than 1 day to more than 1,900 years in the 10 local-scale study wells (13 local-scale study tritium observations) used to calibrate the local-scale model. The simulated mean particle ages for the tritium observations representing selected well depths (shallow, intermediate, and deep) ranged from 2.5 to 15 years. The minimum (youngest) mean particle ages for the selected public-supply well and the Timberhill monitoring wells were for the intermediate well depth, while the youngest mean particle age for the Zarzamora monitoring wells was for the intermediate and deep well depth. The maximum (oldest) mean particle ages for the selected public-supply well and the Zarzamora monitoring wells were for the shallow well depth. The mean of simulated particle ages for tritium observations representing well depths open to the simulated conduit hydrogeologic units was 3.8 years, whereas the mean of simulated particle ages for tritium observations representing well depths not open to the simulated conduit hydrogeologic units was 9.6 years. The effect of short-circuit pathways, for example karst conduits, in the flow system on the movement of young water to the selected public-supply well could greatly alter contaminant arrival times compared to what might be expected from advection in a system without short circuiting. In a forecasting exercise, the simulated concentrations showed rapid initial response at the beginning and end of chemical input, followed by more gradual response as older water moved through the system. The nature of karst groundwater flow, where flow predominantly occurs via conduit flow paths, could lead to relatively rapid water quality responses to land-use changes. Results from the forecasting exercise indicate that timescales for change in the quality of water from the selected public-supply well could be on the order of a few years to decades for land-use changes that occur over days to decades, which has implications for source-water protection strategies that rely on land-use change to achieve water-quality objectives.

Texas↗

USGS2013

Hydrologic data and groundwater flow simulations in the vicinity of Long Lake, Indiana Dunes National Lakeshore, near Gary, Indiana

The U.S. Geological Survey (USGS) collected data and simulated groundwater flow to increase understanding of the hydrology and the effects of drainage alterations to the water table in the vicinity of Long Lake, near Gary, Indiana. East Long Lake and West Long Lake (collectively known as Long Lake) make up one of the largest interdunal lakes within the Indiana Dunes National Lakeshore. The National Park Service is tasked with preservation and restoration of wetlands in the Indiana Dunes National Lakeshore along the southern shoreline of Lake Michigan. Urban development and engineering have modified drainage and caused changes in the distribution of open water, streams and ditches, and groundwater abundance and flow paths. A better understanding of the effects these modifications have on the hydrologic system in the area will help the National Park Service, the Gary Sanitary District (GSD), and local stakeholders manage and protect the resources within the study area. This study used hydrologic data and steady-state groundwater simulations to estimate directions of groundwater flow and the effects of various engineering controls and climatic conditions on the hydrology near Long Lake. Periods of relatively high and low groundwater levels were examined and simulated by using MODFLOW and companion software. Simulated hydrologic modifications examined the effects of (1) removing the beaver dams in US-12 ditch, (2) discontinuing seepage of water from the filtration pond east of East Long Lake, (3) discontinuing discharge from US-12 ditch to the GSD sewer system, (4) decreasing discharge from US-12 ditch to the GSD sewer system, (5) connecting East Long Lake and West Long Lake, (6) deepening County Line Road ditch, and (7) raising and lowering the water level of Lake Michigan. Results from collected hydrologic data indicate that East Long Lake functioned as an area of groundwater recharge during October 2002 and a “flow-through” lake during March 2011, with the groundwater divide south of US-12. Wetlands to the south of West Long Lake act as points of recharge to the surficial aquifer in both dry- and wet-weather conditions. Among the noteworthy results from a dry-weather groundwater flow model simulation are (1) US-12 ditch does not receive water from East Long Lake or West Long Lake, (2) the filtration pond at the east end of East Long Lake, when active, contributed approximately 10 percent of the total water entering East Long Lake, and (3) County Line Road ditch has little effect on simulated water level. Among the noteworthy results from a wet-weather groundwater flow simulation are (1) US-12 ditch does not receive water from East Long Lake or West Long Lake, (2) when the seepage from the filtration pond to the surficial aquifer is not active, sources of inflow to East Long Lake are restricted to only precipitation (46 percent of total) and inflow from the surficial aquifer (54 percent of total), and (3) County Line Road ditch bisects the groundwater divide and creates two water-table mounds south of US-12. The results from a series of model scenarios simulating certain engineering controls and changes in Lake Michigan levels include the following: (1) The simulated removal of beaver dams in US-12 ditch during a wet-weather simulation increased discharge from the ditch to the Gary Sanitary system by 13 percent. (2) Discontinuation of seepage from the filtration pond east of East Long Lake decreased discharge from US-12 ditch to the Gary Sanitary system by 2.3 percent. (3) Simulated discontinuation of discharge from the US-12 ditch to the GSD sewer system increased the area where the water table was estimated to be above the land surface beyond the inundated area in the initial wet-weather simulation. (4) Simulated modifications to the control structure at the discharge point of US-12 ditch to the GSD sewer system can decrease discharge by as much as 61 percent while increasing the simulated inundated area during dry weather and decrease discharge as much as 6 percent while increasing the simulated inundated area during wet weather. (5) Deepening of County Line Road ditch can decrease the discharge from US-12 ditch by 26 percent during dry weather and 24 percent during wet weather, as well as decrease the extent of flooded areas south and east of the filtration pond near Ogden Dunes. (7) The increase of the Lake Michigan water level to match the historical maximum can increase the discharge from US-12 ditch by 14 percent during dry weather and by 9.6 percent during wet weather. (8) The decrease of the Lake Michigan water level to match the historical minimum can decrease the discharge from US-12 ditch by 7.4 percent during dry weather and by 3.1 percent during wet weather. The results of this study can be used by water-resource managers to understand how surrounding ditches affect water levels in East and West Long Lake and in the surrounding wetlands and residential areas. The groundwater model developed in this study can be applied in the future to answer questions about how alterations to the drainage system in the area will affect water levels in East and West Long Lake and surrounding areas. The modeling methods developed in this study provide a template for other studies of groundwater flow and groundwater/surface-water interactions within the shallow surficial aquifer in northern Indiana, and in similar hydrologic settings that include surficial sand aquifers in coastal settings.

Indiana↗

USGS2023

Prediction of the probability of elevated nitrate concentrations at groundwater depths used for drinking-water supply in the Puget Sound basin, Washington, 2004–19

The Puget Sound basin encompasses the 13,700-square-mile area that drains to the Puget Sound and the adjacent marine waters of Washington State. Well more than 4 million people live within the basin, with numbers continuing to increase, who rely on the basin’s natural resources including groundwater. The Puget Sound Partnership was created by a Washington State statute to implement a science-based recovery of the Puget Sound to help address impacts to these resources. As part of the recovery, the partnership developed the Puget Sound Vital Signs as measures of ecosystem health that guide the assessment of progress toward Puget Sound recovery goals. The Puget Sound Partnership Leadership Council adopted a Drinking Water Vital Sign associated with human health and quality of life, recognizing certain indicators as integral to the sustainability of Puget Sound recovery efforts. One such Vital Sign indicator was the vulnerability of groundwater throughout the aquifers of the Puget Sound basin to elevated nitrate concentrations as defined by the probability of exceeding 2 milligrams/liter (mg/L) at a specific location and well depth. The U.S. Geological Survey (USGS) led the effort to characterize groundwater vulnerability. For this study, groundwater vulnerability refers to a probability with which a contaminant applied at or near the land surface can migrate to the aquifer of interest for a given set of land-use practices. Nitrate concentration data were selected for evaluation because elevated nitrate concentrations are typically caused by anthropogenic activities and have been associated with deleterious impacts on human health. To identify groundwater vulnerability to elevated nitrate concentrations, logistic regression was used to relate anthropogenic (human associated) and natural variables to the occurrence of elevated nitrate concentrations in untreated groundwater from large public water supply system wells found within the Washington State Department of Health Sentry database. Variables that were analyzed included well depth, soil hydraulic conductivity, precipitation, population density, fertilizer application amounts, and land-use types. Statistically significant models that predicted the probabilities of groundwater nitrate concentrations greater than 2 mg/L based on the predictor variables were created for the time periods 2000–04, 2005–09, 2010–14, and 2015–19. For all time periods, well depth and a measure of the abundance of urban and agricultural land over or near the well consistently helped explain the vulnerability of the well to elevated nitrate concentrations defined as a probability of exceeding 2 mg/L of nitrate. Precipitation and (or) soil hydraulic conductivity were also important predictor variables in the models. The models for each time period were used to create maps of groundwater vulnerability at 150- and 300-foot depths throughout the Puget Sound basin. As expected, the most vulnerable locations were associated with shallower well depths and increased agriculture and urban land cover. Across all four time periods, groundwater vulnerability throughout the Puget Sound was low, with probabilities of exceeding 2 mg/L concentrations of nitrate at depths at 150 and 300 feet typically less than 50 percent. Results also found a slight decrease in probabilities of elevated nitrate concentrations throughout the basin over time. More specifically, additional statistical tests found that groundwater with probabilities of less than about 60 percent declined from 2000 to 2019 and represented more than 75 percent of the modeled Puget Sound basin aquifer. Wells with greater than 60 percent probability increased over the same time period but represented only about 25 percent of the aquifer. The maps and statistical analysis presented in the study provide valuable and informative evaluation of the vulnerability of groundwater in the Puget Sound basin to elevated nitrate concentrations. The probability maps do not represent measured nitrate concentrations in groundwater, but rather they present the probability that nitrate concentrations exceed 2 mg/L. The models and predictions from this study are a viable indicator for the Puget Sound Partnership’s Healthy Human Population—Drinking Water Vital Sign. The logistic regression modeling approach presented here benefits water managers by allowing them to assess temporal trends in a range of probabilities, explore vulnerability changes as new regional land cover and anthropogenic data are generated, and distinguish vulnerabilities at different depths within the aquifer.

Washington↗

USGS2018

Drain tiles and groundwater resources: Understanding the relations

Executive Summary Drainage for agricultural production over the past 150 years has been an integral component of human-driven change to Minnesota’s rural landscapes. Benefits of drainage Historically, poorly drained soils across much of the State would often remain saturated or flooded after spring snowmelt, preventing timely farm operations such as tilling and planting crops (Arneman, 1963). Installation of agricultural drainage, both surface ditches and subsurface drainage, accelerated transport of water off farm fields and imparted producers higher crop yields (Beauchamp, 1987; Stoner and others, 1993). Agricultural drainage offered many other benefits such as preventing crop drown out, aerating the soil profile for improved plant growth, limiting surface runoff and soil erosion, and allowing farmers better access to croplands (Fausey and others, 1987). Without agricultural drainage on much of Minnesota’s croplands, it would have been difficult to realize high enough crop yields to remain economically viable. Environmental concerns While drainage of Minnesota croplands provided the benefits mentioned above, several environmental concerns result. These include wetland loss, degradation of downstream water quality, and reduced [potential for] recharge. Early agricultural drainage efforts (pre-20th century) led to the disappearance of much of Minnesota’s natural wetlands. Increased focus on preventing or mitigating wetland loss over the last 50 years has helped curtail further losses, even as agricultural drainage proceeds. Prior to establishment of Minnesota statehood, wetlands accounted for more than 10 million acres in Minnesota, including prairie wetlands, peatlands, and forest wetlands that comprised approximately 19 percent of the total land area (Palmer, 1915; King, 1980). In 2018, only half of Minnesota’s pre-settlement wetlands remain, mostly in parts of the State that have not experienced widespread drainage, such as northern Minnesota. Water-quality monitoring has shown that agricultural drainage, in particular the practice of subsurface drainage, provides a direct flow path for nutrient (nitrogen and soluble phosphorus) losses to surface water resources. The negative consequences of agricultural drainage on surface water quality are well documented (for example, Dinnes and others, 2002; Kladivko and others, 2004; Richards and others, 2008; Rozemeijer and others, 2010; Schottler and others, 2013). Agricultural basins with a high percentage of agricultural drainage have been implicated as part of the cause of the Gulf of Mexico hypoxia zone due to excessive nitrogen export (Goolsby and Battaglin, 2001; Randall and Mulla, 2001). The connection of hydrological effects of agricultural subsurface drainage on groundwater recharge and aquifers, on the other hand, has not been well-established. Agricultural subsurface drainage intercepts infiltrating water below croplands and directly discharges the water to nearby surface waters. However, the size of the water balance shift from drained water that would have evapotranspired or run off the land to drained water that would recharge underlying aquifers has been poorly characterized (Schuh, 2008). Drain Tiles and Groundwater Given the poor accounting of subsurface drainage effects on groundwater resources, the Minnesota Ground Water Association (MGWA) deemed it imperative that we document these effects so that groundwater resources in agricultural regions with substantial drainage can be effectively managed. This white paper documents the relations of drain tiles and groundwater resources and discusses the historical significance of agricultural drainage practices, the recognized positive benefits and potential negative consequences of agricultural drainage practices, and the gaps in understanding of the connections between agricultural drainage and groundwater resources. The major messages emerged from the findings of this white paper are: Complex history. Minnesota has a long history of agricultural drainage, spanning over 150 years. Agricultural drainage, and the eventual widespread usage of subsurface drainage, can be separated into at least four distinct periods of time: (a) early drainage to get water off the land, pre-20th century; (b) the boom and bust era (1900-1945); (c) postwar resurgence of subsurface drainage and early conservation efforts (1945-1960); and, (d) emergence of the environmental movement (1960 to present). The State’s regulatory framework that both allowed for drainage and controlled its usage during these periods is complicated and has been governed by a patchwork of both State and Federal statutes. Drainage Provinces. This white paper advances the concept of tile drainage provinces to aid in the discussion of regional differences in subsurface drainage and its overall effect on groundwater resources. Built upon the concept of groundwater provinces (Minnesota Department of Natural Resources, 2001), three distinct tile drainage provinces were conceptualized: (1) the Southeastern Province; (2) the SouthCentral Province; and, (3) the Western Province. The distinct geology and the soils that developed in these regions have implications for each region’s subsurface drainage density and the potential implications for groundwater. Knowledge gaps. Several critical knowledge gaps are identified in this paper, creating opportunities for further research to improve our understanding for better managed water resources: Extent of drainage is unknown. Direct estimates of the extent of subsurface drainage do not exist in Minnesota. However, several indirect methods have been utilized to estimate subsurface drainage, from the field-scale to countylevel through the use of geographic information system (GIS) analysis and aerial photography. Based on a 2012 U.S. Geological Survey estimate of subsurface drainage extent (Nakagaki and Wieczorek, 2016), about 21% of the land area in Minnesota has some density of subsurface drainage. Effect of drainage on underlying aquifers is unknown. A basic understanding of unconfined and confined aquifers and their recharge is necessary to connect any hydrological effects from agricultural drainage to groundwater. The basic goal of subsurface drainage to efficiently drain saturated soils clearly alters the water balance in croplands. However, its overall effect on groundwater resources has been poorly characterized, and is in large part determined by the geology below drained areas and the arrangement of underlying aquifers. Water balance shifts. An improved understanding of historical water balance shifts from pre- to post-drainage periods is necessary to understand long-term implications on net groundwater recharge. Also, more direct field-scale studies and indirect modeling studies are needed to characterize water budgets for fields with subsurface drainage.

Minnesota↗

USGS2011

Hydrogeologic setting and simulation of groundwater flow near the Canterbury and Leadville Mine Drainage Tunnels, Leadville, Colorado

The Leadville mining district is historically one of the most heavily mined regions in the world producing large quantities of gold, silver, lead, zinc, copper, and manganese since the 1860s. A multidisciplinary investigation was conducted by the U.S. Geological Survey, in cooperation with the Colorado Department of Public Health and Environment, to characterize large-scale groundwater flow in a 13 square-kilometer region encompassing the Canterbury Tunnel and the Leadville Mine Drainage Tunnel near Leadville, Colorado. The primary objective of the investigation was to evaluate whether a substantial hydraulic connection is present between the Canterbury Tunnel and Leadville Mine Drainage Tunnel for current (2008) hydrologic conditions. Altitude in the Leadville area ranges from about 3,018 m (9,900 ft) along the Arkansas River valley to about 4,270 m (14,000 ft) along the Continental Divide east of Leadville, and the high altitude of the area results in a moderate subpolar climate. Winter precipitation as snow was about three times greater than summer precipitation as rain, and in general, both winter and summer precipitation were greatest at higher altitudes. Winter and summer precipitation have increased since 2002 coinciding with the observed water-level rise near the Leadville Mine Drainage Tunnel that began in 2003. The weather patterns and hydrology exhibit strong seasonality with an annual cycle of cold winters with large snowfall, followed by spring snowmelt, runoff, and recharge (high-flow) conditions, and then base-flow (low-flow) conditions in the fall prior to the next winter. Groundwater occurs in the Paleozoic and Precambrian fractured-rock aquifers and in a Quaternary alluvial aquifer along the East Fork Arkansas River, and groundwater levels also exhibit seasonal, although delayed, patterns in response to the annual hydrologic cycle. A three-dimensional digital representation of the extensively faulted bedrock was developed and a geophysical direct-current resistivity field survey was performed to evaluate the geologic structure of the study area. The results show that the Canterbury Tunnel is located in a downthrown structural block that is not in direct physical connection with the Leadville Mine Drainage Tunnel. The presence of this structural discontinuity implies there is no direct groundwater pathway between the tunnels along a laterally continuous bedrock unit. Water-quality results for pH and major-ion concentrations near the Canterbury Tunnel showed that acid mine drainage has not affected groundwater quality. Stable-isotope ratios of hydrogen and oxygen in water indicate that snowmelt is the primary source of groundwater recharge. On the basis of chlorofluorocarbon and tritium concentrations and mixing ratios for groundwater samples, young groundwater (groundwater recharged after 1953) was indicated at well locations upgradient from and in a fault block separate from the Canterbury Tunnel. Samples from sites downgradient from the Canterbury Tunnel were mixtures of young and old (pre-1953) groundwater and likely represent snowmelt recharge mixed with older regional groundwater that discharges from the bedrock units to the Arkansas River valley. Discharge from the Canterbury Tunnel contained the greatest percentage of old (pre-1953) groundwater with a mixture of about 25 percent young water and about 75 percent old water. A calibrated three-dimensional groundwater model representing high-flow conditions was used to evaluate large-scale flow characteristics of the groundwater and to assess whether a substantial hydraulic connection was present between the Canterbury Tunnel and Leadville Mine Drainage Tunnel. As simulated, the faults restrict local flow in many areas, but the fracture-damage zones adjacent to the faults allow groundwater to move along faults. Water-budget results indicate that groundwater flow across the lateral edges of the model controlled the majority of flow in and out of the aquifer (79 percent and 63 percent of the total water budget, respectively). The largest contributions to the water budget were groundwater entering from the upper reaches of the watershed and the hydrologic interaction of the groundwater with the East Fork Arkansas River. Potentiometric surface maps of the simulated model results were generated for depths of 50, 100, and 250 m. The surfaces revealed a positive trend in hydraulic head with land-surface altitude and evidence of increased control on fluid movement by the fault network structure at progressively greater depths in the aquifer. Results of advective particle-tracking simulations indicate that the sets of simulated flow paths for the Canterbury Tunnel and the Leadville Mine Drainage Tunnel were mutually exclusive of one another, which also suggested that no major hydraulic connection was present between the tunnels. Particle-tracking simulations also revealed that although the fault network generally restricted groundwater movement locally, hydrologic conditions were such that groundwater did cross the fault network at many locations. This cross-fault movement indicates that the fault network controls regional groundwater flow to some degree but is not a complete barrier to flow. The cumulative distributions of adjusted age results for the watershed indicate that approximately 30 percent of the flow pathways transmit groundwater that was younger than 68 years old (post-1941) and that about 70 percent of the flow pathways transmit old groundwater. The particle-tracking results are consistent with the apparent ages and mixing ratios developed from the chlorofluorocarbon and tritium results. The model simulations also indicate that approximately 50 percent of the groundwater flowing through the study area was less than 200 years old and about 50 percent of the groundwater flowing through the study area is old water stored in low-permeability geologic units and fault blocks. As a final examination of model response, the conductance parameters of the Canterbury Tunnel and Leadville Mine Drainage Tunnel were manually adjusted from the calibrated values to determine if altering the flow discharge in one tunnel affects the hydraulic behavior in the other tunnel. The examination showed no substantial hydraulic connection. The multidisciplinary investigation yielded an improved understanding of groundwater characteristics near the Canterbury Tunnel and the Leadville Mine Drainage Tunnel. Movement of groundwater between the Canterbury Tunnel and Leadville Mine Drainage Tunnel that was central to this investigation could not be evaluated with strong certainty owing to the structural complexity of the region, study simplifications, and the absence of observation data within the upper sections of the Canterbury Tunnel and between the Canterbury Tunnel and the Leadville Mine Drainage Tunnel. There was, however, collaborative agreement between all of the analyses performed during this investigation that a substantial hydraulic connection did not exist between the Canterbury Tunnel and the Leadville Mine Drainage Tunnel under natural flow conditions near the time of this investigation.

Colorado↗

USGS2017

Hydrogeologic framework and selected components of the groundwater budget for the upper Umatilla River Basin, Oregon

Executive Summary This report presents a summary of the hydrogeology of the upper Umatilla River Basin, Oregon, based on characterization of the hydrogeologic framework, horizontal and vertical directions of groundwater flow, trends in groundwater levels, and components of the groundwater budget. The conceptual model of the groundwater flow system integrates available data and information on the groundwater resources of the upper Umatilla River Basin and provides insights regarding key hydrologic processes, such as the interaction between the groundwater and surface water systems and the hydrologic budget. The conceptual groundwater model developed for the study area divides the groundwater flow system into five hydrogeologic units: a sedimentary unit, three Columbia River basalt units, and a basement rock unit. The sedimentary unit, which is not widely used as a source of groundwater in the upper basin, is present primarily in the lowlands and consists of conglomerate, loess, silt and sand deposits, and recent alluvium. The Columbia River Basalt Group is a series of Miocene flood basalts that are present throughout the study area. The basalt is uplifted in the southeastern half of the study area, and either underlies the sedimentary unit, or is exposed at the surface. The interflow zones of the flood basalts are the primary aquifers in the study area. Beneath the flood basalts are basement rocks composed of Paleogene to Pre-Tertiary sedimentary, volcanic, igneous, and metamorphic rocks that are not used as a source of groundwater in the upper Umatilla River Basin. The major components of the groundwater budget in the upper Umatilla River Basin are (1) groundwater recharge, (2) groundwater discharge to surface water and wells, (3) subsurface flow into and out of the basin, and (4) changes in groundwater storage. Recharge from precipitation occurs primarily in the upland areas of the Blue Mountains. Mean annual recharge from infiltration of precipitation for the upper Umatilla River Basin during 1951–2010 is about 9.6 inches per year (in/yr). Annual recharge from precipitation for water year 2010 ranged from 3 in. in the lowland area to about 30 in. in the Blue Mountains. Using Kahle and others (2011) data and methods from the Columbia Plateau regional model, average annual recharge from irrigation is estimated to be about 2.2 in/yr for the 13 square miles of irrigated land in the upper Umatilla River Basin. Groundwater discharges to streams throughout the year and is a large component of annual streamflow in the upper Umatilla River Basin. Upward vertical hydraulic gradients near the Umatilla River indicate the potential for groundwater discharge. Groundwater discharge to the Umatilla River generally occurs in the upper part of the basin, upstream from the main stem. Groundwater development in the upper Umatilla River Basin began sometime after 1950 (Davies-Smith and others, 1988; Gonthier and Bolke, 1991). By water year 2010, groundwater use in the upper Umatilla River Basin was approximately 11,214 acre-feet (acre-ft). Total groundwater withdrawals for the study area were estimated at 7,575 acre-ft for irrigation, 3,173 acre-ft for municipal use, and 466 acre-ft for domestic use. Total groundwater flow into or from the study area depends locally on geology and hydraulic head distribution. Estimates of subsurface flow were calculated using the U.S. Geological Survey Columbia Plateau regional groundwater flow model. Net flux values range from 25,000 to 27,700 acre-ft per year and indicate that groundwater is moving out of the upper Umatilla River Basin into the lower Umatilla River Basin. Water level changes depend on storage changes within an aquifer, and storage changes depend on the storage properties of the aquifer, as well as recharge to or discharge from the aquifer. Groundwater level data in the upper Umatilla River Basin are mostly available from wells in Columbia River basalt units, which indicate areas of long-term water level declines in the Grande Ronde basalt unit near Pendleton and Athena, Oregon. Groundwater levels in the Wanapum basalt unit do not show long-term declines in the upper Umatilla River Basin. Because of pumping, some areas in the upper Umatilla River Basin have shown a decrease, or reversal, in the upward vertical head gradient. Key data needs are improvement of the spatial and temporal distribution of water-level data collection and continued monitoring of streamflow gaging sites. Additionally, refinement of recharge estimates would enhance understanding of the processes that provide the groundwater resources in the upper Umatilla River Basin.

Oregon↗

USGS2009

Groundwater, surface–water, and water-chemistry data, Black Mesa area, northeastern Arizona—2007-2008

The N aquifer is an extensive aquifer and the primary source of groundwater in the 5,400-square-mile Black Mesa area in northeastern Arizona. Availability of water is an important issue in northeastern Arizona because of continued water requirements for industrial and municipal use by a growing population and because of low precipitation in the arid climate of the Black Mesa area, which is typically about 6 to 14 inches per year. The U.S. Geological Survey water-monitoring program in the Black Mesa area began in 1971 and provides information about the long-term effects of groundwater withdrawals from the N aquifer for industrial and municipal uses. This report presents results of data collected as part of the monitoring program in the Black Mesa area from January 2007 to September 2008. The monitoring program includes measurements of (1) groundwater withdrawals, (2) groundwater levels, (3) spring discharge, (4) surface-water discharge, and (5) groundwater chemistry. In 2007, total groundwater withdrawals were 4,270 acre-feet, industrial withdrawals were 1,170 acre-ft, and municipal withdrawals were 3,100 acre-ft. Total withdrawals during 2007 were about 41 percent less than total withdrawals in 2005. From 2006 to 2007, however, total withdrawals increased by 4 percent, industrial withdrawals decreased by approximately 2 percent, and total municipal withdrawals increased by 7 percent. From 2007 to 2008, annually measured water levels in the Black Mesa area declined in 6 of 11 wells measured in the unconfined areas of the N aquifer, and the median change was -0.2 feet. Water levels declined in 9 of 18 wells measured in the confined area of the aquifer. The median change for the confined area of the aquifer was -0.2 feet. From the prestress period (prior to 1965) to 2008, the median water-level change for 33 wells in both the confined and unconfined area was -12.9 feet. Median water-level changes were -1.0 feet for 15 wells measured in the unconfined areas and -33.2 feet for 18 wells measured in the confined area. Spring flow was measured at two springs in 2008. Flow decreased at both Moenkopi School Spring and Pasture Canyon Spring from previous years. Flow fluctuated during the period of record, but a decreasing trend was apparent. Continuous records of surface-water discharge in the Black Mesa area were collected from streamflow-gaging stations at the following sites: Moenkopi Wash at Moenkopi 09401260 (1976 to 2007), Dinnebito Wash near Sand Springs 09401110 (1993 to 2007), Polacca Wash near Second Mesa 09400568 (1994 to 2007), and Pasture Canyon Springs 09401265 (August 2004 to 2007). Median winter flows (November through February) of each water year were used as an index of the amount of groundwater discharge at the above-named sites. For the period of record of each streamflow-gaging station, the median winter flows have generally remained constant, which suggests no change in groundwater. The period of record is too short to determine if there is a trend at Pasture Canyon Spring. In 2008, water samples collected from 6 wells and 2 springs in the Black Mesa area were analyzed for selected chemical constituents and the results compared with previous analyses. Concentrations of dissolved solids, chloride, and sulfate have varied at all 6 wells for the period of record, but neither increasing nor decreasing trends over time were found. Dissolved-solids, chloride, and sulfate concentrations increased at Moenkopi School Spring during the more than 12 years of record at that site. Concentrations of dissolved solids, chloride, and sulfate at Pasture Canyon Spring have not varied much since the early 1980s, and there is no trend in those data.

Arizona↗