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 91 records · Page 5Linked to original sources

USGS2019

Groundwater availability in the Ozark Plateaus aquifer system

Executive Summary The study described in this report, initiated by the U.S. Geological Survey in 2014, was designed to evaluate fresh groundwater resources within the Ozark Plateaus, central United States, as an area within a broader national assessment of groundwater availability. The goals of the Ozark study were to evaluate historical effects of human activities on water levels and groundwater availability, quantify groundwater resources now and under probable future pumping and climate conditions, and evaluate existing monitoring networks for their value in making better predictions of future groundwater resources. Previous studies include simulation of local-scale groundwater flow under varying temporal scales, or simulation of the regional system under steady-state conditions. While these studies are useful, particularly for the problem for which they were designed, there is a need to look at the larger regional system under transient conditions to fully evaluate the water resource over time. This study focused on multiple spatial and temporal scales to examine changes in groundwater pumping, storage, and water-level declines. The regional scale provides a broad view of the sources and demands on the system with time. The study area covers approximately 68,000 square miles in the central United States in parts of Missouri, Arkansas, Kansas, and Oklahoma and encompasses the Ozark Plateaus Physiographic Province (Ozark Plateaus), including the Salem Plateau, Springfield Plateau, and Boston Mountains. Groundwater is withdrawn from the Ozark Plateaus aquifer system (Ozark system) for public supply and for domestic, agriculture (including irrigation and aquaculture), livestock, and non-agricultural use (including industrial, thermoelectric power generation, mining, and commercial). The Ozark system provides an important drinking-water supply for people living in the Ozark Plateaus because public supply and domestic use combined constitute the largest groundwater use. Precipitation is the ultimate source of freshwater to the Ozark system; most rainfall occurs during April, May, and June, and precipitation increases generally from north to south across the study area. Groundwater use currently accounts for only 10 percent of the total water use in the areas overlying the Ozark system, but provides a critical drinking-water resource because public supply and domestic groundwater withdrawals are largely from groundwater resources. The 380 million gallons per day of groundwater withdrawn from the Ozark system in 2010 accounts for approximately 2 percent of recharge. Although groundwater use represents a small component of the hydrologic budget, because of low storage in aquifer units, cones of depression with steep water-level gradients can develop quickly around pumping centers. The amount of water entering and leaving the aquifer system from 1900 to about 1965 was relatively constant at a rate of about 13 billion gallons per day (Bgal/d). Much of this inflow of water is discharged through streams in the system to balance the hydrologic budget. Changes in storage over time (from outflows to inflows) reflect the large variability in recharge: if recharge decreases, water levels will decrease, resulting in less groundwater discharge to streams and more water released from aquifer storage. Conversely, when recharge increases, water levels increase, more groundwater discharges to streams, and aquifer storage is replenished. Although pumping generally increased from 1900 to 2016, it does not appear to correlate with the change in storage over the same time period. Regionally, simulated change in groundwater storage corresponds with changes in recharge, more so than with increases in pumping. Average recharge was 11.6 Bgal/d for the period 1900 to 2016. Recharge was generally above average from predevelopment to 1965, followed by a period of below-average recharge from 1965 to about 1980. Recharge remained consistently above average from 1980 to about 1988, after which there was a period of average or below-average recharge, reflected by a decline through the mid-2000s. The implications and potential effects of increased pumping and long-term climate change on the Ozark Plateaus hydrologic system and groundwater availability are a concern for communities and resource managers in the area. Pumping varies from year to year, but is generally expected to moderately increase with population, industrial, and agricultural needs. Most climate models predict warmer minimum and maximum air temperatures by midcentury in the Ozark Plateaus area, especially from midspring through early fall. Three scenarios were developed to simulate possible future conditions from 2016 to 2060 and assess the potential effects on the hydrologic system and availability of water resources. For each scenario, changes in water levels and hydrologic budget components were evaluated from predevelopment (1900) to present (2016) and 45 years into the future (2060). The baseline scenario represents an extension of the average (1996 to 2016) seasonal pumping and recharge values. The pumping scenario is an extension of the average (1996 to 2016) seasonal recharge values with increases in pumping following the historical trend for the period 2016–2060 of up to 120 percent of the 1996 to 2016 average seasonal pumping values. The general circulation model (GCM) scenario is an extension of the average (1996 to 2016) seasonal pumping values and variable recharge based on seasonal averages of soil water storage from a water-balance model using temperature and precipitation from multiple GCMs. The general patterns of water-level decline are similar for each scenario. The areas of water-level decline in southwest Missouri and northeast Oklahoma are only marginally different by 2060 from those of 2009. In one area south of Springfield, Mo., water-level declines are less in the baseline and GCM scenarios than in 2009. This may be the result of a transition from groundwater use to surface-water supplies for a larger percentage of the demand in the area. For all three scenarios, forecasted pumping, recharge, and aquifer properties play an important role in determining the uncertainty of water-level forecasts at 94 real-time observation wells. Simulated aquifer properties in the productive middle and lower Ozark aquifers and the St. Francois confining unit of the Ozark system contribute most to predictive uncertainty in water levels at approximately 35 percent of the real-time observation wells. Out of the 94 real-time observation wells, 82 are developed in the lower Ozark aquifer.

Arkansas, Kansas, Missouri, Oklahoma↗

USGS2016

Simulated effects of groundwater withdrawals from aquifers in Ocean County and vicinity, New Jersey

Rapid population growth since the 1930s in Ocean County and vicinity, New Jersey, has placed increasing demands upon the area’s freshwater resources. To examine effects of groundwater withdrawals, a three-dimensional groundwater-flow model was developed to simulate the groundwater-flow systems of five area aquifers: the unconfined Kirkwood-Cohansey aquifer system and Vincentown aquifer, and three confined aquifers— the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer. The influence of withdrawals is evaluated by using transient groundwater-flow model simulations that incorporate three withdrawal schemes. These are (1) no-withdrawal conditions; (2) 2000–03 withdrawal conditions, using reported monthly withdrawals at all production wells from January 2000 through December 2003; and (3) maximum-allocation withdrawal conditions using the maximum withdrawal allowed by New Jersey Department of Environmental Protection permits at each well. Particle tracking analysis, using results from model simulations, delineated particle flow paths from production wells to the point of recharge, and estimated particle travel times. Compared with no-withdrawal conditions, 2000–03 withdrawal conditions reduced the amount of groundwater flow out of the Kirkwood-Cohansey aquifer system into streams, increased the net flow of water into other layers, reduced net flow into or out of storage, and reduced flow from the Kirkwood-Cohansey aquifer system to constant head cells. Freshwater discharging to the Barnegat Bay-Little Egg Harbor estuary from streams and groundwater is essential to maintaining the ecology of the bay. Examination of selected stress periods indicates that simulated base flow in streams flowing into the Barnegat Bay-Little Egg Harbor estuary is reduced by as much as 49 cubic feet per second for 2000 to 2003 withdrawal conditions when compared with no-withdrawal conditions. In the three confined aquifers, water levels during periods of low recharge and high withdrawals, and high recharge and low withdrawals, were examined to determine seasonal effects on the confined flow systems. The simulated potentiometric surface of the Rio Grande water-bearing zone and the Atlantic City 800-foot sand during selected stress periods indicates substantial declines from no-withdrawal conditions to 2000–03 conditions as a result of groundwater withdrawals. Cones of depression in Toms River Township, Seaside Heights and Seaside Park Boroughs, and Barnegat Light Borough developed in the potentiometric surface of the Piney Point aquifer in response to withdrawals. Maximum-allocation withdrawals decreased flow out of the Kirkwood-Cohansey aquifer system to constant head cells, increased flow out of the aquifer system to adjacent and lower layers, and reduced groundwater discharge to streams when compared with 2000–03 withdrawal conditions. Increases in withdrawals from the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer result in an increase in simulated net groundwater flow into these aquifers. Base-flow reduction from 2000–03 conditions to maximum-allocation conditions of 25 to 29 cubic feet per second in all streams draining to the Barnegat Bay-Little Egg Harbor also is indicated. Potentiometric surfaces of the Rio Grande water-bearing zone, Atlantic City 800-foot sand, and the Piney Point aquifer during two stress periods of simulated maximum-allocation withdrawal conditions indicated the expansion of several cones of depression developed during 2000–03 withdrawals. Simulation of average 2000–03 withdrawal conditions indicated the extent to which the groundwater-flow system is susceptible to potential saltwater intrusion into near-shore wells. Travel time from recharge to discharge location ranged from 11 to approximately 50,700 years in near-shore Kirkwood-Cohansey aquifer system wells. Those in Seaside Heights Borough, in Island Beach State Park (Berkeley Township), and in Ship Bottom Borough have particle travel times from 140 to 12,000 years and flow paths that originated under Barnegat Bay or the Atlantic Ocean from the simulation of average maximum-allocation withdrawal conditions. Travel time along flow paths to wells screened in the Rio Grande water-bearing zone and the Atlantic City 800-foot sand from recharge to discharge point ranged from nearly 530 years to greater than 3.73 million years from the simulation of average 2000–03 withdrawal conditions. Particle tracking indicated that most wells screened in these aquifers derived a large part of their recharge from the Oswego River Basin, with a small portion of flow originating either beneath Barnegat Bay or to the east beneath the Atlantic Ocean. Travel time along flow paths that start beneath either Barnegat Bay or the Atlantic Ocean ranged from 2,300 to approximately 134,000 years from the simulation of average maximum-allocation withdrawal conditions."

