SEARCH · USGS Science

Results for “groundwater aquifer”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

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

At least 883 records · Page 49Linked to original sources

USGS2019

Controls on spatial and temporal variations of brine discharge to the Dolores River in the Paradox Valley, Colorado, 2016–18

The Paradox Valley in southwestern Colorado is a collapsed anticline formed by movement of the salt-rich Paradox Formation at the core of the anticline. The salinity of the Dolores River, a tributary of the Colorado River, increases substantially as it crosses the valley because of discharge of brine-rich groundwater derived from the underlying salts. Although the brine is naturally occurring, it increases the salinity of the Colorado River, which is a major concern to downstream agricultural, municipal, and industrial water users. The U.S. Geological Survey in cooperation with the Bureau of Reclamation conducted a study to improve the characterization of processes controlling spatial and temporal variations in brine discharge to the Dolores River. For the study, three geophysical surveys were conducted in March, May, and September 2017, and water levels were monitored in selected ponds and groundwater wells from November 2016 to May 2018. The study also utilized streamflow and specific conductance data from two U.S. Geological Survey streamflow-gaging stations on the Dolores River to estimate salt load to the river. River-based continuous resistivity profiling and frequency domain electromagnetic induction surveys made during low-flow conditions indicated a zone of brine-rich groundwater close to the riverbed along an approximately 4-kilometer reach of the river. Under high-flow conditions, the brine was depressed as much as 2 meters below the riverbed, and brine discharge to the river was reduced to a minimum. Direct current electrical resistivity surveys show that the freshwater lens overlying the brine is much thicker (up to 10 meters) on the west bank than on the east bank (less than 5 meters). A large low-conductivity anomaly at river distance 6,800 meters was observed in all surveys and may represent a freshwater discharge zone or a losing reach of the river. Filling and draining of the wildlife ponds on the west side of the river had a negligible effect on salt loads in the river during the study period. Groundwater monitoring showed there was active exchange of water between the river and the adjacent alluvial aquifer. When river stage was low, groundwater flowed towards the river, and brine discharge to the river increased. When the river stage was high, the gradient was reversed, and fresh surface water recharged the alluvial aquifer minimizing brine discharge. Most of the salt load to the river occurred during the winter and appeared to be enhanced by diurnal stage fluctuations. A conceptual model of brine discharge to the river is presented at three scales. Groundwater at the regional scale drives dissolution of salt in the Paradox Formation and flow of brine into the base of the alluvial aquifer. Surface water–groundwater interactions at the scale of the alluvial aquifer control brine discharge to the river seasonally and interannually. At the finest scale, diurnal fluctuations in river stage drive exchange of freshwater with saltier pore water in the hyporheic zone, which appears to increase brine discharge to the river during winter.

Colorado↗

USGS2013

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

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

Oklahoma↗

USGS2016

Hydrological conditions and evaluation of sustainable groundwater use in the Sierra Vista Subwatershed, Upper San Pedro Basin, southeastern Arizona

