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 757 records · Page 42Linked to original sources

Groundwater, surface-water, and water-chemistry data, Black Mesa area, northeastern Arizona—2012–2013

The Navajo (N) aquifer is an extensive aquifer and the primary source of groundwater in the 5,400-square-mile Black Mesa area in northeastern Arizona. Availability of water is an important issue in northeastern Arizona because of continued water requirements for industrial and municipal use by a growing population and because of low precipitation in the arid climate of the Black Mesa area. Precipitation in the area typically is between 6 and 14 inches per year. The U.S. Geological Survey water-monitoring program in the Black Mesa area began in 1971 and provides information about the long-term effects of groundwater withdrawals from the N aquifer for industrial and municipal uses. This report presents results of data collected as part of the monitoring program in the Black Mesa area from January 2012 to September 2013. The monitoring program includes measurements of (1) groundwater withdrawals, (2) groundwater levels, (3) spring discharge, (4) surface-water discharge, and (5) groundwater chemistry. In calendar year 2012, total groundwater withdrawals were 4,010 acre-ft, industrial withdrawals were 1,370 acre-ft, and municipal withdrawals were 2,640 acre-ft. Total withdrawals during 2012 were about 45 percent less than total withdrawals in 2005 because of Peabody Western Coal Company’s discontinued use of water to transport coal in a coal slurry pipeline. From 2011 to 2012 total withdrawals decreased by 10 percent; industrial withdrawals decreased by approximately 1 percent, and total municipal withdrawals decreased by 15 percent. From 2012 to 2013, annually measured water levels in the Black Mesa area declined in 6 of 16 wells that were available for comparison in the unconfined areas of the N aquifer, and the median change was 0.8 feet. Water levels declined in 5 of 16 wells measured in the confined area of the aquifer. The median change for the confined area of the aquifer was 0.3 feet. From the prestress period (prior to 1965) to 2013, the median water-level change for 34 wells in both the confined and unconfined areas was -13.5 feet; the median water-level changes were -0.8 feet for 16 wells measured in the unconfined areas and -51.0 feet for 16 wells measured in the confined area. Spring flow was measured at four springs in 2013; Burro, Unnamed Spring near Dennehotso, Moenkopi School, and Pasture Canyon Springs. Flow fluctuated during the period of record for Burro and Unnamed Springs near Dennehotso, but a decreasing trend was apparent at Moenkopi School Spring and Pasture Canyon Spring. Discharge at Burro Spring has remained relatively constant since it was first measured in the 1980s and discharge at Unnamed Spring near Dennehotso has fluctuated for the period of record at each spring. Trend analysis for discharge at Moenkopi School and Pasture Canyon Springs showed a decreasing trend. Continuous records of surface-water discharge in the Black Mesa area were collected from streamflow-gaging stations at the following sites: Moenkopi Wash at Moenkopi 09401260 (1976 to 2013), Dinnebito Wash near Sand Springs 09401110 (1993 to 2013), Polacca Wash near Second Mesa 09400568 (1994 to 2013), and Pasture Canyon Springs 09401265 (2004 to 2013). Median winter flows (November through February) from these sites for each water year were used as an index of the amount of groundwater discharge. For the period of record of each streamflow-gaging station, the median winter flows have generally remained constant, which suggests no change in groundwater discharge. In 2013, water samples collected from 12 wells and 4 springs in the Black Mesa area were analyzed for selected chemical constituents, and the results were compared with previous analyses. Concentrations of dissolved solids, chloride, and sulfate have varied at all 12 wells for the period of record, but neither increasing nor decreasing trends over time were found. Dissolved solids, chloride, and sulfate concentrations increased at Moenkopi School Spring during the more than 13 years of record at that site. Concentrations of dissolved solids, chloride, and sulfate at Pasture Canyon Spring have not varied significantly since the early 1980s. Concentrations of dissolved solids, chloride, and sulfate at Burro Spring and Unnamed Spring near Dennehotso have varied for the period of record with no increasing or decreasing trend in the data.

Arizona↗

Preliminary investigation of groundwater flow and trichloroethene transport in the Surficial Aquifer System, Naval Industrial Reserve Ordnance Plant, Fridley, Minnesota

Industrial practices at the Naval Industrial Reserve Ordnance Plant, in Fridley, Minnesota, caused soil and groundwater contamination. Some volatile organic compounds from the plant might have discharged to the Mississippi River, forced by the natural hydraulic gradient in the surficial aquifer system. The U.S. Environmental Protection Agency included the Naval Industrial Reserve Ordnance Plant on the Superfund National Priorities List in 1989. This report describes a preliminary characterization of trichloroethene transport in the surficial and Cambrian-Ordovician aquifer systems at the Naval Industrial Reserve Ordnance Plant. The characterization first involved simulation of 2001 conditions using a model, followed by an application of this 2001 simulator to 2011 conditions. The U.S. Geological Survey, in cooperation with the U.S. Department of the Navy, used a steady-state, uniform-density groundwater flow model to simulate measured potentiometric heads in aquifer systems on August 20, 2001, and a single-phase, conservative, non-reactive, miscible transport model to simulate trichloroethene concentrations in aquifer systems measured in 2001. The U.S. Department of the Navy furnished trichloroethene source areas and trichloroethene source area concentrations to the U.S. Geological Survey for this model simulation. Furnished delineations were postulated and informed by data collected from 1995 to 2011. The groundwater flow simulation of August 20, 2001, was superior to the trichloroethene transport simulation at replicating measurements; simulated potentiometric heads matched 90 percent of measured potentiometric heads on August 20, within 2 feet at selected locations whereas simulated trichloroethene concentration contours of 3, 10, 100, 1000, and 10,000 micrograms per liter (µg/L) correctly bounded 52 percent of measured concentrations in 2001 at selected locations. The degree to which the simulated trichloroethene plume does not match trichloroethene measurements in the surficial aquifer system during the 2001 simulation may suggest that furnished trichloroethene source areas and trichloroethene source area concentrations did not accurately represent all trichloroethene sources in the hydrogeologic system. During the model simulation of 2001, trichloroethene discharged to the Mississippi River. A simulated 900-foot-long zone of benthic trichloroethene discharge flux existed in the shallow flow zone, across which simulated trichloroethene discharged from the surficial aquifer system to the Mississippi River at simulated trichloroethene concentrations that ranged from 3 µg/L to more than 100 µg/L. The Mississippi River was not sampled for volatile organic compounds in Fridley, Minn., from 1999 to 2016 (the publication of this report). Trichloroethene concentrations were measured in wells close to the Mississippi River in the surficial aquifer system on the downgradient side of the Naval Industrial Reserve Ordnance Plant groundwater flow field; for example, at well MS–43 in the shallow flow zone of the surficial aquifer system 280 feet east of the Mississippi River between December 1999 and August 2012, trichloroethene concentrations ranged from 130 to 220 µg/L. The 220-µg/L maximum concentration was reached in March 2003 and October 2006. The August 2012 concentration was 140 µg/L. The August 20, 2001, groundwater flow model simulator and the 2001 trichloroethene transport simulator were applied to a groundwater extraction and treatment system that existed in 2011. Furnished trichloroethene source areas and concentrations in the 2001 simulator were replaced with different, furnished, hypothetical source areas and concentrations. Forcing in 2001 was replaced with forcing in 2011. No trichloroethene concentrations greater than 3 µg/L were simulated as discharging to the Mississippi River during applications of the 2001 simulator to the 2011 groundwater extraction and treatment system. These applications were not intended to represent historical conditions. Differences between furnished and actual trichloroethene sources may explain differences between measurements and simulation results for the 2001 trichloroethene transport simulator. Causes of differences between furnished and actual trichloroethene sources may cause differences between hypothetical application results and the performance of the actual U.S. Department of the Navy groundwater extraction and treatment system at the Naval Industrial Reserve Ordnance Plant. Other limitations may also cause differences between application results and performance.

