USGS Science⌕ Search

SEARCH · USGS Science

Results for “Groundwater”

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

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

At least 181 records · Page 10Linked to original sources

Simulation of groundwater flow and analysis of the effects of water-management options in the North Platte Natural Resources District, Nebraska

The North Platte Natural Resources District (NPNRD) has been actively collecting data and studying groundwater resources because of concerns about the future availability of the highly inter-connected surface-water and groundwater resources. This report, prepared by the U.S. Geological Survey in cooperation with the North Platte Natural Resources District, describes a groundwater-flow model of the North Platte River valley from Bridgeport, Nebraska, extending west to 6 miles into Wyoming. The model was built to improve the understanding of the interaction of surface-water and groundwater resources, and as an optimization tool, the model is able to analyze the effects of water-management options on the simulated stream base flow of the North Platte River. The groundwater system and related sources and sinks of water were simulated using a newton formulation of the U.S. Geological Survey modular three-dimensional groundwater model, referred to as MODFLOW–NWT, which provided an improved ability to solve nonlinear unconfined aquifer simulations with wetting and drying of cells. Using previously published aquifer-base-altitude contours in conjunction with newer test-hole and geophysical data, a new base-of-aquifer altitude map was generated because of the strong effect of the aquifer-base topography on groundwater-flow direction and magnitude. The largest inflow to groundwater is recharge originating from water leaking from canals, which is much larger than recharge originating from infiltration of precipitation. The largest component of groundwater discharge from the study area is to the North Platte River and its tributaries, with smaller amounts of discharge to evapotranspiration and groundwater withdrawals for irrigation. Recharge from infiltration of precipitation was estimated with a daily soil-water-balance model. Annual recharge from canal seepage was estimated using available records from the Bureau of Reclamation and then modified with canal-seepage potentials estimated using geophysical data. Groundwater withdrawals were estimated using land-cover data, precipitation data, and published crop water-use data. For fields irrigated with surface water and groundwater, surface-water deliveries were subtracted from the estimated net irrigation requirement, and groundwater withdrawal was assumed to be equal to any demand unmet by surface water. The groundwater-flow model was calibrated to measured groundwater levels and stream base flows estimated using the base-flow index method. The model was calibrated through automated adjustments using statistical techniques through parameter estimation using the parameter estimation suite of software (PEST). PEST was used to adjust 273 parameters, grouped as hydraulic conductivity of the aquifer, spatial multipliers to recharge, temporal multipliers to recharge, and two specific recharge parameters. Base flow of the North Platte River at Bridgeport, Nebraska, streamgage near the eastern, downstream end of the model was one of the primary calibration targets. Simulated base flow reasonably matched estimated base flow for this streamgage during 1950–2008, with an average difference of 15 percent. Overall, 1950–2008 simulated base flow followed the trend of the estimated base flow reasonably well, in cases with generally increasing or decreasing base flow from the start of the simulation to the end. Simulated base flow also matched estimated base flow reasonably well for most of the North Platte River tributaries with estimated base flow. Average simulated groundwater budgets during 1989–2008 were nearly three times larger for irrigation seasons than for non-irrigation seasons. The calibrated groundwater-flow model was used with the Groundwater-Management Process for the 2005 version of the U.S. Geological Survey modular three-dimensional groundwater model, MODFLOW–2005, to provide a tool for the NPNRD to better understand how water-management decisions could affect stream base flows of the North Platte River at Bridgeport, Nebr., streamgage in a future period from 2008 to 2019 under varying climatic conditions. The simulation-optimization model was constructed to analyze the maximum increase in simulated stream base flow that could be obtained with the minimum amount of reductions in groundwater withdrawals for irrigation. A second analysis extended the first to analyze the simulated base-flow benefit of groundwater withdrawals along with application of intentional recharge, that is, water from canals being released into rangeland areas with sandy soils. With optimized groundwater withdrawals and intentional recharge, the maximum simulated stream base flow was 15–23 cubic feet per second (ft 3 /s) greater than with no management at all, or 10–15 ft 3 /s larger than with managed groundwater withdrawals only. These results indicate not only the amount that simulated stream base flow can be increased by these management options, but also the locations where the management options provide the most or least benefit to the simulated stream base flow. For the analyses in this report, simulated base flow was best optimized by reductions in groundwater withdrawals north of the North Platte River and in the western half of the area. Intentional recharge sites selected by the optimization had a complex distribution but were more likely to be closer to the North Platte River or its tributaries. Future users of the simulation-optimization model will be able to modify the input files as to type, location, and timing of constraints, decision variables of groundwater withdrawals by zone, and other variables to explore other feasible management scenarios that may yield different increases in simulated future base flow of the North Platte River.

Nebraska↗

Hydrogeologic characterization of the Cahuilla Valley and Terwilliger Valley Groundwater Basins, Riverside County, California

The relation between the groundwater and the amount of natural recharge to the Cahuilla Valley and Terwilliger Valley groundwater basins is not well understood. During the 20th century, the reliance on groundwater near Anza, California, used for agricultural, domestic, and municipal reasons has increased, and there is the potential for changes in groundwater availability related to climate change. Several types of existing data were evaluated, and new data were collected for this study, with the goal of characterizing the region’s hydrogeology. The study’s scope included constructing a geologic framework model to show where the groundwater-bearing units are present and their relation to each other, estimating the major components of the groundwater budget, and understanding local short-term and regional long-term groundwater flow and how that has changed since the early 1900s. Two electrical resistivity tomography surveys were done in the Durasno Valley about 2,150 feet apart to identify the thickness of the alluvium, its horizontal extent, and the depth-to-basement along two profiles perpendicular to Cahuilla Creek. The subsurface sediments were mostly horizontally layered and the transitional boundary between the alluvium and basement was thinner and shallower along the upgradient profile where the depth-to-basement was about 70 feet below land surface; the depth-to-basement at the downgradient profile was more than about 140 feet below land surface. The results from the surveys were used to place four monitoring wells at two sites along the survey profiles. Artesian flow from the deepest well at the downgradient site indicated that the decomposed and competent basement likely contributed some groundwater to the overlying alluvium, laterally, from below, or both. A digital three-dimensional geologic framework model was constructed using EarthVision software to represent the subsurface geometry of the alluvium, decomposed basement, and competent basement. Maps and cross sections of the modeled thicknesses of the alluvium and decomposed basement, and the modeled elevation of the top of the competent basement, were made to show the subsurface geometry of vertical faults, selected wells, and the groundwater-bearing units. Because natural recharge is related to the variable cycles of precipitation, estimates are difficult to quantify. Recharge and runoff have extreme interannual variability in the study area; recharge and runoff can be sporadic, and a substantive amount may not occur in some years. Estimates of recharge from a previous study and the regional-scale Basin Characterization Model for California for four different periods ranged from 3,800 acre-feet/year for 1897–1947 to 5,900 acre-feet/year for 1971–2000. Potential recharge from the disposal of domestic septic systems may have been as much as 500 acre-feet in 2020. It was estimated that between about 400 and 2,400 acre-feet/year of groundwater is lost through evapotranspiration by vegetation and evaporation from open water bodies, but the main source of discharge is through pumpage, mainly used for agriculture from the alluvium in the Cahuilla Valley and Terwilliger Valley groundwater basins. The estimated total pumpage for 1991–2021 ranged from about 1,140 acre-feet in 2019 to about 3,450 acre-feet in 1994. When summed, the cumulative amount of estimated pumpage between 1991 and 2021 was about 81,400 acre-feet. The general direction of groundwater flow is from the northeast along the San Jacinto fault zone at the headwaters of Cahuilla and Hamilton Creeks, to the surface-water outlets at the west and southeast parts of the study area. Groundwater-level data from the 1950s and earlier indicate that there was a natural groundwater divide between the Cahuilla Valley and Terwilliger Valley groundwater basins, but the changing magnitude and extent of the groundwater depressions caused by pumping since about 1950 indicate that the location of the natural groundwater boundary between the Cahuilla Valley and Terwilliger Valley groundwater basins has migrated over time. Flow from the upper to the lower parts of the Cahuilla Valley groundwater basin roughly follows the course of Cahuilla Creek through the narrow Durasno Valley where an estimated volume of flow in April 2019 was about 10–150 acre-feet/year. Short-term trends in groundwater levels, particularly in wells where groundwater is shallow and in the basement unit, show how some areas respond quickly to recharge and discharge. Wells located further to the east within the Cahuilla Valley groundwater basin in the alluvium show much less of a response to recharge events; areas of sustained pumpage from the alluvium, primarily for agriculture, show long-term declines in groundwater levels and generally do not show the effects of storm events or recent runoff. Groundwater levels in wells that are farthest from where most of the recharge occurs and where pumping has been the greatest, had some of the largest long-term groundwater-level declines at a rate of about 0.8 foot/year between 1971 and 2021.

California↗

Status and understanding of groundwater quality in the San Diego Drainages Hydrogeologic Province, 2004: California GAMA Priority Basin Project