This study assessed progress toward achieving sustainable groundwater use in the Sierra Vista Subwatershed of the Upper San Pedro Basin, Arizona, through evaluation of 14 indicators of sustainable use. Sustainable use of groundwater in the Sierra Vista Subwatershed requires, at a minimum, a stable rate of groundwater discharge to, and thus base flow in, the San Pedro River. Many of the 14 indicators are therefore related to long-term or short-term effects on base flow and provide us with a means to evaluate groundwater discharge to and base flow in the San Pedro River. The indicators were based primarily on 10 to 20 years of data monitoring in the subwatershed, ending in 2012, and included subwatershedwide indicators, riparian-system indicators, San Pedro River indicators, and springs indicators. Groundwater management actions including voluntary retirement of irrigation pumping in the subwatershed resulted in about a 5,100 acre-feet (acre-ft) reduction in net human use from 2002 to 2012. Subwatershed population increased more than 10,000 during the same period. Most of the reduction occurred during 2002–07 and included reductions in groundwater pumping and increases in managed recharge; net human use varied annually by a few hundred acre-ft during 2007–12. The groundwater budget for 2012 showed a deficit of about 5,000 acre-ft, although the total water-budget uncertainty was about 5,500 acre-ft. In the vicinity of the U.S. Army’s Fort Huachuca, regional-aquifer water levels were in steady decline beginning in at least the mid-1990s (in older wells since at least the early-1970s), as the cone of depression centered on the Sierra Vista and Fort Huachuca pumping centers continued to deepen. This was evident in the individual water levels on Fort Huachuca, as well as from the horizontal hydraulic gradients that extend from the pumping centers toward the San Pedro and Babocomari Rivers. Basin water levels in wells southeast of Sierra Vista, away from the river, were also experiencing declines, while some water levels closer to the river were rising. Near-stream vertical gradients along the San Pedro River showed no clear increasing or decreasing trends that would indicate a shift in the direction of subsurface flow between the riverbed and the alluvial aquifer, or a trend in the magnitude of groundwater/surface-water exchange. Annual streamflow permanence data showed no clear change in streamflow permanence trends in any of the river reaches, other than those related to precipitation trends. Similarly, the single-day, dry-season, wet-dry streamflow analysis of all subwatershed river reaches indicated no change in condition over the past 14 years, with the exception of the Hereford reach, which has seen a statistically significant increase in wetted length. Dry-season, alluvial-aquifer water levels in the Hereford reach also showed a statistically significant increase. These improvements are attributed to the end of irrigation pumping in the area. Although data indicate that the length of the Fairbank North wetted reach may be in decline, it is not yet statistically significant. Stable-isotope data indicated reduced groundwater discharge to the Babocomari River in the vicinity of the Babocomari River near Tombstone gaging station and to the San Pedro River near the San Pedro River at Palominas gaging station and near the Lewis Springs DCP stage recorder. The Babocomari River near Tombstone gaging station is downgradient of the major pumping centers. The change in isotopic signature at the Lewis Springs stage recorder could have been the result of alterations in groundwater/surface-water interactions there caused by beaver damming of the river. Base flow in the San Pedro River declined over the periods of record at the three San Pedro River gaging stations in the subwatershed (Palominas, Charleston, and Tombstone), as well as at the Babocomari River near Tombstone gaging station. Precipitation declined slightly from the 1990s to the 2000s, although there is no statistically significant trend in subwatershed precipitation from 1991 to 2012. The occurrence of large winter discharge events appeared to decline and that of large summer discharge events appeared to increase over this same period. Data for physical parameters, general chemistry, nutrient species, select trace elements, and suspended sediment were collected at San Pedro River at Charleston stream-gaging station. These data were summarized over time and analyzed in relation to discharge and season as a means to assess trends over the period of analysis. Federal and State of Arizona drinking-water and human-contact standards were all met and few exceedances occurred for the ecological thresholds investigated. Several constituents showed a significant trend over the period of analysis, but only concentration and flux data for total phosphate, orthophosphate, total nitrogen, suspended sediment, and sulfate were suitable to be used in a weighted regression analysis that statistically accounted for time, discharge, and season. Sulfate concentrations and flux showed a significant downward trend over the period of analysis, whereas total phosphorus and ortho-phosphate showed a relatively small magnitude upward trend relative to standards. Suspended sediment concentrations and flux both showed a significant downward trend in the 1980s, an effect attributed to reduction of cattle in the subwatershed at about this time, and (or) increased cottonwood ( Populus fremontii ) and willow ( Salix goodingii ) recruitment, and (or) the curtailment of sand and gravel mining adjacent to the San Pedro River with the designation of the San Pedro Riparian National Conservation Area in 1988. A spike in sediment flux in 2006 may be attributable to the more than 100 debris flows in the Huachuca Mountains during the summer monsoon of that year. Spring discharge along the San Pedro River generally increased at three sites proximate to the Sierra Vista treated effluent recharge facility and varied somewhat with climate at two other sites. Median annual discharge at the recharge facility peaked in 2006, and at Murray Springs and Horsethief Spring, downgradient of the recharge facility, in 2009. Sampling for trace organic compounds in flow from springs was carried out using both discrete sampling and passive sampling methods. Spring samples thus collected showed the presence of trace-organic compounds. Lewis Springs (background site) had the least number of detections, whereas Murray Springs, located directly downgradient of the City of Sierra Vista’s treated effluent recharge facility, had the greatest number of detections of all the springs. Discrete samples from the recharge facility had more than twice the detections found in discrete samples from Murray Spring and at much higher concentrations. Few similar trace-organic compounds were detected at both the springs and the treated effluent recharge facility, and the number of detections did not increase during the collection period. Limitations of the study prevented the determination of trace-organic concentration in passive samplers and also prevented linking trace organic compounds detected at the treated effluent recharge facility with compounds detected from the springs. In particular, trace organic compounds could also derive from other sources such as septic systems. Looking at the subwatershed as a whole, base flow was in decline along the entire river reach, but determination of the specific cause of the decline was beyond the scope of this report. Conditions in the area from the municipal pumping center of Sierra Vista and Fort Huachuca northeast to the river (from about the Charleston to Tombstone gaging stations) were more commonly in decline than in regions further south. Both long-term indicators, such as regional aquifer groundwater levels and horizontal gradients, and the isotope analysis indicated that groundwater discharge to the river and thus base flow may continue to decline in that area. South of Charleston, indicators were more mixed. Some indicators in the Hereford reach suggest groundwater discharge to the San Pedro River may be increasing there, whereas some indicators in the Palominas reach suggest groundwater discharge to the river there may be declining.

Arizona↗

USGS2020

Use of real-time sensors to temporally characterize water quality in groundwater and surface water in Mason County, Illinois, 2017–19

The persistence of high nitrate concentrations in shallow groundwater has been well documented in the shallow glacial aquifer of Mason County, Illinois. Nitrates in groundwater can be a concern when concentrations exceed 10 milligrams per liter in drinking water. Additionally, nitrate in groundwater can contribute to surface water nitrogen loads that can cause increased algal growth. Algal growth increases oxygen consumption causing anoxic conditions as observed in the Gulf of Mexico Hypoxic Zone. From March 8, 2017, to March 31, 2019, groundwater level, continuous nitrate, dissolved oxygen, specific conductance, water temperature, and pH data were collected in a monitoring well to temporally assess changes in water quality using high frequency data. During this same period, instantaneous field measurements of water quality and groundwater levels were made in surface water and groundwater in and near Quiver Creek in the presumed groundwater flow path about 0.6 mile from the continuous monitoring well. Groundwater nitrate concentrations continuously measured in the aquifer during this time ranged from 14.7 to 23.2 milligrams per liter, whereas instantaneously measured nitrate concentrations in Quiver Creek ranged from 0.9 to 6.4 milligrams per liter. Nitrate concentrations measured by piezometer varied laterally and vertically in the Quiver Creek floodplain and beneath the stream. Irrigation and fertigation for agriculture is widely practiced in Mason County. This may seasonally affect the groundwater flow and movement as well as the persistence of nitrate in this area. Continuously and instantaneously measured nitrate concentrations and groundwater levels indicate that during the irrigation season, discharge to Quiver Creek from the shallow groundwater system may be limited. During and following periods when estimated irrigation use is highest, the low-nitrate deeper groundwater may be the dominant contributor to the Quiver Creek surface water, whereas during recharge events and when the system is not under the stress of irrigation, there is more contribution from the local and higher-nitrate shallow groundwater.