New Jersey↗

USGS2014

Groundwater availability in the Crouch Branch and McQueen Branch aquifers, Chesterfield County, South Carolina, 1900-2012

Chesterfield County is located in the northeastern part of South Carolina along the southern border of North Carolina and is primarily underlain by unconsolidated sediments of Late Cretaceous age and younger of the Atlantic Coastal Plain. Approximately 20 percent of Chesterfield County is in the Piedmont Physiographic Province, and this area of the county is not included in this study. These Atlantic Coastal Plain sediments compose two productive aquifers: the Crouch Branch aquifer that is present at land surface across most of the county and the deeper, semi-confined McQueen Branch aquifer. Most of the potable water supplied to residents of Chesterfield County is produced from the Crouch Branch and McQueen Branch aquifers by a well field located near McBee, South Carolina, in the southwestern part of the county. Overall, groundwater availability is good to very good in most of Chesterfield County, especially the area around and to the south of McBee, South Carolina. The eastern part of Chesterfield County does not have as abundant groundwater resources but resources are generally adequate for domestic purposes. The primary purpose of this study was to determine groundwater-flow rates, flow directions, and changes in water budgets over time for the Crouch Branch and McQueen Branch aquifers in the Chesterfield County area. This goal was accomplished by using the U.S. Geological Survey finite-difference MODFLOW groundwater-flow code to construct and calibrate a groundwater-flow model of the Atlantic Coastal Plain of Chesterfield County. The model was created with a uniform grid size of 300 by 300 feet to facilitate a more accurate simulation of groundwater-surface-water interactions. The model consists of 617 rows from north to south extending about 35 miles and 884 columns from west to east extending about 50 miles, yielding a total area of about 1,750 square miles. However, the active part of the modeled area, or the part where groundwater flow is simulated, totaled about 1,117 square miles. Major types of data used as input to the model included groundwater levels, groundwater-use data, and hydrostratigraphic data, along with estimates and measurements of stream base flows made specifically for this study. The groundwater-flow model was calibrated to groundwater-level and stream base-flow conditions from 1900 to 2012 using 39 stress periods. The model was calibrated with an automated parameter-estimation approach using the computer program PEST, and the model used regularized inversion and pilot points. The groundwater-flow model was calibrated using field data that included groundwater levels that had been collected between 1940 and 2012 from 239 wells and base-flow measurements from 44 locations distributed within the study area. To better understand recharge and inter-aquifer interactions, seven wells were equipped with continuous groundwater-level recording equipment during the course of the study, between 2008 and 2012. These water levels were included in the model calibration process. The observed groundwater levels were compared to the simulated ones, and acceptable calibration fits were achieved. Root mean square error for the simulated groundwater levels compared to all observed groundwater levels was 9.3 feet for the Crouch Branch aquifer and 8.6 feet for the McQueen Branch aquifer. The calibrated groundwater-flow model was then used to calculate groundwater budgets for the entire study area and for two sub-areas. The sub-areas are the Alligator Rural Water and Sewer Company well field near McBee, South Carolina, and the Carolina Sandhills National Wildlife Refuge acquisition boundary area. For the overall model area, recharge rates vary from 56 to 1,679 million gallons per day (Mgal/d) with a mean of 737 Mgal/d over the simulation period (1900–2012). The simulated water budget for the streams and rivers varies from 653 to 1,127 Mgal/d with a mean of 944 Mgal/d. The simulated “storage-in term” ranges from 0 to 565 Mgal/d with a mean of 276 Mgal/d. The simulated “storage-out term” has a range of 0 to 552 Mgal/d with a mean of 77 Mgal/d. Groundwater budgets for the McBee, South Carolina, area and the Carolina Sandhills National Wildlife Refuge acquisition area had similar results. An analysis of the effects of past and current groundwater withdrawals on base flows in the McBee area indicated a negligible effect of pumping from the Alligator Rural Water and Sewer well field on local stream base flows. Simulate base flows for 2012 for selected streams in and around the McBee area were similar with and without simulated groundwater withdrawals from the well field. Removing all pumping from the model for the entire simulation period (1900–2012) produces a negligible difference in increased base flow for the selected streams. The 2012 flow for Lower Alligator Creek was 5.04 Mgal/d with the wells pumping and 5.08 Mgal/d without the wells pumping; this represents the largest difference in simulated flows for the six streams.

South Carolina↗

USGS2022

Mapping areas of groundwater susceptible to transient contamination events from rapid infiltration into shallow fractured-rock aquifers in agricultural regions of the conterminous United States

Current time-invariant groundwater vulnerability assessments may not capture intermittent contamination events in landscape areas that experience rapid infiltration following precipitation or snowmelt. Occurrences of rapid infiltration and intermittent degradation of groundwater quality are frequently reported in fractured-rock aquifers. This investigation identifies landscape areas underlain by fractured rock within the conterminous United States (CONUS) that may be susceptible to rapid infiltration and where groundwater is a principal source of water supply to the population. Our analysis shows that approximately 27 percent of the CONUS, corresponding to a population of approximately 150 million people, is both underlain by fractured rock and denoted as an area of significant groundwater use. The results of this survey identified shallow fractured-rock aquifers underlying glacial sediments in the upper Midwest and northeastern United States as areas that may be subject to rapid infiltration events. Additionally, aquifers associated with the early Mesozoic basins located in the northeastern and mid-Atlantic United States and bands of carbonate aquifers in the southeastern United States show high susceptibility to rapid infiltration. Index values used in this investigation indicate isolated areas in the western half of the United States also show high susceptibility to rapid infiltration. The isolated areas in Oklahoma, Texas, Arkansas, and southwestern Missouri correspond to karst regions of carbonate aquifers. The isolated areas showing high susceptibility to rapid infiltration and contamination from agricultural sources are locations where more detailed investigations of transient contamination events are warranted. This survey also addresses the potential for contaminant longevity in fractured-rock aquifers stemming from intermittent contamination events. Contaminants that can dissolve into the groundwater following infiltration may be introduced into fractures, and the dissolved constituents can diffuse from fractures into the porosity of the adjacent rock matrix. These constituents can then diffuse back into permeable fractures and adversely affect groundwater quality at downgradient locations over an extended time frame. Rock types with larger matrix porosities have the capacity to retain and then release larger quantities of dissolved constituents, resulting in longer residence times for dissolved groundwater contaminants. The magnitude of the dissolved contaminant concentration infiltrating to the water table will also dictate whether the contaminant concentration in the groundwater exceeds limits for human consumption over the duration of a contamination event. In general, sedimentary- and carbonate-rock aquifers have larger matrix porosities in comparison to igneous- and metamorphic-rock aquifers, and thus, they are more susceptible to longer contaminant residence times. Aquifers composed of sedimentary or carbonate rock constitute approximately 51 percent of the CONUS, and 19 percent of the CONUS is associated with sedimentary- or carbonate-rock aquifers that are of significance for groundwater use. Depending on the contaminants of concern and the concentration of the contaminants introduced into the groundwater from infiltrating water, it would be beneficial for investigations of susceptibility to rapid infiltration to also consider the potential for contaminant longevity. This investigation identifies areas of rapid infiltration into fractured rock using index values applied to the attributes (1) depth to the water table, (2) depth to bedrock, and (3) percentage of sand in soil, where larger index values indicate a greater susceptibility to rapid infiltration. These attributes are selected as the most likely factors that affect rapid infiltration to the water table. The combination of depth to water table and depth to bedrock highlight those aquifer settings that are characterized as shallow fractured-rock aquifers, where the water table may reside either in the bedrock or in overlying unconsolidated geologic materials. In addition, we consider the percentage of agricultural use as a land-use attribute when formulating an index of susceptibility to rapid infiltration and contamination. Agricultural areas are well recognized as nonpoint sources of contaminants that can affect groundwater quality because of seasonal amendments applied to the land surface. Rural agricultural areas are also characterized by septic tanks and leach fields for onsite treatment of wastewater, which may also be a source of contamination that may be introduced into the groundwater following precipitation or snowmelt events.

Conterminous United States↗

USGS2012