Groundwater quality in the approximately 3,900-square-mile (mi 2 ) San Diego Drainages Hydrogeologic Province (hereinafter San Diego) study unit was investigated from May through July 2004 as part of the Priority Basin Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is located in southwestern California in the counties of San Diego, Riverside, and Orange. The GAMA Priority Basin Project is being conducted by the California State Water Resources Control Board in collaboration with the U.S. Geological Survey (USGS) and the Lawrence Livermore National Laboratory. The GAMA San Diego study was designed to provide a statistically robust assessment of untreated-groundwater quality within the primary aquifer systems. The assessment is based on water-quality and ancillary data collected by the USGS from 58 wells in 2004 and water-quality data from the California Department of Public Health (CDPH) database. The primary aquifer systems (hereinafter referred to as the primary aquifers) were defined by the depth interval of the wells listed in the California Department of Public Health (CDPH) database for the San Diego study unit. The San Diego study unit consisted of four study areas: Temecula Valley (140 mi 2 ), Warner Valley (34 mi 2 ), Alluvial Basins (166 mi 2 ), and Hard Rock (850 mi 2 ). The quality of groundwater in shallow or deep water-bearing zones may differ from that in the primary aquifers. For example, shallow groundwater may be more vulnerable to surficial contamination than groundwater in deep water-bearing zones. This study had two components: the status assessment and the understanding assessment. The first component of this study-the status assessment of the current quality of the groundwater resource-was assessed by using data from samples analyzed for volatile organic compounds (VOC), pesticides, and naturally occurring inorganic constituents, such as major ions and trace elements. The status assessment is intended to characterize the quality of groundwater resources within the primary aquifers of the San Diego study unit, not the treated drinking water delivered to consumers by water purveyors. The second component of this study-the understanding assessment-identified the natural and human factors that affect groundwater quality by evaluating land use, well construction, and geochemical conditions of the aquifer. Results from these evaluations were used to help explain the occurrence and distribution of selected constituents in the study unit. Relative-concentrations (sample concentration divided by benchmark concentration) were used as the primary metric for relating concentrations of constituents in groundwater samples to water-quality benchmarks for those constituents that have Federal and (or) California benchmarks. For organic and special-interest constituents, relative-concentrations were classified as high (> 1.0), moderate (> 0.1 and &le;1.0), and low (&le;0.1). For inorganic constituents, relative concentrations were classified as high (> 1.0), moderate (> 0.5 and &le;1.0), and low (&le;0.5). Grid-based and spatially weighted approaches were then used to evaluate the proportion of the primary aquifers (aquifer-scale proportions) with high, moderate, and low relative-concentrations for individual compounds and classes of constituents. One or more of the inorganic constituents with health-based benchmarks were high (relative to those benchmarks) in 17.6 percent of the primary aquifers in the Temecula Valley, Warner Valley, and Alluvial Basins study areas (hereinafter also collectively referred to as the Alluvial Fill study areas because they are composed of alluvial fill aquifers), and in 25.0 percent of the Hard Rock study area. Inorganic constituents with health-based benchmarks that were frequently detected at high relative-concentrations included vanadium (V), arsenic (As), and boron (B). Vanadium and As concentrations were not significantly correlated to either urban or agricultural land use indicating natural sources as the primary contributors of these constituents to groundwater. The positive correlation of B concentration to urban land-use was significant which indicates that anthropogenic activities are a contributing source of B to groundwater. The correlation of V, As and B concentrations to pH was positive, indicating that in alkaline groundwater these constituents are being desorbed from, or being inhibited from adsorbing to, particle surfaces. Inorganic constituents with aesthetic benchmarks that were detected at high relative-concentrations include manganese (Mn), iron (Fe), and total dissolved solids (TDS). In the Alluvial Fill study areas, Mn and TDS were detected at high relative-concentrations in 13.7 percent of the primary aquifers, and Fe in 6.9 percent. In the Hard Rock study area, Mn was detected at high relative-concentrations in 33.3 percent of the primary aquifers, and TDS in 16.7 percent; Fe was not detected at high relative-concentrations. Total dissolved solids concentrations were significantly correlated to agricultural land use suggesting that agricultural practices are a contributing source of TDS to groundwater. Manganese and Fe concentrations were highest in groundwater with low dissolved oxygen and pH indicating that the reductive dissolution of oxyhydroxides may be an important mechanism for the mobilization of Mn and Fe in groundwater. TDS concentrations were highest in shallow wells and in modern (< 50 yrs) groundwater which indicates anthropogenic activities as a source of TDS concentrations in groundwater. The relative-concentrations of organic constituents with health-based benchmarks were high in 3.0 percent of the primary aquifers in the Alluvial Fill study areas. A single detection in the Alluvial Basins study area of the discontinued gasoline oxygenate methyl tert-butyl ether (MTBE) was the only organic constituent detected at a high relative-concentration; high relative-concentrations of these constituents were not detected in the Hard Rock study area. Twelve of 88 VOCs and 14 of 123 pesticides and pesticide degradates analyzed in grid wells were detected. Chloroform was the only VOC detected in more than 10 percent of the grid wells. The herbicides simazine, atrazine, and prometon were each detected in greater than 10 percent of the grid wells. Perchlorate was detected in 22 percent of the grid wells sampled. The understanding assessment showed a significant correlation of trihalomethanes (THMs) and solvents to urban land-use, indicating that detections of these constituents are more likely to occur in groundwater underlying urbanized areas of the study unit. MTBE concentrations were negatively correlated to the distance from the nearest leaking underground fuel tank, indicating that point sources are the most significant contributing factor for MTBE concentrations to groundwater in the study unit. The positive correlation of THM and herbicide concentrations to modern groundwater was significant, as was the negative correlation of herbicide concentrations to pH and anoxic groundwater. The negative correlation of herbicides to pH and anoxic groundwater was likely due to the fact that these constituents were detected more frequently in shallow wells where groundwater conditions tend to be oxic with relatively low pH.

California↗

Groundwater and surface-water interactions near White Bear Lake, Minnesota, through 2011

The U.S. Geological Survey, in cooperation with the White Bear Lake Conservation District, the Minnesota Pollution Control Agency, the Minnesota Department of Natural Resources, and other State, county, municipal, and regional planning agencies, watershed organizations, and private organizations, conducted a study to characterize groundwater and surface-water interactions near White Bear Lake through 2011. During 2010 and 2011, White Bear Lake and other lakes in the northeastern part of the Twin Cities Metropolitan Area were at historically low levels. Previous periods of lower water levels in White Bear Lake correlate with periods of lower precipitation; however, recent urban expansion and increased pumping from the Prairie du Chien-Jordan aquifer have raised the question of whether a decline in precipitation is the primary cause for the recent water-level decline in White Bear Lake. Understanding and quantifying the amount of groundwater inflow to a lake and water discharge from a lake to aquifers is commonly difficult but is important in the management of lake levels. Three methods were used in the study to assess groundwater and surface-water interactions on White Bear Lake: (1) a historical assessment (1978-2011) of levels in White Bear Lake, local groundwater levels, and their relation to historical precipitation and groundwater withdrawals in the White Bear Lake area; (2) recent (2010-11) hydrologic and water-quality data collected from White Bear Lake, other lakes, and wells; and (3) water-balance assessments for White Bear Lake in March and August 2011. An analysis of covariance between average annual lake-level change and annual precipitation indicated the relation between the two variables was significantly different from 2003 through 2011 compared with 1978 through 2002, requiring an average of 4 more inches of precipitation per year to maintain the lake level. This shift in the linear relation between annual lake-level change and annual precipitation indicated the net effect of the non-precipitation terms on the water balance has changed relative to precipitation. The average amount of precipitation required each year to maintain the lake level has increased from 33 inches per year during 1978-2002 to 37 inches per year during 2003-11. The combination of lower precipitation and an increase in groundwater withdrawals can explain the change in the lake-level response to precipitation. Annual and summer groundwater withdrawals from the Prairie du Chien-Jordan aquifer have more than doubled from 1980 through 2010. Results from a regression model constructed with annual lake-level change, annual precipitation minus evaporation, and annual volume of groundwater withdrawn from the Prairie du Chien-Jordan aquifer indicated groundwater withdrawals had a greater effect than precipitation minus evaporation on water levels in the White Bear Lake area for all years since 2003. The recent (2003-11) decline in White Bear Lake reflects the declining water levels in the Prairie du Chien-Jordan aquifer; increases in groundwater withdrawals from this aquifer are a likely cause for declines in groundwater levels and lake levels. Synoptic, static groundwater-level and lake-level measurements in March/April and August 2011 indicated groundwater was potentially flowing into White Bear Lake from glacial aquifers to the northeast and south, and lake water was potentially discharging from White Bear Lake to the underlying glacial and Prairie du Chien-Jordan aquifers and glacial aquifers to the northwest. Groundwater levels in the Prairie du Chien-Jordan aquifer below White Bear Lake are approximately 0 to 19 feet lower than surface-water levels in the lake, indicating groundwater from the aquifer likely does not flow into White Bear Lake, but lake water may discharge into the aquifer. Groundwater levels from March/April to August 2011 declined more than 10 feet in the Prairie du Chien-Jordan aquifer south of White Bear Lake and to the north in Hugo, Minnesota. Water-quality analyses of pore water from nearshore lake-sediment and well-water samples, seepage-meter measurements, and hydraulic-head differences measured in White Bear Lake also indicated groundwater was potentially flowing into White Bear Lake from shallow glacial aquifers to the east and south. Negative temperature anomalies determined in shallow waters in the water-quality survey conducted in White Bear Lake indicated several shallow-water areas where groundwater may be flowing into the lake from glacial aquifers below the lake. Cool lake-sediment temperatures (less than 18 degrees Celsius) were measured in eight areas along the northeast, east, south, and southwest shores of White Bear Lake, indicating potential areas where groundwater may flow into the lake. Stable isotope analyses of well-water, precipitation, and lake-water samples indicated wells downgradient from White Bear Lake screened in the glacial buried aquifer or open to the Prairie du Chien-Jordan aquifer receive a mixture of surface water and groundwater; the largest surface-water contributions are in wells closer to White Bear Lake. A wide range in oxygen-18/oxygen-16 and deuterium/protium ratios was measured in well-water samples, indicating different sources of water are supplying water to the wells. Well water with oxygen-18/oxygen-16 and deuterium/protium ratios that plot close to the meteoric water line consisted mostly of groundwater because deuterium/protium ratios for most groundwater usually are similar to ratios for rainwater and snow, plotting close to meteoric water lines. Well water with oxygen-18/oxygen-16 and deuterium/protium ratios that plot between the meteoric water line and ratios for the surface-water samples from White Bear Lake consists of a mixture of surface water and groundwater; the percentage of each source varies relative to its ratios. White Bear Lake is the likely source of the surface water to the wells that have a mixture of surface water and groundwater because (1) it is the only large, deep lake near these wells; (2) these wells are near and downgradient from White Bear Lake; and (3) these wells obtain their water from relatively deep depths, and White Bear Lake is the deepest lake in that area. The percentages of surface-water contribution to the three wells screened in the glacial buried aquifer receiving surface water were 16, 48, and 83 percent. The percentages of surface-water contribution ranged from 5 to 79 percent for the five wells open to the Prairie du Chien-Jordan aquifer receiving surface water; wells closest to White Bear Lake had the largest percentages of surface-water contribution. Water-balance analysis of White Bear Lake in March and August 2011 indicated a potential discharge of 2.8 and 4.5 inches per month, respectively, over the area of the lake from the lake to local aquifers. Most of the sediments from a 12.4-foot lake core collected at the deepest part of White Bear Lake consisted of silts, sands, and gravels likely slumped from shallower waters, with a very low amount of low-permeability, organic material.

Minnesota↗

Preliminary hydrogeologic assessment near the boundary of the Antelope Valley and El Mirage Valley groundwater basins, California