Illinois↗

USGS1990

Sources and distribution of nitrate in ground water at a farmed field irrigated with sewage treatment-plant effluent, Tallahassee, Florida

The city of Tallahassee, Florida began applying sewage treatment-plant effluent to a sprayfield southeast of the city in 1980. Fertilizers containing inorganic nitrogen were also applied in conjunction with the operation of a commercial farm at this site. Analysis of groundwater in the surficial aquifer and the Upper Floridan aquifer have indicated that nitrate concentrations in some wells exceed the prescribed drinking water maximum contaminant level of 10 mg/L (nitrate as nitrogen). Nitrate concentrations greater than the maximum contaminant level were not detected in samples from monitoring wells outside the sprayfield boundary. Analyses of water from the unsaturated zone indicated that conversion of organic nitrogen and ammonia to nitrate was complete before the nitrogen- enriched water reached the water table. Groundwater samples from wells in the surficial and Upper Floridan aquifers less than 100 ft deep located inside sprayed areas had mean concentrations of nitrate much higher than samples from similar wells located outside sprayed areas at the southeast sprayfield. These shallow wells inside the sprayed areas were the only wells in which the maximum contaminant level for nitrate was exceeded. Analyses of the nitrogen isotope ratios in groundwater were used to determine whether the major source of nitrogen was treated sewage or fertilizers. The nitrogen isotope ratios in contaminated groundwater at the southeast sprayfield were compared to those at another sprayfield southwest of the city, where treated sewage was the sole source of nitrogen. Statistical analyses indicated a significant difference in the nitrogen isotope ratios at the two sites, indicating that both nitrogen sources are significant at the southeast sprayfield. (USGS)

Water-Resources Investigations Report↗

USGS1990

Preliminary investigation of soil and ground-water contamination at a U.S. Army Petroleum Training Facility, Fort Lee, Virginia, September-October 1989

Fuel-oil constituents in the soil and groundwater at the Fort Lee Petroleum Training Facility near Petersburg, Virginia, were studied by the U.S. Geological Survey (USGS) in cooperation with the Department of Defense, U.S. Army. The study included installation of 25 groundwater monitoring wells and description of groundwater flow patterns of the shallow-aquifer system underlying the facility. Soil and groundwater samples were collected to determine the concentrations of fuel-oil constituents and to determine the potential for off-site migration of the constituents. Total petroleum hydrocarbon concentrations up to 18,400 mg/km were reported in soil samples. Concentrations of benzene in water from wells at the facility were up to 130 micrograms per liter (ug/L), and concentrations of ethylbenzene and xylene were up to 54 and 120 ug/L, respectively. Potential exists for off-site migration of the contaminants and migration of contaminants downward to deeper aquifers. Further investigations of these potential contamination-migration pathways are warranted. Risk identification at the Petroleum Training Facility cannot be properly addressed because the distribution of the fuel-oil constituents has not been fully characterized. Preliminary identification of risk, however is presented by an examination of toxicity data for the chemical constituents reported in the groundwater at the facility. Concentrations of constituents were compared to the maximum contaminant levels (MCLs) for drinking water established by the U.S. Environmental Protection Agency (USEPA). Concentrations of benzene in water from wells at the facility exceed the USEPA 's 5 ug/L MCL by as much as 26 times. Sufficient data are not available to fully design the remedial-action plan for the facility; however, general responses to contamination of the type associated with the facility include no-action, monitoring, institutional controls, removal, and treatment.

Virginia↗

USGS2013

Appraising options to reduce shallow groundwater tables and enhance flow conditions over regional scales in an irrigated alluvial aquifer system

Some of the world’s key agricultural production systems face big challenges to both water quantity and quality due to shallow groundwater that results from long-term intensive irrigation, namely waterlogging and salinity, water losses, and environmental problems. This paper focuses on water quantity issues, presenting finite-difference groundwater models developed to describe shallow water table levels, non-beneficial groundwater consumptive use, and return flows to streams across two regions within an irrigated alluvial river valley in southeastern Colorado, USA. The models are calibrated and applied to simulate current baseline conditions in the alluvial aquifer system and to examine actions for potentially improving these conditions. The models provide a detailed description of regional-scale subsurface unsaturated and saturated flow processes, thereby enabling detailed spatiotemporal description of groundwater levels, recharge to infiltration ratios, partitioning of ET originating from the unsaturated and saturated zones, and groundwater flows, among other variables. Hybrid automated and manual calibration of the models is achieved using extensive observations of groundwater hydraulic head, groundwater return flow to streams, aquifer stratigraphy, canal seepage, total evapotranspiration, the portion of evapotranspiration supplied by upflux from the shallow water table, and irrigation flows. Baseline results from the two regional-scale models are compared to model predictions under variations of four alternative management schemes: (1) reduced seepage from earthen canals, (2) reduced irrigation applications, (3) rotational lease fallowing (irrigation water leased to municipalities, resulting in temporary dry-up of fields), and (4) combinations of these. The potential for increasing the average water table depth by up to 1.1 and 0.7 m in the two respective modeled regions, thereby reducing the threat of waterlogging and lowering non-beneficial consumptive use from adjacent fallow and naturally-vegetated lands, is demonstrated for the alternative management intervention scenarios considered. Net annual average savings of up to about 9.9 million m 3 (8000 ac ft) and 2.3 million m 3 (1900 ac ft) of non-beneficial groundwater consumptive use is demonstrated for the study periods in each of the two respective study regions. Alternative water management interventions achieve varying degrees of benefits in each of the two regions, suggesting a need to adopt region-specific interventions and avoid a ‘one-size-fits-all’ approach. Impacts of the considered interventions on return flows to the river were predicted to be significant, highlighting the need for flow augmentation to comply with an interstate river compact and portending beneficial impacts on solute loading.