Simulated effects of alternative withdrawal strategies on groundwater flow in the unconfined Kirkwood-Cohansey aquifer system, the Rio Grande water-bearing zone, and the Atlantic City 800-foot sand in the Great Egg Harbor and Mullica River Basins, New Jersey

Groundwater is essential for water supply and plays a critical role in maintaining the environmental health of freshwater and estuarine ecosystems in the Atlantic Coastal basins of New Jersey. The unconfined Kirkwood-Cohansey aquifer system and the confined Atlantic City 800-foot sand are major sources of groundwater in the area, and each faces different water-supply concerns. The U.S. Geological Survey (USGS), in cooperation with the New Jersey Department of Environmental Protection (NJDEP), conducted a study to simulate the effects of withdrawals in the Kirkwood-Cohansey aquifer system, the Atlantic City 800-foot sand, and the Rio Grande water-bearing zone and to evaluate potential scenarios. The study area encompasses Atlantic County and parts of Burlington, Camden, Gloucester, Ocean, Cape May, and Cumberland Counties. The major hydrogeologic units affecting water supply in the study area are the surficial Kirkwood-Cohansey aquifer system, a thick diatomaceous clay confining unit in the upper part of Kirkwood Formation; the Rio Grande water-bearing zone; and the Atlantic City 800-foot sand of the Kirkwood Formation. Hydrogeologic data from 18 aquifer tests and specific capacity data from 230 wells were analyzed to provide horizontal hydraulic conductivity of the aquifers. Groundwater withdrawals are greatest from the Kirkwood-Cohansey aquifer system, and 65 percent of the water is used for public supply. Groundwater withdrawals from the Atlantic City 800-foot sand are about half those from the Kirkwood-Cohansey aquifer system. Ninety-five percent of the withdrawals from the Atlantic City 800-foot sand is used for public supply. Data from six streamgaging stations and 51 low-flow partial record sites were used to estimate base flow in the area. Base flow ranges from 60 to 92 percent of streamflow. A groundwater flow model of the Kirkwood-Cohansey aquifer system, the Rio Grande water-bearing zone, and the Atlantic City 800-foot sand was developed and calibrated using water-level data from 148 wells and base-flow data from 22 gaging or low-flow partial record stations. The Kirkwood-Cohansey aquifer system within the Great Egg Harbor River and the Mullica River Basins was simulated on a monthly basis from 1998 through 2006. An existing regional model of the New Jersey Coastal Plain was revised to provide boundary conditions for the Great Egg Harbor and Mullica River Basin model (referred to as the Great Egg-Mullica model). In the Great Egg-Mullica model, monthly groundwater recharge rates used in the model ranged from 10-15 inches per year in 2001 to 20-25 inches per year in 2005. The mean-absolute error for 10 of the 14 long-term hydrographs used in model calibration was less than 5 ft. Groundwater flow budgets for the Great Egg-Mullica model calibration periods, May 2005 and September 2006, and for the entire model calibration period 1998 to 2006, showed that nearly 70 percent of the water entering the Atlantic City 800-foot sand came from the horizontal connection with the Kirkwood-Cohansey aquifer system in updip areas. The groundwater flow model was used to simulate scenarios under three possible conditions: average 1998 to 2006 withdrawals (Average scenario), full-allocation withdrawals (Full Allocation scenario), and projected 2050-demand withdrawals (2050 Demand scenario). Withdrawals in the Full Allocation scenario are nearly twice the withdrawals from the Average scenario, primarily because of the potential for large agricultural withdrawals if all allocations are used. Withdrawals for the 2050 Demand scenario are about 50 percent greater than those for the Average scenario, primarily due to expected increases in withdrawals for public supply. Monthly base-flow depletion criteria were determined using the Low-Flow Margin method, currently under consideration by NJDEP, to estimate available water on an annual basis at the Hydrologic Unit Code 11 (HUC11) level and to determine whether a water-supply deficit exists. Simulations of various groundwater-withdrawal scenarios were made using the calibrated model, and results were compared with baseline conditions (no withdrawals) to determine where and when base-flow deficits may be occurring and may be expected to occur in the future. Scenarios were simulated to assess base-flow depletion that could occur from different groundwater-withdrawal situations. In the Average scenario, deficits occurred in 7 of the 14 subbasins. In the Full Allocation scenario, deficits occurred in 11 of the subbasins. In the 2050 Demand scenario, deficits occurred in 9 of the 14 subbasins. The largest deficits occurred in the Absecon Creek subbasin because the base-flow depletion criteria for this subbasin is small due to the surface-water diversions that are already occurring there and because existing groundwater withdrawals in the subbasin have resulted in base-flow depletion under current (1998-2006) conditions. Three adjusted scenarios, variations of the Average, Full Allocation, and 2050 Demand scenarios, were simulated; for the adjusted scenarios, the withdrawals were modified in stages with the intent to successively eliminate or minimize the base-flow deficits. Modifications included shifting withdrawals to a deeper part of the Kirkwood-Cohansey aquifer system, implementing seasonal conjunctive use of shallow and deep aquifers, and specifying reductions in withdrawals within a HUC11 subbasin in deficit. The adjusted scenarios are intended to show the relative effectiveness of each of the three approaches in reducing the deficits. Most of the deficits under the Average, Full Allocation, and 2050 Demand scenarios were eliminated by reductions in withdrawals or allocations. Shifting withdrawals to a deeper part of the Kirkwood-Cohansey aquifer system or seasonal conjunctive use did not eliminate deficits for any subbasin. Reductions in withdrawals accounted for more than 95 percent of the total reduction of deficits in all but one subbasin.

New Jersey↗

USGS2017

Groundwater quality in the shallow aquifers of the Tulare, Kaweah, and Tule Groundwater Basins and adjacent highlands areas, Southern San Joaquin Valley, California

Groundwater provides more than 40 percent of California’s drinking water. To protect this vital resource, the State of California created the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The Priority Basin Project of the GAMA Program provides a comprehensive assessment of the State’s groundwater quality and increases public access to groundwater-quality information. The shallow aquifers of the Tulare, Kaweah, and Tule groundwater basins and adjacent highlands areas of the southern San Joaquin Valley constitute one of the study units being evaluated.

California↗

USGS2019

Hydraulic and water-quality indicators of aquifer zones contributing groundwater flow to wells in the Santa Fe Group aquifer system near southeast Albuquerque, New Mexico, 2013–16

An ethylene dibromide (EDB) plume extends approximately 5,880 feet northeast from the Bulk Fuels Facility on Kirtland Air Force Base. The leading edge of the EDB plume is about 3,700 feet upgradient from several water-supply wells. The water-supply wells are screened in the upper Santa Fe Group aquifer system. Within the upper Santa Fe Group, two primary clay-rich layers, the A1 and A2 units, separate water-producing zones. The U.S. Geological Survey, in cooperation with the Albuquerque Bernalillo County Water Utility Authority and the U.S. Air Force, installed four sentinel well nests and two aquifer-test pumping wells between the EDB plume and the water-supply wells. The purpose of the sentinel wells is to provide early warning of EDB plume migration towards water-supply wells. The sentinel well nests include at least three wells that are screened above, in between, and below the A1 and A2 units. The two aquifer-test pumping wells, installed for performing hydraulic tests on the aquifer system, are screened across both clay layers, as are the nearby water-supply wells. Well-bore flow logging indicated that greater than 60 percent of groundwater flow to the wells was through the deepest interval below the lowermost clay layer (A1 unit). The interval between the A1 and A2 units is the second most productive interval. Water-quality data also indicated that water drawn from the aquifer-test pumping wells and previously studied nearby water-supply wells is most similar in composition to water from the sentinel wells screened in the middle and deep intervals.

New Mexico↗

USGS2013

Enhancements to the Mississippi Embayment Regional Aquifer Study (MERAS) groundwater-flow model and simulations of sustainable water-level scenarios

