Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Nature Scientific Data”

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 1,225 records · Page 68Linked to original sources

Potentiometric surfaces and water-level trends in the Cockfield (upper Claiborne) aquifer in southern Arkansas and the Wilcox (lower Wilcox) aquifer of northeastern and southern Arkansas, 2012

The Cockfield aquifer, located in southern Arkansas, is composed of Eocene-age sand beds found near the base of the Cockfield Formation of Claiborne Group. The Wilcox aquifer, located in northeastern and southern Arkansas, is composed of Paleocene-age sand beds found in the middle to lower part of the Wilcox Group. The Cockfield and Wilcox aquifers are primary sources of groundwater. In 2010, withdrawals from the Cockfield aquifer in Arkansas totaled 19.2 million gallons per day (Mgal/d), and withdrawals from the Wilcox aquifer totaled 36.5 Mgal/d. A study was conducted by the U.S. Geological Survey in cooperation with the Arkansas Natural Resources Commission and the Arkansas Geological Survey to measure water levels associated with the Cockfield aquifer and the Wilcox aquifer in northeastern and southern Arkansas. Water levels were measured at 43 wells completed in the Cockfield aquifer and 47 wells completed in the Wilcox aquifer in February and March 2012. Measurements from 2012 are presented as potentiometric-surface maps and in combination with measurements from 2006 as water-level difference maps. Trends in water-level change over time within the Cockfield and Wilcox aquifers were determined using the water-level difference maps and selected well hydrographs. The Cockfield aquifer study area in southern Arkansas is bounded on the east by the Mississippi River and on the west by the area that contains outcrops and subcrops of the Cockfield Formation. The northern boundary of the Cockfield aquifer study area is defined by the area that contains observation wells completed in the Cockfield aquifer and the southern boundary is the Louisiana State line. The Wilcox aquifer study area in northeastern Arkansas is bounded on the east by the Mississippi River and on the north by the Missouri State line. The southern and western boundaries are defined by areas containing observation wells completed in the Wilcox aquifer or by outcrop areas on or near Crowleys Ridge. The Wilcox aquifer study area in southern Arkansas is defined by observation wells completed in the Wilcox aquifer or by areas that contain outcrops of the Wilcox Group, or both. The potentiometric-surface map of the Cockfield aquifer shows the regional direction of groundwater flow was generally toward the east-southeast, except in areas of intense groundwater withdrawals such as southwestern Ashley County, where groundwater flows toward the town of Crossett. The highest water-level altitude measured was 350 feet (ft) above National Geodetic Vertical Datum of 1929 (NGVD 29) in central Columbia County. The lowest water-level altitude measured was 40 ft above NGVD 29 in southeastern Lincoln County. The water-level difference map for the Cockfield aquifer in Arkansas was constructed using 42 water-level measurements made during 2006 and 2012. The difference in water levels for the Cockfield aquifer ranged from 27.4 ft to -10.4 ft. The largest water-level rise was in Calhoun County, and the largest water-level decline was 10.4 ft in Union County. Of the 42 wells, 13 wells had a rise in water level, and the remaining 29 wells had a decline in water level. Hydrographs for 32 wells in the Cockfield aquifer with historical water-level data were evaluated using linear regression to calculate the annual rise or decline for each well. These data were aggregated by county and statistically evaluated for the range, mean, and median of water-level change in each county. Hydrographs for Bradley, Calhoun, Chicot, Columbia, and Union Counties indicated both rising and declining water levels. The mean annual water-level rise or decline for Calhoun County was 0.00 foot per year (ft/yr) or unchanged. The mean annual water-level for Ashley, Bradley, Chicot, Cleveland, Columbia, Lincoln, and Union Counties show declines ranging from -0.02 to -1.10 ft/yr. Two potentiometric-surface maps, one for the southern area and one for the northeastern area, were constructed to show the altitude of the water surface in the Wilcox aquifer. The direction of groundwater flow in the northeastern area was generally towards the south-southwest except for some areas immediately adjacent to the Mississippi River where the flow was more eastward towards the river. The highest water-level altitude was 219 ft in northern Mississippi County, and the lowest water-level altitude was 123 ft near West Memphis in Crittenden County. The direction of groundwater flow in the northern part of the southern area was generally towards the southwest. The direction of groundwater flow in the southern part was in all directions because of two cones of depression and two water-level mounds. The highest water-level altitude measured was 394 ft at the center of a water-level mound in eastern Hot Spring County and a water-level mound in southwestern Hempstead County. The lowest water-level altitude measured was 145 ft at the center of the cone of depression in Clark County. Water-level difference maps for the Wilcox aquifer in Arkansas were constructed using 47 water-level measurements made during 2006 and 2012. The difference in water levels for the Wilcox aquifer in the northeastern area ranged from 22.0 ft to -17.9 ft. The largest rise in water level occurred in Crittenden County, and the largest decline occurred in Lee County. Twenty-one wells had rising water levels, and 10 wells had declining water levels. The difference in water levels for the Wilcox aquifer in the southern area ranged from 18.1 ft to -4.2 ft. The largest rise and the largest decline in water level occurred in Nevada County. Twelve wells had rising water levels, and 4 wells had declining water levels. Linear regression analysis of long-term hydrographs was used to determine the mean annual water-level rise and decline in the Wilcox aquifer in the northeastern and southern areas of Arkansas. In the northeastern area, the mean annual water level declined in all seven counties. The mean annual declines ranged from -0.55 ft/yr in Craighead County to -1.46 ft/yr in St. Francis County. In the southern area, the annual rise and decline calculations for wells with over 20 years of records indicate rising and declining water levels in Clark, Hot Spring, and Nevada Counties. The mean annual water level declined in all counties except Hot Spring County.

Arkansas↗

Groundwater quality for 75 domestic wells in Lycoming County, Pennsylvania, 2014

Groundwater is a major source of drinking water in Lycoming County and adjacent counties in north-central and northeastern Pennsylvania, which are largely forested and rural and are currently undergoing development for hydrocarbon gases. Water-quality data are needed for assessing the natural characteristics of the groundwater resource and the potential effects from energy and mineral extraction, timber harvesting, agriculture, sewage and septic systems, and other human influences. This report, prepared in cooperation with Lycoming County, presents analytical data for groundwater samples from 75 domestic wells sampled throughout Lycoming County in June, July, and August 2014. The samples were collected using existing pumps and plumbing prior to any treatment and analyzed for physical and chemical characteristics, including nutrients, major ions, metals and trace elements, volatile organic compounds, gross-alpha particle and gross beta-particle activity, uranium, and dissolved gases, including methane and radon-222. Results indicate groundwater quality generally met most drinking-water standards, but that some samples exceeded primary or secondary maximum contaminant levels (MCLs) for arsenic, iron, manganese, total dissolved solids (TDS), chloride, pH, bacteria, or radon-222. Arsenic concentrations were higher than the MCL of 10 micrograms per liter (µg/L) in 9 of the 75 (12 percent) well-water samples, with concentrations as high as 23.6 μg/L; arsenic concentrations were higher than the health advisory level (HAL) of 2 μg/L in 23 samples (31 percent). Total iron concentrations exceeded the secondary maximum contaminant level (SMCL) of 300 μg/L in 20 of the 75 samples. Total manganese concentrations exceeded the SMCL of 50 μg/L in 20 samples and the HAL of 300 μg/L in 2 of those samples. Three samples had chloride concentrations that exceeded the SMCL of 250 milligrams per liter (mg/L); two of those samples exceeded the SMCL of 500 mg/L for TDS. The pH ranged from 5.3 to 9.15 and did not meet the SMCL range of 6.5 to 8.5 in 22 samples, with 17 samples having a pH less than 6.5 and 8 samples having pH greater than 8.5. Generally, the samples that had elevated TDS, chloride, or arsenic concentrations had high pH. Total coliform bacteria were detected in 39 of 75 samples (52 percent), with Escherichia coli detected in 10 of those 39 samples. Radon-222 activities ranged from non-detect to 7,420 picocuries per liter (pCi/L), with a median of 863 pCi/L, and exceeded the proposed drinking-water standard of 300 pCi/L in 50 (67 percent) of the 75 samples; radon-222 activities were higher than the alternative proposed standard of 4,000 pCi/L in 3 samples. Water from 15 of 75 (20 percent) wells had concentrations of methane greater than the reporting level of 0.01 mg/L; detectable methane concentrations ranged from 0.04 to 16.8 mg/L. Two samples had methane concentrations (13.1 and 16.8 mg/L) exceeding the action level of 7 mg/L. Low levels of ethane (up to 0.12 mg/L) were present in the five samples with the highest methane concentrations (near or above 1 mg/L) that were analyzed for hydrocarbon compounds and isotopic composition. The isotopic composition of methane in four of these groundwater samples, from the Catskill and Lock Haven Formations and the Hamilton Group, have sample carbon isotopic ratio delta values (carbon-13/carbon-12) ranging from –42.36 to –36.08 parts per thousand (‰) and hydrogen isotopic ratio delta values (deuterium/protium) ranging from –212.0 to –188.4 ‰, which are consistent with the isotopic compositions reported for mud-gas logging samples from these geologic units and a thermogenic source of the methane. However, the isotopic composition and ratios of methane to ethane in a fifth sample indicate the methane in that sample may be of microbial origin that subsequently underwent oxidation. The fifth sample had the highest concentration of methane, 16.8 mg/L, with an carbon isotopic ratio delta values of -50.59 ‰ and a hydrogen isotopic ratio delta values of -209.7 ‰. The six well-water samples with the highest methane concentrations also had among the highest pH values (8.25 to 9.15) and elevated concentrations of sodium, lithium, boron, fluoride, arsenic, and bromide. Relatively elevated concentrations of some other constituents, such as barium, strontium, and chloride, commonly were present in, but not limited to, those well-water samples with elevated methane. Three of the six groundwater samples with the highest methane concentrations had chloride/bromide ratios that indicate mixing with a small amount of brine (0.02 percent or less) similar in composition to those reported at undetermined depth below the freshwater aquifer and for gas and oil well brines in Pennsylvania. The sample with the highest methane concentration and most other samples with low methane concentrations (less than about 1 mg/L) have chloride/bromide ratios that indicate predominantly anthropogenic sources of chloride, such as road-deicing salt, septic systems, and (or) animal waste. Brines that are naturally present may originate from deeper parts of the aquifer system, while anthropogenic sources are more likely to affect shallow groundwater because they occur on or near the land-surface. The spatial distribution of groundwater compositions generally indicate that (1) uplands along the western border of Lycoming County usually have dilute, slightly acidic oxygenated, calcium-bicarbonate type waters; (2) intermediate altitudes or areas of carbonate bedrock usually have water of near neutral pH, with highest amounts of hardness (calcium and magnesium); (3) stream valleys, low elevations where groundwater may be discharging, and deep wells in uplands usually have water with pH values greater than 8 and highest arsenic, sodium, lithium, bromide concentrations. Geochemical modeling indicated that for samples with elevated pH, sodium, lithium, bromide, and alkalinity, the water chemistry could have resulted by dissolution of calcite (calcium carbonate) combined with cation-exchange and mixing with a small amount of brine. Through cation-exchange reactions between water and bedrock, which are equivalent to processes in a water softener, calcium ions released by calcite dissolution are exchanged for sodium ions on clay minerals. Thus, the assessment of groundwater quality in Lycoming County indicates groundwater is generally of good quality, but various parts of Lycoming County can have groundwater with low to moderate concentrations of methane and other constituents that appear in naturally present brine and produced waters from gas and oil wells at high concentrations."

Pennsylvania↗

Computing daily mean streamflow at ungaged locations in Iowa by using the Flow Anywhere and Flow Duration Curve Transfer statistical methods