Colorado↗

USGS2018

Hydrogeology of, simulation of groundwater flow in, and potential effects of sea-level rise on the Kirkwood-Cohansey aquifer system in the vicinity of Edwin B. Forsythe National Wildlife Refuge, New Jersey

The Edwin B. Forsythe National Wildlife Refuge encompasses more than 47,000 acres of New Jersey coastal habitats, including salt marshes, freshwater wetlands, tidal wetlands, barrier beaches, woodlands, and swamps. The refuge is along the Atlantic Flyway and provides breeding habitat for fish, migratory birds, and other wildlife species. The refuge area may be threatened by global climate change, including sea-level rise (SLR). The Kirkwood-Cohansey aquifer system underlies the Edwin B. Forsythe National Wildlife Refuge. Groundwater is an important source of freshwater flow into the refuge, but information about the interaction of surface water and groundwater in the refuge area and the potential effects of SLR on the underlying aquifer system is limited. The U.S. Geological Survey (USGS), in cooperation with the U.S. Fish and Wildlife Service (USFWS), conducted a hydrologic assessment of the refuge in New Jersey and developed a groundwater flow model to improve understanding of the geohydrology of the refuge area and to serve as a tool to evaluate changes in groundwater-level altitudes that may result from a rise in sea level. Groundwater flow simulations completed for this study include a calibrated baseline simulation that represents 2005–15 hydraulic conditions and three SLR scenarios―20, 40, and 60 centimeters (cm) (0.656, 1.312, and 1.968 feet, respectively). Results of the three SLR simulations indicate that the water table in the unconfined Kirkwood-Cohansey aquifer system in the refuge area will rise, resulting in increased discharge of fresh groundwater to freshwater wetlands and streams. As sea level rises, simulated groundwater discharge to the salt marsh, bay, and ocean is projected to decrease. Flow from the salt marsh, bay, and ocean to the overlying surface water is projected to increase as sea level rises. The simulated movement of the freshwater-seawater interface as sea level rises depends on the hydraulic-head gradient. In the center of the Forsythe model area, topographic relief is 23 feet (ft) and the hydraulic-head gradient is 0.0033. In the center of the Forsythe model area, the simulated interface moved inland about 600 ft and downward about 15 ft from the baseline simulation to scenario 3 as a result of a SLR of 60 cm. In the southern part of the Forsythe model area, the topography is flatter (relief of 8 ft) and the hydraulic-head gradient is smaller (0.001). In the southern part of the Forsythe model study area, the simulated interface in this area is projected to move inland about 200 ft from the baseline simulation to scenario 3 and does not move downward.

New Jersey↗

USGS1987

Simulation of flow in the Edwards Aquifer, San Antonio region, Texas, and refinement of storage and flow concepts

The Edwards aquifer is a complexly faulted, carbonate aquifer lying within the Balcones fault zone of south-central Texas. The aquifer is recharged mainly by streamflow losses in the outcrop area of the Edwards aquifer and is discharged by major springs located at considerable distances, as much as 150 mi, from the areas of recharge, and by wells. Groundwater flow within the Edwards aquifer of the San Antonio region was simulated to investigate concepts relating to the storage and flow characteristics. A general purpose, finite difference model, modified to provide the capability of representing barrier faults, was used to simulate groundwater flow and storage in the aquifer. The simulations investigated the effects of complex geologic structures and significant changes in transmissivity, anisotropy, and storage coefficient, with initial values based on concepts developed in previous studies. Results of the simulations confirmed the original estimates of transmissivity values (> 100 sq ft/sec) in the confined zone of the aquifer between San Antonio and Comal Springs. A storage coefficient of 0.05 in the unconfined zone of the aquifer produced the best simulation of water levels and springflow. A major interpretation resulting from the simulations is that two essentially independent areas of regional flow were identified in the west and central part of the study area. Flow from the two areas converges at Comal Springs. The directions of computed flux vectors reflected the presence of major barrier faults which locally deflect patterns of groundwater movement. The most noticeable deflection is the convergence of flow through the geologic structural opening, the Knippa gap, in eastern Uvalde County. A second significant interpretation is that groundwater flow in northeastern Bexar, Comal, and Hays Counties is diverted by barrier faults toward San Marcos Springs, a regional discharge point. (Lantz-PTT)

Texas↗

USGS2024

Inset groundwater-flow models for the Cache and Grand Prairie Critical Groundwater Areas, northeastern Arkansas