Arkansas continues to be one of the largest users of groundwater in the Nation. As such, long-term planning and management are essential to ensure continued availability of groundwater and surface water for years to come. The Mississippi Embayment Regional Aquifer Study (MERAS) model was developed previously as a tool to evaluate groundwater availability within the Mississippi embayment, which encompasses much of eastern Arkansas where the majority of groundwater is used. The Arkansas Water Plan is being updated for the first time since 1990 and serves as the State’s primary, comprehensive water-resources planning and guidance document. The MERAS model was selected as the best available tool for evaluation of specific water-use pumping scenarios that are currently being considered by the State of Arkansas. The model, developed as part of the U.S. Geological Survey Groundwater Resources Program’s assessment of the Nation’s groundwater availability, is proving to be invaluable to the State as it works toward development of a sustained yield pumping strategy. One aspect of this investigation was to evaluate multiple methods to improve the match of observed to simulated groundwater levels within the Mississippi River Valley alluvial and middle Claiborne (Sparta) aquifers in the MERAS model. Five primary methods were evaluated: (1) explicit simulation of evapotranspiration (ET), (2) upgrade of the Multi-Node Well (MNW2) Package, (3) geometry improvement within the Streamflow Routing (SFR) Package, (4) parameter estimation of select aquifer properties with pilot points, and (5) modification of water-use estimates. For the planning purposes of the Arkansas Water Plan, three scenarios were developed to evaluate potential future conditions: (1) simulation of previously optimized pumping values within the Mississippi River Valley alluvial and the middle Claiborne aquifers, (2) simulated prolonged effects of pumping at average recent (2000–5) rates, and (3) simulation of drawdown constraints on most pumping wells. The explicit simulation of ET indicated little, if any, improvement of model fit at the expense of much longer simulation time and was not included in further simulations. Numerous attempts to fully utilize the MNW2 Package were unsuccessful in achieving model stability, though modifications made to the water-use dataset remained intact. Final improvements in the residual statistics may be attributed to a single method, or a cumulative effect of all other methods (geometry improvement with the SFR Package, parameter estimation with pilot points, and modification of water-use estimates) attempted. The root mean squared error (RMSE) for all observations in the model is 22.65 feet (ft) over a range in observed hydraulic head of 741.66 ft. The RMSE for water-level observations in the Mississippi River Valley alluvial aquifer is 14.14 ft (an improvement of almost 3 ft) over a range in observed hydraulic head of 297.25 ft. The RMSE for the Sparta aquifer is 32.02 ft (an improvement of approximately 3 ft) over a range in observed hydraulic head of 634.94 ft. Three scenarios were developed to utilize a steady-state version of the MERAS model. Scenario 1 was developed to use pumping values resulting from the optimization of baseline rates (typically 1997 pumping rates) from previous optimization modeling of the alluvial aquifer and the Sparta aquifer. Scenario 2 was developed to evaluate the prolonged effects of pumping from the alluvial aquifer at recent pumping rates. Scenario 3A was designed to evaluate withdrawal limits from the alluvial aquifer by utilizing drawdown constraints equal to an altitude of approximately 50 percent of the predevelopment saturated thickness of the alluvial aquifer or 30 ft above the bottom of the alluvial aquifer, whichever was greater. The results of scenario 1 indicate large water-level declines throughout the area of the alluvial aquifer, regardless of the substitution of the optimized pumping values from earlier model simulations. The results of scenario 2 also indicate large areas of water-level decline, as compared to half of the saturated thickness, throughout the alluvial aquifer. The results of scenario 3A reveal some effects from the inclusion of multiple aquifers in a single simulation. The initial configuration of scenario 3A resulted in water levels well below the defined drawdown constraint, and some areas of depleted aquifer (water levels that are near or below the bottom of the aquifer) in east-central Arkansas. A fourth simulation (scenario 3B) was configured to apply the same drawdown constraints from the alluvial aquifer wells to the Sparta aquifer wells in the depleted area. These drawdown constraints reduce leakage from the alluvial aquifer to the underlying Sparta aquifer. This configuration did not produce depleted areas within the alluvial aquifer. Scenarios 3A and 3B indicate that even when pumping is limited in the alluvial aquifer, water levels in the alluvial aquifer may continue to decline in some areas because of pumping in the underlying Sparta aquifer.

Arkansas↗

USGS2016

The Maryland Coastal Plain Aquifer Information System: A GIS-based tool for assessing groundwater resources

Groundwater is the source of drinking water for &sim;1.4 million people in the Coastal Plain Province of Maryland (USA). In addition, groundwater is essential for commercial, industrial, and agricultural uses. Approximately 0.757 × 109 L d –1 (200 million gallons/d) were withdrawn in 2010. As a result of decades of withdrawals from the coastal plain confined aquifers, groundwater levels have declined by as much as 70 m (230 ft) from estimated prepumping levels. Other issues posing challenges to long-term groundwater sustainability include degraded water quality from both man-made and natural sources, reduced stream base flow, land subsidence, and changing recharge patterns (drought) caused by climate change. In Maryland, groundwater supply is managed primarily by the Maryland Department of the Environment, which seeks to balance reasonable use of the resource with long-term sustainability. The chief goal of groundwater management in Maryland is to ensure safe and adequate supplies for all current and future users through the implementation of appropriate usage, planning, and conservation policies. To assist in that effort, the geographic information system (GIS)–based Maryland Coastal Plain Aquifer Information System was developed as a tool to help water managers access and visualize groundwater data for use in the evaluation of groundwater allocation and use permits. The system, contained within an ESRI ArcMap desktop environment, includes both interpreted and basic data for 16 aquifers and 14 confining units. Data map layers include aquifer and confining unit layer surfaces, aquifer extents, borehole information, hydraulic properties, time-series groundwater-level data, well records, and geophysical and lithologic logs. The aquifer and confining unit layer surfaces were generated specifically for the GIS system. The system also contains select groundwater-quality data and map layers that quantify groundwater and surface-water withdrawals. The aquifer information system can serve as a pre- and postprocessing environment for groundwater-flow models for use in water-supply planning, development, and management. The system also can be expanded to include features that evaluate constraints to groundwater development, such as insufficient available drawdown, degraded groundwater quality, insufficient aquifer yields, and well-field interference. Ultimately, the aquifer information system is intended to function as an interactive Web-based utility that provides a broad array of information related to groundwater resources in Maryland’s coastal plain to a wide-ranging audience, including well drillers, consultants, academia, and the general public.

Special Paper of the Geological Society of America↗

USGS2015

Conceptual and numerical models of groundwater flow in the Ogallala and Arikaree aquifers, Pine Ridge Indian Reservation area, South Dakota, water years 1980-2009