The U.S. Geological Survey (USGS) maintains approximately 148 real-time streamgages in Iowa for which daily mean streamflow information is available, but daily mean streamflow data commonly are needed at locations where no streamgages are present. Therefore, the USGS conducted a study as part of a larger project in cooperation with the Iowa Department of Natural Resources to develop methods to estimate daily mean streamflow at locations in ungaged watersheds in Iowa by using two regression-based statistical methods. The regression equations for the statistical methods were developed from historical daily mean streamflow and basin characteristics from streamgages within the study area, which includes the entire State of Iowa and adjacent areas within a 50-mile buffer of Iowa in neighboring states. Results of this study can be used with other techniques to determine the best method for application in Iowa and can be used to produce a Web-based geographic information system tool to compute streamflow estimates automatically. The Flow Anywhere statistical method is a variation of the drainage-area-ratio method, which transfers same-day streamflow information from a reference streamgage to another location by using the daily mean streamflow at the reference streamgage and the drainage-area ratio of the two locations. The Flow Anywhere method modifies the drainage-area-ratio method in order to regionalize the equations for Iowa and determine the best reference streamgage from which to transfer same-day streamflow information to an ungaged location. Data used for the Flow Anywhere method were retrieved for 123 continuous-record streamgages located in Iowa and within a 50-mile buffer of Iowa. The final regression equations were computed by using either left-censored regression techniques with a low limit threshold set at 0.1 cubic feet per second (ft3/s) and the daily mean streamflow for the 15th day of every other month, or by using an ordinary-least-squares multiple linear regression method and the daily mean streamflow for the 15th day of every other month. The Flow Duration Curve Transfer method was used to estimate unregulated daily mean streamflow from the physical and climatic characteristics of gaged basins. For the Flow Duration Curve Transfer method, daily mean streamflow quantiles at the ungaged site were estimated with the parameter-based regression model, which results in a continuous daily flow-duration curve (the relation between exceedance probability and streamflow for each day of observed streamflow) at the ungaged site. By the use of a reference streamgage, the Flow Duration Curve Transfer is converted to a time series. Data used in the Flow Duration Curve Transfer method were retrieved for 113 continuous-record streamgages in Iowa and within a 50-mile buffer of Iowa. The final statewide regression equations for Iowa were computed by using a weighted-least-squares multiple linear regression method and were computed for the 0.01-, 0.05-, 0.10-, 0.15-, 0.20-, 0.30-, 0.40-, 0.50-, 0.60-, 0.70-, 0.80-, 0.85-, 0.90-, and 0.95-exceedance probability statistics determined from the daily mean streamflow with a reporting limit set at 0.1 ft 3 /s. The final statewide regression equation for Iowa computed by using left-censored regression techniques was computed for the 0.99-exceedance probability statistic determined from the daily mean streamflow with a low limit threshold and a reporting limit set at 0.1 ft 3 /s. For the Flow Anywhere method, results of the validation study conducted by using six streamgages show that differences between the root-mean-square error and the mean absolute error ranged from 1,016 to 138 ft 3 /s, with the larger value signifying a greater occurrence of outliers between observed and estimated streamflows. Root-mean-square-error values ranged from 1,690 to 237 ft 3 /s. Values of the percent root-mean-square error ranged from 115 percent to 26.2 percent. The logarithm (base 10) streamflow percent root-mean-square error ranged from 13.0 to 5.3 percent. Root-mean-square-error observations standard-deviation-ratio values ranged from 0.80 to 0.40. Percent-bias values ranged from 25.4 to 4.0 percent. Untransformed streamflow Nash-Sutcliffe efficiency values ranged from 0.84 to 0.35. The logarithm (base 10) streamflow Nash-Sutcliffe efficiency values ranged from 0.86 to 0.56. For the streamgage with the best agreement between observed and estimated streamflow, higher streamflows appear to be underestimated. For the streamgage with the worst agreement between observed and estimated streamflow, low flows appear to be overestimated whereas higher flows seem to be underestimated. Estimated cumulative streamflows for the period October 1, 2004, to September 30, 2009, are underestimated by -25.8 and -7.4 percent for the closest and poorest comparisons, respectively. For the Flow Duration Curve Transfer method, results of the validation study conducted by using the same six streamgages show that differences between the root-mean-square error and the mean absolute error ranged from 437 to 93.9 ft 3 /s, with the larger value signifying a greater occurrence of outliers between observed and estimated streamflows. Root-mean-square-error values ranged from 906 to 169 ft 3 /s. Values of the percent root-mean-square-error ranged from 67.0 to 25.6 percent. The logarithm (base 10) streamflow percent root-mean-square error ranged from 12.5 to 4.4 percent. Root-mean-square-error observations standard-deviation-ratio values ranged from 0.79 to 0.40. Percent-bias values ranged from 22.7 to 0.94 percent. Untransformed streamflow Nash-Sutcliffe efficiency values ranged from 0.84 to 0.38. The logarithm (base 10) streamflow Nash-Sutcliffe efficiency values ranged from 0.89 to 0.48. For the streamgage with the closest agreement between observed and estimated streamflow, there is relatively good agreement between observed and estimated streamflows. For the streamgage with the poorest agreement between observed and estimated streamflow, streamflows appear to be substantially underestimated for much of the time period. Estimated cumulative streamflow for the period October 1, 2004, to September 30, 2009, are underestimated by -9.3 and -22.7 percent for the closest and poorest comparisons, respectively.

Illinois;Iowa;Minnesota;Missouri;Nebraska;Wisconsi↗

Status and understanding of groundwater quality in the Bear Valley and Lake Arrowhead Watershed Study Unit, 2010: California GAMA Priority Basin Project

Groundwater quality in the 112-square-mile Bear Valley and Lake Arrowhead Watershed (BEAR) study unit was investigated as part of the Priority Basin Project (PBP) of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit comprises two study areas (Bear Valley and Lake Arrowhead Watershed) in southern California in San Bernardino County. The GAMA-PBP is conducted by the California State Water Resources Control Board (SWRCB) in cooperation with the U.S. Geological Survey (USGS) and the Lawrence Livermore National Laboratory. The GAMA BEAR study was designed to provide a spatially balanced, robust assessment of the quality of untreated (raw) groundwater from the primary aquifer systems in the two study areas of the BEAR study unit. The assessment is based on water-quality collected by the USGS from 38 sites (27 grid and 11 understanding) during 2010 and on water-quality data from the SWRCB-Division of Drinking Water (DDW) database. The primary aquifer system is defined by springs and the perforation intervals of wells listed in the SWRCB-DDW water-quality database for the BEAR study unit. This study included two types of assessments: (1) a status assessment , which characterized the status of the quality of the groundwater resource as of 2010 by using data from samples analyzed for volatile organic compounds, pesticides, and naturally present inorganic constituents, such as major ions and trace elements, and (2) an understanding assessment , which evaluated the natural and human factors potentially affecting the groundwater quality. The assessments were intended to characterize the quality of groundwater resources in the primary aquifer system of the BEAR study unit, not the treated drinking water delivered to consumers. Bear Valley study area and the Lake Arrowhead Watershed study area were also compared statistically on the basis of water-quality results and factors potentially affecting the groundwater quality. Relative concentrations (RCs), which are sample concentration of a particular constituent divided by its associated health- or aesthetic-based benchmark concentrations, were used for evaluating the groundwater quality for those constituents that have Federal or California regulatory or non-regulatory benchmarks for drinking-water quality. An RC greater than 1.0 indicates a concentration greater than a benchmark. Organic (volatile organic compounds and pesticides) and special-interest (perchlorate) constituent RCs were classified as “high” (RC greater than 1.0), “moderate” (RC less than or equal to 1.0 and greater than 0.1), or “low” (RC less than or equal to 0.1). For inorganic (radioactive, trace element, major ion, and nutrient) constituents, the boundary between low and moderate RCs was set at 0.5. Aquifer-scale proportion was used as the primary metric in the status assessment for evaluating groundwater quality at the study-unit scale or for its component areas. High aquifer-scale proportion was defined as the percentage of the area of the primary aquifer system with a RC greater than 1.0 for a particular constituent or class of constituents; the percentage is based on area rather than volume. Moderate and low aquifer-scale proportions were defined as the percentage of the primary aquifer system with moderate and low RCs, respectively. A spatially weighted statistical approach was used to evaluate aquifer-scale proportions for individual constituents and classes of constituents. The status assessment for the Bear Valley study area found that inorganic constituents with health-based benchmarks were detected at high RCs in 9.0 percent of the primary aquifer system and at moderate RCs in 13 percent. The high RCs of inorganic constituents primarily reflected high aquifer-scale proportions of fluoride (in 5.4 percent of the primary aquifer system) and arsenic (3.6 percent). The RCs of organic constituents with health-based benchmarks were high in 1.0 percent of the primary aquifer system, moderate in 8.1 percent, and low in 70 percent. Organic constituents were detected in 79 percent of the primary aquifer system. Two groups of organic constituents and two individual organic constituents were detected at frequencies greater than 10 percent of samples from the USGS grid sites: trihalomethanes (THMs), solvents, methyl tert -butyl ether (MTBE), and simazine. The special-interest constituent perchlorate was detected in 93 percent of the primary aquifer system; it was detected at moderate RCs in 7.1 percent and at low RCs in 86 percent. The status assessment in the Lake Arrowhead Watershed study area showed that inorganic constituents with human-health benchmarks were detected at high RCs in 25 percent of the primary aquifer system and at moderate RCs in 41 percent. The high aquifer-scale proportion of inorganic constituents primarily reflected high aquifer-scale proportions of radon‑222 (in 62 percent of the primary aquifer system) and uranium (26 percent). RCs of organic constituents with health-based benchmarks were moderate in 7.7 percent of the primary aquifer system and low in 46 percent. Organic constituents were detected in 54 percent of the primary aquifer system. The only organic constituents that were detected at frequencies greater than 10 percent of samples from the USGS grid sites were THMs. Perchlorate was detected in 62 percent of the primary aquifer system at uniformly low RCs. The second component of this study, the understanding assessment , identified the natural and human factors that could have affected the groundwater quality in the BEAR study unit by evaluating statistical correlations between water-quality constituents and potential explanatory factors. The potential explanatory factors evaluated were land use (including density of septic tanks and leaking or formerly leaking underground fuel tanks), site type, aquifer lithology, well construction (well depth and depth to the top-of-perforated interval), elevation, aridity index, groundwater-age distribution, and oxidation-reduction condition (including pH and dissolved oxygen concentration). Results of the statistical evaluations were used to explain the distribution of constituents in groundwater of the BEAR study unit. In the Bear Valley study area, high and moderate RCs of fluoride were found in sites known to be influenced by hydrothermic conditions or that had high concentrations of fluoride historically. The high RC of arsenic can likely be attributed to desorption of arsenic from aquifer sediments saturated in old groundwater with high pH under reducing conditions. The THMs were detected more frequently at USGS grid sites that were wells, part of a large urban water system, and surrounded by urban land use. Solvents, MTBE, and simazine were all detected more frequently at USGS grid sites that were wells with a greater urban percentage of surrounding land use and that accessed older groundwater than other USGS grid sites. Comparison between the observed and predicted detection frequencies of perchlorate at USGS grid sites indicated that anthropogenic sources could have contributed to low levels of perchlorate in the groundwater of the Bear Valley study area. In the Lake Arrowhead Watershed study area, high and moderate RCs of radon-222 and uranium can be attributed to older groundwater from the granitic fractured-rock primary aquifer system. Low RCs of THMs were detected at USGS grid sites that were wells and part of small water systems. The similarities between the observed and predicted detection frequencies of perchlorate in samples from USGS grid sites indicated that the source and distribution of perchlorate were most likely attributable to precipitation (rain and snow), with minimal, if any, contribution from anthropogenic sources.

California↗