The increasing demands on groundwater for water supply in desert areas in California and the western United States have resulted in the need to better understand groundwater sources, availability, and sustainability. This is true for a 650-square-mile area that encompasses the Antelope Valley, El Mirage Valley, and Upper Mojave River Valley groundwater basins, about 50 miles northeast of Los Angeles, California, in the western part of the Mojave Desert. These basins have been adjudicated to ensure that groundwater rights are allocated according to legal judgments. In an effort to assess if the boundary between the Antelope Valley and El Mirage Valley groundwater basins could be better defined, the U.S. Geological Survey began a cooperative study in 2014 with the Mojave Water Agency to better understand the hydrogeology in the area and investigate potential controls on groundwater flow and availability, including basement topography. Recharge is sporadic and primarily from small ephemeral washes and streams that originate in the San Gabriel Mountains to the south; estimates range from about 400 to 1,940 acre-feet per year. Lateral underflow from adjacent basins has been considered minor in previous studies; underflow from the Antelope Valley to the El Mirage Valley groundwater basin has been estimated to be between 100 and 1,900 acre-feet per year. Groundwater discharge is primarily from pumping, mostly by municipal supply wells. Between October 2013 and September 2014, the municipal pumpage in the Antelope Valley and El Mirage Valley groundwater basins was reported to be about 800 and 2,080 acre-feet, respectively. This study was motivated by the results from a previously completed regional gravity study, which suggested a northeast-trending subsurface basement ridge and saddle approximately 3.5 miles west of the boundary between the Antelope Valley and El Mirage Valley groundwater basins that might influence groundwater flow. To better define potential basement structures that could affect groundwater flow between the groundwater basins in the study area, gravity data were collected using more closely spaced measurements in September 2014. Groundwater-level data was gathered and collected from March 2014 through March 2015 to determine depth to water and direction of groundwater flow. The gravity and groundwater-level data showed that the saturated thickness of the alluvium was about 2,000 feet thick to the east and about 130 feet thick above the northward-trending basement ridge near Llano, California. Although it was uncertain whether the basement ridge affects the groundwater system, a potential barrier to groundwater flow could be created if the water table fell below the altitude of the basement ridge, effectively causing the area to the west of the basement ridge to become hydraulically isolated from the area to the east. In addition, the direction of regional-groundwater flow likely will be influenced by future changes in the number and distribution of pumping wells and the thickness of the saturated alluvium from which water is withdrawn. Three-dimensional animations were created to help visualize the relation between the basins’ basement topography and the groundwater system in the area. Further studies that could help to more accurately define the basins and evaluate the groundwater-flow system include exploratory drilling of multi-depth monitoring wells; collection of depth-dependent water-quality samples; and linking together existing, but separate, groundwater-flow models from the Antelope Valley and El Mirage Valley groundwater basins into a single, calibrated groundwater-flow model.

California↗

Hydrogeologic and geochemical characterization of groundwater resources in Pine and Wah Wah Valleys, Iron, Beaver, and Millard Counties, Utah

Pine and Wah Wah Valleys are neighboring structural basins that encompass about 1,330 square miles in Beaver, Iron, and Millard Counties in Utah, approximately 50 miles northwest of Cedar City, Utah, and 50 miles southeast of Baker, Nevada. Perennial streamflow is limited and only exists in higher-altitude reaches of small mountain streams in both basins. Groundwater is in unconsolidated basin-fill aquifers and bedrock mountain aquifers. Groundwater in Pine and Wah Wah Valleys is being targeted for large-scale groundwater extraction and export to provide municipal supply to the growing population in Iron County, Utah. Concern about declining groundwater levels and spring flows from proposed groundwater withdrawals has increased interest in an improved understanding of the groundwater system. Previous studies have indicated that an average of 28,000 acre-feet per year of recharge occurs mostly as infiltration of precipitation in high-altitude regions in the two basins. Groundwater discharge in the mountain hydrologic systems was estimated to average 8,500 acre-feet per year and is assumed to be consumed before subsequently recharging the valley basin-fill aquifers. Subsurface groundwater outflow moves from basin-fill aquifers in Pine and Wah Wah Valleys northward to adjacent regional basins and was estimated to average 19,500 acre-feet per year. An updated water-level map for the basin-fill aquifers in Pine and Wah Wah Valleys indicates that groundwater moves northward along the lengths of both valleys toward adjacent basins. Measured depths to water range from about 210 to 750 feet below land surface in Wah Wah Valley, and from about 300 to 620 feet below land surface in Pine Valley. Long-term water levels at seven wells completed in the basin-fill aquifers of Pine and Wah Wah Valleys with records spanning more than 40 years are generally stable with observed fluctuations of less than 5 feet. Observed discharge from two springs monitored between 2013 and 2016 also is generally stable. Groundwater leaving Pine and Wah Wah Valleys through the subsurface moves northward, converges with regional groundwater flow, and discharges by evapotranspiration at regional groundwater discharge areas, likely Tule Valley, Utah. In this study, basin-scale groundwater discharge was estimated by (1) mapping the groundwater discharge areas in each valley; (2) evaluating the 2005–11 summer multispectral satellite images against the Basin and Range carbonate-rock aquifer system study evapotranspiration measurements to select scenes broadly representative of average conditions in the study area and partitioning the groundwater discharge areas into evapotranspiration units using the selected satellite images and field reconnaissance; and (3) scaling evapotranspiration to the evapotranspiration units using evapotranspiration-rate estimates from several studies in the Great Basin. The resulting updated estimates of average annual groundwater evapotranspiration in the Tule Valley and Sevier Lake groundwater discharge areas were 35,000 and 10,500 acre-feet per year, respectively, with a likely uncertainty of plus or minus 35 percent. Groundwater samples from 13 sites in Pine Valley and 11 sites in Wah Wah Valley were analyzed for major ions and nutrients, to characterize geochemistry and water quality. Groundwater samples also were analyzed for the stable isotopes of oxygen, hydrogen, and carbon, the radioactive isotopes of carbon and hydrogen, and dissolved noble gases including helium-3, helium-4, neon, argon, krypton and xenon. Groundwater sampling sites included 12 wells and 12 springs. Carbon-14 and tritium/helium groundwater age dating indicate that groundwater in the basin-fill aquifers is typically thousands to tens of thousands of years older than groundwater in the shallow mountain aquifers. Dissolved-solids concentrations are lower and noble-gas temperatures are warmer in the valley wells compared to almost all groundwater sampled from wells and springs in the surrounding mountains. These results indicate a hydraulic discontinuity between the mountain and valley aquifers throughout much of the study area, and that much of the valley recharge is not derived from direct infiltration of precipitation in the mountains.

Utah↗

Long-term groundwater availability in the Waihe‘e, ‘Īao, and Waikapū aquifer systems, Maui, Hawai‘i

Groundwater levels have declined since the 1940s in the Wailuku area of central Maui, Hawai‘i, on the eastern flank of West Maui volcano, mainly in response to increased groundwater withdrawals. Available data since the 1980s also indicate a thinning of the freshwater lens and an increase in chloride concentrations of pumped water from production wells. These trends, combined with projected increases in demand for groundwater in central Maui, have led to concerns over groundwater availability and have highlighted a need to improve general understanding of the hydrologic effects of proposed groundwater withdrawals in the Waihe‘e, ‘Īao, and Waikapū areas of central Maui. A numerical groundwater model was constructed to simulate the flow and salinity of groundwater in central Maui. The model simulates the effects of changes in groundwater withdrawals and recharge on water levels, freshwater-lens thicknesses, and chloride concentrations of pumped water from production wells. The model incorporates updated water-budget estimates of groundwater recharge from infiltration and direct recharge, seepage in stream channels, and inflow from inland areas. Mean annual groundwater recharge from infiltration and direct recharge was estimated using a daily water-budget model and the most current data, including the distributions of monthly rainfall and potential evapotranspiration, for the study area for nine historical periods from 1926 through 2012: 1926–69, 1970–79, 1980–84, 1985–89, 1990–94, 1995–99, 2000–04, 2005–09, and 2010–12. The water-budget model also estimated groundwater recharge based on one hypothetical scenario that used 1980–2010 rainfall and 2017 land cover. For the nine historical periods, estimated recharge from infiltration and direct recharge within the area of the groundwater model ranged from 30.4 million gallons per day (Mgal/d) during 2010–12 to 98.7 Mgal/d during 1926–69. Variability in recharge during these periods mainly reflects changes in rainfall and irrigation over time. Between 2010 and 2014, streamflow restoration in previously diverted streams resulted in an estimated increase in recharge from seepage in stream channels of about 12.5 Mgal/d. Average groundwater inflow of about 39.6 Mgal/d from inland, dike-intruded areas to the main area of interest was estimated from an existing island-wide numerical groundwater-flow model, which is at a larger scale and incorporates a greater number of simplifying assumptions. The numerical groundwater model developed for this study was calibrated to 1926–2012 transient water levels, vertical salinity profiles, and chloride concentrations of water pumped by production wells in the study area. The model was then used to evaluate one future recharge and six selected withdrawal scenarios, developed in consultation with the Maui Department of Water Supply, in terms of long-term changes in water level and 50-percent ocean-water salinity surface. The groundwater model was also used to simulate the future salinity of water withdrawn by existing and proposed production wells. The simulations were run to steady-state conditions, providing an estimate of the long-term effects of changes in withdrawal and recharge on the groundwater resource. Results of the simulated future withdrawal scenarios indicate that, relative to 2017–18 rates, the scenarios’ long-term effect of increased withdrawals ultimately leads to lower water levels and a higher 50-percent ocean-water salinity surface indicating a thinning of the freshwater lens. Results also indicate that the increased withdrawals produce some groundwater with chloride concentration below 250 milligrams per liter and some groundwater with higher chloride concentration. The amount of drawdown near production wells and the quality of water withdrawn from production wells is dependent on the rate and spatial distribution of the withdrawals. The model was also used to evaluate how groundwater availability may be affected for a drier recharge scenario based on a published study of future climate. Model results of the future recharge scenario indicate that the rate of groundwater recharge is a controlling factor for (1) water levels, (2) the 50-percent ocean-water salinity surface, and (3) the quality of water withdrawn from production wells in the Wailuku area. Coupled with reduced groundwater recharge (with all other factors remaining equal), the modeled future withdrawals in the scenario would tend to cause lower water levels, a higher 50-percent ocean-water salinity surface, and increased salinity of water withdrawn from production wells. The three-dimensional numerical groundwater model developed for this study utilizes the latest available hydrologic and geologic information and is a useful tool for understanding the long-term hydrologic effects of additional groundwater withdrawals in central Maui. The model has several limitations, including its non-uniqueness and inability to account for local-scale heterogeneities. Short-term effects of changes in recharge and withdrawals—and optimization of pumping rates to meet increased demand for water with acceptable salinity—are possible conditions for future simulation analyses.