The Ogallala and Arikaree aquifers are the largest sources of groundwater on the Pine Ridge Indian Reservation and are used extensively for irrigation and public and domestic water supplies. To assess the potential for decreased water levels and discharge to streams in the Pine Ridge Indian Reservation, conceptual and numerical models of groundwater flow in the Ogallala and Arikaree aquifers in southwestern South Dakota were developed by the U.S. Geological Survey in cooperation with the Oglala Sioux Tribe. The study area includes most of the Pine Ridge Reservation in Jackson and Shannon Counties and Indian trust lands in Bennett County in southwestern South Dakota. The High Plains aquifer, which includes the Ogallala and Arikaree aquifers, generally is less developed in South Dakota compared with other areas underlain by this aquifer; therefore, water levels in the High Plains aquifer in South Dakota generally fluctuated by less than 5 feet (ft) from 1980 to 1999. Despite minimal water-level changes in the High Plains aquifer in South Dakota, extensive withdrawals of groundwater for irrigation have caused water-level declines in many areas and increased concerns about the long-term sustainability of the aquifer; therefore, continued or increased withdrawals from the aquifer or prolonged drought may have the potential to affect water levels within the aquifer and discharge to important streams in the area. The Ogallala and Arikaree aquifers generally consist of poorly consolidated claystones, siltstones, sandstones, and shale deposited in fluvial and lacustrine environments. Saturated thicknesses ranged from 10 to 314 ft for the Ogllala aquifer and from 10 to 862 ft for the Arikaree aquifer. Previous hydraulic conductivity estimates ranged from less than 1 to 180 feet per day (ft/d) for the Ogallala aquifer and from less than 1 to 13 ft/d for the Arikaree aquifer. Recharge to the Ogallala and Arikaree aquifers is from precipitation on the outcrop areas, and discharge occurs through evapotranspiration, discharge to streams, and well withdrawals. Evapotranspiration generally occurs in topographically low areas along streams, and maximum evapotranspiration occurs when the water level is at the land surface. The generalized groundwater-flow direction is to the northeast with local flow towards streams. Precipitation for water years 1980–2009 ranged from about 11 to 39 inches per year (in/yr) and averaged about 19 in/yr. Estimated mean recharge for water years 1980–2009 was about 17.3 percent of precipitation for the Ogallala aquifer and 7.9 percent of precipitation for the Arikaree aquifer. The estimated mean maximum evapotranspiration for water years 1980–2009 was about 35 in/yr. Estimated mean base flow for gaged streams was about 0.06 cubic foot per second (ft 3 /s) per square mile of drainage area. Estimated mean total water use for water years 1980–2009 was 5.4 ft 3 /s from the Ogallala aquifer and 7.1 ft 3 /s from the Arikaree aquifer. A two-layer numerical groundwater-flow model was constructed using MODFLOW–NWT with a uniformly spaced grid consisting of 166 rows and 288 columns with cells 1,640 ft on a side. The numerical model of the Ogallala and Arikaree aquifers was used to simulate steady-state and transient conditions for water years 1980–2009. Model calibration was accomplished using the Parameter ESTimation (PEST) program that adjusted individual model input parameters and assessed the difference between estimated and model-simulated values of hydraulic head and base flow. Aquifer boundaries were no-flow on the northern and western sides and constant-head on the southern and eastern sides. The mean arithmetic difference was 1.4 ft between the 731 simulated and observed hydraulic heads in the Ogallala aquifer and 9.8 ft between the 2,754 simulated and observed hydraulic heads in the Arikaree aquifer. Simulated mean discharge from the Ogallala and Arikaree aquifers to selected stream reaches was 92.1 ft 3 /s compared to estimated discharge of 88.7 ft 3 /s. Calibrated recharge for the transient simulation averaged 3.3 in/yr for the Ogallala aquifer and 1.1 in/yr for the Arikaree aquifer. The mean maximum potential evapotranspiration rate was 35.4 in/yr. Streambed conductance for perennial stream reaches averaged 530 feet squared per day. Horizontal hydraulic conductivity averaged 27 ft/d for the Ogallala aquifer and 1.0 ft/d for the Arikaree aquifer. The vertical hydraulic conductivity averaged 1.4 ft/d for the Ogallala aquifer and 0.004 ft/d for the Arikaree aquifer. Specific yield for the Ogallala aquifer was 0.15 (dimensionless) and averaged 0.02 for the Arikaree aquifer. Specific storage for the Arikaree aquifer was 1.7x10 -6 per foot. Simulated steady-state model inflow and outflow was 459 ft 3 /s. The percentages of inflows were 17 percent from constant-head boundaries, 9 percent from streams, and 74 percent from recharge. Percentages of outflow were 8 percent to constant-head boundaries, 1 percent to wells, 31 percent to streams, and 59 percent to evapotranspiration. Simulated net inflow from the Ogallala aquifer to the Arikaree aquifer ranged from about 22 ft 3 /s in dry years to about 37 ft 3 /s in wet years. Two hypothetical future stress scenarios were simulated using input from the 30-year calibrated simulation of water years 1980–2009. The first hypothetical scenario represented an increase in groundwater withdrawals from 50 hypothetical production wells completed in the Arikaree aquifer. At the end of the 30-year hypothetical increased pumping simulation, water levels declined as much as 66 ft in the Arikaree aquifer, decreased discharge to streams accounted for about 26 percent (2.6 ft 3 /s) of increased withdrawals, and decreased evapotranspiration accounted for about 53 (5.3 ft 3 /s) percent of increased withdrawals. The second hypothetical scenario represented a 30-year period of decreased recharge (drought) by decreasing recharge 0.2 inch (24 ft 3 /s) for each water year. At the end of the hypothetical drought simulation, water levels declined as much as 10.9 ft in the Arikaree aquifer, decreased discharge to streams accounted for about 23 percent (5.5 ft 3 /s) of decreased recharge, and decreased evapotranspiration accounted for about 72 percent (17.3 ft 3 /s) of decreased recharge. The numerical model is a tool that could be used to better understand the flow system of the Ogallala and Arikaree aquifers, to approximate hydraulic heads in the aquifer, and to estimate discharge to rivers, springs, and seeps in the Pine Ridge Reservation area in Bennett, Jackson, and Shannon Counties. The model also is useful to help assess the response of the aquifer to additional stress, including potential increased well withdrawals and potential drought conditions.

South Dakota↗

USGS2017

Geophysics- and geochemistry-based assessment of the geochemical characteristics and groundwater-flow system of the U.S. part of the Mesilla Basin/Conejos-Médanos aquifer system in Doña Ana County, New Mexico, and El Paso County, Texas, 2010–12

One of the largest rechargeable groundwater systems by total available volume in the Rio Grande/Río Bravo Basin (hereinafter referred to as the “Rio Grande”) region of the United States and Mexico, the Mesilla Basin/Conejos-Médanos aquifer system, supplies water for irrigation as well as for cities of El Paso, Texas; Las Cruces, New Mexico; and Ciudad Juárez, Chihuahua, Mexico. The U.S. Geological Survey in cooperation with the Bureau of Reclamation assessed the groundwater resources in the Mesilla Basin and surrounding areas in Doña Ana County, N. Mex., and El Paso County, Tex., by using a combination of geophysical and geochemical methods. The study area consists of approximately 1,400 square miles in Doña Ana County, N. Mex., and 100 square miles in El Paso County, Tex. The Mesilla Basin composes most of the study area and can be divided into three parts: the Mesilla Valley, the West Mesa, and the East Bench. The Mesilla Valley is the part of the Mesilla Basin that was incised by the Rio Grande between Selden Canyon to the north and by a narrow valley (about 4 miles wide) to the southeast near El Paso, Tex., named the Paso del Norte, which is sometimes referred to in the literature as the “El Paso Narrows.” Previously published geophysical data for the study area were compiled and these data were augmented by collecting additional geophysical and geochemical data. Geophysical resistivity measurements from previously published helicopter frequency domain electromagnetic data, previously published direct-current resistivity soundings, and newly collected (2012) time-domain electromagnetic soundings were used in the study to detect spatial changes in the electrical properties of the subsurface, which reflect changes that occur within the hydrogeology. The geochemistry of the groundwater system was evaluated by analyzing groundwater samples collected in November 2010 for physicochemical properties, major ions, trace elements, nutrients, pesticides (reported but not used in the assessment), and environmental tracers. The data obtained from these samples (with the exception of the pesticide data) were used to gain insights into processes controlling the groundwater movement through the groundwater system in the study area. Results from the geophysical and geochemical assessments facilitated the interpretation of the geochemical characteristics of the groundwater sources and geochemical groups within the groundwater system. The groundwater-flow system in the study area consists primarily of the Mesilla Basin aquifer system, which can be divided into four hydrogeologic units by using an informal classification scheme based on basin-fill stratigraphy and sedimentology with an emphasis on aquifer characteristics. The four hydrogeologic units are (1) the Rio Grande alluvium, which is the shallow aquifer of the Mesilla Basin within the confines of the Mesilla Valley, and the three hydrogeologic units that compose the Santa Fe Group: (2) the lower part of the Santa Fe Group, which is the least productive zone, (3) the middle part of the Santa Fe Group, which is the primary water-bearing hydrogeologic unit in the basin and is generally saturated, and (4) the upper part of the Santa Fe Group, which is the most productive water-bearing unit within the Santa Fe Group but is only partially saturated in the north and largely unsaturated in the south and western parts of the Mesilla Basin. The helicopter frequency domain electromagnetic survey results indicated that approximately half of the resistivity values were less than 10 ohm-meters at depths of 50 and 100 feet with a transition where the resistivity values changed from relatively high values (greater than 20 ohm-meters) to relatively low resistivity values (less than 10 ohm-meters) near Vado, New Mexico. Slightly more than 25 percent of the gridded resistivity values from the three-dimensional grid of the combined inverse modeling results of the direct-current resistivity and time-domain electromagnetic soundings were equal to or less than 10 ohm-meters with large regions of low resistivity becoming apparent in the southernmost part of the study area near the Paso Del Norte where these low resistivity features are spatially the widest at or below the top of the bedrock. These low resistivity values might represent clayey deposits, sediments composed largely of sand and gravel saturated with saline water, or both. Historical dissolved-solids-concentration data within the surface geophysical subset area of the study area were compiled and compared to the inverse modeling results of the combined direct-current resistivity and time-domain soundings; this comparison was done to strengthen the interpretation made from the combined inverse modeling results that the low resistivity features were representative of sand and gravel deposits saturated with saline water and not clayey deposits. Water-level altitudes within the Rio Grande alluvium generally decreased from north to south, with a west to east decrease in water-level altitudes near Las Cruces, New Mexico, as a result of groundwater pumping. Groundwater flow within the Santa Fe Group is more complex than the groundwater flow within the Rio Grande alluvium because of the larger lateral and vertical extent of the Santa Fe Group compared to the Rio Grande alluvium. Groundwater from the Organ Mountains flows directly south towards the Paso del Norte. Groundwater from the Robledo Mountains, the Rough and Ready Hills, and the Sleeping Lady Hills generally flows to the southeast. Groundwater flowing near the north end of the midbasin uplift generally continues east towards the Rio Grande and then flows south on the east side of the midbasin uplift. Groundwater flowing near the west side of the midbasin uplift generally continues south parallel to the faults that make up the midbasin uplift and then flows east towards the Paso del Norte when it reaches the south end of the midbasin uplift. Groundwater from the Aden Hills and the East and West Potrillo Mountains flows to the south end of the midbasin uplift and then continues east towards the Paso del Norte. Throughout most of the Mesilla Valley, the vertical hydraulic gradient was downward because the water-level altitude in the Rio Grande alluvium was higher than it was in the Santa Fe Group, but in some areas (typically in the middle and southern parts of the Mesilla Valley), the vertical hydraulic gradient was substantially reduced or even reversed to an upward hydraulic gradient. The geochemistry data indicate that there was a complex system of multiple geochemical endmembers and mixing between these endmembers with recharge to the Rio Grande alluvium and Santa Fe Group composed mostly of seepage from the Rio Grande, inflows from deeper or neighboring water systems, and mountain-front recharge. Five distinct geochemical groups were identified in the Mesilla Basin study area: (1) ancestral Rio Grande (pre-Pleistocene) geochemical group, (2) modern Rio Grande (Pleistocene to present) geochemical group, (3) mountain-front geochemical group, (4) deep groundwater upwelling geochemical group, and (5) unknown freshwater geochemical group. The ancestral Rio Grande groundwater was water that recharged into the system as seepage losses from the ancestral Rio Grande; this groundwater generally flows from north to south-southeast towards the Paso del Norte. Groundwater on the west side of the midbasin uplift generally flows south until it reaches the southern part of the study area; from the southern part of the study area, the groundwater flows east towards the Paso del Norte. Groundwater on the east side of the midbasin uplift flows south-southeast towards the Paso del Norte where it mixes with groundwater from the modern Rio Grande, uplifted areas in the west, and the deep saline source. The water type of the modern Rio Grande geochemical group ranged from calcium-sulfate water type in the northern part of the study area to sodium-chloride-sulfate water type in the southern part of the study area; from north to south there was a substantial increase in specific conductance, strontium-87/strontium-86 ratio, potassium, and the trace metals of iron and lithium, changing the water chemistry such that it became similar to the water chemistry of the deep groundwater upwelling geochemical group. From age-dating results, water in the modern Rio Grande geochemical group was recharged to the Rio Grande alluvium within the past 10 years. The mountain-front geochemical group was generally old water (apparent age was greater than 10,000 carbon-14 years before present) that was somewhat mineralized and has relatively high concentrations of fluoride and silica, which might indicate longer exposure to volcanic and siliciclastic rocks or aluminosilicate minerals. There were five different locations of recharge determined from the groundwater geochemistry within the mountain-front geochemical group, all having a slightly different geochemical signature: (1) the Rough and Ready Hills, Robledo Mountains, and the Sleeping Lady Hills, (2) the Doña Ana Mountains, (3) the Aden Hills and West Potrillo Mountains, (4) the East Potrillo Mountains, and (5) the Sierra Juárez in Mexico. The groundwater from the Rough and Ready Hills, Robledo Mountains, the Sleeping Lady Hills, and the Doña Ana Mountains generally flows toward the Rio Grande and eventually mixes together and with the modern Rio Grande groundwater. The groundwater originating from the Aden Hills and East and West Potrillo Mountains generally flows east to southeast at a slow rate and eventually mixes and continues east, where it mixes with groundwater from the ancestral Rio Grande geochemical group and with the groundwater from the Sierra Juárez. The groundwater from the Sierra Juárez flows north and then east towards the Paso del Norte where it mixes with groundwater from the uplifted areas in the west, ancestral and modern Rio Grande groundwater, and the upwelling groundwater from a deep saline source. The deep groundwater upwelling geochemical group had the highest concentrations of bicarbonate, potassium, silica, aluminum, iron, and lithium within the study area, indicating that it had been in contact with carbonate and siliciclastic rocks for a much longer period of time and at higher temperatures compared to the other geochemical groups, and was most likely ancient marine groundwater originating from the Paleozoic and Cretaceous carbonate rocks which was upwelling into the Mesilla Basin aquifer system in the southeastern part of the study area through the extensive fault systems. Direct-current resistivity and time-domain electromagnetic soundings support the interpretation of ancient marine groundwater upwelling into the Mesilla Basin aquifer system, as do the analytical results from wells, and the helicopter frequency domain electromagnetic data collected along the Rio Grande. The hydrogen-2/hydrogen-1 ratio and oxygen-18/oxygen-16 ratio isotopic results for samples in the unknown freshwater geochemical group did not plot on the Rio Grande evaporation line, indicating this group did not have a Rio Grande signature (that is, there was no isotopic evidence of a component of Rio Grande water) and it also had the lowest mineralized content of any geochemical group in the study area.