Water movement through the unsaturated zone of the High Plains Aquifer in the Central Platte Natural Resources District, Nebraska, 2008-12

Uncertainty about the effects of land use and climate on water movement in the unsaturated zone and on groundwater recharge rates can lead to uncertainty in water budgets used for groundwater-flow models. To better understand these effects, a cooperative study between the U.S. Geological Survey and the Central Platte Natural Resources District was initiated in 2007 to determine field-based estimates of recharge rates in selected land-use areas of the Central Platte Natural Resources District in Nebraska. Measured total water potential and unsaturated-zone profiles of tritium, chloride, nitrate as nitrogen, and bromide, along with groundwater-age dates, were used to evaluate water movement in the unsaturated zone and groundwater recharge rates in the central Platte River study area. Eight study sites represented an east-west precipitation contrast across the study area—four beneath groundwater-irrigated cropland (sites 2, 5, and 6 were irrigated corn and site 7 was irrigated alfalfa/corn rotation), three beneath rangeland (sites 1, 4, and 8), and one beneath nonirrigated cropland, or dryland (site 3). Measurements of transient vertical gradients in total water potential indicated that periodic wetting fronts reached greater mean maximum depths beneath the irrigated sites than beneath the rangeland sites, in part, because of the presence of greater and constant antecedent moisture. Beneath the rangeland sites, greater temporal variation in antecedent moisture and total water potential existed and was, in part, likely a result of local precipitation and evapotranspiration. Moreover, greater variability was noticed in the total water potential profiles beneath the western sites than the corresponding eastern sites, which was attributed to less mean annual precipitation in the west. The depth of the peak post-bomb tritium concentration or the interface between the pre-bomb/post-bomb tritium, along with a tritium mass balance, within sampled soil profiles were used to estimate water fluxes in the unsaturated zone at three of the eight study sites: site 2 (irrigated), site 3 (dryland), and site 8 (rangeland). Estimates for recharge were about 68 millimeters per year [(mm/yr), post-bomb peak], 133 to 159 mm/yr (tritium interface), and 137 mm/yr (mass balance) at site 2 (irrigated); about 63 mm/yr (tritium interface) and 12 mm/yr (mass balance) at site 3 (dryland); and about 53 mm/yr (tritium interface) and 10 mm/yr (mass balance) at site 8 (rangeland). Recharge values from the mass balance at site 2 were more than an order of magnitude greater than recharge values at site 3, suggesting irrigation is an important control on water movement through the unsaturated zone. For the remaining five sites, the post-bomb tritium had flushed through the system and recharge was considered modern (within 10 years of sampling). The chloride mass-balance method was used to determine water fluxes below the root zone (less than 2 meters below land surface) at the rangeland sites: sites 1, 4, and 8. At these rangeland sites, water fluxes ranged from 1.8 to 96 mm/yr at site 1, 1.1 to 9.6 mm/yr at site 4, and 1.1 to 68 mm/yr at site 8, with mean rates of 21, 4.3, and 13 mm/yr, respectively. Site 1 had a greater mean water flux, which was consistent with the greater precipitation in the east than at site 8 in the west. Chloride mass balance was not calculated at the irrigated and dryland sites because of uncertainty about additional sources of chloride. Concentrations of nitrate as nitrogen in pore water in the unsaturated zone were larger beneath the irrigated and dryland (agricultural) sites compared with the rangeland sites. The larger concentrations at the agricultural sites are consistent with the application of nitrogen fertilizer at the agricultural sites and no substantial accumulation at the rangeland sites. The shape of the nitrate as nitrogen and chloride concentration profiles at site 1 (rangeland) indicate a reasonably larger and more consistent water flux in the UZ than beneath the other two rangeland sites (sites 4 or 8). Excluding site 7, the general shape of the nitrate as nitrogen profiles was similar beneath the agricultural sites and supports the estimates of water movement and recharge rates determined from the tritium and chloride methods. Movement of bromide through the unsaturated zone indicated greater water fluxes are found beneath irrigated lands than beneath rangeland. Bromide profiles in the unsaturated zone, determined from center of mass and peak displacement methods, document water fluxes ranged from 58 to 394 mm/yr beneath irrigated sites and 9 to 201 mm/yr beneath rangeland sites. Water-flux estimates from the potassium bromide tests at most sites did not represent overall recharge rates because the bromide remained primarily in the root zone. Apparent groundwater age was used to determine the groundwater residence time at the eight sites and to estimate recharge rates. Groundwater ages in the study area ranged from old water (defined here as groundwater that was recharged more than 50 years ago) to modern (defined here as groundwater that has recharged within the past 10 years). Groundwater ages indicated that the shallow monitoring wells generally had younger residence times, whereas the deeper monitoring wells generally had the older residence times. Groundwater dates from the shallowest monitoring wells were used to determine recharge rates at the water table. These rates generally were similar to recharge rates determined from tritium and chloride mass-balance methods. Groundwater recharge rates generally increased with well depth, and the deeper monitoring wells likely do not represent local recharge conditions but recharge from a regional flow system that receives recharge from distant sources. Overall, these data generally indicate that water movement within the unsaturated zone primarily is affected by spatial contrasts in mean annual precipitation and by the land use or land cover. The eight unsaturated-zone sites each generated unique, valuable datasets that likely will improve the understanding of water movement and recharge rates in the central Platte River valley.

Nebraska↗

Response of Green Lake, Wisconsin, to changes in phosphorus loading, with special emphasis on near-surface total phosphorus concentrations and metalimnetic dissolved oxygen minima

Green Lake is the deepest natural inland lake in Wisconsin, with a maximum depth of about 72 meters. In the early 1900s, the lake was believed to have very good water quality (low nutrient concentrations and good water clarity) with low dissolved oxygen (DO) concentrations occurring in only the deepest part of the lake. Because of increased phosphorus (P) inputs from anthropogenic activities in its watershed, total phosphorus (TP) concentrations in the lake have increased; these changes have led to increased algal production and low DO concentrations not only in the deepest areas but also in the middle of the water column (metalimnion). The U.S. Geological Survey has routinely monitored the lake since 2004 and its tributaries since 1988. Results from this monitoring led the Wisconsin Department of Natural Resources (WDNR) to list the lake as impaired because of low DO concentrations in the metalimnion, and they identified elevated TP concentrations as the cause of impairment. As part of this study by the U.S. Geological Survey, in cooperation with the Green Lake Sanitary District, the lake and its tributaries were comprehensively sampled in 2017–18 to augment ongoing monitoring that would further describe the low DO concentrations in the lake (especially in the metalimnion). Empirical and process-driven water-quality models were then used to determine the causes of the low DO concentrations and the magnitudes of P-load reductions needed to improve the water quality of the lake enough to meet multiple water-quality goals, including the WDNR’s criteria for TP and DO. Data from previous studies showed that DO concentrations in the metalimnion decreased slightly as summer progressed in the early 1900s but, since the late 1970s, have typically dropped below 5 milligrams per liter (mg/L), which is the WDNR criterion for impairment. During 2014–18 (the baseline period for this study), the near-surface geometric mean TP concentration during June–September in the east side of the lake was 0.020 mg/L and in the west side was 0.016 mg/L (both were above the 0.015-mg/L WDNR criterion for the lake), and the metalimnetic DO minimum concentrations (MOMs) measured in August ranged from 1.0 to 4.7 mg/L. The degradation in water quality was assumed to have been caused by excessive P inputs to the lake; therefore, the TP inputs to the lake were estimated. The mean annual external P load during 2014–18 was estimated to be 8,980 kilograms per year (kg/yr), of which monitored and unmonitored tributary inputs contributed 84 percent, atmospheric inputs contributed 8 percent, waterfowl contributed 7 percent, and septic systems contributed 1 percent. During fall turnover, internal sediment recycling contributed an additional 7,040 kilograms that increased TP concentrations in shallow areas of the lake by about 0.020 mg/L. The elevated TP concentrations then persisted until the following spring. On an annual basis, however, there was a net deposition of P to the bottom sediments. Empirical models were used to describe how the near-surface water quality of Green Lake would be expected to respond to changes in external P loading. Predictions from the models showed a relatively linear response between P loading and TP and chlorophyll- a (Chl- a ) concentrations in the lake, with the changes in TP and Chl- a concentrations being less on a percentage basis (50–60 percent for TP and 30–70 percent for Chl- a ) than the changes in P loading. Mean summer water clarity, quantified by Secchi disk depths, had a greater response to decreases in P loading than to increases in P loading. Based on these relations, external P loading to the lake would need to be decreased from 8,980 kg/yr to about 5,460 kg/yr for the geometric mean June–September TP concentration in the east side of the lake, with higher TP concentrations than in the west side, to reach the WDNR criterion of 0.015 mg/L. This reduction of 3,520 kg/yr is equivalent to a 46-percent reduction in the potentially controllable external P sources (all external sources except for precipitation, atmospheric deposition, and waterfowl) from those measured during water years 2014–18. The total external P loading would need to decrease to 7,680 kg/yr (a 17-percent reduction in potentially controllable external P sources) for near-surface June–September TP concentrations in the west side of the lake to reach 0.015 mg/L. Total external P loading would need to decrease to 3,870–5,320 kg/yr for the lake to be classified as oligotrophic, with a near-surface June–September TP concentration of 0.012 mg/L. Results from the hydrodynamic water-quality model GLM–AED (General Lake Model coupled to the Aquatic Ecodynamics modeling library) indicated that MOMs are driven by external P loading and internal sediment recycling that lead to high TP concentrations during spring and early summer, which in turn lead to high phytoplankton production, high metabolism and respiration, and ultimately DO consumption in the upper, warmer areas of the metalimnion. GLM–AED results indicated that settling of organic material during summer might be slowed by the colder, denser, and more viscous water in the metalimnion and thus increase DO consumption. Based on empirical evidence from a comparison of MOMs with various meteorological, hydrologic, water quality, and in-lake physical factors, MOMs were lower during summers, when metalimnetic water temperatures were warmer, near-surface Chl- a and TP concentrations were higher, and Secchi depths were lower. GLM–AED results indicated that the external P load would need to be reduced to about 4,060 kg/yr, a 57-percent reduction from that measured in 2014–18, to eliminate the occurrence of MOMs less than 5 mg/L during more than 75 percent of the years (the target provided by the WDNR). Large reductions in external P loading are expected to have an immediate effect on the near-surface TP concentrations and metalimnetic DO concentrations in Green Lake; however, it may take several years for the full effects of the external-load reduction to be observed because internal sediment recycling is an important source of P for the following spring.

Wisconsin↗

Altitude, age, and quality of groundwater, Papio-Missouri River Natural Resources District, eastern Nebraska, 1992 to 2009