Hawaii↗

Water availability and use pilot: Methods development for a regional assessment of groundwater availability, southwest alluvial basins, Arizona

Executive Summary: Arizona is located in an arid to semiarid region in the southwestern United States and is one of the fastest growing States in the country. Population in Arizona surpassed 6.5 million people in 2008, an increase of 140 percent since 1980, when the last regional U.S. Geological Survey (USGS) groundwater study was done as part of the Regional Aquifer System Analysis (RASA) program. The alluvial basins of Arizona are part of the Basin and Range Physiographic Province and cover more than 73,000 mi 2 , 65 percent of the State's total land area. More than 85 percent of the State's population resides within this area, accounting for more than 95 percent of the State's groundwater use. Groundwater supplies in the area are expected to undergo further stress as an increasing population vies with the State's important agricultural sector for access to these limited resources. To provide updated information to stakeholders addressing issues surrounding limited groundwater supplies and projected increases in groundwater use, the USGS Groundwater Resources Program instituted the Southwest Alluvial Basins Groundwater Availability and Use Pilot Program to evaluate the availability of groundwater resources in the alluvial basins of Arizona. The principal products of this evaluation of groundwater resources are updated groundwater budget information for the study area and a proof-of-concept groundwater-flow model incorporating several interconnected groundwater basins. This effort builds on previous research on the assessment and mapping of groundwater conditions in the alluvial basins of Arizona, also supported by the USGS Groundwater Resources Program. Regional Groundwater Budget: The Southwest Alluvial Basins-Regional Aquifer System Analysis (SWAB-RASA) study produced semiquantitative groundwater budgets for each of the alluvial basins in the SWAB-RASA study area. The pilot program documented in this report developed new quantitative estimates of groundwater budget components using recent (2000-2007) data and methods of data analysis. Estimates of inflow components, including mountain-front recharge, incidental recharge from irrigation of agriculture, managed recharge from recharge facilities, interbasin underflow from upgradient basins, and streamflow losses, are quantified for recent time periods. Mountain-front recharge is the greatest inflow component to the groundwater system and was estimated using two methods: a basin characteristic model and new precipitation information used in a previously developed regression equation. Annual mountain-front recharge for the study area for 1940-2007 estimated by the two methods is 730,000 acre-ft for the basin characteristic model and 643,000 acre-ft for the regression equation, representing 1.5 percent and 1.3 percent of precipitation, respectively. Outflow components, including groundwater withdrawals, evapotranspiration, and interbasin flow to downgradient basins, are also presented for recent time periods. Groundwater withdrawals accounted for the largest share of the water budget, with nearly 2.4 million acre-ft per year withdrawn from the study area in recent years. Evapotranspiration from groundwater was estimated at nearly 1.3 million acre-ft per year for the study area using a newly developed method incorporating vegetation indices from satellite images and land cover information. For water-budget components with temporal variation that could be assessed from available data, estimates for intervening time periods since before development were also developed. An estimate of aquifer storage change, representing both gains to and losses from the groundwater system since before development, was derived for the most developed basins in the study area using available estimates of groundwater-level changes and storage coefficients. An overall storage loss of 74.5 million acre-ft was estimated for these basins within the study area. Demonstration

Arizona↗

Hydrogeologic framework and occurrence, movement, and chemical characterization of groundwater in Dixie Valley, west-central Nevada

Dixie Valley, a primarily undeveloped basin in west-central Nevada, is being considered for groundwater exportation. Proposed pumping would occur from the basin-fill aquifer. In response to proposed exportation, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation and Churchill County, conducted a study to improve the understanding of groundwater resources in Dixie Valley. The objective of this report is to characterize the hydrogeologic framework, the occurrence and movement of groundwater, the general water quality of the basin-fill aquifer, and the potential mixing between basin-fill and geothermal aquifers in Dixie Valley. Various types of geologic, hydrologic, and geochemical data were compiled from previous studies and collected in support of this study. Hydrogeologic units in Dixie Valley were defined to characterize rocks and sediments with similar lithologies and hydraulic properties influencing groundwater flow. Hydraulic properties of the basin-fill deposits were characterized by transmissivity estimated from aquifer tests and specific-capacity tests. Groundwater-level measurements and hydrogeologic-unit data were combined to create a potentiometric surface map and to characterize groundwater occurrence and movement. Subsurface inflow from adjacent valleys into Dixie Valley through the basin-fill aquifer was evaluated using hydraulic gradients and Darcy flux computations. The chemical signature and groundwater quality of the Dixie Valley basin-fill aquifer, and potential mixing between basin-fill and geothermal aquifers, were evaluated using chemical data collected from wells and springs during the current study and from previous investigations. Dixie Valley is the terminus of the Dixie Valley flow system, which includes Pleasant, Jersey, Fairview, Stingaree, Cowkick, and Eastgate Valleys. The freshwater aquifer in the study area is composed of unconsolidated basin-fill deposits of Quaternary age. The basin-fill hydrogeologic unit can be several orders of magnitude more transmissive than surrounding and underlying consolidated rocks and Dixie Valley playa deposits. Transmissivity estimates in the basin fill throughout Dixie Valley ranged from 30 to 45,500 feet squared per day; however, a single transmissivity value of 0.1 foot squared per day was estimated for playa deposits. Groundwater generally flows from the mountain range uplands toward the central valley lowlands and eventually discharges near the playa edge. Potentiometric contours east and west of the playa indicate that groundwater is moving eastward from the Stillwater Range and westward from the Clan Alpine Mountains toward the playa. Similarly, groundwater flows from the southern and northern basin boundaries toward the basin center. Subsurface groundwater flow likely enters Dixie Valley from Fairview and Stingaree Valleys in the south and from Jersey and Pleasant Valleys in the north, but groundwater connections through basin-fill deposits were present only across the Fairview and Jersey Valley divides. Annual subsurface inflow from Fairview and Jersey Valleys ranges from 700 to 1,300 acre-feet per year and from 1,800 to 2,300 acre-feet per year, respectively. Groundwater flow between Dixie, Stingaree, and Pleasant Valleys could occur through less transmissive consolidated rocks, but only flow through basin fill was estimated in this study. Groundwater in the playa is distinct from the freshwater, basin-fill aquifer. Groundwater mixing between basin-fill and playa groundwater systems is physically limited by transmissivity contrasts of about four orders of magnitude. Total dissolved solids in playa deposit groundwater are nearly 440 times greater than total dissolved solids in the basin-fill groundwater. These distinctive physical and chemical flow restrictions indicate that groundwater interaction between the basin fill and playa sediments was minimal during this study period (water years 2009–11). Groundwater in Dixie Valley generally can be characterized as a sodium bicarbonate type, with greater proportions of chloride north of the Dixie Valley playa, and greater proportions of sulfate south of the playa. Analysis of major ion water chemistry data sampled during the study period indicates that groundwater north and south of Township 22N differ chemically. Dixie Valley groundwater quality is marginal when compared with national primary and secondary drinking-water standards. Arsenic and fluoride concentrations exceed primary drinking water standards, and total dissolved solids and manganese concentrations exceed secondary drinking water standards in samples collected during this study. High concentrations of boron and tungsten also were observed. Chemical comparisons between basin-fill and geothermal aquifer water indicate that most basin-fill groundwater sampled could contain 10–20 percent geothermal water. Geothermal indicators such as high temperature, lithium, boron, chloride, and silica suggest that mixing occurs in many wells that tap the basin-fill aquifer, particularly on the north, south, and west sides of the basin. Magnesium-lithium geothermometers indicate that some basin-fill aquifer water sampled for the current study likely originates from water that was heated above background mountain-block recharge temperatures (between 3 and 15 degrees Celsius), highlighting the influence of mixing with warm water that was possibly derived from geothermal sources.

Nevada↗

Simulation of groundwater withdrawal scenarios for the Redwall-Muav and Coconino Aquifer Systems of northern and central Arizona

The Northern Arizona Regional Groundwater Flow Model was used to estimate the hydrologic changes, including water-level change and groundwater discharge to streams and springs, that may result from future changes in groundwater withdrawals in and near the Coconino Plateau Water Advisory Council study area, Coconino and Navajo Counties, Arizona. Three future groundwater withdrawal scenarios for tribal and nontribal uses were developed by the Coconino Plateau Water Advisory Council and were simulated for the period representing the years from 2006 through 2105. Scenario 1 assumes no major changes in groundwater use except for increased demand based on population projections. Scenario 2 assumes that a pipeline will provide a source of surface water from Lake Powell to areas near Cameron and Moenkopi that would replace local groundwater withdrawals. Scenario 3 assumes that the pipeline extends to the Flagstaff and Williams areas, and would replace groundwater demands for water in the area. The Coconino Plateau Water Advisory Council withdrawal scenarios primarily influence water levels and groundwater discharge in the Coconino Plateau basin, near the western margin of the Little Colorado River Plateau basin, and the Verde Valley subbasin. Simulated effects of the withdrawal scenarios are superimposed on effects of previous variations in groundwater withdrawals and artificial and incidental recharge. Pre-scenario variations include changes in water-levels in wells; groundwater storage; discharge to streams and springs; and evapotranspiration by plants that use groundwater. Future variations in groundwater discharge and water-levels in wells will continue to occur as a result of both the past and any future changes. Water-level variations resulting from post-2005 stresses, including groundwater withdrawals and incidental and artificial recharge, in the area of the withdrawal scenarios are primarily localized and superimposed on the regional changes caused by variations in stresses that occurred since the beginning of the initial stresses in the early 1900s through 2005. Withdrawal scenario 1 produced a broad region on the Coconino Plateau where water-levels declined 3–5 feet by 2105, and local areas with water-level declines of 100 feet or more where groundwater withdrawals are concentrated, near the City of Flagstaff Woody Mountain and Lake Mary well fields, and the towns of Tusayan, Williams, and Moenkopi. Water-level rises of 100 feet or more were simulated at areas of incidental recharge near wastewater treatment facilities near Flagstaff, Tusayan, Grand Canyon South Rim, Williams, and Munds Park. Simulated water-level change from 2006 through 2105 for scenarios 2 and 3 is mostly different from water-level change simulated for scenario 1 at the local level. For scenarios 2 and 3, water levels near Cameron in 2105 where 1–3 feet higher than simulated for scenario 1. Water levels at Moenkopi are more than 100 feet higher due to the elimination of a proposed withdrawal well that was simulated in scenario 1. Scenario 3 eliminates more groundwater withdrawals in the Flagstaff and Williams areas, simulates 1–3 feet less water-level decline than scenario 1 across much of the Coconino Plateau, and water levels that are as much as 50 feet higher than simulated by scenario 1 near withdrawal wells in the Williams and Flagstaff areas. Scenario 1 simulated the most change in groundwater discharge for the Little Colorado River below Cameron and for Oak Creek above Page Springs where declines in discharge of about 1.3 and 0.9 cubic feet per second (ft 3 /s), respectively, were simulated. Other simulated changes in discharge through 2105 in scenario 1 are losses of less than 0.4 ft 3 /s at the Upper Verde River, losses of less than 0.3 ft 3 /s at Havasu Creek and at Colorado River below Havasu Creek, losses of less than 0.1 ft 3 /s at Clear Creek, and increases in flow at the south rim springs and Chevelon Creek of less than 0.1 and 0.3 ft 3 /s, respectively. Simulated changes in discharge for scenarios 2 and 3 are less than for scenario 1 because of lower rates of groundwater withdrawal. Scenario 3 resulted in greater groundwater discharge than scenarios 1 and 2 at all major groundwater discharge features from 2006 through 2105 except for Clear and Chevelon Creeks, where the same groundwater discharge was simulated by each of the three scenarios. Changes in groundwater discharge are expected to occur after 2105 to all major surface features that discharge from the Redwall-Muav and Coconino aquifers because change in aquifer storage was occurring at the end of the simulation in 2105. The accuracy of simulated changes resulting from the Coconino Plateau Water Advisory Council groundwater withdrawal scenarios is dependent on the persistence of several hydrologic assumptions that are inherent in the Northern Arizona Regional Groundwater Flow Model including, but not limited to, the reasonably accurate simulation of (1) transmissivity distributions, (2) distributions of vertical hydraulic properties, (3) distributions of spatial rates of withdrawal and incidental recharge, (4) aquifer extents, and (5) hydrologic barriers and conduits.