The water resources in the Mississippi alluvial plain, located in parts of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas, supports a multibillion-dollar agricultural industry that relies heavily on pumping of groundwater for irrigation of crops and aquaculture. The primary source of groundwater for agricultural-related pumping is the Mississippi River Valley alluvial aquifer, which has declined in storage for decades; secondary groundwater sources include the middle Claiborne aquifer and Wilcox aquifer system. Two areas in northeastern Arkansas that lie within the Mississippi alluvial plain, part of the Cache and Grand Prairie regions, have been designated as Critical Groundwater Areas owing to decades of groundwater declines that resulted from past and current water use. The multidisciplinary Mississippi Alluvial Plain project, led by the U.S. Geological Survey, and funded by their Water Availability and Use Science Program, included objectives to develop numerical groundwater models in focus regions, including the part of the Cache and Grand Prairie regions of northeastern Arkansas. Two inset models were developed using the child model capabilities of MODFLOW 6, the U.S. Geological Survey’s Modular Hydrologic Model simulation software. Both models, called the Cache model and Grand Prairie model, simulated the groundwater system and surface-water/groundwater interactions for the Mississippi River Valley alluvial aquifer and underlying Tertiary-age aquifers and confining units to the Midway confining unit. Each model was spatially discretized into 500-meter x 500-meter orthogonal cells on a grid with 5-meter constant-thickness vertical layers that represented the Mississippi River Valley alluvial aquifer and increasing thickness layers for the aquifers and confining units below the alluvial aquifer. The Cache and Grand Prairie models were calibrated with the PEST++ iterative ensemble smoother Version 5 and employed high dimensional parameterization schemes of 13,740 and 30,436 parameters, respectively. The Cache mean absolute residual for groundwater-level observations within each model domain for the priority well was 1.58 meters. Grand Prairie mean absolute residuals for the alluvial aquifer and middle Claiborne aquifer groundwater-level observations were 2.71 and 10.78 meters, respectively. The groundwater budgets for the Cache and Grand Prairie models were characterized by substantial outflows to irrigation wells, which constituted about 52 and 54 percent of all outflows, with the primary source of water to those wells being releases from unconfined aquifer storage.

Arkansas↗

USGS2022

Groundwater resources of the Harney Basin, southeastern Oregon

Groundwater development has increased substantially in southeastern Oregon’s Harney Basin since 2010, mainly for the purpose of large-scale irrigation. Concurrently, some areas of the basin experienced groundwater-level declines of more than 100 feet, and some shallow wells have gone dry. The Oregon Water Resources Department has limited new groundwater development in the basin until an improved understanding of the groundwater-flow system is available. This report describes the results of a hydrologic investigation undertaken to provide that understanding. The investigation encompasses the groundwater hydrology of the entire 5,240-square-mile Harney Basin. Most of the precipitation in the Harney Basin falls in the higher-elevation areas of the Blue Mountains and Steens Mountain. Although considerable groundwater recharge occurs in these upland areas, most (83 percent) re-emerges as streams and springs in the uplands. Groundwater recharge in the lowlands is provided through infiltration of surface water flowing onto the lowlands from rivers and streams leaving the uplands and as groundwater flow from the surrounding upland rocks. Water-balance calculations indicate that the rate of groundwater recharge to the Harney Basin lowlands (where most groundwater is withdrawn) averages 173,000 acre-feet per year (acre-ft/yr). Groundwater in the Harney Basin lowlands mainly discharges through evapotranspiration from groundwater-irrigated (supplied from wells) crops or from natural vegetation drawing groundwater from the shallow water table and capillary fringe. Groundwater discharge in the lowlands is estimated to be about 283,000 acre-ft/yr, which exceeds the estimated groundwater recharge to the lowlands by about 110,000 acre-ft/yr. This imbalance results in removal of groundwater from storage in the aquifer system and is evidenced by the large declines observed in groundwater levels in the areas of greatest groundwater pumpage. To a large degree, the location and depth of pumpage dictate the timing and distribution of the effects of groundwater use in the Harney Basin. Pumpage is commonly greatest in the areas where higher-permeability geologic units allow for higher well yields. However, many of these higher-permeability units are bounded by lower-permeability units that cannot supply groundwater at a sufficient rate to replenish the areas of greatest pumpage, resulting in groundwater-level declines. Three Harney Basin areas with a combined area exceeding 140 square miles have experienced groundwater-level declines exceeding 40 feet compared to pre-development conditions: near the Weaver Spring/Dog Mountain area, in the northeastern floodplains along Highway 20, and near Crane. Areas of more modest groundwater-level decline (about 10 feet) were identified in the Virginia Valley area and the Silver Creek floodplain north of Riley. Smaller localized areas of groundwater-level depression have also formed around individual wells or groups of wells throughout the Harney Basin lowlands. Most groundwater being pumped from the Harney Basin lowlands, including all three areas experiencing large groundwater-level declines, was recharged more than 12,000 years ago, near the end of the last glacial period when the climate in the basin was cooler and wetter than today. Geochemical evidence indicates that modern recharge generally circulates to a depth no greater than 100 feet below the floodplains of major rivers and streams in the lowlands. Away from the major river and stream corridors, pre-modern water commonly is found at the water table. Recharge to groundwater and recovery of groundwater levels in the most heavily pumped areas in the Harney Basin lowlands are restricted by the limited spatial extent and depth of modern recharge in the Harney Basin lowlands and the relatively fine-grained deposits underlying most of the lowland areas.

Oregon↗

USGS2026

Groundwater source water assessment area delineation of public supply wells in the Treasure Valley, western Idaho