The U.S. Geological Survey, in cooperation with the Papio-Missouri River Natural Resources District (PMRNRD), conducted this study to map the water-level altitude of 2009 within the Elkhorn River Valley, Missouri River Valley, and Platte River Valley alluvial aquifers; to present the predevelopment potentiometric-surface altitude within the Dakota aquifer; and to describe the age and quality of groundwater in the five principal aquifers of the PMRNRD in eastern Nebraska using data collected from 1992 to 2009. In addition, implications of alternatives to the current PMRNRD groundwater-quality monitoring approach are discussed. In the PMRNRD, groundwater altitude, relative to National Geodetic Vertical Datum of 1929, ranged from about 1,080 feet (ft) to 1,180 ft in the Elkhorn River Valley alluvial aquifer and from about 960 ft to 1,080 ft in the Missouri River Valley and Platte River Valley alluvial aquifers. In the PMRNRD, the estimated altitude of the potentiometric surface of the Dakota aquifer, predevelopment, ranged from about 1,100 ft to 1,200 ft. To assess groundwater age and quality, groundwater samples were collected from a total of 217 wells from 1992 to 2009 for analysis of various analytes. Groundwater samples collected in the PMRNRD from 1992 to 2009 and interpreted in this report were analyzed for age-dating analytes (chlorofluorocarbons), dissolved gases, major ions, trace elements, nutrients, stable isotope ratios, pesticides and pesticide degradates, volatile organic compounds, explosives, and 222radon. Apparent groundwater age was estimated from concentrations of chlorofluorocarbons measured in samples collected in 2000. Apparent groundwater-recharge dates ranged from older than 1940 in samples from wells screened in the Missouri River Valley alluvial aquifer to the early 1980s in samples from wells screened in the Dakota aquifer. Concentrations of major ions in the most recent sample per well collected from 1992 to 2009 indicate that the predominant water type was calcium bicarbonate. Samples from 4 wells exceeded the U.S. Environmental Protection Agency (USEPA) Secondary Drinking Water Regulation (SDWR) for sulfate [250 milligrams per liter (mg/L)], and samples from 4 wells exceeded the USEPA Drinking Water Advisory Table for sodium (30-60 mg/L). Eighteen of the 21 trace elements analyzed in samples from PMRNRD wells have USEPA drinking-water standards. Sixteen of the trace elements with USEPA standards were detected in the selected samples. In the samples selected for trace-element analysis, the only trace-element concentration that exceeded an enforceable USEPA drinking-water standard, the Maximum Contaminant Level (MCL), was for arsenic; arsenic concentration exceeded the USEPA MCL of 10 micrograms per liter (μg/L) in 4 percent of the samples. Trace-element concentrations that exceeded the USEPA SDWR or Lifetime Health Advisory level were iron (46 percent of the samples were greater than USEPA SDWR of 300 μg/L), manganese (70 percent of the samples were greater than USEPA SDWR of 50 μg/L), and strontium (4 percent of the samples were greater than USEPA Lifetime Health Advisory level of 4,000 μ/L). The concentration of nitrate plus nitrite as nitrogen (nitrate-N) in the most recent nutrient samples collected from the network wells and from one randomly selected well in the well nests from 1992 to 2009 for most wells (80 percent) ranged from less than 0.06 to 8.55 mg/L, with a median value of 0.12 mg/L. Concentrations of nitrate-N in 13 (7 percent) nutrient samples, 1992 to 2009, were greater than or equal to the USEPA MCL and Nebraska Title-118 standard of 10 mg/L, and concentrations of nitrate-N in 35 (18 percent) nutrient samples, 1992 to 2009, were greater than or equal to 5 mg/L, which is the PMRNRD action level for possible management implementation to reduce nitrate concentrations in groundwater. Of the 61 pesticides or pesticide degradates analyzed from 2007 to 2009, 21 were detected. Three of the 21 pesticides detected (alachlor, atrazine, and metolachlor) have established health-based criteria; all detections of these compounds were at concentrations less than their USEPA standards. From 2007 to 2009, 1 or more pesticide compounds were detected in 16 of the 82 network wells and in 18 of the 26 wells in well nests. From 2007 to 2009, the individual pesticide compounds that were detected most frequently were alachlor ethane sulfonic acid, a degradate of alachlor; deethylcyanazine acid, a degradate of cyanazine; and atrazine. Analytes with concentrations that exceeded 30 percent of the applicable Nebraska Title-118 standard were identified so that the PMRNRD can plan to monitor groundwater in the area and consider possible actions should the analyte concentrations continue to rise. The analytical results from the most recent samples collected in the network wells and all the wells in well nests from 1992 to 2009 indicate that, in at least 1 sample, there was a concentration that exceeded 30 percent of the Nebraska Title-118 standard for at least 1 of 3 major ions (chloride, fluoride, and sulfate), 1 nutrient (nitrate-N), 1 pesticide (atrazine), or 3 trace elements (arsenic, iron, and manganese). In addition, 30 percent of the USEPA MCL or Nebraska Title-118 standard for gross alpha activity likely was exceeded in samples from three wells screened in the Dakota aquifer. Study findings indicate that some alternatives to the current PMRNRD groundwater-sampling approach that could be considered are to collect fewer samples for nutrient analysis and to collect samples periodically for determining concentrations of additional analytes, particularly the analytes with concentrations that were at least 30 percent or more than the Nebraska Title-118 standard.

Nebraska↗

Trends in streamflow and concentrations and flux of nutrients and total suspended solids in the Upper White River at Muncie, near Nora, and near Centerton, Indiana

The U.S. Geological Survey (USGS), in cooperation with The Nature Conservancy, completed a study to estimate and assess trends in streamflow and annual mean concentrations and flux of nutrients (nitrate plus nitrite, total Kjeldahl nitrogen, and total phosphorus) and total suspended solids at three USGS streamgages (hereafter referred to as “study gages”) on the Upper White River at Muncie (USGS station 03347000), near Nora (USGS station 03351000), and near Centerton (USGS station 03354000), Indiana. Water-quality data used in the analyses were collected by several agencies between calendar years 1991 and 2017, and streamflow (discharge) data were collected by the USGS. For most of the water-quality constituents, there were suitable data to facilitate an analysis of the 26-year period extending from calendar years 1991 to 2017 (water years 1992 to 2017); however, shorter analytical periods were necessary for total Kjeldahl nitrogen for the study gages at Muncie and near Centerton and for total suspended solids for the study gage near Centerton. Temporal trends in streamflows at the study gages for the period extending from water years 1978 to 2017 were assessed using Exploration and Graphics for RivEr Trends (EGRET) and Mann-Kendall and Pettitt tests. With just one exception, the annual maximum and mean daily streamflows and the annual minimum 7-day mean streamflows at the study gages demonstrated upward trends (increasing streamflows) in the EGRET analyses. The exception was the annual 7-day minimum streamflow at the study gage near Nora, which indicated no trend. Mann-Kendall tests also indicated that the average trend for the annual maximum daily, annual mean daily, and annual 7-day minimum streamflow statistics between water years 1978 and 2017 was upward at each of the study gages; however, only the trends in the annual mean daily streamflows at the study gage at Muncie and the annual maximum daily streamflows at the study gages near Nora and near Centerton were statistically significant at a 0.05 probability level. The Pettitt tests indicated that a statistically significant step trend (abrupt change) in annual mean daily streamflows occurred at each of the study gages around water year 2001. The seasonal distributions of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen concentrations at the study gages were evaluated to identify patterns and other distinguishing characteristics by examining boxplots of concentrations as a function of month of the year. Seasonal distributions of nitrate plus nitrite concentrations and total suspended solids concentrations differed from each other but were generally similar among the three study gages for a given constituent. Median concentrations of nitrate plus nitrite were highest during the January–June months, whereas median concentrations of total suspended solids were highest during June and July. Seasonal distributions of total phosphorus concentrations were similar at the study gages near Nora and near Centerton, but the seasonal distribution was noticeably different at the study gage at Muncie, which had monthly median concentrations that were substantially lower than at the two downstream study gages (near Nora and near Centerton). The seasonal distribution of total Kjeldahl nitrogen concentrations differed in pattern among the three study gages; however, in general, some of the higher monthly median total Kjeldahl nitrogen concentrations at each study gage were associated with the late spring and summer periods. The Weighted Regressions on Time, Discharge, and Season (WRTDS) method implemented in EGRET was used to estimate water-year annual mean daily concentrations and flux of nutrients and total suspended solids, as well as estimates of concentrations and flux that were “normalized” to remove the effect of year-to-year variation in streamflow. The approximate coefficients of determination for the WRTDS regression models ranged from a high of 0.82 for total phosphorus for the study gage near Centerton to a low of 0.19 for nitrate plus nitrite for the study gage near Nora. Loads and yields of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen were estimated for analytical periods consisting of the longest periods of concurrent record at the three study gages. Loads of each of the constituents increased sequentially from the most upstream study gage to the most downstream study gage; however, the same was not true for yields. The highest yields of total suspended solids, total phosphorus, and total Kjeldahl nitrogen occurred at the most upstream study gage (at Muncie); however, the highest yield of nitrate plus nitrite occurred at the most downstream study gage (near Centerton). WRTDS bootstrap tests were used to assess the magnitude, direction, and likelihood of changes in annual flow-normalized mean daily concentrations and flux of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen at the study gages between water years 1997 and 2017. Changes in flow-normalized concentrations and flux of the constituents between water years 1997 and 2017 were mostly downward (decreasing). The exceptions were likely to highly likely upward (increasing) changes in (1) flow-normalized annual mean daily concentration and annual flux for total suspended solids and total phosphorus at the study gage at Muncie, (2) flow-normalized annual mean daily total phosphorus concentration at the study gage near Centerton, (3) flow-normalized annual flux of total phosphorus at the study gage near Centerton, and (4) flow-normalized annual mean daily nitrate plus nitrite concentration at the study gage near Centerton. Although an upward change in flow-normalized nitrate plus nitrite concentrations was likely at the study gage near Centerton, flow-normalized annual flux of nitrate plus nitrite at that study gage was determined to have a highly likely downward change. EGRET and Exploration and Graphics for RivEr Trends Confidence Intervals (EGRETci) analyses can be used to improve our understanding of how concentrations and flux change as functions of time and streamflow, as well as provide information on how the relations between streamflow and constituent concentrations have changed within the calendar year between any 2 years included in the analyses. Examples of those uses, illustrating changes between calendar years 1992 and 2017, were given for total suspended solids concentrations at the study gage near Nora and for nitrate plus nitrite concentrations at the study gage near Centerton.

Indiana↗

Evaluation of salinity and nutrient conditions in the Heart River Basin, North Dakota, 1970–2020