Arizona↗

Simulation of groundwater flow and pumping scenarios for 1900–2050 near Mount Pleasant, South Carolina

Groundwater withdrawals from the Upper Cretaceous-age Middendorf aquifer in South Carolina have created a large, regional cone of depression in the potentiometric surface of the Middendorf aquifer in Charleston and Berkeley Counties, South Carolina. Groundwater-level declines of as much as 249 feet have been observed in wells over the past 125 years and are a result of groundwater use for public water supply, irrigation, and private industry. To address the concerns of users of the Middendorf aquifer, the U.S. Geological Survey, in cooperation with Mount Pleasant Waterworks (MPW), recalibrated an existing groundwater-flow model to incorporate additional groundwater-use and water-level data since 2008. This recalibration process consisted of a technique of parameter estimation that uses regularized inversion and employs “pilot points” for spatial hydraulic property characterization. The groundwater-flow system of the Coastal Plain physiographic province of South Carolina and parts of Georgia and North Carolina was simulated using the U.S. Geological Survey finite-difference computer code MODFLOW-2000. After the model recalibration, the following six predictive water-management scenarios were created to simulate potential changes in groundwater flow and groundwater-level conditions in the Mount Pleasant, South Carolina, area: Scenario 1—maximize MPW reverse-osmosis plant capacity by increasing groundwater withdrawals from the Middendorf aquifer from 3.9 million gallons per day (Mgal/d), which was the amount withdrawn in 2015, to 8.58 Mgal/d; Scenario 2—same as Scenario 1, but with the addition of a 0.5 Mgal/d supply well in the Middendorf aquifer near Moncks Corner, South Carolina; Scenario 3—same as Scenario 1, but with the addition of a 1.5 Mgal/d supply well in the Middendorf aquifer near Moncks Corner, South Carolina; Scenario 4—maximize MPW well capacity by increasing withdrawals from the Middendorf aquifer from 3.9 Mgal/d (in 2015) to 10.16 Mgal/d; Scenario 5—minimize MPW surface-water purchase from the Charleston Water System by adding supply wells and increasing withdrawals from the Middendorf aquifer from 3.9 Mgal/d (in 2015) to 12.16 Mgal/d; and Scenario 6—same as Scenario 1, but with he addition of quarterly model stress periods to simulate seasonal variations in the groundwater withdrawals. Results from the simulations indicated further decline of groundwater levels creating cones of depressions near pumping wells in the Middendorf aquifer in the Mount Pleasant, South Carolina, area between 2015 and 2050 for all six scenarios. Simulation results from Scenario 1 showed an average decline of about 150 feet in the groundwater levels of the MPW production wells. Simulated hydrographs for two area observation wells illustrate the gradual decline in groundwater levels with overall changes in water-level altitudes of –92 and –33 feet, respectively. Simulated groundwater altitudes at a hypothetical observation well located in the MPW well field declined 121 feet between 2015 and 2050. Scenarios 2 and 3 have the same pumping rates as Scenario 1 for the MPW production wells; however, a single hypothetical pumping well was added in the Middendorf aquifer near the town of Moncks Corner, South Carolina. This hypothetical pumping well has a withdrawal rate of 0.5 Mgal/d for Scenario 2 and 1.5 Mgal/d for Scenario 3. A comparison to the 2050 Scenario 1 simulation indicates groundwater altitudes for Scenarios 2 and Scenario 3 are 3 feet and 8 feet lower, respectively, at the MPW production wells. Scenario 4 simulates the maximum pumping capacity of 10.16 Mgal/d for the MPW network of production wells. Simulated 2050 groundwater altitudes for this simulation declined to –359 feet. Simulated hydrographs for two observation wells show groundwater-level declines of 116 and 41 feet, respectively. Simulated differences in groundwater altitudes at a hypothetical observation well located in the MPW well field indicate a water-level decline of 164 feet between 2015 and 2050. Scenario 5 is a modification of Scenario 4 with the addition of two new MPW production wells. For this scenario, the MPW network of production wells were simulated the same as in Scenario 4, but withdrawals from the two new production wells were added in 2020. Simulated 2050 groundwater altitudes for this simulation declined to – 405 feet. Simulated hydrographs for two observation wells show groundwater-level declines of 143 and 51 feet, respectively. Simulated groundwater altitudes at a hypothetical observation well located in the MPW well field declined 199 feet between 2015 and 2050. Scenario 6 is a modification of Scenario 1, in which 140 additional quarterly stress periods were added to simulate MPW seasonal demands. Simulated groundwater altitudes for Scenario 6 declined to –353 feet during 2050. For Scenario 6, simulated hydrographs for two observation wells and the hypothetical observation well show similar groundwater-level declines as seen in Scenario 1, but with seasonal fluctuations of as much as 56 feet in the hypothetical observation well. Water budgets for the model area immediately surrounding Mount Pleasant, South Carolina, were calculated for 2015 and for 2050. The water budget for 2015 is equal for all of the scenarios because it represents the year prior to the hypothetical pumping beginning in 2016. The largest flow component in the 2015 water budget for the Mount Pleasant area is discharge to wells at a rate of 4.17 Mgal/d. Additionally, 0.23 Mgal/d flows laterally out of the Middendorf aquifer in this area of the model due to the regional horizontal hydraulic gradient. Flow into this zone consists predominantly of lateral flow within the Middendorf aquifer at 4.08 Mgal/d. Additionally, 0.02 Mgal/d is released into this zone from aquifer storage. Vertically, 0.06 Mgal/d flows down from the Middendorf confining unit located above the Middendorf aquifer, and 0.25 Mgal/d flows up from the Cape Fear confining unit below. The largest flow component in the 2050 water budget for all six scenarios is discharge to wells in the Mount Pleasant area at rates between 8.89 and 12.47 Mgal/d. Flow into this zone consists mostly of lateral flow between 8.47 and 11.77 Mgal/d within the Middendorf aquifer. Between 0.003 and 0.46 Mgal/d is released into this zone from aquifer storage. Between 0.004 and 0.15 Mgal/d flows laterally out of this zone into adjacent areas of the Middendorf aquifer due to the regional horizontal hydraulic gradient. Finally, between 0.15 and 0.22 Mgal/d flows vertically into this zone from confining units above and below the Middendorf aquifer.

South Carolina↗

Groundwater contributions to excessive algal growth in the East Fork Carson River, Carson Valley, west-central Nevada, 2010 and 2012

Excessive algal growth and low dissolved oxygen concentrations were observed during low streamflow conditions during summer months along a 5,800-foot reach of the East Fork Carson River in Carson Valley, west-central Nevada. Algal growth from nutrient enrichment of a stream reduces aquatic diversity, threatens fish ecology and stream health, and can be a recreational nuisance. In response to concerns that groundwater discharging to the 5,800-foot reach of the East Fork Carson River may be a source of nutrients to the stream, the U.S. Geological Survey, in cooperation with the Carson Water Subconservancy District and the Nevada Division of Environmental Protection, conducted studies during the summers of 2010 and 2012 to gain an improved understanding of the contributions of nutrients to the stream from groundwater, characterize algal conditions and algal effects on water quality, assess potential sources of nitrate in groundwater discharging to the stream, and evaluate nitrate reduction in groundwater from denitrification. A reconnaissance study in the summer of 2010 along the 5,800-foot study reach located a subreach with clear evidence of nutrient-rich groundwater discharging to the stream. At the subreach, nitrate plus nitrite (referred to hereafter as nitrate) concentrations in groundwater discharging to the stream were high (average 2.75 milligrams per liter as nitrogen) along the right bank. The stream at this location had the highest stream nitrate concentrations (average 0.056 milligrams per liter as nitrogen) compared to other locations upstream and downstream of the subreach. As a result, the 2012 study focused on a 405-foot subreach of the East Fork Carson River centered where results from the 2010 study found the highest stream and groundwater concentrations of nitrate, as well as the greatest observed contributions of groundwater discharge to the stream. Groundwater nutrient concentrations were much higher than stream nutrient concentrations during the summer of 2012 during low streamflow conditions at the 405-foot subreach of the East Fork Carson River. Average groundwater nitrate and orthophosphate concentrations along the right bank of the 405‑foot subreach were 9 and 12 times higher, respectively, than in the stream at this subreach. Groundwater discharge rates to the study reach based on different methods varied from 0.09 to 1.2 cubic feet per second per mile. Estimated groundwater discharge rates to the right bank of the study subreach were used to calculate groundwater nutrient load estimates to the subreach right bank, which were found to be low when compared to stream nutrient loads. Elevated algal biomass levels above nuisance thresholds were observed during the summers of 2010 and 2012. The study reach was characterized as mesotrophic-eutrophic during the 2010 study and eutrophic during the 2012 study. The presence of algae caused daily dissolved oxygen and pH fluctuations in the stream, resulting in exceedances of the State of Nevada water-quality standards owing to low dissolved oxygen concentrations and high pH concentrations, although the standards might not have been applicable during 2012 because of extremely low streamflow. The addition of nutrients to the stream from the constant supply in groundwater discharge sustains the growth of algae during low streamflow conditions. In the summer when streamflow is low or very low, nutrient-rich groundwater discharge enters the stream through the sediment-water interface at the streambed. Because the attached algae is thick and stream velocity is low, the nutrient-rich water pools at the sediment-water interface. Higher nutrient concentrations at the streambed create a favorable microenvironment for algae attached to the substrate to consume available nutrients from the groundwater before the groundwater mixes with overlying stream water. The source of nitrate in groundwater in this subreach is anthropogenic because nitrate concentrations are greater than background groundwater nitrate concentrations in Douglas County, high groundwater nitrate concentrations are only found at the right bank of the stream near a housing development, and organic wastewater compounds indicative of human-derived sources were also detected in groundwater wells on the right bank of the stream. Nitrogen and oxygen isotope concentrations of nitrate in shallow groundwater were used to determine the specific source of the nitrate, but the isotopic values indicated denitrification was occurring. Further investigation is needed to determine the specific anthropogenic source of the nitrate in the groundwater because the denitrification present in all samples obscures the original source of nitrogen.