The Idaho Department of Environmental Quality (IDEQ) is responsible for assessing the source water for public water supply systems including the delineation of source water assessment areas for wells, which are defined as the subsurface areas contributing water to a well during a period of time. Part of the IDEQ assessment plan is to update source water assessments over time. The U.S. Geological Survey (USGS) developed a new groundwater-flow model for part of the western Snake River Plain aquifer in southwestern Idaho. This groundwater- flow model provided a new opportunity for IDEQ to update source water assessment for public water supply wells within the model boundary. The USGS, in cooperation with IDEQ, used the USGS Treasure Valley Groundwater Flow Model along with a MODPATH particle tracking post-processing program to reassess public-water supply source water assessment areas. The USGS created an automated script that requests a well location, depth, and pumpage rate, and then incorporates that information into the Treasure Valley Groundwater Flow Model combined with the MODPATH model to create 3-, 6-, and 10-year time of travel zones that are output to a geographic information system shapefile. The TVSWA_Delineator script provides a fast, transparent tool for IDEQ to reassess or create new public-water supply source water assessment area delineations based on new and updated information and is available as a Python script and a Windows executable with the compiled script and all necessary components. TVSWA_Delineator utilizes a quadtree approach to rediscretize the Treasure Valley Groundwater Flow Model grid around the pumping well location from 1-mile cell lengths telescoped down to 1/16 mile. This rediscretization increases the source water assessment area delineation resolution near the pumping well and minimizes the effect of weak sinks in the MODPATH model. The stress inputs for each month of the final 10 years (2006 to 2015) of the Treasure Valley Groundwater Flow model historical simulation were averaged for the particle tracking analysis. The MODPATH model used backwards tracking of particles released at each stress period to create the source water assessment area delineations. This report documents the creation of the TVSWA_Delineator Python script, including adaptations made to the Treasure Valley Groundwater Flow Model, creation of a companion MODPATH model, and post-processing results to provide source water assessment areas.

Idaho↗

USGS2017

Archaeal diversity and CO2 fixers in carbonate-/siliciclastic-rock groundwater ecosystems

Groundwater environments provide habitats for diverse microbial communities, and although Archaea usually represent a minor fraction of communities, they are involved in key biogeochemical cycles. We analysed the archaeal diversity within a mixed carbonate-rock/siliciclastic-rock aquifer system, vertically from surface soils to subsurface groundwater including aquifer and aquitard rocks. Archaeal diversity was also characterized along a monitoring well transect that spanned surface land uses from forest/woodland to grassland and cropland. Sequencing of 16S rRNA genes showed that only a few surface soil-inhabiting Archaea were present in the groundwater suggesting a restricted input from the surface. Dominant groups in the groundwater belonged to the marine group I (MG-I) Thaumarchaeota and the Woesearchaeota. Most of the groups detected in the aquitard and aquifer rock samples belonged to either cultured or predicted lithoautotrophs (e.g., Thaumarchaeota or Hadesarchaea). Furthermore, to target autotrophs, a series of 13 CO 2 stable isotope-probing experiments were conducted using filter pieces obtained after filtration of 10,000 L of groundwater to concentrate cells. These incubations identified the SAGMCG Thaumarchaeota and Bathyarchaeota as groundwater autotrophs. Overall, the results suggest that the majority of Archaea on rocks are fixing CO 2 , while archaeal autotrophy seems to be limited in the groundwater.

Archaea↗

USGS1996

Hydrogeology of the Hawaii Scientific Drilling Project borehole KP-1 1. Hydraulic conditions adjacent to the well bore

Temperature and formation resistivity logs obtained in borehole KP-1 of the Hawaii Scientific Drilling Project indicate that the adjacent formation is characterized by several zones of distinctly different average temperature and water salinity. A series of hydraulic analyses and water sampling programs were conducted to rule out the possibility of local hydraulic effects associated with the presence of the borehole in the generation of these apparent groundwater zones. Hydraulic tests and sampling with the borehole cased to a depth of 710 m and open below that depth indicate that the deep aquifer contains seawater at a temperature nearly identical to that of the open ocean at the same depth. Various analyses give estimates of aquifer transmissivity of about 10-3 m2/s in the vicinity of the borehole. Isolation of this deeper aquifer from the overlying groundwater zones was investigated by perforating the casing at six locations and then measuring the changes in water level in the borehole, in the salinity of the fluid column, in the temperature profile of the fluid column, and in the rate of flow in the fluid column induced by the perforations. These results positively confirm that the zones of distinctly different formation properties indicated on the temperature and resistivity logs are not caused by flow in or around casing. Flow and fluid column salinity induced by the perforations also confirm significant differences between the hydraulic heads and geochemistry of the different groundwater zones inferred from the well logs.

Journal of Geophysical Research B: Solid Earth↗

USGS1980

Hydrogeochemistry and simulated solute transport, Piceance Basin, northwestern Colorado

Oil-shale mining activities in Piceance basin in northwestern Colorado could adversely affect the ground- and surface-water quality in the basin. This study of the hydrology and geochemistry of the area used groundwater solute-transport-modeling techniques to investigate the possible impact of the mines on water quality. Maps of the extent and structure of the aquifer were prepared and show that a saturated thickness of 2,000 feet occurs in the northeast part of the basin. Ground-water recharge in the upland areas in the east, south, and west parts of the basin moves down into deeper zones in the aquifer and laterally to the discharge areas along Piceance and Yellow Creeks. The saline zone and the unsaturated zone provide the majority of the dissolved solids found in the ground water. Precipitation, ion-exchange, and oxidation-reduction reactions are also occurring in the aquifer. Model simulations of groundwater pumpage in tracts C-a and C-b indicate that the altered direction of groundwater movement near the pumped mines will cause an improvement in groundwater quality near the mines and a degradation of water quality downgradient from the tracts. Model simulations of mine leaching in tract C-a and C-b indicate that equal rates of mine leaching in the tracts will produce much different effects on the water quality in the basin. Tract C-a, by virtue of its remote location from perennial streams, will primarily degrade the groundwater quality over a large area to the northeast of the tract. Tract C-b, by contrast, will primarily degrade the surface-water quality in Piceance Creek, with only localized effects on the groundwater quality. (USGS)