The Heart River Basin is predominantly an agricultural basin in western North Dakota and is approximately 3,350 square miles. The U.S. Geological Survey, in cooperation with the U.S. Department of Agriculture Natural Resources Conservation Service and the Grant County Soil Conservation District, completed a study to assess spatial and temporal patterns of water quality in the Heart River Basin. The purpose of this report is to describe the methods and results of a study to evaluate salinity and nutrients in the Heart River Basin in western North Dakota. Water-quality and streamflow data used in the study were compiled from 1970 to 2020 using the National Water Quality Monitoring Council Water Quality Portal and National Water Information System. Changes in streamflow characteristics were investigated at three sites from 1970 to 2020, and changes in water quality were investigated at four sites from 1974 to 2019. Streamflow analysis indicated decreasing streamflow from 1970 until the late 1990s followed by increasing streamflow through 2020, with the largest increase in the 7-day minimum streamflow or base flow. For the historical water-quality trend period (1974–2019), total dissolved solids, sulfate, sodium, chloride, and sodium adsorption ratio concentrations have increased since the mid-1970s through 2019. Potassium concentrations during the historical period remained mostly constant with some small fluctuations. Calcium and magnesium concentrations increased since the mid-1970s at all sites, except for a decrease at one site between 1974 and 1999. During the recent trend period (1999–2019), increasing concentrations in total dissolved solids, sulfate, sodium, chloride, calcium, magnesium, and sodium adsorption ratios were observed across the Heart River Basin. The magnitude of the increases was smaller at tributary sites compared to main-stem sites. During the recent period, potassium was mostly constant, although small (−0.9 milligram per liter or less) decreases on tributaries and minor (1.3 milligrams per liter) increases on the main-stem sites were detected. Unlike dissolved ion concentrations, significant increases in nutrient concentrations were not detected from 1999 to 2019, but nitrate plus nitrite concentrations most likely decreased upstream from Lake Tschida. Inverse modeling for period 1 (1974–99) in model zone 1 (Heart River reach from site 5 to site 6) had eight reasonable models that indicated the clay mineral-water interactions and dissolution of evaporites control the geochemistry. Results of the inverse modeling for period 2 (1999–2019) in model zone 1 also had eight reasonable models that indicated that the dissolution of evaporites was the major geochemical control. Results of the geochemical modeling for period 1 (1974–99) in model zone 2 (Heart River and Sweetbriar Creek reach from sites 20 and 21 to site 22) produced seven reasonable models, and the geochemical control of the system was the dissolution of sulfate evaporite minerals. Geochemical modeling results for period 2 (1999–2019) in model zone 2 produced 11 reasonable models and was also controlled by the dissolution of sulfate evaporite minerals. Differences between the two model zones indicated that geology controls some of the water-quality changes in the Heart River Basin. Loads were estimated for total dissolved solids, sulfate, sodium, and chloride and total phosphorus. Annual loads estimated for the Heart River from 2013 through 2020 at the Heart River site upstream from Lake Tschida (site 5) and near Mandan (site 22) were generally greatest in 2014 and least in 2016 for total dissolved solids, sulfate, sodium, and chloride. Most of the annual loads of total dissolved solids, sulfate, sodium, and chloride are delivered in March through July in the Heart River at these sites and are likely from snowmelt and spring and summer rains. The mean annual yields of total dissolved solids and sodium from 2013 to 2020 generally were largest in Big Muddy Creek (site 18), whereas yields of sulfate and chloride were largest at Sweetbriar Creek (site 21) compared to the other selected sites in the Heart River Basin. Larger yields of total dissolved solids, sulfate, sodium, and chloride at sites located on Big Muddy Creek and Sweet Briar Creek in the lower Heart River Basin were likely a result of differences in geology and soils upstream from the selected sites. A mass balance of total dissolved solids, sulfate, sodium, and chloride was estimated for the lower Heart River Basin, specifically the reach below Lake Tschida to Mandan (site 7 to site 22). Intervening flow was the largest contributor to the dissolved ion loads in the lower Heart River Basin and is an important part of understanding the transport of dissolved ions in the basin. The intervening load can include groundwater discharge, irrigation return flow, local runoff, and input from smaller ephemeral tributaries. Tributaries in the lower Heart River Basin contributed portions of the total dissolved solids, sulfate, sodium, and chloride loads at the Heart River near Mandan (site 22) that generally were proportional to the streamflow contributions. Annual loads for total phosphorus between 2013 and 2020 at the Heart River site upstream from Lake Tschida (site 5) and near Mandan (site 22) generally were largest in 2019 and smallest in 2016. Most of the total phosphorus loads for main-stem sites 5 and 22 were transported in March, April, and June, likely from snowmelt and early summer rains. The mean annual yields of total phosphorus for 2013–20 were largest on the main-stem site upstream from Lake Tschida (site 5) and Sweetbriar Creek (site 21), whereas the smallest yields were in Big Muddy Creek (site 18). Much of the phosphorus that enters Lake Tschida from the upper basin does not get transported downstream to the lower basin, and much of the phosphorus in the lower basin was attributed to intervening flow.

North Dakota↗

Quality of water from crystalline rock aquifers in New England, New Jersey, and New York, 1995-2007

Crystalline bedrock aquifers in New England and parts of New Jersey and New York (NECR aquifers) are a major source of drinking water. Because the quality of water in these aquifers is highly variable, the U.S. Geological Survey (USGS) statistically analyzed chemical data on samples of untreated groundwater collected from 117 domestic bedrock wells in New England, New York, and New Jersey, and from 4,775 public-supply bedrock wells in New England to characterize the quality of the groundwater. The domestic-well data were from samples collected by the USGS National Water-Quality Assessment (NAWQA) Program from 1995 through 2007. The public-supply-well data were from samples collected for the U.S. Environmental Protection Agency (USEPA) Safe Drinking Water Act (SDWA) Program from 1997 through 2007. Chemical data compiled from the domestic wells include pH, specific conductance, dissolved oxygen, alkalinity, and turbidity; 6 nitrogen and phosphorus compounds, 14 major ions, 23 trace elements, 222 radon gas (radon), 48 pesticide compounds, and 82 volatile organic compounds (VOCs). Additional samples were collected from the domestic wells for the analysis of gross alpha- and gross beta-particle radioactivity, radium isotopes, chlorofluorocarbon isotopes, and the dissolved gases methane, carbon dioxide, nitrogen, and argon. Chemical data compiled from the public-supply wells include pH, specific conductance, nitrate, iron, manganese, sodium, chloride, fluoride, arsenic, uranium, radon, combined radium ( 226 radium plus 228 radium), gross alpha-particle radioactivity, and methyl tert -butyl ether (M t BE). Patterns in fluoride, arsenic, uranium, and radon distributions were discernable when the data were compared to lithology groupings of the bedrock, indicating that the type of bedrock has an effect on the quality of groundwater from NECR aquifers. Fluoride concentrations were significantly higher in groundwater samples from the alkali granite, peraluminous granite, and metaluminous granite lithology groups than from samples in the other lithology groups. Water samples from 1.4 percent of 2,167 studied wells had fluoride concentrations that were equal to or greater than the maximum contaminant level (MCL) of 4 milligrams per liter (mg/L) and 7.5 percent of the wells had fluoride concentrations that were equal to or greater than the secondary MCL of 2 mg/L. For arsenic, groundwater samples from the calcareous metasedimentary rocks in the New Hampshire-Maine geologic province, peraluminous granite, and pelitic rocks lithology groups had higher concentrations than did samples from the other lithology groups. Water samples from 13.3 percent of 2,054 studied wells had arsenic concentrations that were equal to or greater than the MCL of 10 micrograms per liter (μg/L), about double the national rate of occurrence in community-supply systems and in domestic wells of the United States. Uranium concentrations were significantly higher in groundwater samples from the peraluminous granite, alkali granite, and calcareous metasedimentary rocks in the New Hampshire-Maine geologic province lithology groups than from samples in the other lithology groups. Water samples from 14.2 percent of 556 studied wells had uranium concentrations equal to or greater than the MCL of 30 μg/L. Radon activities were equal to or greater than the proposed MCL of 300 picocuries per liter (pCi/L) in 95 percent of 943 studied wells, and 33 percent of the wells had radon activities were equal to or greater than the proposed alternative maximum contaminant level (AMCL) of 4,000 pCi/L. Radon activities exceeded the proposed AMCL in 20 percent or more of groundwater samples in each of the studied lithology groups with a minimum of 9 samples, but radon activities were significantly higher in groundwater samples from the alkali granite, peraluminous granite, and Narragansett basin metasedimentary rocks lithology groups. Water samples from 3.2 percent of 564 studied wells had combined radium activities equal to or greater than the MCL of 5 pCi/L; however, combined radium activities were not significantly different among the studied lithology groups. Land use and population density also were evaluated to explain patterns in water quality. Concentrations of nitrate, sodium, chloride, and MtBE from the studied wells were significantly greater in areas of high population density (≥50 persons per square kilometer) than in areas of low population density (<50 persons per square kilometer). Concentrations of sodium, chloride, and M t BE from the studied wells were significantly greater in areas classified as developed (urban lands) than in areas classified as undeveloped (forested), agricultural, or mixed (no dominant land use). Nitrate concentrations from the public-supply wells were not significantly different among the four land use categories, but nitrate concentrations from the domestic wells were significantly greater in areas classified as developed than in areas classified as undeveloped, agricultural, or mixed. Chloride to bromide mass ratios in the domestic well samples indicate that the groundwater was probably affected by at least three halogen sources: local precipitation and recharge waters, remnant seawater and connate waters evolved from seawater, and recharge waters affected by road salt. The groundwater in the NECR aquifers generally contained low concentrations of nitrate, VOCs, and pesticides. Less than 1 percent of water samples from 4,781 studied wells had concentrations of nitrate greater than the MCL of 10 mg/L. Less than 1 percent of water samples from 1,299 studied wells exceeded the USEPA advisory level of 20 to 40 μg/L for M t BE. None of the other studied VOCs exceeded a human health benchmark. M t BE (36 percent frequency detection) and chloroform (32.9 percent frequency detection) were the most frequently detected (>0.02 μg/L) VOCs in the domestic wells. M t BE was detected more often in water samples with apparent ages of less than 25 years than in water samples with apparent ages greater than 25 years. This finding is consistent with the time period of high M t BE use in areas in the United States where reformulated gasoline was mandated. The largest pesticide concentration was an estimated concentration of 0.06 μg/L for the herbicide metolachlor. Deethylatrazine, a degradate of atrazine, (18 percent frequency detection) and atrazine (8 percent frequency detection) were the only pesticide compounds detected (>0.001 μg/L) in more than 3 percent of the domestic wells. None of the detected pesticide compounds exceeded human health benchmarks. Concentrations of nitrate and gross alpha-particle activities were significantly greater in the water samples from the domestic wells than in samples from the public-supply wells. Concentrations of sodium, chloride, iron, manganese, and uranium were significantly greater in the water samples from the public-supply wells than in the samples from the domestic wells. One possible explanation may be related to differences in field processing (filtered samples from the domestic wells compared to unfiltered samples from the public-supply wells). The high frequency of detections for a wide variety of manmade and naturally occurring contaminants in both domestic and public-supply wells shows the vulnerability of NECR aquifers to contamination. The highly variable water quality and the association with highly variable lithology of crystalline bedrock underscores the importance of testing individual wells to determine if concentrations for the most commonly detected contaminants exceed human health benchmarks.

Connecticut, Maine, Massachusetts, New Hampshire, ↗

Occurrence and hydrogeochemistry of radiochemical constituents in groundwater of Jefferson County and surrounding areas, southwestern Montana, 2007 through 2010