Nevada↗

Arsenic, chromium, uranium, and vanadium in rock, alluvium, and groundwater, Mojave River and Morongo Areas, western Mojave Desert, southern California

Trace elements within groundwater that originate from aquifer materials and pose potential public-health hazards if consumed are known as geogenic contaminants. The geogenic contaminants arsenic, chromium, and vanadium can form negatively charged ions with oxygen known as oxyanions. Uranium complexes with bicarbonate and carbonate to form negatively charged ions having aqueous chemistry similar to oxyanions. The concentrations of arsenic, chromium, uranium, and vanadium in groundwater result from the combined effects of (1) geologic abundance within aquifer materials; (2) the fraction of these elements that have weathered from and sorbed to the surfaces of mineral grains and are potentially available to groundwater; and (3) the aqueous chemistry of dissolved oxyanions in groundwater during different redox conditions and pH, both of which are affected by hydrogeology, including the length of time groundwater has been in contact with aquifer materials. Concentrations of arsenic, chromium, uranium, and vanadium were measured in samples of (1) rock, surficial alluvium, and drill cuttings using portable (handheld) X-ray fluorescence (pXRF); (2) operationally defined fractions extractable from these materials; and (3) water from wells sampled between 2000 and 2018 within the 3,500 square mile Mojave River area and Morongo area of the western Mojave Desert, southern California. Regionally, rock and surficial alluvium in the Mojave River and Morongo areas are high in arsenic, low in chromium and uranium, and near the average bulk continental crust concentration for vanadium. Locally, high chromium concentrations are present in mafic rock within the San Gabriel Mountains; high uranium concentrations are present in felsic rock within the San Bernardino Mountains; and high arsenic, uranium, and vanadium concentrations are present in extrusive (volcanic) felsic rock within uplands surrounding groundwater basins along the Mojave River downstream from Barstow, California. Elemental assemblages identified using principal component analyses (PCA) of pXRF data were used to characterize felsic, mafic, and felsic volcanic source terranes in rock, surficial alluvium, and in geologic material penetrated by selected monitoring wells drilled between 1994 and 2018. Highly felsic alluvium associated with recent deposition from the Mojave River was identified along the 90-mile length of the floodplain aquifer along the river. The thickness of these highly felsic alluvial deposits ranged from 200 feet (ft) near Victorville and near Barstow to a thin veneer about 30 ft thick downstream from Victorville and downstream portions of the floodplain aquifer within the Mojave Valley. Groundwater in the Mojave River and Morongo areas was generally oxic and alkaline (pH≥7.5). Maximum concentrations of arsenic, hexavalent chromium [Cr(VI)], uranium, and vanadium in water from as many as 498 wells sampled between 2000 and 2018 were 360, 140, 1,470, and 690 micrograms per liter (μg/L), respectively. Water from 22 percent of sampled wells exceeded the U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL) for arsenic of 10 μg/L, with arsenic concentrations commonly exceeding the MCL in water from wells east of Barstow, deep wells in the Victorville fan, and in suboxic or reduced groundwater within the floodplain aquifer. Water from about 1 percent of sampled wells had Cr(VI) concentrations greater than the California MCL for total chromium of 50 μg/L, whereas 13 percent of sampled wells had Cr(VI) concentrations greater than the former California MCL of 10 μg/L. Hexavalent chromium concentrations were highest in water from wells in the Sheep Creek alluvial fan, eroded from mafic rock in the San Gabriel Mountains, although Cr(VI) concentrations greater than the former California MCL also were present elsewhere in the study area where mafic materials or older groundwater were present. Water from about 9 percent of sampled wells exceeded the EPA MCL for uranium of 30 μg/L, with concentrations exceeding the MCL commonly associated with irrigation return from agricultural land overlying the floodplain aquifer. Water from about 7 percent of sampled wells had vanadium concentrations greater than the California notification level of 50 μg/L; most of these wells were in the Victorville fan within the Mojave River area. In general, arsenic concentrations were higher in suboxic or reduced water; chromium concentrations were higher in oxic, alkaline (pH≥7.5) water; uranium concentrations were higher in circumneutral to slightly alkaline water (pH≤7.4); and vanadium concentrations were higher in highly alkaline (pH≥8.0) water, independent of redox status. Concentrations within geologic source terranes are not the sole factor controlling the concentrations of geogenic elements in groundwater. Differences in mineral weathering, pH-dependent sorption to surface-exchange sites on mineral grains, and aqueous geochemistry (especially redox status and pH) affect geogenic element concentrations in groundwater. Consequently, the relative abundances of arsenic, Cr(VI), uranium, and vanadium in groundwater differ from their relative abundances in the average bulk continental crust and their regional abundances in rock and surficial alluvium within groundwater basins of the western Mojave Desert. Processes that control the concentrations of arsenic, chromium, uranium, and vanadium in groundwater operate at the mineral-grain and aquifer scale. At the mineral-grain scale, sequential chemical extraction data show arsenic and uranium are more available to groundwater (under specific geochemical conditions) than chromium or vanadium, which largely are unavailable within unweathered mineral grains. Additionally, chromium and vanadium form few aqueous complexes and bind tightly with iron minerals within surface coatings on mineral grains making them less available to groundwater, whereas complexation with other dissolved ions enhances the solubility of uranium and, to a lesser extent, arsenic. Complexation also increases the valence (less negative charge) and increases the size of dissolved oxyanions and uranium complexes with bicarbonate and carbonate making them less readily sorbed to aquifer materials. At the aquifer scale, hydrogeology (including isolation of water in aquifers from surface sources of recharge, older groundwater age, and long contact times between groundwater and aquifer materials) combined with geochemical processes (such as silicate weathering) to produce alkaline groundwater. Desorption from sorption sites on the surfaces of mineral grains with increasing pH increases arsenic, chromium, and vanadium concentrations in water from wells and increases Cr(VI) concentrations as long as water remains oxic. Aqueous geochemistry and concentrations of geogenic contaminants also are affected by anthropogenic activities including (1) discharge of treated municipal wastewater, which may change the redox status of groundwater; (2) return from irrigated agriculture, which may alter the chemistry of groundwater and increase the solubility of trace elements such as uranium; and (3) groundwater pumping and subsequent water-level declines, which may change the source of water yielded by wells. The quality of water imported from northern California and infiltrated from ponds for groundwater recharge may be altered by naturally present trace elements, especially uranium in areas of agricultural land use or chromium within mafic alluvium.

California↗

Groundwater availability of the Mississippi embayment