New Mexico, Texas↗

USGS2007

New isotopic evidence for the origin of groundwater from the Nubian Sandstone Aquifer in the Negev, Israel

The geochemistry and isotopic composition (H, O, S, O sulfate , C, Sr) of groundwater from the Nubian Sandstone (Kurnub Group) aquifer in the Negev, Israel, were investigated in an attempt to reconstruct the origin of the water and solutes, evaluate modes of water–rock interactions, and determine mean residence times of the water. The results indicate multiple recharge events into the Nubian sandstone aquifer characterized by distinctive isotope signatures and deuterium excess values. In the northeastern Negev, groundwater was identified with deuterium excess values of ∼16‰, which suggests local recharge via unconfined areas of the aquifer in the Negev anticline systems. The δ 18 O H2O and δ 2 H values (−6.5‰ and −35.4‰) of this groundwater are higher than those of groundwater in the Sinai Peninsula and southern Arava valley (−7.5‰ and −48.3‰) that likewise have lower deuterium excess values of ∼10‰. Based on the geochemical differences between groundwater in the unconfined and confined zones of the aquifer, a conceptual geochemical model for the evolution of the groundwater in the Nubian sandstone aquifer has been reconstructed. The isotopic composition of shallow groundwater from the unconfined zone indicates that during recharge oxidation of pyrite to SO 4 ( δ 34 S SO4 ∼−13‰; δ 18 O SO4 ∼+7.7‰) and dissolution of CaCO 3 ( 87 Sr/ 86 Sr ∼0.70787; δ 13 C DIC = −3.7‰) occur. In the confined zone of the aquifer, bacterial SO 4 reduction removes a significant part of dissolved SO 4 2 - "> SO42- , thereby modifying its isotopic composition ( δ 34 S SO4 ∼−2‰; δ 18 O SO4 ∼+8.5‰) and liberating dissolved inorganic C that contains little or no radiocarbon ( 14 C-free) with low δ 13 C DIC values (<−12‰). In addition to local recharge, the Sr and S isotopic data revealed contribution of external groundwater sources to the Nubian Sandstone aquifer, resulting in further modifications of the groundwater chemical and isotopic signatures. In the northeastern Negev, it is shown that SO 4 -rich groundwater from the underlying Jurassic aquifer contributes significantly to the salt budget of the Nubian Sandstone aquifer. The unique chemical and isotopic composition of the Jurassic groundwater ( δ 34 S SO4 ∼ +14‰; δ 18 O SO4 ∼ 14‰; 87 Sr/ 86 Sr ∼0.70764) is interpreted as reflecting dissolution of Late Triassic marine gypsum deposits. In the southern Arava Valley the authors postulate that SO 4 -rich groundwater with distinctively high Br/Cl (3 × 10 −3 ) low 87 Sr/ 86 Sr (0.70734), and high δ 34 S SO4 values (+15‰) is derived from mixing with underlying brines from the Paleozoic units. The radiocarbon measurements reveal low 14 C activities (0.2–5.8 pmc) in both the northeastern Negev and southern Arava Valley. Taking into account dissolution of carbonate rocks and bacterial SO 4 reduction in the unconfined area, estimated mean residence times of groundwater in the confined zone in the northeastern Negev are on the order of 21–38 ka, which suggests recharge predominantly during the last glacial period. The 14 C signal in groundwater from the southern Arava Valley is equally low but due to evidence for mixing with external water sources the residence time estimates are questionable.

Applied Geochemistry↗

USGS2025

An inset groundwater-flow model to evaluate the effects of layering configuration on model calibration and assess managed aquifer recharge near Shellmound, Mississippi