The U.S. Geological Survey, in cooperation with Jefferson County and the Jefferson Valley Conservation District, sampled groundwater in southwestern Montana to evaluate the occurrence and concentration of naturally-occurring radioactive constituents and to identify geologic settings and environmental conditions in which elevated concentrations occur. A total of 168 samples were collected from 128 wells within Broadwater, Deer Lodge, Jefferson, Lewis and Clark, Madison, Powell, and Silver Bow Counties from 2007 through 2010. Most wells were used for domestic purposes and were primary sources of drinking water for individual households. Water-quality samples were collected from wells completed within six generalized geologic units, and analyzed for constituents including uranium, radon, gross alpha-particle activity, and gross beta-particle activity. Thirty-eight wells with elevated concentrations or activities were sampled a second time to examine variability in water quality throughout time. These water-quality samples were analyzed for an expanded list of radioactive constituents including the following: three isotopes of uranium (uranium-234, uranium-235, and uranium-238), three isotopes of radium (radium-224, radium-226, and radium-228), and polonium-210. Existing U.S. Geological Survey and Montana Bureau of Mines and Geology uranium and radon water-quality data collected as part of other investigations through 2011 from wells within the study area were compiled as part of this investigation. Water-quality data from this study were compared to data collected nationwide by the U.S. Geological Survey through 2011. Radionuclide samples for this study typically were analyzed within a few days after collection, and therefore data for this study may closely represent the concentrations and activities of water being consumed locally from domestic wells. Radioactive constituents were detected in water from every well sampled during this study regardless of location or geologic unit. Nearly 41 percent of sampled wells had at least one radioactive constituent concentration that exceeded U.S. Environmental Protection Agency drinking-water standards or screening levels. Uranium concentrations were higher than the U.S. Environmental Protection Agency maximum contaminant level (MCL) of 30 micrograms per liter in samples from 14 percent of the wells. Radon concentrations exceeded a proposed MCL of 4,000 picocuries per liter in 27 percent of the wells. Combined radium (radium-226 and radium-228) exceeded the MCL of 5 picocuries per liter in samples from 10 of 47 wells. About 40 percent (42 of 104 wells) of the wells had gross alpha-particle activities (72-hour count) at or greater than a screening level of 15 pCi/L. Gross beta-particle activity exceeded the U.S. Environmental Protection Agency 50 picocuries per liter screening level in samples from 5 of 104 wells. Maximum radium-224 and polonium-210 activities in study wells were 16.1 and 3.08 picocuries per liter, respectively; these isotopes are constituents of human-health concern, but the U.S. Environmental Protection Agency has not established MCLs for them. Radioactive constituent concentrations or activities exceeded at least one established drinking-water standard, proposed drinking-water standard, or screening level in groundwater samples from five of six generalized geologic units assessed during this study. Radioactive constituent concentrations or activities were variable not only within each geologic unit, but also among wells that were completed in the same geologic unit and in close proximity to one another. Established or proposed drinking-water standards were exceeded most frequently in water from wells completed in the generalized geologic unit that includes rocks of the Boulder batholith and other Tertiary through Cretaceous igneous intrusive rocks (commonly described as granite). Specifically, of the wells completed in the Boulder batholith and related rocks sampled as part of this study, 24 percent exceeded the MCL of 30 micrograms per liter for uranium, 50 percent exceeded the proposed alternative MCL of 4,000 picocuries per liter for radon, and 27 percent exceeded the MCL of 5 micrograms per liter for combined radium-226 and radium-228. Elevated radioactive constituent values were detected in samples representing a large range of field properties and water types. Correlations between radioactive constituents and pH, dissolved oxygen, and most major ions were not statistically significant (p-value > 0.05) or were weakly correlated with Spearman correlation coefficients (rho) ranging from -0.5 to 0.5. Moderate correlations did exist between gross beta-particle activity and potassium (rho = 0.72 to 0.82), likely because one potassium isotope (potassium-40) is a beta-particle emitter. Total dissolved solids and specific conductance also were moderately correlated (rho = 0.62 to 0.71) with gross alpha-particle and gross beta-particle activity, indicating that higher radioactivity values can be associated with higher total dissolved solids. Correlations were evaluated among radioactive constituents. Moderate to strong correlations occurred between gross alpha-particle and beta-particle activities (rho = 0.77 to 0.96) and radium isotopes (rho = 0.78 to 0.92). Correlations between gross alpha-particle activity (72-hour count) and all analyzed radioactive constituents were statistically significant (p-value < 0.05), and therefore, gross alpha-particle activity likely may be used as a screening tool for determining the presence of radionuclides in area waters. In this study, gross alpha-particle activities of 7 picocuries per liter or greater were associated with all radioactive constituents whose concentrations exceeded drinking-water standards or screening levels. Radiochemical results varied temporally in samples from several of the thirty-eight wells sampled at least twice during the study. The time between successive sampling events ranged from about 1 to 10 months for 29 wells to about 3 years for the other 9 wells. Radiochemical constituents that varied by greater than 30 percent between sampling events included uranium (29 percent of the resampled wells), and radon (11 percent of the resampled wells), gross alpha-particle activity (38 percent of the resampled wells), and gross beta-particle activity (15 percent of the resampled wells). Variability in uranium concentrations from two wells was sufficiently large that concentrations were less than the MCL in the first set of samples and greater than the MCL in the second. Sample holding times affect analytical results in this study. Gross alpha-particle and gross beta-particle activities were measured twice, 72 hours and 30 days after sample collection. Gross alpha-particle activity decreased an average of 37 percent between measurements, indicating the presence of short-lived alpha-emitting radionuclides in these samples. Gross beta-particle activity increased an average of 31 percent between measurements, indicating ingrowth of longer-lived beta-emitting radionuclides.

Montana↗

Status and understanding of groundwater quality in the Cascade Range and Modoc Plateau study unit, 2010: California GAMA Priority Basin Project

Groundwater quality in the Cascade Range and Modoc Plateau study unit was investigated as part of the California State Water Resources Control Board&rsquo;s Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project. The study was designed to provide a statistically unbiased assessment of untreated groundwater quality in the primary aquifer system. The depth of the primary aquifer system for the Cascade Range and Modoc Plateau study unit was delineated by the depths of the screened or open intervals of wells in the State of California&rsquo;s database of public-supply wells. Two types of assessments were made: a status assessment that described the current quality of the groundwater resource, and an understanding assessment that made evaluations of relations between groundwater quality and potential explanatory factors representing characteristics of the primary aquifer system. The assessments characterize the quality of untreated groundwater, not the quality of treated drinking water delivered to consumers by water distributors. The status assessment was based on water-quality data collected in 2010 by the U.S. Geological Survey from 90 wells and springs (USGS-grid wells) and on water-quality data compiled from the State of California&rsquo;s regulatory compliance database for samples collected from 240 public-supply wells between September 2007 and September 2010. To provide context, the water-quality data discussed in this report were compared to California and Federal drinking-water regulatory and non-regulatory benchmarks for treated drinking water. Groundwater quality is defined in terms of relative concentrations (RCs), which are calculated by dividing the concentration of a constituent in groundwater by the concentration of the benchmark for that constituent. The RCs for inorganic constituents (major ions, trace elements, nutrients, and radioactive constituents) were classified as &ldquo;high&rdquo; (the RC is greater than 1.0, indicating that the concentration is above the benchmark), &ldquo;moderate&rdquo; (the RC is from 1.0 to greater than 0.5), or &ldquo;low&rdquo; (the RC is less than or equal to 0.5). For organic constituents (volatile organic compounds and pesticides) and special-interest constituents (perchlorate), the boundary between moderate and low RCs was set at 0.1. All benchmarks used for organic constituents were health-based. For inorganic constituents, health-based and aesthetic-based benchmarks were used. Constituents without benchmarks were not considered in the status assessment . The primary metric used for quantifying regional-scale groundwater quality was the aquifer-scale proportion&mdash;the areal percentages of the primary aquifer system with high, moderate, and low RCs for a given constituent or class of constituents. The study unit was divided into six study areas on the basis of geologic differences (Eastside Sacramento Valley, Honey Lake Valley groundwater basin, Cascade Range and Modoc Plateau Low Use Basins, Quaternary Volcanic Areas, Shasta Valley and Mount Shasta Volcanic Area, and Tertiary Volcanic Areas), and each study area was divided into equal-area grid cells. Aquifer-scale proportions were calculated for individual constituents and constituent classes for each of the six study areas and for the study unit as a whole by using grid-based (one well per cell) and spatially weighted (many wells per cell) statistical methods. The status assessment showed that inorganic constituents were present at high and moderate RCs in greater proportions of the Cascade Range and Modoc Plateau study unit than were organic constituents. One or more inorganic constituents with health-based benchmarks were present at high RCs in 9.4 percent, and at moderate RCs in 14.7 percent of the primary aquifer system. Arsenic was present at high RCs in approximately 3 percent of the primary aquifer system; boron, molybdenum, uranium, and vanadium each were present at high RCs in approximately 2 percent of the primary aquifer system. One or more inorganic constituents with aesthetic-based benchmarks were present at high RCs in 15.1 percent of the primary aquifer system and at moderate RCs in 4.9 percent. Manganese, iron, and total dissolved solids were present at high RCs in approximately 12 percent, 5 percent, and 2 percent, respectively, of the primary aquifer system. Organic constituents were not detected at high or moderate RCs in the primary aquifer system, and one or more organic constituents were detected at low RCs in approximately 40 percent of the primary aquifer system. Two classes of organic constituents were detected in more than 10 percent of the primary aquifer system: trihalomethanes (chloroform only) and herbicides. The special interest constituent perchlorate was not detected at high RCs, but was detected at moderate RCs in approximately 2 percent of the primary aquifer system. The understanding assessment relied on statistical tests to evaluate relations between concentrations of constituents and values of potential explanatory factors representing geology, land use, well construction, hydrologic conditions, groundwater age, and geochemical conditions. The majority of the high and moderate RCs of arsenic, boron, molybdenum, uranium, and total dissolved solids were in samples from the Honey Lake Valley groundwater basin study area. Groundwater mixing with hydrothermal fluids present in the study area, evaporative concentration of groundwater in the Honey Lake playa, presence of uranium-bearing sediment derived from the adjacent Sierra Nevada, and release of arsenic and other trace elements from sediments under high pH and low dissolved oxygen conditions all appeared to contribute to these elevated concentrations. Thermal springs are in many parts of the Cascade Range and Modoc Plateau study unit and could account for locally elevated concentrations of arsenic, boron, molybdenum, and total dissolved solids in samples from the other study areas. Vanadium concentrations were greater in oxic samples than in anoxic samples, but were not correlated with pH, contrary to expectations from previous studies. Organic constituents were not detected at high or moderate RCs, and the occurrence of low organic constituents at low RCs ranged from 27 percent to 73 percent of the primary aquifers system in the six study areas. The Shasta Valley and Mount Shasta Volcanic study area had significantly greater occurrence of low RCs of herbicides compared to all of the other study areas, which could reflect the greater prevalence of modern groundwater in the Shasta Valley and Mount Shasta Volcanic study area and the presence of potential sources of herbicides, including applications to timberlands and roadside rights-of-way. The Eastside Sacramento Valley study area had the greatest occurrence of low concentrations of chloroform, and chloroform occurrence was most strongly associated with the combination of septic-tank density greater than two tanks per square kilometer and urban land use greater than 10 percent within a radius of 500 meters of the well. These conditions were most prevalent in the Eastside Sacramento Valley study area. The detection frequency of low concentrations of perchlorate was consistent with the probability of occurrence expected under natural conditions, except in the Eastside Sacramento Valley study area, where detection frequencies were much higher than expected and could not be explained by known anthropogenic sources of perchlorate.

California↗

Evaluation of the expected moments algorithm and a multiple low-outlier test for flood frequency analysis at streamgaging stations in Arizona