Groundwater is an important resource for agricultural and municipal uses in the Mississippi embayment. Arkansas ranks first in the Nation for rice and third for cotton production, with both crops dependent on groundwater as a major source of irrigation requirements. Multiple municipalities rely on the groundwater resources to provide water for industrial and public use, which includes the city of Memphis, Tennessee. The demand for the groundwater resource has resulted in groundwater availability issues in the Mississippi embayment including: (1) declining groundwater levels of 50 feet or more in the Mississippi River Valley alluvial aquifer in parts of eastern Arkansas from agricultural pumping, (2) declining groundwater levels of over 360 feet over the last 90 years in the confined middle Claiborne aquifer in southern Arkansas and northern Louisiana from municipal pumping, and (3) litigation between the State of Mississippi and a Memphis water utility over water rights in the middle Claiborne aquifer. To provide information to stakeholders addressing the groundwater-availability issues, the U.S. Geological Survey Groundwater Resources Program supported a detailed assessment of groundwater availability through the Mississippi Embayment Regional Aquifer Study (MERAS). This assessment included (1) an evaluation of how these resources have changed over time through the use of groundwater budgets, (2) development of a numerical modeling tool to assess system responses to stresses from future human uses and climate trends, and (3) application of statistical tools to evaluate the importance of individual observations within a groundwater-monitoring network. An estimated 12 million acre-feet per year (11 billion gallons per day) of groundwater was pumped in 2005 from aquifers in the Mississippi embayment. Irrigation constitutes the largest groundwater use, accounting for approximately 10 million acre-feet per year (9 billion gallons per day) in 2000 from the Mississippi River Valley alluvial aquifer in Arkansas, Louisiana, Mississippi, and Missouri, and to a lesser extent in Illinois, Kentucky, and Tennessee. Predevelopment groundwater flow is represented in the MERAS model as a steady-state stress period, assumed to be prior to 1870. The simulated groundwater-flow budget indicates the largest predevelopment inflow to the system is net recharge to the alluvial aquifer. This inflow is balanced by outflow to gaining streams. Overall, water enters as net recharge to the alluvial aquifer or through outcrop areas of the various hydrogeologic units. Away from the outcrop areas, groundwater flow in the deeper formations is primarily upward into overlying units, ultimately discharging to streams through the alluvial aquifer. Total net recharge and discharge (sum of inflows or outflows) for the model ranged from about 0.66 million acre-feet per year during predevelopment to 20.16 million acre-feet per year by the end of the simulation (final simulated irrigation period in summer of 2006). This change in the model budget reflects increases in withdrawals compared to predevelopment conditions. Cumulative storage within aquifers simulated in the MERAS model indicates overall depletion of 140 million acre-feet (equivalent to 2.8 feet of water covering the entire study area). Postdevelopment inflow to the system is still through net recharge to the alluvial aquifer and the outcrop areas of the several hydrogeologic units, however, the flow between each unit is no longer upward to the alluvial aquifer. Groundwater flow during the summer of 2006 was primarily downward to offset demand from pumping. Early in the model simulation (1870-1920s), the primary components of the water budget were simulated as outflow from stream leakage and inflow from net recharge. As pumpage increased through time, water that would otherwise flow to streams reversed, and net stream leakage became an inflow to the system. The largest reversals began in the mid-1980s, but indications of the reversal began in the early 1960s with a trend in loss of streamflow leakage coupled with the first consistent inflow from storage. While groundwater pumped out of the alluvial aquifer was derived primarily from storage, pumpage out of the middle Claiborne aquifer was derived primarily from other aquifers (up to 15 percent from the alluvial aquifer), followed by flow from storage and net recharge. The potential consequences of climate change have been identified as a major concern facing the sustainability of the Nation's groundwater resources. To address this concern, two climate simulations were developed through the use of the MERAS model by extending the simulation period by 30 years to the year 2038 using extrapolated precipitation based on frequency analysis of historic climate cycles. There is little difference between the dry and wet scenarios in terms of percent water-level change. Both scenarios resulted in 14.6 to 13.9 percent of the area containing more than 100 feet of decline, 14.5 to 13.8 percent containing between 75 and 100 feet of decline, and 15.8 to 15.7 percent containing 51 to 75 feet of decline in the alluvial aquifer. The middle Claiborne aquifer water-level changes also were similar between the two scenarios. These scenarios indicate that even with a 25-percent increase in precipitation from that of the dry scenario, there is little difference in the resultant water levels. This is in large part because of the magnitude of differences between changes in net recharge and changes in pumping. When compared to the volume of water pumped out of the system, the effect of this change in net recharge is negligible. The groundwater-level monitoring network used to construct the 2007 middle Claiborne aquifer potentiometric surface was used as an example case to demonstrate statistical technique and to evaluate the importance of individual groundwater-level observations. To calculate the importance of each water-level observation to a prediction, predictions were specified as water-level altitudes near the end of the dry scenario simulation. These predictions were located near the center of cones of depression. Many of the observations that have a high importance are in close proximity to stressed areas of the aquifer.

Arkansas, Kentucky, Louisiana, Mississippi, Missou↗

Status and understanding of groundwater quality in the Tahoe-Martis, Central Sierra, and Southern Sierra study units, 2006-2007--California GAMA Priority Basin Project

Groundwater quality in the Tahoe-Martis, Central Sierra, and Southern Sierra study units was investigated as part of the Priority Basin Project of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program. The three study units are located in the Sierra Nevada region of California in parts of Nevada, Placer, El Dorado, Madera, Tulare, and Kern Counties. The GAMA Priority Basin Project is being conducted by the California State Water Resources Control Board, in collaboration with the U.S. Geological Survey (USGS) and the Lawrence Livermore National Laboratory. The project was designed to provide statistically robust assessments of untreated groundwater quality within the primary aquifer systems used for drinking water. The primary aquifer systems (hereinafter, primary aquifers) for each study unit are defined by the depth of the screened or open intervals of the wells listed in the California Department of Public Health (CDPH) database of wells used for municipal and community drinking-water supply. The quality of groundwater in shallower or deeper water-bearing zones may differ from that in the primary aquifers; shallower groundwater may be more vulnerable to contamination from the surface. The assessments for the Tahoe-Martis, Central Sierra, and Southern Sierra study units were based on water-quality and ancillary data collected by the USGS from 132 wells in the three study units during 2006 and 2007 and water-quality data reported in the CDPH database. Two types of assessments were made: (1) status, assessment of the current quality of the groundwater resource, and (2) understanding, identification of the natural and human factors affecting groundwater quality. The assessments characterize untreated groundwater quality, not the quality of treated drinking water delivered to consumers by water purveyors. Relative-concentrations (sample concentrations divided by benchmark concentrations) were used for evaluating groundwater quality for those constituents that have Federal or California regulatory or non-regulatory benchmarks for drinking-water quality. A relative-concentration (RC) greater than (>) 1.0 indicates a concentration above a benchmark. RCs for organic constituents (volatile organic compounds and pesticides) and special-interest constituents were classified as "high" (RC > 1.0), "moderate" (1.0 &ge; RC > 0.1), or "low" (RC &le; 0.1). For inorganic constituents (major ions, trace elements, nutrients, and radioactive constituents), the boundary between low and moderate RCs was set at 0.5. A new metric, aquifer-scale proportion, was used in the status assessment as the primary metric for evaluating regional-scale groundwater quality. High aquifer-scale proportion is defined as the percentage of the area of the primary aquifers with RC > 1.0 for a particular constituent or class of constituents; moderate and low aquifer-scale proportions are defined as the percentages of the area of the primary aquifer with moderate and low RCs, respectively. Percentages are based on an areal rather than a volumetric basis. Two statistical approaches&mdash;grid-based, which used one value per grid cell, and spatially weighted, which used multiple values per grid cell&mdash;were used to calculate aquifer-scale proportions for individual constituents and classes of constituents. The spatially weighted estimates of high aquifer-scale proportions were within the 90-percent (%) confidence intervals of the grid-based estimates in all cases. The status assessment showed that inorganic constituents had greater high and moderate aquifer-scale proportions than did organic constituents in all three study units. In the Tahoe-Martis study unit, RCs for inorganic constituents with health-based benchmarks (primarily arsenic) were high in 20% of the primary aquifer, moderate in 13%, and low in 67%. In the Central Sierra study unit, aquifer-scale proportions for inorganic constituents with health-based benchmarks (primarily arsenic, uranium, fluoride, and molybdenum) were 41% high, 36% moderate, and 23% low. In the Southern Sierra study unit, 32, 34, and 34% of the primary aquifer had high, moderate, and low RCs of inorganic constituents with health-based benchmarks (primarily arsenic, uranium, fluoride, boron, and nitrate). The high aquifer-scale proportions for inorganic constituents with non-health-based benchmarks were 14, 34, and 24% for the Tahoe-Martis, Central Sierra, and Southern Sierra study units, respectively, and the primary constituent was manganese for all three study units. Organic constituents with health-based benchmarks were not present at high RCs in the primary aquifers of the Central Sierra and Southern Sierra study units, and were present at high RCs in only 1% of the Tahoe-Martis study unit. Moderate aquifer-scale proportions for organic constituents were < 5% in all three study units. Of the 173 organic constituents analyzed, 22 were detected, and of those 22, 17 have health-based benchmarks. Organic constituents were detected in 20, 27, and 40% of the primary aquifers in the Tahoe-Martis, Central Sierra, and Southern Sierra study units, respectively. Four organic constituents had study-unit detection frequencies of > 10%: the trihalomethane chloroform in the Tahoe-Martis study unit; chloroform and the herbicide simazine in the Central Sierra study unit; and chloroform, simazine, the herbicide atrazine, and the solvent perchloroethene in the Southern Sierra study unit. The second component of this study, the understanding assessment, identified the natural and human factors that may have affected groundwater quality in the three study units by evaluating statistical correlations between water-quality constituents and potential explanatory factors. The potential explanatory factors evaluated were land use, septic tank density, climate, relative position in the regional flow system, aquifer lithology, geographic location, well depth and depth to the top of the screened or open interval in the well, groundwater age distribution, pH, and dissolved oxygen concentration. Results of the statistical evaluations were used to explain the occurrence and distribution of constituents in the study units. Aquifer lithology (granitic, metamorphic, sedimentary, or volcanic rocks), groundwater age distribution [modern (recharged since 1952), pre-modern (recharged before 1952), or mixed (containing both modern and pre-modern recharge)], geographic location, pH, and dissolved oxygen were the most significant factors explaining the occurrence patterns of most inorganic constituents. High and moderate RCs of arsenic were associated with pre-modern and mixed-age groundwater and two distinct sets of geochemical conditions: (1) oxic, high-pH conditions, particularly in volcanic rocks, and (2) low-oxygen to anoxic conditions and low- to neutral-pH conditions, particularly in granitic rocks. In granitic and metamorphic rocks, high and moderate RCs of uranium were associated with pre-modern and mixed-age groundwater, low-oxygen to anoxic conditions, and location within parts of the Central Sierra and Southern Sierra study units known to have rocks with anomalously high uranium content compared to other parts of the Sierra Nevada. High and moderate RCs of uranium in sedimentary rocks were associated with pre-modern-age groundwater, oxic and high-pH conditions, and location in the Tahoe Valley South subbasin within the Tahoe-Martis study unit. Land use within 500 meters of the well and groundwater age were the most significant factors explaining occurrence patterns of organic constituents. Herbicide detections were most strongly associated with modern- and mixed-age groundwater from wells with agricultural land use. Trihalomethane detections were most strongly associated with modern- and mixed-age groundwater from wells with > 10% urban land use and (or) septic tank density > 7 tanks per square kilometer. Solvent detections were not significantly related to groundwater age. Eighty-three percent of the wells with modern- or mixed-age groundwater, and 86% of wells with detections of herbicides and (or) THMs had depths to the top of the screened or open interval of < 170 feet. These observations suggest that modern groundwater has infiltrated to depths of approximately 170 feet below land surface. Land use and occurrence of herbicides and solvents were the most significant factors explaining the occurrence of nitrate. Wells with > 5% agricultural land use and detection of a herbicide or solvent had the highest nitrate concentrations. Comparison between observed and predicted detection frequencies of perchlorate suggests that the perchlorate detected at concentrations < 1 microgram per liter likely reflects the distribution of perchlorate under natural conditions, and that the perchlorate detected at higher concentrations may reflect redistribution of originally natural perchlorate salts by irrigation in the agricultural areas of the Southern Sierra study unit.

California↗

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&ndash;2010 is part of the U.S. Geological Survey&rsquo;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 &ldquo;irrigation-recharge layer&rdquo; 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&mdash;Maximum Contaminant Levels (MCLs), and Health-Based Screening Levels (HBSLs)&mdash;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&rsquo; -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, ↗

Simulation of groundwater flow, effects of artificial recharge, and storage volume changes in the Equus Beds aquifer near the city of Wichita, Kansas well field, 1935–2008