The U.S. Geological Survey has developed a high-resolution inset groundwater-flow model in the Mississippi Delta as part of an interdisciplinary collaboration coordinated by the Mississippi Alluvial Plain project to provide a tool that stakeholders can use to support water-resource management decisions. Groundwater withdrawals from the Mississippi River Valley alluvial (MRVA) aquifer have been vital to support agricultural production in the region, but substantial groundwater-level declines near Shellmound, Mississippi, have caused concerns for long-term sustainability of the aquifer. To better understand the subsurface and try to mitigate the long-term groundwater-level declines, stakeholders have undertaken actions including a Groundwater Transfer and Injection Pilot (GTIP) project using a riverbank filtration-based managed aquifer recharge approach. The pilot project consisted of extracting groundwater near the Tallahatchie River and reinjecting it into the aquifer 3 kilometers west where water levels have substantially declined. A high-resolution airborne electromagnetic (AEM) survey was also completed to collect electrical resistivity data to support the GTIP project and the development of the groundwater model. The inset groundwater-flow model was developed to (1) integrate the AEM data into the optimal layering configuration of the MRVA aquifer that the available observation data can support through calibration, and (2) assess the potential effect of the GTIP project on the groundwater levels. The AEM data were processed into three different layering configurations leading to the development of model A (18 layers), model B (16 layers), and model C (8 layers), all at a 100- x 100-meter cell spatial resolution using the U.S. Geological Survey modular finite-difference flow model 6 code with Newton-Raphson formulation. The model development process integrated recent advances in modeling, such as the incorporation of AEM data, the use of outputs from the soil-water-balance (SWB) model, and the Aquaculture and Irrigation Water-Use Model, and was facilitated by robust automation using the open-source python packages Modflow-setup and SFRmaker. Using Parameter Estimation ++ Iterative Ensemble Smoother, the three numerical groundwater-flow models (models A, B, and C) were calibrated against a set of observations, which included aquifer groundwater levels, streamflows, stream stage, and aquifer transmissivity. Results indicate that the detailed representation of MRVA aquifer layers in model A produced the best calibrated model by history matching, and the integration of data representing surficial connectivity played a key role in improving groundwater recharge and enhancing the ability of the model to match groundwater levels in the cone of depression. A forecast model simulated the managed aquifer recharge approach, and the results indicated that, given average irrigation and recharge conditions (2010–15), the GTIP project has the potential to induce groundwater-level increases of as much as 3 meters around the injection site, but a sustained increase would require repetition in subsequent years of water transfer at 2022 rates or above.

Mississippi↗

USGS2012

Groundwater quality in the Columbia Plateau, Snake River Plain, and Oahu basaltic-rock and basin-fill aquifers in the Northwestern United States and Hawaii, 1992-2010

This assessment of groundwater-quality conditions of the Columbia Plateau, Snake River Plain, and Oahu for the period 1992–2010 is part of the U.S. Geological Survey’s National Water Quality Assessment (NAWQA) program. It shows where, when, why, and how specific water-quality conditions occur in groundwater of the three study areas and yields science-based implications for assessing and managing the quality of these water resources. The primary aquifers in the Columbia Plateau, Snake River Plain, and Oahu are mostly composed of fractured basalt, which makes their hydrology and geochemistry similar. In spite of the hydrogeologic similarities, there are climatic differences that affect the agricultural practices overlying the aquifers, which in turn affect the groundwater quality. Understanding groundwater-quality conditions and the natural and human factors that control groundwater quality is important because of the implications to human health, the sustainability of rural agricultural economies, and the substantial costs associated with land and water management, conservation, and regulation. The principal regional aquifers of the Columbia Plateau, Snake River Plain, and Oahu are highly vulnerable to contamination by chemicals applied at the land surface; essentially, they are as vulnerable as many shallow surficial aquifers elsewhere. The permeable and largely unconfined character of principal aquifers in the Columbia Plateau, Snake River Plain, and Oahu allow water and chemicals to infiltrate to the water table despite depths to water commonly in the hundreds of feet. The aquifers are essentially unconfined over large areas, having few extensive clay layers to impede infiltration through permeable volcanic rock and alluvial sediments. Agriculture is intensive in all three study areas, and heavy irrigation has imposed large artificial flows of irrigation recharge that rival or exceed natural recharge rates. Fertilizers and pesticides applied at land surface are leached from soil and transported to deep water tables with the infiltrating irrigation recharge, resulting in a layer of degraded water quality overlying better quality regional groundwater beneath. This “irrigation-recharge layer” is best known on Oahu, where it has been studied since the 1960s; however, the extent of nitrate and pesticide contamination in the Columbia Plateau and Snake River Plain indicate that the same situation exists in those areas. Contamination from agricultural and urban activities is present not only at shallow depths in surficial materials of the three areas, but extends regionally in the deep, principal bedrock aquifers that are tapped for drinking water by domestic and public-supply wells. Naturally occurring constituents and nitrate concentrations above human-health benchmarks—Maximum Contaminant Levels (MCLs), and Health-Based Screening Levels (HBSLs)—were more common in the Columbia Plateau and the Snake River Plain than in Oahu. Concentrations of anthropogenic constituents (constituents related to human activities) above human-health benchmarks were more common in Oahu. Naturally occurring contaminants, such as arsenic and radon, may be present in groundwater at concentrations of potential concern for human health in relatively undeveloped settings that otherwise may not be perceived as susceptible to contamination. Even though the median depth to groundwater in Oahu is more than 300 feet, the common occurrence of anthropogenic compounds in groundwater indicates that Oahu has a high susceptibility to contamination. Nitrate concentrations in groundwater were above the national background concentrations of 1 milligram per liter (mg/L) in all three study areas. In the Columbia Plateau, nitrate exceeded the human-health benchmark of 10 mg/L in 20 percent of the wells sampled. In the Snake River Plain, nitrate exceeded the human-health benchmark of 10 mg/L in 3 percent of the wells sampled. Nitrate can persist in groundwater for years and even decades in the oxygen-rich groundwater of the Columbia Plateau and the Snake River Plain, so prudent groundwater protection measures are critical to protect drinking water resources by reducing nitrate leaching from the land surface. Nitrate logistic regression models indicated that areas with a high percentage of land in crops (such as potatoes or sugarcane) and soils with low amounts of organic matter are most likely to have elevated nitrate concentrations in the groundwater. Areas where agricultural activities were absent had much lower probabilities of detecting elevated nitrate concentrations. The Columbia Plateau had a much higher probability of having elevated nitrate concentrations, with most of the land area having greater than a 50 percent probability of elevated nitrate concentrations. Oahu and the Snake River Plain had a much lower probability of having elevated nitrate concentrations because of their lower percentage of agricultural land. Pesticides were detected at many sites in groundwater of the Columbia Plateau, Snake River Plain, and Oahu but generally at low concentrations below human-health benchmarks. Atrazine and its degradate (a compound produced from the breakdown of a parent pesticide), deethylatrazine, were the most commonly detected pesticides in groundwater sampled in the Columbia Plateau and Snake River Plain. Bromacil was the most commonly detected pesticide on Oahu. The other pesticides most commonly detected in the study areas include simazine, hexazinone, metribuzin, diuron, prometon, metolachlor, p,p’ -DDE, dieldrin, 2-4-D, and alachlor. DDE (a degradate of DDT) and dieldrin are still being detected in groundwater despite having been banned for more than 30 years. Codetection of multiple pesticides in water from a single well was common. The widespread occurrence of pesticides in groundwater in the study areas indicates that the groundwater is highly susceptible to pesticide contamination. Some pesticides were detected in groundwater samples from all three study areas, but other pesticides were detected only in samples from Oahu, or only in samples from the Columbia Plateau and Snake River Plain. This is because some pesticides (such as atrazine) are broad-spectrum pesticides that are used on many crops in many different areas of the United States. Other pesticides (such as simazine, metribuzin, and metolachlor) are used on row crops (such as potatoes, barley, and alfalfa) grown in the Columbia Plateau and Snake River Plain, but not on pineapple or sugarcane grown in Oahu. Atrazine logistic-regression models indicate that areas with a high percentage of land in crops (such as potatoes or sugarcane), a low percentage of fallow land, and highly permeable soils with low amounts of organic matter are most likely to have atrazine detected in the groundwater. Areas where agricultural activities were absent had much lower probabilities of atrazine being detected. The Snake River Plain had a much higher probability of atrazine detections, with more than 50 percent of the land area having greater than a 50 percent probability of atrazine contamination. Oahu had a much lower probability of atrazine contamination, with only 24 percent of the land area having greater than a 50 percent probability of atrazine contamination. Oahu and the Columbia Plateau had some of the highest percentages of soil fumigant detections in groundwater in the United States. Soil fumigants are volatile organic compounds (VOCs) used as pesticides, which are applied to soils to reduce populations of plant parasitic nematodes (harmful rootworms), weeds, fungal pathogens, and other soil-borne microorganisms. They are used in Oahu and the Columbia Plateau on crops such as pineapple and potatoes. All three areas (Columbia Plateau, Snake River Plain, and Oahu) had fumigant concentrations exceeding human-health benchmarks for drinking water.

Hawai'i, Idaho, Nevada, Oregon, Utah, Washington, ↗

USGS2009

Hydrogeology and Simulation of Groundwater Flow in the Plymouth-Carver-Kingston-Duxbury Aquifer System, Southeastern Massachusetts