Minnesota↗

Groundwater, surface-water, and water-chemistry data, Black Mesa area, northeastern Arizona—2013–2015

The Navajo (N) aquifer is an extensive aquifer and the primary source of groundwater in the 5,400-square-mile Black Mesa area in northeastern Arizona. Availability of water is an important issue in northeastern Arizona because of continued water requirements for industrial and municipal use by a growing population and because of low precipitation in the arid climate of the Black Mesa area. Precipitation in the area typically is between 6 and 16 inches per year. The U.S. Geological Survey water-monitoring program in the Black Mesa area began in 1971 and provides information about the long-term effects of groundwater withdrawals from the N aquifer for industrial and municipal uses. This report presents results of data collected as part of the monitoring program in the Black Mesa area from January 2013 to December 2015. The monitoring program includes measurements of (1) groundwater withdrawals (pumping), (2) groundwater levels, (3) spring discharge, (4) surface-water discharge, and (5) groundwater chemistry. In 2013, total groundwater withdrawals were 3,980 acre-feet (ft), in 2014 total withdrawals were 4,170 acre-ft, and in 2015 total withdrawals were 3,970 acre-ft. From 2013 to 2015 total withdrawals varied by less than 5 percent. From 2014 to 2015, annually measured water levels in the Black Mesa area declined in 9 of 15 wells that were available for comparison in the unconfined areas of the N aquifer, and the median change was -0.1 feet. Water levels declined in 3 of 16 wells measured in the confined area of the aquifer. The median change for the confined area of the aquifer was 0.6 feet. From the prestress period (prior to 1965) to 2015, the median water-level change for 34 wells in both the confined and unconfined areas was -13.2 feet; the median water-level changes were -1.7 feet for 16 wells measured in the unconfined areas and -42.3 feet for 18 wells measured in the confined area. Spring flow was measured at four springs in 2014. Flow fluctuated during the period of record for Burro Spring and Unnamed Spring near Dennehotso, but a decreasing trend was statistically significant (p<0.05) at Moenkopi School Spring and Pasture Canyon Spring. Discharge at Burro Spring has remained relatively constant since it was first measured in the 1980s and discharge at Unnamed Spring near Dennehotso has fluctuated for the period of record. Trend analysis for discharge at Moenkopi and Pasture Canyon Springs yielded a slope significantly different (p<0.05) from zero. Continuous records of surface-water discharge in the Black Mesa area were collected from streamflow-gaging stations at the following sites: Moenkopi Wash at Moenkopi 09401260 (1976 to 2015), Dinnebito Wash near Sand Springs 09401110 (1993 to 2015), Polacca Wash near Second Mesa 09400568 (1994 to 2015), and Pasture Canyon Springs 09401265 (2004 to 2015). Median winter flows (November through February) of each water year were used as an index of the amount of groundwater discharge at the above-named sites. For the period of record of each streamflow-gaging station, the median winter flows have generally remained constant, which suggests no change in groundwater discharge. In 2014, water samples collected from four springs in the Black Mesa area were analyzed for selected chemical constituents, and the results were compared with previous analyses. Dissolved solids, chloride, and sulfate concentrations increased at Moenkopi School Spring during the more than 25 years of record at that site. Concentrations of dissolved solids, chloride, and sulfate at Pasture Canyon Spring have not varied significantly (p>0.05) since the early 1980s, and there is no increasing or decreasing trend in those data. Concentrations of dissolved solids, chloride, and sulfate at Burro Spring and Unnamed Spring near Dennehotso have varied for the period of record, but there is no increasing or decreasing statistical trend in the data.

Arizona↗

Interactive tool to estimate groundwater elevations in central and eastern North Dakota

This report describes an interactive tool (NDakGWtool) in which a statistical model is developed using locally weighted regression to estimate monthly mean groundwater elevations for a specified latitude and longitude, referred to as the “user-specified location.” For each user-specified location, seven models are developed for each month from April through October. Localized, high spatial-resolution maps of estimated monthly mean groundwater surface elevations are produced from the models. The tool was evaluated for glacial drift aquifers of the 32-county study area in central and eastern North Dakota. Although groundwater elevations from 1960 to 2017 were available to develop the tool, groundwater elevations from 1995 to 2015 were used for model testing and development of the model domain. There are 413 grid cells of 0.1-degree latitude by 0.1-degree longitude size in the model domain, and the tool produces maps of estimated monthly mean groundwater surface elevations for the cell containing the user-specified location. Additionally, the NDakGWtool produces maps of estimated groundwater depth below land surface and ArcGIS files of estimated groundwater surface elevations and groundwater depth below land surface. The tool is composed of four main components: data input, statistical model, output, and user-interactive process.

North Dakota↗

Effect of groundwater withdrawals, river stage, and precipitation on water-table elevations in the Iowa River alluvial aquifer near Tama, Iowa, 2017–20

The Sac and Fox Tribe of the Mississippi in Iowa is the only federally recognized Tribe in the State of Iowa and is commonly known as the Meskwaki Nation. The Tribe owns more than 8,100 acres, referred to as the “Meskwaki Settlement.” The Meskwaki Settlement uses a well field that withdraws water from the Iowa River alluvial aquifer (IRAA) to supply drinking water to members of the Tribe. Increased severity and timing of flooding and drought conditions, coupled with water-quality concerns in the Iowa River, have prompted the Meskwaki Nation to start identifying tools to provide a better understanding of how extreme climate events (changes in streamflow, flood frequency, and magnitude and persistence of drought conditions), increasing water-supply demands, and groundwater storage depletion will affect water availability in the IRAA. From June 2017 through September 2020, the U.S. Geological Survey, in cooperation with the Meskwaki Nation, collected continuous and discrete groundwater level data from 11 wells in a U.S. Geological Survey monitoring-well network. Groundwater level data collected at these wells were assessed with daily precipitation data and compared to changes in stream level elevations and daily groundwater withdrawals to determine how these changes affect groundwater-table elevations. Results from this study could be used to guide the development of a conceptual model for groundwater flow and a groundwater flow model for the IRAA to quantify and forecast the effect of groundwater withdrawals, Iowa River streamflow, and local precipitation on the water table in the IRAA.