Open-File Report↗

USGS2026

Multiple machine-learning estimation of groundwater levels and trends for the regional Mississippi River Valley alluvial aquifer

The Mississippi River Valley alluvial aquifer provides irrigation, public, and domestic water supplies across the south-central United States. Declining groundwater levels require improved characterization of changing conditions. Traditional potentiometric-surface mapping does not use all available water-level data or quantify uncertainty. To address these limitations, we developed a data-driven multiple machine-learning (MML) framework delivered through two open-source R packages. The covMRVAgen1 software assembles covariates to 155,960 monthly groundwater levels from 57,695 wells; the mmlMRVAgen1 software trains Cubist and Random Forest models, blends them, and makes 1-kilometer gridded predictions of monthly potentiometric surfaces for the period January 1980–December 2022. The MML approach provides a methodological foundation for region-scale spatiotemporal groundwater prediction and uncertainty quantification, generating 90-percent prediction limits with appropriate empirical coverage. Model performance is acceptable, with a root-mean-square error of about 4.2 feet, standard deviation of 24.82 feet, and a normalized Nash–Sutcliffe efficiency of 0.973.

Arkansas, Illinois, Louisiana, Mississippi, Missou↗

USGS2024

Simulation of groundwater flow and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River affects agricultural, municipal, and industrial water users and is an important concern in the Western United States. In the Paradox Valley of southwestern Colorado, natural discharge of sodium-chloride brine to the Dolores River from the underlying core of a salt-valley anticline accounts for about 6 percent of the salinity load to the Colorado River. Formation of the Paradox Valley began during the Miocene, and subsequent erosion exposed the Pennsylvania Paradox Formation in the core of the anticline where a cap rock, collapse features, breccia, and sodium-chloride saturated brine developed at the top of the exposed salt diapir. The discharge of brine to the Dolores River is affected by these dissolution features, along with seasonal hydrologic conditions and density-dependent flow between older dense brine and the younger fresh groundwater in the overlying alluvial aquifer. To reduce TDS concentrations in the Dolores River through the Paradox Valley, the Bureau of Reclamation has pumped brine from a series of shallow wells adjacent to the river since July 1996. The pumped brine is collected and piped to a deep disposal well where it is injected into the Mississippian Leadville Limestone at a depth of about 4,570-meters below land surface. The pumping and injection operation is collectively known as the Paradox Valley Unit (PVU), and by 2015, the PVU had substantially reduced TDS concentrations in the Dolores River by about 70 percent. Since 2019, injection-pressure limits and related seismic activity have constrained deep-well injection and thus brine pumping at the PVU. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey developed a MODFLOW-6 three-dimensional, variable-density groundwater flow and TDS transport model of the Paradox Valley to evaluate the effects of PVU pumping operations on brine discharge to the Dolores River and to guide additional research. The finite-difference model grid consists of 76 rows and 48 columns oriented from northwest to southeast in alignment with valley topography and groundwater-flow directions in the near-surface freshwater alluvial aquifer. A 7-layer hydrogeologic framework was developed from existing datasets to represent the alluvial aquifer, cap rock, collapse breccia, and groundwater flow and TDS transport from the underlying Paradox Formation salt to the Dolores River. The model represents a 33-year transient calibration period from 1987 through 2020 that includes pre-PVU conditions from 1987 through June 1996 and post-PVU conditions from July 1996 through 2020. A 1,000-year simulation of groundwater flow and coupled TDS transport computed the initial conditions for the subsequent 33-year transient simulation. Observations of precipitation, streamflow, evaporation, agricultural land use, and PVU brine pumping rates were used to specify appropriate boundary conditions to the model representing time-varying recharge, tributary streamflow, groundwater underflow, evapotranspiration (ET), and PVU pumping. Values for average monthly streamflow and TDS concentration at the upstream streamgage, the Dolores River at Bedrock (USGS streamgage 09169500), were specified as model input where the Dolores River enters Paradox Valley. Observed pumping from the PVU, water levels and TDS concentrations in groundwater, and streamflow and estimated TDS concentrations at the downstream streamgage, the Dolores River near Bedrock (USGS streamgage 09171100), were calibration targets that constrained the manual calibration of model parameters representing aquifer hydraulic conductivity, storage, streambed conductance, recharge, and (ET). Two primary model-calibration targets were the match between observed and simulated TDS mass flux from PVU pumping wells and the match between estimated and simulated TDS mass flux to the Dolores River. The simulated TDS mass withdrawn by pumping wells is calculated by the model as the product of the assigned pumping rate and simulated groundwater TDS concentrations. Because actual pumping rates were assigned as simulated values, the total simulated PVU pumping for the 33-year calibration is within 0.5 percent of the observed values. However, simulated concentrations and thus mass flux of TDS withdrawn by the PVU pumping wells were consistently about 26 percent less than observed values for all the simulated time periods (33-year simulation, pre-PVU, and post-PVU). The representation of brine inflow was explored through additional modeling to evaluate the effect of the simulated brine source on groundwater TDS concentrations. Results indicated that a saturated-salt constant-flux brine source best replicated the magnitude and transient pattern observed for TDS mass flux from PVU pumping wells. The simulated TDS mass flux to the Dolores River is compared to estimates based on observed streamflow and specific conductance (SC) data for the downstream streamgage. The calibrated model provided a close fit of simulated to measured streamflow at the downstream streamgage, and the calibrated model fit to estimated TDS concentrations at the downstream streamgage was reasonable. The greatest differences between simulated and estimated values occurred during drought periods from June 2000 to March 2003, May 2012 to June 2013, and October 2013 to October 2014, when simulated TDS concentrations in the river were greater than estimated concentrations. In general, simulated TDS mass flux to the river for the pre-PVU period is in good agreement with estimated values (2-percent difference), but the model overestimated TDS mass flux to the river by about 41 percent during the post-PVU period. The model uncertainty with respect to TDS mass flux to the river indicates other processes or model parameters not well represented by the model are affecting the system, especially during drought. During model calibration, the most sensitive parameters were identified as vertical hydraulic conductivity of the alluvial aquifer, conductance of the Dolores River streambed, ET extinction depth and rate, and recharge rate. Five 5-year scenarios of conditions for 2021–25 were simulated to assist evaluation of alternative strategies to manage the discharge of brine into the Dolores River. The first scenario simulates no PVU pumping and serves as a base case for comparison to the other scenarios. Two scenarios simulate the effects of varying withdrawal timing at an annual rate about one-third less than during 2010 through 2018. During high-flow spring snowmelt runoff periods when brine discharge is naturally minimized, PVU pumping does not substantially affect salinity in the Dolores River, and comparison of these two scenarios indicates that scheduling brine withdrawals during times of low river stage is nearly as effective at reducing TDS mass flux to the river as pumping brine year-round. Cessation of pumping during periods of high river stage may be advantageous for system maintenance, brine injection, and seismic-risk reduction. The fourth scenario tested the effect of reducing irrigation-return flow on brine discharge and predicted a slight reduction of TDS mass flux to the Dolores River, but not as great a reduction as that of using the PVU to remove brine. The fifth scenario simulated 5 years of drought conditions without PVU pumping and indicates brine discharge during drought about 15 percent greater than during average hydrologic conditions. Results from scenario 5 are consistent with the calibrated model results and indicate that aquifer properties and ET processes and parameters may be affecting simulation results during drought. The Paradox Valley groundwater model provides a reasonable overall match to observed conditions in the Dolores River. The model is useful for evaluating relative differences between brine management scenarios to inform PVU operational decisions and to identify gaps in data and process understanding. Representation of the brine source, hydraulic-conductivity parameters, and recharge and ET processes were identified as potential areas for additional field and modeling research. Additional research in the Paradox Valley might include field-data collection that provides additional information on the hydrogeologic framework, groundwater levels, groundwater TDS concentrations, stream characteristics, and aquifer properties. Additional modeling efforts could benefit from applying advanced tools for model development, calibration, and visualization including parameter-estimation and sensitivity analysis. Statistical evaluation of known model uncertainties such as hydraulic conductivity, streambed conductance, representations of the brine source, recharge, and ET could improve the match between simulated and estimated TDS mass flux from PVU pumping wells and to the Dolores River further informing model predictions and system understanding for the Paradox Valley.