The Equus Beds aquifer is a primary water-supply source for Wichita, Kansas and the surrounding area because of shallow depth to water, large saturated thickness, and generally good water quality. Substantial water-level declines in the Equus Beds aquifer have resulted from pumping groundwater for agricultural and municipal needs, as well as periodic drought conditions. In March 2006, the city of Wichita began construction of the Equus Beds Aquifer Storage and Recovery project to store and later recover groundwater, and to form a hydraulic barrier to the known chloride-brine plume near Burrton, Kansas. In October 2009, the U.S. Geological Survey, in cooperation with the city of Wichita, began a study to determine groundwater flow in the area of the Wichita well field, and chloride transport from the Arkansas River and Burrton oilfield to the Wichita well field. Groundwater flow was simulated for the Equus Beds aquifer using the three-dimensional finite-difference groundwater-flow model MODFLOW-2000. The model simulates steady-state and transient conditions. The groundwater-flow model was calibrated by adjusting model input data and model geometry until model results matched field observations within an acceptable level of accuracy. The root mean square (RMS) error for water-level observations for the steady-state calibration simulation is 9.82 feet. The ratio of the RMS error to the total head loss in the model area is 0.049 and the mean error for water-level observations is 3.86 feet. The difference between flow into the model and flow out of the model across all model boundaries is -0.08 percent of total flow for the steady-state calibration. The RMS error for water-level observations for the transient calibration simulation is 2.48 feet, the ratio of the RMS error to the total head loss in the model area is 0.0124, and the mean error for water-level observations is 0.03 feet. The RMS error calculated for observed and simulated base flow gains or losses for the Arkansas River for the transient simulation is 7,916,564 cubic feet per day (91.6 cubic feet per second) and the RMS error divided by (/) the total range in streamflow (7,916,564/37,461,669 cubic feet per day) is 22 percent. The RMS error calculated for observed and simulated streamflow gains or losses for the Little Arkansas River for the transient simulation is 5,610,089 cubic feet per day(64.9 cubic feet per second) and the RMS error divided by the total range in streamflow (5,612,918/41,791,091 cubic feet per day) is 13 percent. The mean error between observed and simulated base flow gains or losses was 29,999 cubic feet per day (0.34 cubic feet per second) for the Arkansas River and -1,369,250 cubic feet per day (-15.8 cubic feet per second) for the Little Arkansas River. Cumulative streamflow gain and loss observations are similar to the cumulative simulated equivalents. Average percent mass balance difference for individual stress periods ranged from -0.46 to 0.51 percent. The cumulative mass balance for the transient calibration was 0.01 percent. Composite scaled sensitivities indicate the simulations are most sensitive to parameters with a large areal distribution. For the steady-state calibration, these parameters include recharge, hydraulic conductivity, and vertical conductance. For the transient simulation, these parameters include evapotranspiration, recharge, and hydraulic conductivity. The ability of the calibrated model to account for the additional groundwater recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project was assessed by using the U.S. Geological Survey subregional water budget program ZONEBUDGET and comparing those results to metered recharge for 2007 and 2008 and previous estimates of artificial recharge. The change in storage between simulations is the volume of water that estimates the recharge credit for the aquifer storage and recovery system. The estimated increase in storage of 1,607 acre-ft in the basin storage area compared to metered recharge of 1,796 acre-ft indicates some loss of metered recharge. Increased storage outside of the basin storage area of 183 acre-ft accounts for all but 6 acre-ft or 0.33 percent of the total. Previously estimated recharge credits for 2007 and 2008 are 1,018 and 600 acre-ft, respectively, and a total estimated recharge credit of 1,618 acre-ft. Storage changes calculated for this study are 4.42 percent less for 2007 and 5.67 percent more for 2008 than previous estimates. Total storage change for 2007 and 2008 is 0.68 percent less than previous estimates. The small difference between the increase in storage from artificial recharge estimated with the groundwater-flow model and metered recharge indicates the groundwater model correctly accounts for the additional water recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project. Small percent differences between inflows and outflows for all stress periods and all index cells in the basin storage area, improved calibration compared to the previous model, and a reasonable match between simulated and measured long-term base flow indicates the groundwater model accurately simulates groundwater flow in the study area. The change in groundwater level through recent years compared to the August 1940 groundwater level map has been documented and used to assess the change of storage volume of the Equus Beds aquifer in and near the Wichita well field for three different areas. Two methods were used to estimate changes in storage from simulation results using simulated change in groundwater levels in layer 1 between stress periods, and using ZONEBUDGET to calculate the change in storage in the same way the effects of artificial recharge were estimated within the basin storage area. The three methods indicate similar trends although the magnitude of storage changes differ. Information about the change in storage in response to hydrologic stresses is important for managing groundwater resources in the study area. The comparison between the three methods indicates similar storage change trends are estimated and each could be used to determine relative increases or decreases in storage. Use of groundwater level changes that do not include storage changes that occur in confined or semi-confined parts of the aquifer will slightly underestimate storage changes; however, use of specific yield and groundwater level changes to estimate storage change in confined or semi-confined parts of the aquifer will overestimate storage changes. Using only changes in shallow groundwater levels would provide more accurate storage change estimates for the measured groundwater levels method. The value used for specific yield is also an important consideration when estimating storage. For the Equus Beds aquifer the reported specific yield ranges between 0.08 and 0.35 and the storage coefficient (for confined conditions) ranges between 0.0004 and 0.16. Considering the importance of the value of specific yield and storage coefficient to estimates of storage change over time, and the wide range and substantial overlap for the reported values for specific yield and storage coefficient in the study area, further information on the distribution of specific yield and storage coefficient within the Equus Beds aquifer in the study area would greatly enhance the accuracy of estimated storage changes using both simulated groundwater level, simulated groundwater budget, or measured groundwater level methods.

Kansas↗

Groundwater and surface-water interaction within the upper Smith River Watershed, Montana 2006-2010

The 125-mile long Smith River, a tributary of the Missouri River, is highly valued as an agricultural resource and for its many recreational uses. During a drought starting in about 1999, streamflow was insufficient to meet all of the irrigation demands, much less maintain streamflow needed for boating and viable fish habitat. In 2006, the U.S. Geological Survey, in cooperation with the Meagher County Conservation District, initiated a multi-year hydrologic investigation of the Smith River watershed. This investigation was designed to increase understanding of the water resources of the upper Smith River watershed and develop a detailed description of groundwater and surface-water interactions. A combination of methods, including miscellaneous and continuous groundwater-level, stream-stage, water-temperature, and streamflow monitoring was used to assess the hydrologic system and the spatial and temporal variability of groundwater and surface-water interactions. Collectively, data are in agreement and show: (1) the hydraulic connectedness of groundwater and surface water, (2) the presence of both losing and gaining stream reaches, (3) dynamic changes in direction and magnitude of water flow between the stream and groundwater with time, (4) the effects of local flood irrigation on groundwater levels and gradients in the watershed, and (5) evidence and timing of irrigation return flows to area streams. Groundwater flow within the alluvium and older (Tertiary) basin-fill sediments generally followed land-surface topography from the uplands to the axis of alluvial valleys of the Smith River and its tributaries. Groundwater levels were typically highest in the monitoring wells located within and adjacent to streams in late spring or early summer, likely affected by recharge from snowmelt and local precipitation, leakage from losing streams and canals, and recharge from local flood irrigation. The effects of flood irrigation resulted in increased hydraulic gradients (increased groundwater levels relative to stream stage) or even reversed gradient direction at several monitoring sites coincident with the onset of nearby flood irrigation. Groundwater-level declines in mid-summer were due to groundwater withdrawals and reduced recharge from decreased precipitation, increased evapotranspiration, and reduced leakage in some area streams during periods of low flow. Groundwater levels typically rebounded in late summer, a result of decreased evapotranspiration, decreased groundwater use for irrigation, increased flow in losing streams, and the onset of late-season flood irrigation at some sites. The effect of groundwater and surface-water interactions is most apparent along the North and South Forks of the Smith River where the magnitude of streamflow losses and gains can be greater than the magnitude of flow within the stream. Net gains consistently occurred over the lower 15 miles of the South Fork Smith River. A monitoring site near the mouth of the South Fork Smith River gained (flow from the groundwater to the stream) during all seasons, with head gradients towards the stream. Two upstream sites on the South Fork Smith River exhibited variable conditions that ranged from gaining during the spring, losing (flowing from the stream to the groundwater) during most of the summer as groundwater levels declined, and then approached or returned to gaining conditions in late summer. Parts of the South Fork Smith River became dry during periods of losing conditions, thus classifying this tributary as intermittent. The North Fork Smith River is highly managed at times through reservoir releases. The North Fork Smith River was perennial throughout the study period although irrigation diversions removed a large percentage of streamflow at times and losing conditions persisted along a lower reach. The lowermost reach of the North Fork Smith River near its mouth transitioned from a losing reach to a gaining reach throughout the study period. Groundwater and surface-water interactions occur downstream from the confluence of the North and South Fork Smith Rivers, but are less discernible compared to the overall magnitude of the main-stem streamflow. The Smith River was perennial throughout the study. Monitoring sites along the Smith River generally displayed small head gradients between the stream and the groundwater, while one site consistently showed strongly gaining conditions. Synoptic streamflow measurements during periods of limited irrigation diversion in 2007 and 2008 consistently showed gains over the upper 41.4 river miles of the main stem Smith River where net gains ranged from 13.0 to 28.9 cubic feet per second. Continuous streamflow data indicated net groundwater discharge and small-scale tributary inflow contributions of around 25 cubic feet per second along the upper 10-mile reach of the Smith River for most of the 2010 record. A period of intense irrigation withdrawal during the last two weeks in May was followed by a period (early June 2010 to mid-July 2010) with the largest net increase (an average of 71.1 cubic feet per second) in streamflow along this reach of the Smith River. This observation is likely due to increased groundwater discharge to the Smith River resulting from irrigation return flow. By late July, the apparent effects of return flows receded, and the net increase in streamflow returned to about 25 cubic feet per second. Two-dimensional heat and solute transport VS2DH models representing selected stream cross sections were used to constrain the hydraulic properties of the Quaternary alluvium and estimate temporal water-flux values through model boundaries. Hydraulic conductivity of the Quaternary alluvium of the modeled sections ranged from 3x10-6 to 4x10-5 feet per second. The models showed reasonable approximations of the streambed and shallow aquifer environment, and the dynamic changes in water flux between the stream and the groundwater through different model boundaries.

Montana↗