Iowa↗

Groundwater, surface-water, and water-chemistry data, Black Mesa area, northeastern Arizona—2016–2018

The Navajo (N) aquifer is the primary source of groundwater in the 5,400-square-mile Black Mesa area in northeastern Arizona. Availability of water is an important issue in the Black Mesa area because of continued water requirements for industrial and municipal use by a growing population and because of its arid climate. Precipitation in the area typically ranges from less than 6 to more than 16 inches per year depending on location. The U.S. Geological Survey water-monitoring program in the Black Mesa area began in 1971 and provides information about the long-term effects of groundwater withdrawals from the N aquifer for industrial and municipal uses. This report presents results of data collected as part of the monitoring program in the Black Mesa area from November 2016 to December 2018. The monitoring program includes measurements of (1) groundwater withdrawals (pumping), (2) groundwater levels, (3) spring discharge, (4) surface-water discharge, and (5) groundwater and surface-water chemistry. In calendar year 2017, total groundwater withdrawals were 3,710 acre-feet (acre-ft), industrial withdrawals were 1,110 acre-ft, and municipal withdrawals were 2,600 acre-ft. In calendar year 2018, total groundwater withdrawals were 3,670 acre-ft, industrial withdrawals were 1,170 acre-ft, and municipal withdrawals were 2,500 acre-ft. Total withdrawals during 2017 and 2018 were about 49 percent less than total withdrawals in 2005 because of Peabody Western Coal Company’s discontinued use of water to transport coal in a coal slurry pipeline. From the prestress period (prior to 1965) to 2018, measured water levels available for comparison in wells completed in the unconfined areas of the N aquifer within the Black Mesa area declined in 8 of 14 wells, the changes ranged from +12.1 feet to −39.4 feet, and the median change was -0.6 feet. Water levels also declined in 15 of 18 wells measured in the confined area of the aquifer. The median change for the confined area of the aquifer was −40.2 feet (ft), with changes ranging from +14.2 ft to −189.0 ft. From the prestress period to 2018, the median water-level change for all 32 wells in both the confined and unconfined areas was −9.4 ft. Spring flow was measured at four springs in 2017 and 2018. Flow fluctuated during the period of record for Burro Spring and Pasture Canyon Spring, but a decreasing trend was statistically significant (p<0.05) at Moenkopi School Spring and Unnamed Spring near Dennehotso. Discharge at Burro Spring has remained relatively constant since it was first measured in the 1980s and discharge at Pasture Canyon Spring has fluctuated for the period of record. Continuous records of surface-water discharge in the Black Mesa area were collected from streamflow-gaging stations at the following sites: Moenkopi Wash at Moenkopi 09401260 (1976 to 2018), Dinnebito Wash near Sand Springs 09401110 (1993 to 2018), Polacca Wash near Second Mesa 09400568 (1994 to 2018), and Pasture Canyon Springs 09401265 (2004 to 2018). Median winter flows (November through February) of each water year were used as an index of the amount of groundwater discharge at the above-named sites. For the period of record, the median winter flows have generally remained constant at Dinnebito Wash and Polacca Wash, whereas a decreasing trend was indicated at Moenkopi Wash and Pasture Canyon Springs. In 2017 and 2018, water samples collected from two wells, four springs, and three streams in the Black Mesa area were analyzed for selected chemical constituents. The results from wells and springs were compared with previous analyses from the same wells and springs. At the Peabody 2 well, a significant (p<0.05) decreasing trend in dissolved solids over time was found, while concentrations of dissolved solids have not varied significantly (p>0.05) at the Kykotsmovi PM2 well. Dissolved solids, chloride, and sulfate concentrations increased at Moenkopi School Spring during the more than 30 years of record at that site. Concentrations of dissolved solids, chloride, and sulfate at Pasture Canyon Spring have not varied significantly (p>0.05) since the early 1980s, and there is no increasing or decreasing trend in those data. Concentrations of dissolved solids, chloride, and sulfate at Burro Spring and Unnamed Spring near Dennehotso have varied for the period of record, but there is no statistical trend in the data. Baseflow water chemistry samples were collected from Moenkopi, Dinnebito, and Polacca washes in 2017. Samples from all three washes had total-dissolved solids concentrations higher than is typically found in the N aquifer water.

Arizona↗

Description and evaluation of numerical groundwater flow models for the Edwards Aquifer, south-central Texas

A substantial number of public water system wells in south-central Texas withdraw groundwater from the karstic, highly productive Edwards aquifer. However, the use of numerical groundwater flow models to aid in the delineation of contributing areas for public water system wells in the Edwards aquifer is problematic because of the complex hydrogeologic framework and the presence of conduit-dominated flow paths in the aquifer. The U.S. Geological Survey, in cooperation with the Texas Commission on Environmental Quality, evaluated six published numerical groundwater flow models (all deterministic) that have been developed for the Edwards aquifer San Antonio segment or Barton Springs segment, or both. This report describes the models developed and evaluates each with respect to accessibility and ease of use, range of conditions simulated, accuracy of simulations, agreement with dye-tracer tests, and limitations of the models. These models are (1) GWSIM model of the San Antonio segment, a FORTRAN computer-model code that pre-dates the development of MODFLOW; (2) MODFLOW conduit-flow model of San Antonio and Barton Springs segments; (3) MODFLOW diffuse-flow model of San Antonio and Barton Springs segments; (4) MODFLOW Groundwater Availability Modeling [GAM] model of the Barton Springs segment; (5) MODFLOW recalibrated GAM model of the Barton Springs segment; and (6) MODFLOW–DCM (dual conductivity model) conduit model of the Barton Springs segment. The GWSIM model code is not commercially available, is limited in its application to the San Antonio segment of the Edwards aquifer, and lacks the ability of MODFLOW to easily incorporate newly developed processes and packages to better simulate hydrologic processes. MODFLOW is a widely used and tested code for numerical modeling of groundwater flow, is well documented, and is in the public domain. These attributes make MODFLOW a preferred code with regard to accessibility and ease of use. The MODFLOW conduit-flow model incorporates improvements over previous models by using (1) a user-friendly interface, (2) updated computer codes (MODFLOW–96 and MODFLOW–2000), (3) a finer grid resolution, (4) less-restrictive boundary conditions, (5) an improved discretization of hydraulic conductivity, (6) more accurate estimates of pumping stresses, (7) a long transient simulation period (54 years, 1947–2000), and (8) a refined representation of high-permeability zones or conduits. All of the models except the MODFLOW–DCM conduit model have limitations resulting from the use of Darcy’s law to simulate groundwater flow in a karst aquifer system where non-Darcian, turbulent flow might actually dominate. The MODFLOW–DCM conduit model is an improvement in the ability to simulate karst-like flow conditions in conjunction with porous-media-type matrix flow. However, the MODFLOW–DCM conduit model has had limited application and testing and currently (2008) lacks commercially available pre- and post-processors. The MODFLOW conduit-flow and diffuse-flow Edwards aquifer models are limited by the lack of calibration for the northern part of the Barton Springs segment (Travis County) and their reliance on the use of the calibrated hydraulic conductivity and storativity values from the calibrated Barton Springs segment GAM model. The major limitation of the Barton Springs segment GAM and recalibrated GAM models is that they were calibrated to match measured water levels and springflows for a restrictive range of hydrologic conditions, with each model having different hydraulic conductivity and storativity values appropriate to the hydrologic conditions that were simulated. The need for two different sets of hydraulic conductivity and storativity values increases the uncertainty associated with the accuracy of either set of values, illustrates the non-uniqueness of the model solution, and probably most importantly demonstrates the limitations of using a one-layer model to represent the heterogeneous hydrostratigraphic units composing the Edwards aquifer. In general, the best matches or agreement between groundwater flow directions inferred by numerical model simulation, and by dye-tracer tests, are observed where model outputs accurately reproduce the configuration of the potentiometric surface with regard to the positions of major and minor groundwater troughs and divides. None of the models, with the possible exception of the MODFLOW–DCM conduit model, has a documented capability to accurately simulate travel times for conduit-dominated velocities in the Edwards aquifer. Public water system assessments of wells in the Barton Springs segment of the Edwards aquifer, and elsewhere where conduit-flow conditions are thought to dominate aquifer hydraulic behavior, might be enhanced by use of either the MODFLOW–DCM model or the newly developed U.S. Geological Survey MODFLOW Conduit-Flow Process module for MODFLOW–2005, because each incorporates a type of dual or triple hydraulic conductivity approach and has the capability to explicitly simulate turbulent flow and conduit hydraulic characteristics.