Colorado↗

USGS2023

Determining seasonal recharge, storage changes, and specific yield using repeat microgravity and water-level measurements in the Mesilla Basin alluvial aquifer, New Mexico, 2016–2018

Increasing water demand and multi-year drought conditions within the Mesilla/Conejos-Médanos Basin near the New Mexico-Texas- Chihuahua border have resulted in diminished surface-water supplies and increased groundwater withdrawals. To better understand recharge to the shallow aquifer, the spatial and temporal groundwater storage changes, and the variability of specific yield (Sy) in the aquifer, seasonal groundwater elevation and repeat microgravity measurements were made during the irrigation release and non-release seasons of 2016, 2017, and 2018 at a network of locations near Las Cruces, New Mexico. The data collected during this investigation were able to capture seasonal change in groundwater elevations and storage from various sources of recharge at multiple sites in the shallow aquifer. Seasonal recharge in the study area was attributed to streamflow, the application and conveyance of irrigation water, and large or sustained precipitation events. However, increasing groundwater gradients in recent decades between piezometers close to the river and those more than a kilometer from the river suggests that recharge from river seepage has become localized at the seasonal scale. Overall, there was a net increase in storage of almost 8.4 cubic hectometers in the study reach between the start and end of the study, largely following the increased surface-water availability and above average precipitation in 2017. Specific yield, estimated by comparing the groundwater-level changes and storage changes at six sites in the study area, ranged from 0.14 (+/− 0.05) to 0.30 (+/− 0.06), which is slightly greater than previously reported estimates (0.10 to 0.25), but still within the error of the estimates. Most of the variability in the estimated storage change, that was not well-correlated with groundwater elevation change, is thought to be from soil moisture in the unsaturated zone. This investigation demonstrates the value of adding repeat microgravity measurements to conventional groundwater monitoring to better understand the sources and extent of recharge as well as the variability of Sy in the aquifer. Continued monitoring, under a variety of available surface water and meteorological conditions, could provide a more comprehensive understanding of the water budget and reduce the specific yield estimation uncertainty. Evaluating water-levels and storage conditions prior to, and following, local recharge events may help managers identify threshold conditions for aquifer storage depletions and recoveries.

New Mexico↗