Flooding is among the costliest natural disasters in terms of loss of life and property in Arizona, which is why the accurate estimation of flood frequency and magnitude is crucial for proper structural design and accurate floodplain mapping. Current guidelines for flood frequency analysis in the United States are described in Bulletin 17B (B17B), yet since B17B’s publication in 1982 (Interagency Advisory Committee on Water Data, 1982), several improvements have been proposed as updates for future guidelines. Two proposed updates are the Expected Moments Algorithm (EMA) to accommodate historical and censored data, and a generalized multiple Grubbs-Beck (MGB) low-outlier test. The current guidelines use a standard Grubbs-Beck (GB) method to identify low outliers, changing the determination of the moment estimators because B17B uses a conditional probability adjustment to handle low outliers while EMA censors the low outliers. B17B and EMA estimates are identical if no historical information or censored or low outliers are present in the peak-flow data. EMA with MGB (EMA-MGB) test was compared to the standard B17B (B17B-GB) method for flood frequency analysis at 328 streamgaging stations in Arizona. The methods were compared using the relative percent difference (RPD) between annual exceedance probabilities (AEPs), goodness-of-fit assessments, random resampling procedures, and Monte Carlo simulations. The AEPs were calculated and compared using both station skew and weighted skew. Streamgaging stations were classified by U.S. Geological Survey (USGS) National Water Information System (NWIS) qualification codes, used to denote historical and censored peak-flow data, to better understand the effect that nonstandard flood information has on the flood frequency analysis for each method. Streamgaging stations were also grouped according to geographic flood regions and analyzed separately to better understand regional differences caused by physiography and climate. The B17B-GB and EMA-MGB RPD-boxplot results showed that the median RPDs across all streamgaging stations for the 10-, 1-, and 0.2-percent AEPs, computed using station skew, were approximately zero. As the AEP flow estimates decreased (that is, from 10 to 0.2 percent AEP) the variability in the RPDs increased, indicating that the AEP flow estimate was greater for EMA-MGB when compared to B17B-GB. There was only one RPD greater than 100 percent for the 10- and 1-percent AEP estimates, whereas 19 RPDs exceeded 100 percent for the 0.2-percent AEP. At streamgaging stations with low-outlier data, historical peak-flow data, or both, RPDs ranged from −84 to 262 percent for the 0.2-percent AEP flow estimate. When streamgaging stations were separated by the presence of historical peak-flow data (that is, no low outliers or censored peaks) or by low outlier peak-flow data (no historical data), the results showed that RPD variability was greatest for the 0.2-AEP flow estimates, indicating that the treatment of historical and (or) low-outlier data was different between methods and that method differences were most influential when estimating the less probable AEP flows (1, 0.5, and 0.2 percent). When regional skew information was weighted with the station skew, B17B-GB estimates were generally higher than the EMA-MGB estimates for any given AEP. This was related to the different regional skews and mean square error used in the weighting procedure for each flood frequency analysis. The B17B-GB weighted skew analysis used a more positive regional skew determined in USGS Water Supply Paper 2433 (Thomas and others, 1997), while the EMA-MGB analysis used a more negative regional skew with a lower mean square error determined from a Bayesian generalized least squares analysis. Regional groupings of streamgaging stations reflected differences in physiographic and climatic characteristics. Potentially influential low flows (PILFs) were more prevalent in arid regions of the State, and generally AEP flows were larger with EMA-MGB than with B17B-GB for gaging stations with PILFs. In most cases EMA-MGB curves would fit the largest floods more accurately than B17B-GB. In areas of the State with more baseflow, such as along the Mogollon Rim and the White Mountains, streamgaging stations generally had fewer PILFs and more positive skews, causing estimated AEP flows to be larger with B17B-GB than with EMA-MGB. The effect of including regional skew was similar for all regions, and the observed pattern was increasingly greater B17B-GB flows (more negative RPDs) with each decreasing AEP quantile. A variation on a goodness-of-fit test statistic was used to describe each method’s ability to fit the largest floods. The mean absolute percent difference between the measured peak flows and the log-Pearson Type 3 (LP3)-estimated flows, for each method, was averaged over the 90th, 75th, and 50th percentiles of peak-flow data at each site. In most percentile subsets, EMA-MGB on average had smaller differences (1 to 3 percent) between the observed and fitted value, suggesting that the EMA-MGB-LP3 distribution is fitting the observed peak-flow data more precisely than B17B-GB. The smallest EMA-MGB percent differences occurred for the greatest 10 percent (90th percentile) of the peak-flow data. When stations were analyzed by USGS NWIS peak flow qualification code groups, the stations with historical peak flows and no low outliers had average percent differences as high as 11 percent greater for B17B-GB, indicating that EMA-MGB utilized the historical information to fit the largest observed floods more accurately. A resampling procedure was used in which 1,000 random subsamples were drawn, each comprising one-half of the observed data. An LP3 distribution was fit to each subsample using B17B-GB and EMA-MGB methods, and the predicted 1-percent AEP flows were compared to those generated from distributions fit to the entire dataset. With station skew, the two methods were similar in the median percent difference, but with weighted skew EMA-MGB estimates were generally better. At two gages where B17B-GB appeared to perform better, a large number of peak flows were deemed to be PILFs by the MGB test, although they did not appear to depart significantly from the trend of the data (step or dogleg appearance). At two gages where EMA-MGB performed better, the MGB identified several PILFs that were affecting the fitted distribution of the B17B-GB method. Monte Carlo simulations were run for the LP3 distribution using different skews and with different assumptions about the expected number of historical peaks. The primary benefit of running Monte Carlo simulations is that the underlying distribution statistics are known, meaning that the true 1-percent AEP is known. The results showed that EMA-MGB performed as well or better in situations where the LP3 distribution had a zero or positive skew and historical information. When the skew for the LP3 distribution was negative, EMA-MGB performed significantly better than B17B-GB and EMA-MGB estimates were less biased by more closely estimating the true 1-percent AEP for 1, 2, and 10 historical flood scenarios.

Arizona↗

Borehole geophysical investigation of a formerly used defense site, Machiasport, Maine, 2003-2006

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, collected borehole geophysical logs in 18 boreholes and interpreted the data along with logs from 19 additional boreholes as part of an ongoing, collaborative investigation at three environmental restoration sites in Machiasport, Maine. These sites, located on hilltops overlooking the seacoast, formerly were used for military defense. At each of the sites, chlorinated solvents, used as part of defense-site operations, have contaminated the fractured-rock aquifer. Borehole geophysical techniques and hydraulic methods were used to characterize bedrock lithology, fractures, and hydraulic properties. In addition, each geophysical method was evaluated for effectiveness for site characterization and for potential application for further aquifer characterization and (or) evaluation of remediation efforts. Results of borehole geophysical logging indicate the subsurface is highly fractured, metavolcanic, intrusive, metasedimentary bedrock. Selected geophysical logs were cross-plotted to assess correlations between rock properties. These plots included combinations of gamma, acoustic reflectivity, electromagnetic induction conductivity, normal resistivity, and single-point resistance. The combined use of acoustic televiewer (ATV) imaging and natural gamma logs proved to be effective for delineating rock types. Each of the rock units in the study area could be mapped in the boreholes, on the basis of the gamma and ATV reflectivity signatures. The gamma and mean ATV reflectivity data were used along with the other geophysical logs for an integrated interpretation, yielding a determination of quartz monzonite, rhyolite, metasedimentary units, or diabase/gabbro rock types. The interpretation of rock types on the basis of the geophysical logs compared well to drilling logs and geologic mapping. These results may be helpful for refining the geologic framework at depth. A stereoplot of all fractures intersecting the boreholes indicates numerous fractures, a high proportion of steeply dipping fractures, and considerable variation in fracture orientation. Low-dip-angle fractures associated with unloading and exfoliation are also present, especially at a depth of less than 100 feet below the top of casing. These sub-horizontal fractures help to connect the steeply dipping fractures, making this a highly connected fracture network. The high variability in the fracture orientations also increases the connectivity of the fracture network. A preliminary comparison of all fracture data from all the boreholes suggests fracturing decreases with depth. Because all the boreholes were not drilled to the same depth, however, there is a clear sampling bias. Hence, the deepest boreholes are analyzed separately for fracture density. For the deepest boreholes in the study, the intensity of fracturing does not decline significantly with depth. It is possible the fractures observed in these boreholes become progressively tighter or closed with depth, but this is difficult to verify with the borehole methods used in this investigation. The fact that there are more sealed fractures at depth (observed in optical televiewer logs in some of the boreholes) may indicate less opening of the sealed fractures, less water moving through the rock, and less weathering of the fracture infilling minerals. Although the fracture orientation remained fairly constant with depth, differences in the fracture patterns for the three restoration sites indicate the orientation of fractures varies across the study area. The fractures in boreholes on Miller Mountain predominantly strike northwest-southeast, and to a lesser degree they strike northeast. The fractures on or near the summit of Howard Mountain strike predominantly east-west and dip north and south, and the fractures near the Transmitter Site strike northeast-southwest and dip northwest and southeast. The fracture populations for the boreholes on or near the summit of Howard Mountain show more variation than at the other two sites. This variation may be related to the proximity of the fault, which is northeast of the summit of Howard Mountain. In a side-by-side comparison of stereoplots from selected boreholes, there was no clear correspondence between fracture orientation and proximity to the fault. There is, however, a difference in the total populations of fractures for the boreholes on or near the summit of Howard Mountain and the boreholes near the Transmitter Site. Further to the southwest and further away from the fault, the fractures at the Transmitter Site predominantly strike northeast-southwest and northwest-southeast.Heat-pulse flowmeter (HPFM) logging was used to identify transmissive fractures and to estimate the hydraulic properties along the boreholes. Ambient downflow was measured in 13 boreholes and ambient upflow was measured in 9 boreholes. In nine other bedrock boreholes, the HPFM did not detect measurable vertical flow. The observed direction of vertical flow in the boreholes generally was consistent with the conceptual flow model of downward movement in recharge locations and upward flow in discharge locations or at breaks in the slope of land surface. Under low-rate pumping or injection rates [0.25 to 1 gallon per minute (gal/min)], one to three inflow zones were identified in each borehole. Two limitations of HPFM methods are (1) the HPFM can only identify zones within 1.5 to 2 orders of magnitude of the most transmissive zone in each borehole, and (2) the HPFM cannot detect flow rates less than 0.010 + or - 0.005 gal/min, which corresponds to a transmissivity of about 1 foot squared per day (ft2/d). Consequently, the HPFM is considered an effective tool for identifying the most transmissive fractures in a borehole, down to its detection level. Transmissivities below that cut-off must be measured with another method, such as packer testing or fluid-replacement logging. Where sufficient water-level and flowmeter data were available, HPFM results were numerically modeled. For each borehole model, the fracture location and measured flow rates were specified, and the head and transmissivity of each fracture zone were adjusted until a model fit was achieved with the interpreted ambient and stressed flow profiles. The transmissivities calculated by this method are similar to the results of an open-hole slug test; with the added information from the flowmeter, however, the head and transmissivity of discrete zones also can be determined. The discrete-interval transmissivities ranged from 0.16 to 330 ft2/d. The flowmeter-derived open-hole transmissivity, which is the combined total of each of the transmissive zones, ranged from 1 to 511 ft2/d. The whole-well open-hole transmissivity values determined with HPFM methods were compared to the results of open-hole hydraulic tests. Despite the fact that the flowmeter-derived transmissivities consistently were lower than the estimates derived from open-hole hydraulic tests alone, the correlation was very strong (with a coefficient of determination, R2, of 0.9866), indicating the HPFM method provides a reasonable estimate of transmissivities for the most transmissive fractures in the borehole. Geologic framework, fracture characterization, and estimates of hydraulic properties were interpreted together to characterize the fracture network. The data and interpretation presented in this report should provide information useful for site investigators as the conceptual site groundwater flow model is refined. Collectively, the results and the conceptual site model are important for evaluating remediation options and planning or implementing the design of a well field and borehole completions that will be adequate for monitoring flow, remediation efforts, groundwater levels, and (or) water quality. Similar kinds of borehole geophysical logging (specifically the borehole imaging, gamma, fluid logs, and HPFM) should be conducted in any newly installed boreholes and integrated with interpretations of any nearby boreholes. If boreholes are installed close to existing or other new boreholes, cross-hole flowmeter surveys may be appropriate and may help characterize the aquifer properties and connections between the boreholes.

Maine↗

Nutrient concentrations and loads in the northeastern United States - Status and trends, 1975-2003