The glacial sediments that underlie the Plymouth-Carver-Kingston-Duxbury area of southeastern Massachusetts compose an important aquifer system that is the primary source of water for a region undergoing rapid development. Population increases and land-use changes in this area has led to two primary environmental effects that relate directly to groundwater resources: (1) increases in pumping that can adversely affect environmentally sensitive groundwater-fed surface waters, such as ponds, streams, and wetlands; and (2) adverse effects of land use on the quality of water in the aquifer. In response to these concerns, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, began an investigation in 2005 to improve the understanding of the hydrogeology in the area and to assess the effects of changing pumping and recharge conditions on groundwater flow in the Plymouth-Carver-Kingston-Duxbury aquifer system. A numerical flow model was developed based on the USGS computer program MODFLOW-2000 to assist in the analysis of groundwater flow. Model simulations were used to determine water budgets, flow directions, and the sources of water to pumping wells, ponds, streams, and coastal areas. Model-calculated water budgets indicate that approximately 298 million gallons per day (Mgal/d) of water recharges the Plymouth-Carver-Kingston-Duxbury aquifer system. Most of this water (about 70 percent) moves through the aquifer, discharges to streams, and then reaches the coast as surface-water discharge. Of the remaining 30 percent of flow, about 25 percent of the water that enters the aquifer as recharge discharges directly to coastal areas and 5 percent discharges to pumping wells. Groundwater withdrawals are anticipated to increase from the current (2005) rate of about 14 Mgal/d to about 21 Mgal/d by 2030. Pumping from large-capacity production wells decreases water levels and increases the potential for effects on surface-water bodies, which are affected by pumping and wastewater disposal locations and rates. Pumping wells that are upgradient of surface-water bodies potentially capture water that would otherwise discharge to these surface-water bodies, thereby reducing streamflow and pond levels. The areas most affected by proposed increases in groundwater withdrawals are in the Towns of Plymouth and Wareham where more than half of the proposed increase in pumping will occur. In response to an increase of about 7 Mgal/d of pumping, groundwater discharge to streams is reduced by about 6 cubic feet per second (ft3/s) (about 4 Mgal/d) from a total of about 325 ft3/s. Reduction in streamflow is moderated by an increase of artificial recharge from wastewater returned to the aquifer by onsite domestic septic systems and centralized wastewater treatment facilities. It is anticipated that about 3 Mgal/d of the 7 Mgal/d of increase in pumped water will be returned to the aquifer as wastewater by 2030. Currently (2005) about 3 percent of groundwater discharge to streams is from wastewater return flow to the aquifer during average conditions. During drought conditions, the component of streamflow augmented by wastewater return flow doubles as wastewater recharge remains constant and aquifer recharge rates decrease. Wastewater return flow, whether as direct groundwater discharge to streams or as an additional source of aquifer recharge, increases the height of the water table near streams, thereby moderating the effects of increased groundwater withdrawals on streamflow. An analysis of a simulated drought similar to the 1960s drought of record indicates that the presence of streams moderates the effects on water levels of reduced aquifer recharge. The area where water-table altitudes were least affected by drought was in the Weweantic River watershed in the Town of Carver. Water levels decreased by less than 2 feet from current average conditions compared to decreases of greater than 5

Scientific Investigations Report↗

USGS2011

Regional groundwater-flow model of the Redwall-Muav, Coconino, and alluvial basin aquifer systems of northern and central Arizona

A numerical flow model (MODFLOW) of the groundwater flow system in the primary aquifers in northern Arizona was developed to simulate interactions between the aquifers, perennial streams, and springs for predevelopment and transient conditions during 1910 through 2005. Simulated aquifers include the Redwall-Muav, Coconino, and basin-fill aquifers. Perennial stream reaches and springs that derive base flow from the aquifers were simulated, including the Colorado River, Little Colorado River, Salt River, Verde River, and perennial reaches of tributary streams. Simulated major springs include Blue Spring, Del Rio Springs, Havasu Springs, Verde River headwater springs, several springs that discharge adjacent to major Verde River tributaries, and many springs that discharge to the Colorado River. Estimates of aquifer hydraulic properties and groundwater budgets were developed from published reports and groundwater-flow models. Spatial extents of aquifers and confining units were developed from geologic data, geophysical models, a groundwater-flow model for the Prescott Active Management Area, drill logs, geologic logs, and geophysical logs. Spatial and temporal distributions of natural recharge were developed by using a water-balance model that estimates recharge from direct infiltration. Additional natural recharge from ephemeral channel infiltration was simulated in alluvial basins. Recharge at wastewater treatment facilities and incidental recharge at agricultural fields and golf courses were also simulated. Estimates of predevelopment rates of groundwater discharge to streams, springs, and evapotranspiration by phreatophytes were derived from previous reports and on the basis of streamflow records at gages. Annual estimates of groundwater withdrawals for agriculture, municipal, industrial, and domestic uses were developed from several sources, including reported withdrawals for nonexempt wells, estimated crop requirements for agricultural wells, and estimated per capita water use for exempt wells. Accuracy of the simulated groundwater-flow system was evaluated by using observational control from water levels in wells, estimates of base flow from streamflow records, and estimates of spring discharge. Major results from the simulations include the importance of variations in recharge rates throughout the study area and recharge along ephemeral and losing stream reaches in alluvial basins. Insights about the groundwater-flow systems in individual basins include the hydrologic influence of geologic structures in some areas and that stream-aquifer interactions along the lower part of the Little Colorado River are an effective control on water level distributions throughout the Little Colorado River Plateau basin. Better information on several aspects of the groundwater flow system are needed to reduce uncertainty of the simulated system. Many areas lack documentation of the response of the groundwater system to changes in withdrawals and recharge. Data needed to define groundwater flow between vertically adjacent water-bearing units is lacking in many areas. Distributions of recharge along losing stream reaches are poorly defined. Extents of aquifers and alluvial lithologies are poorly defined in parts of the Big Chino and Verde Valley sub-basins. Aquifer storage properties are poorly defined throughout most of the study area. Little data exist to define the hydrologic importance of geologic structures such as faults and fractures. Discharge of regional groundwater flow to the Verde River is difficult to identify in the Verde Valley sub-basin because of unknown contributions from deep percolation of excess surface water irrigation.

Scientific Investigations Report↗

USGS2023

Modeling groundwater-level responses to multiple stresses using transfer-function models and wavelet analysis in a coastal aquifer system

In coastal aquifers, dynamic stresses such as climate forcings, groundwater withdrawals, and ocean tidal fluctuations cause nonlinear responses to groundwater levels. Such responses to the stresses impact groundwater resources and related flooding and infrastructure risks at multiple scales. We used time-series models such as transfer-function models and wavelet analysis to quantify the relative contribution of these stresses to groundwater-level fluctuation in wells from the unconfined and confined aquifers in an Atlantic coastal aquifer. Climate forcings, such as precipitation and temperature, explained most of the groundwater-level variation for wells in the unconfined aquifer, whereas groundwater withdrawals were the dominant driver of groundwater levels for wells in the confined aquifer. The impact of groundwater withdrawals also was detected in several wells in the unconfined aquifer. Although the influence of ocean tides on groundwater levels commonly is observed in coastal aquifers, we found that daily groundwater withdrawals can obscure the semi-diurnal coherence signal of the two series. The magnitude of groundwater-level fluctuation that could be explained solely by tides was minor compared to that explained by climate or withdrawal stresses. Transfer-function modeling showed seasonal withdrawals from wells in confined aquifers had a significant, yet heterogeneous influence on groundwater levels in coastal aquifers, which highlights climate and withdrawals as key compounding stresses in coastal hydrology. This study demonstrates the value of time-series approaches to advance characterization of groundwater systems in areas with limited hydrogeologic parameter information.

Virginia↗

USGS2022

Composite regional groundwater hydrographs for selected principal aquifers in New Mexico, 1980–2019

Groundwater is an important source of water for New Mexico. An estimated 48 percent of the total water used comes from groundwater sources, and groundwater levels generally are declining over large areas of New Mexico. Groundwater levels are affected by local and regional recharge or discharge processes. Groundwater hydrographs show the history of groundwater-level changes at a well. A single hydrograph is not necessarily representative of the larger regional area; however, individual hydrographs from several wells can be combined into a composite hydrograph to show average groundwater changes for a regional area. The U.S. Geological Survey, in cooperation with the New Mexico Office of the State Engineer, has been measuring groundwater levels in a network of wells since about 1925. Although groundwater levels in the statewide well network have been measured at various frequencies, most wells have been measured in 5-year cycles since about 1980. The composite hydrographs in this report were developed to show groundwater-level changes for selected principal aquifers in New Mexico. Composite hydrographs were developed using wells in the Colorado Plateaus aquifers, the High Plains aquifer, the Pecos River Basin alluvial aquifer, the Rio Grande aquifer system, and the Roswell Basin aquifer system. Statewide, groundwater levels generally have declined or remained steady over the time period in aquifers analyzed for this study. The largest water-level declines occurred in the Colorado Plateaus and High Plains aquifers and in the Rio Grande aquifer system (north-central New Mexico), where median water-level declines ranged from 17 to 40 feet and mean water-level declines ranged from 3.8 to 32 feet. Groundwater-level declines (or rises) were generally smaller in other areas of New Mexico.

New Mexico↗