Texas↗

Review of Trace-Element Field-Blank Data Collected for the California Groundwater Ambient Monitoring and Assessment (GAMA) Program, May 2004-January 2008

Trace-element quality-control samples (for example, source-solution blanks, field blanks, and field replicates) were collected as part of a statewide investigation of groundwater quality in California, known as the Priority Basins Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The GAMA Priority Basins Project is being conducted by the U.S. Geological Survey (USGS) in cooperation with the California State Water Resources Control Board (SWRCB) to assess and monitor the quality of groundwater resources used for drinking-water supply and to improve public knowledge of groundwater quality in California. Trace-element field blanks were collected to evaluate potential bias in the corresponding environmental data. Bias in the environmental data could result from contamination in the field during sample collection, from the groundwater coming into contact with contaminants on equipment surfaces or from other sources, or from processing, shipping, or analyzing the samples. Bias affects the interpretation of environmental data, particularly if any constituents are present solely as a result of extrinsic contamination that would have otherwise been absent from the groundwater that was sampled. Field blanks were collected, analyzed, and reviewed to identify and quantify extrinsic contamination bias. Data derived from source-solution blanks and laboratory quality-control samples also were considered in evaluating potential contamination bias. Eighty-six field-blank samples collected from May 2004 to January 2008 were analyzed for the concentrations of 25 trace elements. Results from these field blanks were used to interpret the data for the 816 samples of untreated groundwater collected over the same period. Constituents analyzed were aluminum (Al), antimony (Sb), arsenic (As), barium (Ba), beryllium (Be), boron (B), cadmium (Cd), chromium (Cr), cobalt (Co), copper (Cu), iron (Fe), lead (Pb), lithium (Li), manganese (Mn), mercury (Hg), molybdenum (Mo), nickel (Ni), selenium (Se), silver (Ag), strontium (Sr), thallium (Tl), tungsten (W), uranium (U), vanadium (V), and zinc (Zn). The detection frequency and the 90th percentile concentration at greater than 90 percent confidence were determined from the field-blank data for each trace element, and these results were compared to each constituent's long-term method detection level (LT-MDL) to determine whether a study reporting level (SRL) was necessary to ensure that no more than 10 percent of the detections in groundwater samples could be attributed solely to contamination bias. Only two of the trace elements analyzed, Li and Se, had zero detections in the 86 field blanks. Ten other trace elements (Sb, As, Be, B, Cd, Co, Mo, Ag, Tl, and U) were detected in fewer than 5 percent of the field blanks. The field-blank results for these constituents did not necessitate establishing SRLs. Of the 13 constituents that were detected in more than 5 percent of the field blanks, six (Al, Ba, Cr, Mn, Hg, and V) had field-blank results that indicated a need for SRLs that were at or below the highest laboratory reporting levels (LRL) used during the sampling period; these SRLs were needed for concentrations between the LT-MDLs and LRLs. The other seven constituents with detection frequencies above 5 percent (Cu, Fe, Pb, Ni, Sr, W, and Zn) had field-blank results that necessitated SRLs greater than the highest LRLs used during the study period. SRLs for these seven constituents, each set at the 90th percentile of their concentrations in the field blanks, were at least an order of magnitude below the regulatory thresholds established for drinking water for health or aesthetic purposes; therefore, reporting values below the SRLs as less than or equal to (=) the measured value would not prevent the identification of values greater than the drinking-water thresholds. The SRLs and drinking-water thresholds, respectively, for these 7 trace elements are Cu (1.7 ?g/L and 1,300

Scientific Investigations Report↗

Hydrogeologic framework, groundwater and surface-water systems, land use, pumpage, and water budget of the Chamokane Creek basin, Stevens County, Washington