The U.S. Geological Survey (USGS) National Water-Quality Assessment Program (NAWQA) began regional studies in 2003 to synthesize information on nutrient concentrations, trends, stream loads, and sources. In the northeastern United States, a study area that extends from Maine to central Virginia, nutrient data were evaluated for 130 USGS water-quality monitoring stations. Nutrient data were analyzed for trends in flow-adjusted concentrations, modeled instream (non-flow-adjusted) concentrations, and stream loads for 32 stations with 22 to 29 years of water-quality and daily mean streamflow record during 1975-2003 (termed the long-term period), and for 46 stations during 1993-2003 (termed the recent period), by using a coupled statistical model of streamflow and water quality developed by the USGS. Recent trends in flow-adjusted concentrations of one or more nutrients also were analyzed for 90 stations by using Tobit regression. Annual stream nutrient loads were estimated, and annual nutrient yields were calculated, for 47 stations for the long-term and recent periods, and for 37 additional stations that did not have a complete streamflow and water-quality record for 1993-2003. Nutrient yield information was incorporated for 9 drainage basins evaluated in a national NAWQA study, for a total of 93 stations evaluated for nutrient yields. Long-term downward trends in flow-adjusted concentrations of total nitrogen and total phosphorus (18 and 19 of 32 stations, respectively) indicate regional improvements in nutrient-related water-quality conditions. Most of the recent trends detected for total phosphorus were upward (17 of 83 stations), indicating possible reversals to the long-term improvements. Concentrations of nutrients in many streams persist at levels that are likely to affect aquatic habitat adversely and promote freshwater or coastal eutrophication. Recent trends for modeled instream concentrations, and modeled reference concentrations, were evaluated relative to ecoregion-based nutrient criteria proposed by the U.S. Environmental Protection Agency. Instream concentrations of total nitrogen and total phosphorus persist at levels higher than proposed criteria at more than one-third and about one-half, respectively, of the 46 stations analyzed. Long-term trends in nutrient loads were primarily downward, with downward trends in total nitrogen and total phosphorus loads detected at 12 and 17 of 32 stations, respectively. Upward trends were rare, with one upward trend for total nitrogen loads and none for total phosphorus. Trends in loads of nitrite-plus-nitrate nitrogen included 7 upward and 8 downward trends among 32 stations. Downward trends in loads of ammonia nitrogen and total Kjeldahl nitrogen were detected at all six stations evaluated. Long-term downward trends detected in four of the five largest drainage basins evaluated include: total nitrogen loads for the Connecticut, Delaware, and James Rivers; total Kjeldahl nitrogen and ammonia nitrogen loads for the Susquehanna River; ammonia nitrogen and nitrite-plus-nitrate nitrogen loads for the James River; and total phosphorus loads for the Connecticut and Delaware Rivers. No trends in load were detected for the Potomac River. Nutrient yields were evaluated relative to the extent of land development in 93 drainage basins. The undeveloped land-use category included forested drainage basins with undeveloped land ranging from 75 to 100 percent of basin area. Median total nitrogen yields for the 27 undeveloped drainage basins evaluated, including 9 basins evaluated in a national NAWQA study, ranged from 290 to 4,800 pounds per square mile per year (lb/mi 2 /yr). Total nitrogen yields even in the most pristine drainage basins may be elevated relative to natural conditions, because of high rates of atmospheric deposition of nitrogen in parts of the northeastern United States. Median total phosphorus yields ranged from 12 to 330 lb/mi 2 /yr for the 26 undeveloped basins evaluated. The undeveloped category includes some large drainage basins with point-source discharges and small percentages of developed land; in these basins, streamflow from undeveloped headwater areas dilutes streamflow in more urbanized reaches, and dampens but does not eliminate the point-source "signal" of higher nutrient loads. Median total nitrogen yields generally do not exceed 1,700 lb/mi 2 /yr, and median total phosphorus yields generally do not exceed 100 lb/mi 2 /yr, in the drainage basins that are least affected by human land-use and waste-disposal practices. Agricultural and urban land use has increased nutrient yields substantially relative to undeveloped drainage basins. Median total nitrogen yields for 24 agricultural basins ranged from 1,700 to 26,000 lb/mi 2 /yr, and median total phosphorus yields ranged from 94 to 1,000 lb/mi 2 /yr. The maximum estimated total nitrogen and total phosphorus yields, 32,000 and 16,000 lb/mi 2 /yr, respectively, for all stations in the region were in small (less than 50 square miles (mi 2 )) agricultural drainage basins. Median total nitrogen yields ranged from 1,400 to 17,000 lb/mi 2 /yr in 26 urbanized drainage basins, and median total phosphorus yields ranged from 43 to 1,900 lb/mi 2 /yr. Urbanized drainage basins with the highest nutrient yields are generally small (less than 300 mi2) and are drained by streams that receive major point-source discharges. Instream nutrient loads were evaluated relative to loads from point-source discharges in four drainage basins: the Quinebaug River Basin in Connecticut, Massachusetts, and Rhode Island; the Raritan River Basin in New Jersey; the Patuxent River Basin in Maryland; and the James River Basin in Virginia. Long-term downward trends in nutrient loads, coupled with similar trends in flow-adjusted nutrient concentrations, indicate long-term reductions in the delivery of most nutrients to these streams. However, the absence of recent downward trends in load for most nutrients, coupled with instream concentrations that exceed proposed nutrient criteria in several of these waste-receiving streams, indicates that challenges remain in reducing delivery of nutrients to streams from point sources. During dry years, the total nutrient load from point sources in some of the drainage basins approached or equaled the nutrient load transported by the stream.

Connecticut;Delaware;Maine;Maryl;Massachusetts;New↗

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

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

California↗

Status and understanding of groundwater quality in the Monterey-Salinas Shallow Aquifer Study Unit, 2012–13: California GAMA Priority Basin Project

Groundwater quality in the approximately 7,820-square-kilometer (km 2 ) Monterey-Salinas Shallow Aquifer (MS-SA) study unit was investigated from October 2012 to May 2013 as part of the second phase of the Priority Basin Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is in the central coast region of California in the counties of Santa Cruz, Monterey, and San Luis Obispo. The GAMA Priority Basin Project is being conducted by the California State Water Resources Control Board in cooperation with the U.S. Geological Survey and the Lawrence Livermore National Laboratory. The MS-SA study was designed to provide a statistically robust assessment of untreated-groundwater quality in the shallow aquifer systems. The assessment was based on water-quality samples collected by the U.S. Geological Survey from 100 groundwater sites and 70 household tap sites, along with ancillary data such as land use and well-construction information. The shallow aquifer systems were defined by the depth interval of wells associated with domestic supply. The MS-SA study unit consisted of four study areas—Santa Cruz (210 km 2 ), Pajaro Valley (360 km 2 ), Salinas Valley (2,000 km 2 ), and Highlands (5,250 km 2 ). This study had two primary components: the status assessment and the understanding assessment . The first primary component of this study—the status assessment —assessed the quality of the groundwater resource indicated by data from samples analyzed for volatile organic compounds (VOCs), pesticides, and naturally present inorganic constituents, such as major ions and trace elements. The status assessment is intended to characterize the quality of groundwater resources in the shallow aquifer system of the MS-SA study unit, not the treated drinking water delivered to consumers by water purveyors. As opposed to the public wells, however, water from private wells, which often tap the shallow aquifer, is usually consumed without any treatment. The second component of this study—the understanding assessment —identified the natural and human factors that potentially affect groundwater quality by evaluating land-use characteristics, measures of location, geologic factors, groundwater age, and geochemical conditions of the shallow aquifer. An additional component of this study was a comparison of MS-SA water-quality results to those of the GAMA Monterey Bay and Salinas Valley Groundwater Basins study unit. This study unit covered much of the same areal extent as the MS-SA, but assessed the deeper, public drinking-water aquifer system. Relative concentrations (sample concentration divided by the benchmark concentration) were used to evaluate concentrations of constituents in groundwater samples relative to water-quality benchmarks for those constituents that have Federal or California benchmarks, such as maximum contaminant levels. For organic and special-interest constituents, relative concentrations were classified as high, greater than 1.0; moderate, greater than 0.1 and less than or equal to 1.0; or low, less than or equal to 0.1. For inorganic constituents, relative concentrations were classified as high, greater than 1.0; moderate, greater than 0.5 and less than or equal to 1.0; or low, less than or equal to 0.5. A relative concentration greater than 1.0 indicates that the concentration was greater than a benchmark. Aquifer-scale proportions were used to quantify regional-scale groundwater quality. The aquifer-scale proportions are the areal percentages of the shallow aquifer system where relative concentrations for a given constituent or class of constituents were high, moderate, or low. Inorganic constituents were measured at high and moderate relative concentrations more frequently than organic constituents. In the MS-SA study unit, inorganic constituents with benchmarks were detected at high relative concentrations in 51 percent of the study unit. The greatest proportions of high relative concentrations of trace elements and radioactive constituents were in the Highlands and Santa Cruz study areas, whereas high relative concentrations of nutrients were most often detected in the Salinas Valley and Pajaro Valley study areas and salinity indicators were most often detected in the Highlands and Salinas Valley study areas. The trace elements detected at high relative concentrations were arsenic, boron, iron, manganese, molybdenum, selenium, and strontium. The radioactive constituents detected at high relative concentrations were adjusted gross alpha radioactivity and uranium. The nutrient detected at high relative concentrations was nitrate plus nitrite. The salinity indicators detected at high relative concentrations were chloride, sulfate, and total dissolved solids. Organic constituents (VOCs and pesticides) were not detected at high relative concentrations in any of the study areas. The fumigant 1,2-dichloropropane was detected at moderate relative concentrations. The VOC chloroform and the pesticide simazine were the only organic constituents detected in more than 10 percent of samples. The constituents of special interest NDMA ( N -nitrosodimethylamine) and perchlorate were detected at high relative concentrations in the MS-SA study unit. Selected constituents were evaluated with explanatory factors to identify potential sources or processes that could explain their presence and distribution. Trace elements and radioactive constituents came from natural sources and were not elevated by anthropogenic sources or processes, except for selenium and the radioactive constituent uranium. Arsenic, manganese, iron, selenium, and uranium concentrations were all influenced by oxidation-reduction conditions. Unlike other trace elements and radioactive constituents, uranium and selenium can be affected by agricultural practices. Uranium and selenium can be released from aquifer sediments as a result of irrigation recharge water interacting with bicarbonate systems. Nitrate can be strongly affected by anthropogenic sources. Nitrate concentrations were significantly higher in modern groundwater, indicating recent inputs of nitrate to the shallow aquifer system. Nitrate was positively correlated with agricultural land use, indicating that irrigation-return water could be leaching nitrogen fertilizer and naturally present nitrate to elevate nitrate concentrations in shallow groundwater. The salinity indicators total dissolved solids, chloride, and sulfate all had natural sources in the MS-SA study unit, primarily marine sediments. Concentrations of the constituents were elevated as a result of evaporative concentration of irrigation water or precipitation. Sulfate concentrations were significantly correlated to agricultural land use, indicating that agricultural land-use practices are a contributing source of sulfate to groundwater. The samples with most of the detections of VOCs were from sites in the Pajaro Valley and northern part of the Salinas Valley. Most of the samples with pesticide detections were from sites in the Salinas Valley study area. The herbicide simazine was positively correlated to the percentage of agricultural land use, and its concentrations were higher in modern groundwater than in pre-modern groundwater. Perchlorate, similar to nitrate, has natural and anthropogenic sources. Correlations of perchlorate to dissolved oxygen, nitrate, and percentage of agricultural land use indicated that the irrigation-return water could be leaching naturally present perchlorate, as well as perchlorate from historical applications of Chilean nitrate fertilizer, to increase perchlorate concentrations in groundwater. The quality of the water in the shallow aquifer system from this study was compared with the quality of water in the public drinking-water aquifer in a previous GAMA (MS-PA) study in the same area. The shallow system was more oxic and had more sites with modern groundwater than the public drinking-water aquifer, which was more anoxic and had sites with more pre-modern groundwater. Arsenic and selenium were found at high relative concentrations in a greater proportion of the shallow system. Manganese and iron were found at high relative concentrations in a greater proportion of the public drinking-water aquifer. Uranium was found at higher relative concentrations in a greater proportion of the shallow system. Concentrations of arsenic, iron, manganese, and molybdenum are not likely to change much as groundwater percolates from the shallow system to the public drinking-water aquifer because there are no anthropogenic sources affecting these constituents. Uranium and selenium concentrations in the public drinking-water aquifer could be affected by the higher concentrations in the shallow system because of irrigation-return water, however. Nitrate and salinity indicators had concentrations that were much higher in the shallow system than the deeper public drinking-water aquifer. High concentrations of these constituents in the shallow system could lead to increased concentrations in the public drinking-water aquifer in parts of the study units because of land-use practices, such as irrigated agriculture. Organic constituents were detected more frequently in the public drinking-water aquifer than in the shallow system, possibly because more of the sites sampled in the public drinking-water aquifer were in urban areas compared to the sites sampled for the shallow system or because sources of contamination have decreased as a result of changes in use at the land surface.

California↗

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↗