A study of the water resources of the unconsolidated groundwater system of the Chamokane Creek basin was conducted to determine the hydrogeologic framework, interactions of shallow and deep parts of the groundwater system with each other and the surface-water system, changes in land use and land cover, and water-use estimates. Chamokane Creek basin is a 179 mi2 area that borders and partially overlaps the Spokane Indian Reservation in southern Stevens County in northeastern Washington State. Aquifers within the Chamokane Creek basin are part of a sequence of glaciofluvial and glaciolacustrine sediment that may reach total thicknesses of about 600 ft. In 1979, most of the water rights in the Chamokane Creek basin were adjudicated by the United States District Court requiring regulation in favor of the Spokane Tribe of Indians' senior water right. The Spokane Tribe, the State of Washington, and the United States are concerned about the effects of additional groundwater development within the basin on Chamokane Creek. Information provided by this study will be used to evaluate the effects of potential increases in groundwater withdrawals on groundwater and surface-water resources within the basin. The hydrogeologic framework consists of six hydrogeologic units: The Upper outwash aquifer, the Landslide Unit, the Valley Confining Unit, the Lower Aquifer, the Basalt Unit, and the Bedrock Unit. The Upper outwash aquifer occurs along the valley floors of the study area and consists of sand, gravel, cobbles, boulders, with minor silt and (or) clay interbeds in places. The Lower aquifer is a confined aquifer consisting of sand and gravel that occurs at depth below the Valley confining unit. Median horizontal hydraulic conductivity values for the Upper outwash aquifer, Valley confining unit, Lower aquifer, and Basalt unit were estimated to be 540, 10, 19, and 3.7 ft/d, respectively. Many low-flow stream discharge measurements at sites on Chamokane Creek and its tributaries were at or near zero flow. The most notable exception is where Chamokane Creek is supported by discharge of large springs from the Upper outwash aquifer in the southern part of the basin. Most high-flow measurements indicated gains in streamflow (groundwater discharging to the stream). Large streamflow losses, however, were recorded near the north end of Walkers Prairie where streamflow directly recharges the Upper outwash aquifer. The similarity in seasonal water-level fluctuations in the Upper outwash aquifer and the Lower aquifer indicate that these systems may be fairly well connected. Land use and land cover change analysis indicates that Chamokane Creek basin has been dominated by forests with some pasture and agricultural lands with sparse residential development from the 1980s to present. Loss in forest cover represents the largest change in land cover in the basin between 1987 and 2009. This appears to be mostly due to forestry activities, especially in the northern part of the basin. Since 1987, more than 18,000 acres of evergreen forest have been logged and are at various stages of regrowth. Estimated average annual total groundwater pumpage in the basin increased from 224 million gallons per year (Mgal/yr) in 1980 to 1,330 Mgal/yr in 2007. The largest withdrawals during 2007 were to supply two fish hatcheries, with a combined total annual pumpage of about 1,150 Mgal. Annual groundwater pumpage values from 1980 through 2007 for the study area ranged from 21.1 to 28.9 Mgal/yr for domestic wells and 0.38 to 23.7 Mgal/yr for public supply. An approximate water budget for a typical year in the Chamokane Creek basin indicates that 19.6 in. of precipitation are balanced by 4.7 in. of streamflow discharge from the basin, and 14.9 in. of evapotranspiration.

Washington↗

Numerical simulation of the groundwater-flow system in the Chambers-Clover Creek Watershed and Vicinity, Pierce County, Washington

A groundwater-flow model was developed to contribute to an improved understanding of water resources in the Chambers-Clover Creek Watershed. The model covers an area of about 491 square miles in western Pierce County, Washington, and is bounded to the northeast by the Puyallup River valley, to the southwest by the Nisqually River valley, and extends northwest to Puget Sound, and southeast to Tanwax Creek. The Puyallup and Nisqually Rivers occupy large, relatively flat alluvial valleys that are separated by a broad, poorly drained, upland area that covers most of the model area. Chambers and Clover Creeks drain much of the central uplands and flow westward to Puget Sound. The model area is underlain by a northwest-thickening sequence of unconsolidated glacial (till and outwash) and interglacial (fluvial and lacustrine) deposits. Ten unconsolidated hydrogeologic units in the model area form the basis of the groundwater-flow model. Groundwater flow in the Chambers-Clover Creek Watershed and vicinity was simulated using the groundwater-flow model, MODFLOW-2000. The finite-difference model grid comprises 146 rows, 132 columns, and 11 layers. Each model cell has a horizontal dimension of 1,000 by 1,000 feet, and the model contains a total of 123,602 active cells. The thickness of model layers varies throughout the model area and ranges from 1.5 feet in the A3 aquifer unit to 1,567 feet in the G undifferentiated unit. Groundwater flow was simulated for both steady-state and transient conditions. Steady-state conditions were simulated using average recharge, discharge, and water levels for the 24-month period September 2006-August 2008. Transient conditions were simulated for the period September 2006-August 2008 using 24 monthly stress periods. Resource managers and local stakeholders intend to use the model to evaluate a range of water resource issues under both steady-state and transient conditions. Initial conditions for the transient model were developed from a 3-year "lead-in" period that used recorded precipitation and river levels, and temporal extrapolations of other boundary conditions. During model calibration, variables were adjusted within probable ranges to minimize differences between measured and simulated groundwater levels and stream baseflows. The model as calibrated to steady-state conditions has a standard deviation for heads and flows of 28.42 feet and 2.12 cubic feet per second, respectively; the model as calibrated to transient conditions has a standard deviation for heads and flows of 23.01 feet and 2.67 cubic feet per second, respectively. Simulated steady-state inflow to the model area from precipitation and secondary recharge was 477,266 acre-feet per year (acre-ft/yr) (79 percent of total simulated inflow), and simulated inflow from stream and lake leakage was 129,778 acre-ft/yr (21 percent of total simulated inflow). Simulated outflow from the model primarily was through discharge to streams, lakes, springs, seeps, and Puget Sound (559,192 acre-ft/yr; 92 percent of total simulated outflow), and withdrawals from wells (47,863 acre-ft/yr; 8 percent of total simulated outflow). Six scenarios were formulated and simulated using the calibrated model to provide representative examples of how the model can be used to evaluate the effects on groundwater levels and stream baseflows of potential changes in groundwater withdrawals, in consumptive use, and in recharge.

Washington↗

A comparison of U.S. Geological Survey three-dimensional model estimates of groundwater source areas and velocities to independently derived estimates, Idaho National Laboratory and vicinity, Idaho

The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy, evaluated a three-dimensional model of groundwater flow in the fractured basalts and interbedded sediments of the eastern Snake River Plain aquifer at and near the Idaho National Laboratory to determine if model-derived estimates of groundwater movement are consistent with (1) results from previous studies on water chemistry type, (2) the geochemical mixing at an example well, and (3) independently derived estimates of the average linear groundwater velocity. Simulated steady-state flow fields were analyzed using backward particle-tracking simulations that were based on a modified version of the particle tracking program MODPATH. Model results were compared to the 5-microgram-per-liter lithium contour interpreted to represent the transition from a water type that is primarily composed of tributary valley underflow and streamflow-infiltration recharge to a water type primarily composed of regional aquifer water. This comparison indicates several shortcomings in the way the model represents flow in the aquifer. The eastward movement of tributary valley underflow and streamflow-infiltration recharge is overestimated in the north-central part of the model area and underestimated in the central part of the model area. Model inconsistencies can be attributed to large contrasts in hydraulic conductivity between hydrogeologic zones. Sources of water at well NPR-W01 were identified using backward particle tracking, and they were compared to the relative percentages of source water chemistry determined using geochemical mass balance and mixing models. The particle tracking results compare reasonably well with the chemistry results for groundwater derived from surface-water sources (-28 percent error), but overpredict the proportion of groundwater derived from regional aquifer water (108 percent error) and underpredict the proportion of groundwater derived from tributary valley underflow from the Little Lost River valley (-74 percent error). These large discrepancies may be attributed to large contrasts in hydraulic conductivity between hydrogeologic zones and (or) a short-circuiting of underflow from the Little Lost River valley to an area of high hydraulic conductivity. Independently derived estimates of the average groundwater velocity at 12 well locations within the upper 100 feet of the aquifer were compared to model-derived estimates. Agreement between velocity estimates was good at wells with travel paths located in areas of sediment-rich rock (root-mean-square error [RMSE] = 5.2 feet per day [ft/d]) and poor in areas of sediment-poor rock (RMSE = 26.2 ft/d); simulated velocities in sediment-poor rock were 2.5 to 4.5 times larger than independently derived estimates at wells USGS 1 (less than 14 ft/d) and USGS 100 (less than 21 ft/d). The models overprediction of groundwater velocities in sediment-poor rock may be attributed to large contrasts in hydraulic conductivity and a very large, model-wide estimate of vertical anisotropy (14,800).

Idaho↗

A comparison of U.S. Geological Survey three-dimensional model estimates of groundwater source areas and velocities to independently derived estimates, Idaho National Laboratory and vicinity, Idaho

The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy, evaluated a three-dimensional model of groundwater flow in the fractured basalts and interbedded sediments of the eastern Snake River Plain aquifer at and near the Idaho National Laboratory to determine if model-derived estimates of groundwater movement are consistent with (1) results from previous studies on water chemistry type, (2) the geochemical mixing at an example well, and (3) independently derived estimates of the average linear groundwater velocity. Simulated steady-state flow fields were analyzed using backward particle-tracking simulations that were based on a modified version of the particle tracking program MODPATH. Model results were compared to the 5-microgram-per-liter lithium contour interpreted to represent the transition from a water type that is primarily composed of tributary valley underflow and streamflow-infiltration recharge to a water type primarily composed of regional aquifer water. This comparison indicates several shortcomings in the way the model represents flow in the aquifer. The eastward movement of tributary valley underflow and streamflow-infiltration recharge is overestimated in the north-central part of the model area and underestimated in the central part of the model area. Model inconsistencies can be attributed to large contrasts in hydraulic conductivity between hydrogeologic zones. Sources of water at well NPR-W01 were identified using backward particle tracking, and they were compared to the relative percentages of source water chemistry determined using geochemical mass balance and mixing models. The particle tracking results compare reasonably well with the chemistry results for groundwater derived from surface-water sources (-28 percent error), but overpredict the proportion of groundwater derived from regional aquifer water (108 percent error) and underpredict the proportion of groundwater derived from tributary valley underflow from the Little Lost River valley (-74 percent error). These large discrepancies may be attributed to large contrasts in hydraulic conductivity between hydrogeologic zones and (or) a short-circuiting of underflow from the Little Lost River valley to an area of high hydraulic conductivity. Independently derived estimates of the average groundwater velocity at 12 well locations within the upper 100 feet of the aquifer were compared to model-derived estimates. Agreement between velocity estimates was good at wells with travel paths located in areas of sediment-rich rock (root-mean-square error [RMSE] = 5.2 feet per day [ft/d]) and poor in areas of sediment-poor rock (RMSE = 26.2 ft/d); simulated velocities in sediment-poor rock were 2.5 to 4.5 times larger than independently derived estimates at wells USGS 1 (less than 14 ft/d) and USGS 100 (less than 21 ft/d). The models overprediction of groundwater velocities in sediment-poor rock may be attributed to large contrasts in hydraulic conductivity and a very large, model-wide estimate of vertical anisotropy (14,800).

Idaho↗

Simulation of groundwater and surface-water interaction and effects of pumping in a complex glacial-sediment aquifer, east central Massachusetts

The effects of groundwater pumping on surface-water features were evaluated by use of a numerical groundwater model developed for a complex glacial-sediment aquifer in northeastern Framingham, Massachusetts, and parts of surrounding towns. The aquifer is composed of sand, gravel, silt, and clay glacial-fill sediments up to 270 feet thick over an irregular fractured bedrock surface. Surface-water bodies, including Cochituate Brook, the Sudbury River, Lake Cochituate, Dudley Pond, and adjoining wetlands, are in hydraulic connection with the aquifer and can be affected by groundwater withdrawals. Groundwater and surface-water interaction was simulated with MODFLOW-NWT under current conditions and a variety of hypothetical pumping conditions. Simulations of hypothetical pumping at reactivated water supply wells indicate that captured groundwater would decrease baseflow to the Sudbury River and induce recharge from Lake Cochituate. Under constant (steady-state) pumping, induced groundwater recharge from Lake Cochituate was equal to about 32 percent of the simulated pumping rate, and flow downstream in the Sudbury River decreased at the same rate as pumping. However, surface water responded quickly to pumping stresses. When pumping was simulated for 1 month and then stopped, streamflow depletions decreased by about 80 percent within 2 months and by about 90 percent within about 4 months. The fast surface water response to groundwater pumping offers the potential to substantially reduce streamflow depletions during periods of low flow, which are of greatest concern to the ecological integrity of the river. Results indicate that streamflow depletion during September, typically the month of lowest flow, can be reduced by 29 percent by lowering the maximum pumping rates to near zero during September. Lowering pumping rates for 3 months (July through September) reduces streamflow depletion during September by 79 percent as compared to constant pumping. These results demonstrate that a seasonal or streamflow-based groundwater pumping schedule can reduce the effects of pumping during periods of low flow.

Massachusetts↗

Evaluation of groundwater quality and selected hydrologic conditions in the South Coast aquifer, Santa Isabel area, Puerto Rico, 2008–09

The source of drinking water in the Santa Isabel and Coamo areas of Puerto Rico (Molina and Gómez-Gómez, 2008) is the South Coast aquifer (hereafter referred to as the aquifer), which supplies about 30,700 cubic meters per day (m³/d) to Puerto Rico Aqueduct and Sewer Authority (PRASA) public-supply wells. In addition, approximately 45 wells provide an estimated 33,700 m³/d of groundwater to irrigate crops in the area. In 1967, baseline nitrate concentrations in groundwater throughout most of the aquifer were generally less than 6 milligrams per liter (mg/L) as nitrogen in collected water samples (U.S. Geological Survey, 2012). In 2007, elevated nitrate concentrations were detected in the aquifer, near Santa Isabel and the foothills north of the coastal plain at Santa Isabel as part of a regional groundwater-quality assessment conducted by the U.S. Geological Survey (USGS) during 2007 (Rodríguez and Gómez-Gómez, 2008). The increase in nitrate concentrations has been of concern to local government agencies because of its potential effect on public supply. To address public-supply concerns, the USGS, in cooperation with the Puerto Rico Department of Natural and Environmental Resources (PRDNER), evaluated groundwater quality in the aquifer near the Santa Isabel area between January 2008 and May 2009. The objectives of the study were to (1) define the groundwater-quality conditions of the aquifer, with emphasis on the distribution of nitrate concentrations; (2) identify potential sources leading to elevated nitrate concentrations; (3) estimate the nitrate loads from major sources identified; and (4) estimate the groundwater withdrawals by principal-use categories in the area. Results of this study will be used by Commonwealth of Puerto Rico and Federal agencies in developing strategies that can result in containment of high nitrate groundwater to minimize degradation of fresh groundwater in the aquifer.

Santa Isabel↗

Pesticides in Wyoming Groundwater, 2008-10

Groundwater samples were collected from 296 wells during 1995-2006 as part of a baseline study of pesticides in Wyoming groundwater. In 2009, a previous report summarized the results of the baseline sampling and the statistical evaluation of the occurrence of pesticides in relation to selected natural and anthropogenic (human-related) characteristics. During 2008-10, the U.S. Geological Survey, in cooperation with the Wyoming Department of Agriculture, resampled a subset (52) of the 296 wells sampled during 1995-2006 baseline study in order to compare detected compounds and respective concentrations between the two sampling periods and to evaluate the detections of new compounds. The 52 wells were distributed similarly to sites used in the 1995-2006 baseline study with respect to geographic area and land use within the geographic area of interest. Because of the use of different types of reporting levels and variability in reporting-level values during both the 1995-2006 baseline study and the 2008-10 resampling study, analytical results received from the laboratory were recensored. Two levels of recensoring were used to compare pesticides—a compound-specific assessment level (CSAL) that differed by compound and a common assessment level (CAL) of 0.07 microgram per liter. The recensoring techniques and values used for both studies, with the exception of the pesticide 2,4-D methyl ester, were the same. Twenty-eight different pesticides were detected in samples from the 52 wells during the 2008-10 resampling study. Pesticide concentrations were compared with several U.S. Environmental Protection Agency drinking-water standards or health advisories for finished (treated) water established under the Safe Drinking Water Act. All detected pesticides were measured at concentrations smaller than U.S. Environmental Protection Agency drinking-water standards or health advisories where applicable (many pesticides did not have standards or advisories). One or more pesticides were detected at concentrations greater than the CAL in water from 16 of 52 wells sampled (about 31 percent) during the resampling study. Detected pesticides were classified into one of six types: herbicides, herbicide degradates, insecticides, insecticide degradates, fungicides, or fungicide degradates. At least 95 percent of detected pesticides were classified as herbicides or herbicide degradates. The number of different pesticides detected in samples from the 52 wells was similar between the 1995-2006 baseline study (30 different pesticides) and 2008-2010 resampling study (28 different pesticides). Thirteen pesticides were detected during both studies. The change in the number of pesticides detected (without regard to which pesticide was detected) in groundwater samples from each of the 52 wells was evaluated and the number of pesticides detected in groundwater did not change for most of the wells (32). Of those that did have a difference between the two studies, 17 wells had more pesticide detections in groundwater during the 1995-2006 baseline study, whereas only 3 wells had more detections during the 2008-2010 resampling study. The difference in pesticide concentrations in groundwater samples from each of the 52 wells was determined. Few changes in concentration between the 1995-2006 baseline study and the 2008-2010 resampling study were seen for most detected pesticides. Seven pesticides had a greater concentration detected in the groundwater from the same well during the baseline sampling compared to the resampling study. Concentrations of prometon, which was detected in 17 wells, were greater in the baseline study sample compared to the resampling study sample from the same well 100 percent of the time. The change in the number of pesticides detected (without regard to which pesticide was detected) in groundwater samples from each of the 52 wells with respect to land use and geographic area was calculated. All wells with land use classified as agricultural had the same or a smaller number of pesticides detected in the resampling study compared to the baseline study. All wells in the Bighorn Basin geographic area also had the same or a smaller number of pesticides detected in the resampling study compared to the baseline study.

Wyoming↗

Geology, water-quality, hydrology, and geomechanics of the Cuyama Valley groundwater basin, California, 2008--12

To assess the water resources of the Cuyama Valley groundwater basin in Santa Barbara County, California, a series of cooperative studies were undertaken by the U.S. Geological Survey and the Santa Barbara County Water Agency. Between 2008 and 2012, geologic, water-quality, hydrologic and geomechanical data were collected from selected sites throughout the Cuyama Valley groundwater basin. Geologic data were collected from three multiple-well groundwater monitoring sites and included lithologic descriptions of the drill cuttings, borehole geophysical logs, temperature logs, as well as bulk density and sonic velocity measurements of whole-core samples. Generalized lithologic characterization from the monitoring sites indicated the water-bearing units in the subsurface consist of unconsolidated to partly consolidated sand, gravel, silt, clay, and occasional cobbles within alluvial fan and stream deposits. Analysis of geophysical logs indicated alternating layers of finer- and coarser-grained material that range from less than 1 foot to more than 20 feet thick. On the basis of the geologic data collected, the principal water-bearing units beneath the monitoring-well sites were found to be composed of younger alluvium of Holocene age, older alluvium of Pleistocene age, and the Tertiary-Quaternary Morales Formation. At all three sites, the contact between the recent fill and younger alluvium is approximately 20 feet below land surface. Water-quality samples were collected from 12 monitoring wells, 27 domestic and supply wells, 2 springs, and 4 surface-water sites and were analyzed for a variety of constituents that differed by site, but, in general, included trace elements; nutrients; dissolved organic carbon; major and minor ions; silica; total dissolved solids; alkalinity; total arsenic and iron; arsenic, chromium, and iron species; and isotopic tracers, including the stable isotopes of hydrogen and oxygen, activities of tritium, and carbon-14 abundance. Of the 39 wells sampled, concentrations of total dissolved solids and sulfate from 38 and 37 well samples, respectively, were greater than the U.S. Environmental Protection Agency’s secondary maximum contaminant levels. Concentrations greater than the maximum contaminant levels for nitrate were observed in five wells and were observed for arsenic in four wells. Differences in the stable-isotopic values of hydrogen and oxygen among groundwater samples indicated that water does not move freely between different formations or between different zones within the Cuyama Valley. Variations in isotopic composition indicated that recharge is derived from several different sources. The age of the groundwater, expressed as time since recharge, was between 600 and 38,000 years before present. Detectable concentrations of tritium indicated that younger water, recharged since the early 1950s, is present in parts of the groundwater basin. Hydrologic data were collected from 12 monitoring wells, 56 domestic and supply wells, 3 surface-water sites, and 4 rainfall-gaging stations. Rainfall in the valley averaged about 8 inches annually, whereas the mountains to the south received between 12 and 19 inches. Stream discharge records showed seasonal variability in surface-water flows ranging from no-flow to over 1,500 cubic feet per second. During periods when inflow to the valley exceeds outflow, there is potential recharge from stream losses to the groundwater system Water-level records included manual quarterly depth-to-water measurements collected from 68 wells, time-series data collected from 20 of those wells, and historic water levels from 16 wells. Hydrographs of the manual measurements showed declining water levels in 16 wells, mostly in the South-Main zone, and rising water levels in 14 wells, mostly in the Southern Ventucopa Uplands. Time-series hydrographs showed daily, seasonal, and longer-term effects associated with local pumping. Water-level data from the multiple-well monitoring sites indicated seasonal fluctuations as great as 80 feet and water-level differences between aquifers as great as 40 feet during peak pumping season. Hydrographs from the multiple-well groundwater monitoring sites showed vertical hydraulic gradients were upward during the winter months and downward during the irrigation season. Historic hydrographs showed water-level declines in the Southern-Main, Western Basin, Caliente Northern-Main, and Southern Sierra Madre zone ranging from 1 to 7 feet per year. Hydrographs of wells in the Southern Ventucopa Uplands zone showed several years with marked increases in water levels that corresponded to increased precipitation in the Cuyama Valley. Investigation of hydraulic properties included hydraulic conductivity and transmissivity estimated from aquifer tests performed on 63 wells. Estimates of horizontal hydraulic conductivity ranged from about 1.5 to 28 feet per day and decreased with depth. The median estimated hydraulic conductivity for the older alluvium was about five times that estimated for the Morales Formation. Estimates of transmissivity ranged from 560 to 163,400 gallons per day per foot and decreased with depth. The median estimated transmissivity for the younger alluvium was about three times that estimated for the older alluvium. Geomechanical analysis included land-surface elevation changes at five continuously operating global positioning systems (GPS) and land-subsidence detection at five interferometric synthetic aperture radar (InSAR) reference points. Analysis of data collected from continuously operating GPS stations showed the mountains to the south and west moved upward about 1 millimeter (mm) annually, whereas the station in the center of the Southern-Main zone moved downward more than 7 mm annually, indicating subsidence. It is likely that this subsidence is inelastic (permanent) deformation and indicates reduced storage capacity in the aquifer sediments. Analysis of InSAR data showed local and regional changes that appeared to be dependent, in part, on the time span of the interferogram, seasonal variations in pumping, and tectonic uplift. Long-term InSAR time series showed a total maximum detected subsidence rate of approximately 12 mm per year at one location and approximately 8 mm per year at a second location, while short-term InSAR time series showed maximum subsidence of about 15 mm at one location and localized maximum uplift of about 10 mm at another location.

California↗

Recharge sources and residence times of groundwater as determined by geochemical tracers in the Mayfield Area, southwestern Idaho, 2011–12

Parties proposing residential development in the area of Mayfield, Idaho are seeking a sustainable groundwater supply. During 2011–12, the U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources, used geochemical tracers in the Mayfield area to evaluate sources of aquifer recharge and differences in groundwater residence time. Fourteen groundwater wells and one surface-water site were sampled for major ion chemistry, metals, stable isotopes, and age tracers; data collected from this study were used to evaluate the sources of groundwater recharge and groundwater residence times in the area. Major ion chemistry varied along a flow path between deeper wells, suggesting an upgradient source of dilute water, and a downgradient source of more concentrated water with the geochemical signature of the Idaho Batholith. Samples from shallow wells had elevated nutrient concentrations, a more positive oxygen-18 signature, and younger carbon-14 dates than deep wells, suggesting that recharge comes from young precipitation and surface-water infiltration. Samples from deep wells generally had higher concentrations of metals typical of geothermal waters, a more negative oxygen-18 signature, and older carbon-14 values than samples from shallow wells, suggesting that recharge comes from both infiltration of meteoric water and another source. The chemistry of groundwater sampled from deep wells is somewhat similar to the chemistry in geothermal waters, suggesting that geothermal water may be a source of recharge to this aquifer. Results of NETPATH mixing models suggest that geothermal water composes 1–23 percent of water in deep wells. Chlorofluorocarbons were detected in every sample, which indicates that all groundwater samples contain at least a component of young recharge, and that groundwater is derived from multiple recharge sources. Conclusions from this study can be used to further refine conceptual hydrological models of the area.

Idaho↗

Numerical simulation of the groundwater-flow system in Chimacum Creek Basin and vicinity, Jefferson County, Washington

A groundwater-flow model was developed to evaluate potential future effects of growth and of water-management strategies on water resources in the Chimacum Creek Basin. The model covers an area of about 64 square miles (mi 2 ) on the Olympic Peninsula in northeastern Jefferson County, Washington. The Chimacum Creek Basin drains an area of about 53 mi 2 and consists of Chimacum Creek and its tributary East Fork Chimacum Creek, which converge near the town of Chimacum and discharge to Port Townsend Bay near the town of Irondale. The topography of the model area consists of north-south oriented, narrow, regularly spaced parallel ridges and valleys that are characteristic of fluted glaciated surfaces. Thick accumulations of peat occur along the axis of East Fork Chimacum Creek and provide rich soils for agricultural use. The study area is underlain by a north-thickening sequence of unconsolidated glacial (till and outwash) and interglacial (fluvial and lacustrine) deposits, and sedimentary and igneous bedrock units that crop out along the margins and the western interior of the model area. Six hydrogeologic units in the model area form the basis of the groundwater-flow model. They are represented by model layers UC (upper confining), UA (upper aquifer), MC (middle confining), LA (lower aquifer), LC (lower confining), and OE (bedrock). Groundwater flow in the Chimacum Creek Basin and vicinity was simulated using the groundwater-flow model, MODFLOW-2005. The finite-difference model grid comprises 245 columns, 313 rows, and 6 layers. Each model cell has a horizontal dimension of 200 × 200 feet (ft). The thickness of model layers varies throughout the model area and ranges from 5 ft in the non-bedrock units to more than 2,400 ft in the bedrock. Groundwater flow was simulated for steady-state conditions, which were simulated for calibration of the model using average recharge, discharge, and water levels for the 180-month period October 1994–September 2009. The model as calibrated has a mean residual of 4.5 ft and a standard error on the mean of 2.1 ft for heads, and 0.64±0.42 cubic feet per second for streamflows. After the model was calibrated, a Current Conditions simulation was developed to reflect current (October 2008–September 2009) hydrologic conditions, with representative pumping, return flows, and “normal” recharge (based on National Weather Service average precipitation for 1981 to 2010). The Current Conditions simulation was used to estimate current flow quantities, and as a basis to compare other simulations.Simulated steady-state inflow to the model area from precipitation and secondary recharge, or “return flow,” was 16,347 acre-feet per year (acre-ft/yr); groundwater inflow from other basins to the north of the model boundary was 1,518 acre-ft/yr (net, 3,114 acre-ft/yr in and 1,596 acre-ft/yr out) and simulated inflow from lake leakage was 613 acre-ft/yr (net, 684 acre-ft/yr in and 71 acre-ft/yr out). Simulated outflow from the model primarily was through discharge to Puget Sound (10,022 acre-ft/yr), streams (5,424 acre-ft/yr ), springs and seeps (1,521 acre-ft/yr), and through withdrawals from wells (1,506 acre-ft/yr). Four simulations were formulated using the calibrated model—one to represent current conditions (2009, the end of the period used for calibration) and three to provide representative examples of how the model can be used to evaluate the relative effects of potential changes in groundwater withdrawals and consumptive use on groundwater levels and stream base flows: Probable Future Use, based on population projections; Full Beneficial Use, based on Jefferson County Public Utility District #1 water rights; Sanitary Sewer, based on eliminating septic return flows in the Urban Growth Area. Particle tracking was used to assess flowpaths from sources and to sinks, and the effects of the presence of irrigation wells and their depths was assessed.

Washington↗