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.

SourceResearch connectionFrom yearThrough year

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

At least 1,171 records · Page 65Linked to original sources

USGS2014

Evaluation of groundwater and surface-water interactions in the Caddo Nation Tribal Jurisdictional Area, Caddo County, Oklahoma, 2010-13

Streamflows, springs, and wetlands are important natural and cultural resources to the Caddo Nation. Consequently, the Caddo Nation is concerned about the vulnerability of the Rush Springs aquifer to overdrafting and whether the aquifer will continue to be a viable source of water to tribal members and other local residents in the future. Interest in the long-term viability of local water resources has resulted in ongoing development of a comprehensive water plan by the Caddo Nation. As part of a multiyear project with the Caddo Nation to provide information and tools to better manage and protect water resources, the U.S. Geological Survey studied the hydraulic connection between the Rush Springs aquifer and springs and streams overlying the aquifer. The Caddo Nation Tribal Jurisdictional Area is located in southwestern Oklahoma, primarily in Caddo County. Underlying the Caddo Nation Tribal Jurisdictional Area is the Permian-age Rush Springs aquifer. Water from the Rush Springs aquifer is used for irrigation, public, livestock and aquaculture, and other supply purposes. Groundwater from the Rush Springs aquifer also is withdrawn by domestic (self-supplied) wells, although domestic use was not included in the water-use summary in this report. Perennial streamflow in many streams and creeks overlying the Rush Springs aquifer, such as Cobb Creek, Lake Creek, and Willow Creek, originates from springs and seeps discharging from the aquifer. This report provides information on the evaluation of groundwater and surface-water resources in the Caddo Nation Jurisdictional Area, and in particular, information that describes the hydraulic connection between the Rush Springs aquifer and springs and streams overlying the aquifer. This report also includes data and analyses of base flow, evidence for groundwater and surface-water interactions, locations of springs and wetland areas, groundwater flows interpreted from potentiometric-surface maps, and hydrographs of water levels monitored in the Caddo Nation Tribal Jurisdictional Area from 2010 to 2013. Flow in streams overlying the Rush Springs aquifer, on average, were composed of 50 percent base flow in most years. Monthly mean base flow appeared to maintain streamflows throughout each year, but periods of zero flow were documented in daily hydrographs at each measured site, typically in the summer months. A pneumatic slug-test technique was used at 15 sites to determine the horizontal hydraulic conductivity of streambed sediments in streams overlying the Rush Springs aquifer. Converting horizontal hydraulic conductivities (Kh) from the slug-test analyses to vertical hydraulic conductivities (Kv) by using a ratio of Kv/Kh = 0.1 resulted in estimates of vertical streambed hydraulic conductivity ranging from 0.1 to 8.6 feet per day. Data obtained from a hydraulic potentiomanometer in streambed sediments and streams in August 2012 indicate that water flow was from the streambed sediments to the stream (gaining) at 6 of 15 sites, and that water flow was from the stream to the streambed sediments (losing) at 9 of 15 sites. The groundwater and surface-water interaction data collected at the Cobb Creek near Eakly, Okla., streamflow gaging station (07325800), indicate that the bedrock groundwater, alluvial groundwater, and surface-water resources are closely connected. Because of this hydrologic connection, large perennial streams in the study area may change from gaining to losing streams in the summer. The timing and severity of this change from a gaining to a losing condition probably is affected by the local or regional withdrawal of groundwater for irrigation in the summer growing season. Wells placed closer to streams have a greater and more immediate effect on alluvial groundwater levels and stream stages than wells placed farther from streams. Large-capacity irrigation wells, even those completed hundreds of feet below land surface in the bedrock aquifer, can induce surface-water flow from nearby streams by lowering alluvial groundwater levels below the stream altitude. Twenty-five new springs visible from public roads and paths were documented during a survey of springs in 2011. Most of the springs are in upland draws on the flanks of topographic ridges. Wetlands primarily were identified by using a combination of data sources including the National Wetlands Inventory, Soil Survey Geographic database frequently flooded soils maps, and aerial photographs. Regional flow directions were determined by analysis of water levels measured in 29 wells completed in the Rush 2 Springs aquifer in Caddo County and the Caddo Nation Tribal Jurisdictional Area. Water levels were monitored every 30 minutes in five wells by using a vented pressure transducer and a data-collection platform with real-time transmitting equipment in each well. Those five wells ranged in depth from 210 to 350 feet. Water levels in these five wells indicate that there was a decrease in water storage in the Rush Springs aquifer from October 2010 to June 2013.

Oklahoma↗

USGS2016

Estimating peak-flow frequency statistics for selected gaged and ungaged sites in naturally flowing streams and rivers in Idaho

The U.S. Geological Survey, in cooperation with the Idaho Transportation Department, updated regional regression equations to estimate peak-flow statistics at ungaged sites on Idaho streams using recent streamflow (flow) data and new statistical techniques. Peak-flow statistics with 80-, 67-, 50-, 43-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities (1.25-, 1.50-, 2.00-, 2.33-, 5.00-, 10.0-, 25.0-, 50.0-, 100-, 200-, and 500-year recurrence intervals, respectively) were estimated for 192 streamgages in Idaho and bordering States with at least 10 years of annual peak-flow record through water year 2013. The streamgages were selected from drainage basins with little or no flow diversion or regulation. The peak-flow statistics were estimated by fitting a log-Pearson type III distribution to records of annual peak flows and applying two additional statistical methods: (1) the Expected Moments Algorithm to help describe uncertainty in annual peak flows and to better represent missing and historical record; and (2) the generalized Multiple Grubbs Beck Test to screen out potentially influential low outliers and to better fit the upper end of the peak-flow distribution. Additionally, a new regional skew was estimated for the Pacific Northwest and used to weight at-station skew at most streamgages. The streamgages were grouped into six regions (numbered 1_2, 3, 4, 5, 6_8, and 7, to maintain consistency in region numbering with a previous study), and the estimated peak-flow statistics were related to basin and climatic characteristics to develop regional regression equations using a generalized least squares procedure. Four out of 24 evaluated basin and climatic characteristics were selected for use in the final regional peak-flow regression equations. Overall, the standard error of prediction for the regional peak-flow regression equations ranged from 22 to 132 percent. Among all regions, regression model fit was best for region 4 in west-central Idaho (average standard error of prediction=46.4 percent; pseudo-R 2 >92 percent) and region 5 in central Idaho (average standard error of prediction=30.3 percent; pseudo-R 2 >95 percent). Regression model fit was poor for region 7 in southern Idaho (average standard error of prediction=103 percent; pseudo-R 2 <78 percent) compared to other regions because few streamgages in region 7 met the criteria for inclusion in the study, and the region’s semi-arid climate and associated variability in precipitation patterns causes substantial variability in peak flows. A drainage area ratio-adjustment method, using ratio exponents estimated using generalized least-squares regression, was presented as an alternative to the regional regression equations if peak-flow estimates are desired at an ungaged site that is close to a streamgage selected for inclusion in this study. The alternative drainage area ratio-adjustment method is appropriate for use when the drainage area ratio between the ungaged and gaged sites is between 0.5 and 1.5. The updated regional peak-flow regression equations had lower total error (standard error of prediction) than all regression equations presented in a 1982 study and in four of six regions presented in 2002 and 2003 studies in Idaho. A more extensive streamgage screening process used in the current study resulted in fewer streamgages used in the current study than in the 1982, 2002, and 2003 studies. Fewer streamgages used and the selection of different explanatory variables were likely causes of increased error in some regions compared to previous studies, but overall, regional peak‑flow regression model fit was generally improved for Idaho. The revised statistical procedures and increased streamgage screening applied in the current study most likely resulted in a more accurate representation of natural peak-flow conditions. The updated, regional peak-flow regression equations will be integrated in the U.S. Geological Survey StreamStats program to allow users to estimate basin and climatic characteristics and peak-flow statistics at ungaged locations of interest. StreamStats estimates peak-flow statistics with quantifiable certainty only when used at sites with basin and climatic characteristics within the range of input variables used to develop the regional regression equations. Both the regional regression equations and StreamStats should be used to estimate peak-flow statistics only in naturally flowing, relatively unregulated streams without substantial local influences to flow, such as large seeps, springs, or other groundwater-surface water interactions that are not widespread or characteristic of the respective region.

Idaho↗

USGS2013

Characterization of water quality and biological communities, Fish Creek, Teton County, Wyoming, 2007-2011

Fish Creek, an approximately 25-kilometer-long tributary to Snake River, is located in Teton County in western Wyoming near the town of Wilson. Fish Creek is an important water body because it is used for irrigation, fishing, and recreation and adds scenic value to the Jackson Hole properties it runs through. Public concern about nuisance growths of aquatic plants in Fish Creek has been increasing since the early 2000s. To address these concerns, the U.S. Geological Survey conducted a study in cooperation with the Teton Conservation District to characterize the hydrology, water quality, and biologic communities of Fish Creek during 2007–11. The hydrology of Fish Creek is strongly affected by groundwater contributions from the area known as the Snake River west bank, which lies east of Fish Creek and west of Snake River. Because of this continuous groundwater discharge to the creek, land-use activities in the west bank area can affect the groundwater quality. Evaluation of nitrate isotopes and dissolved-nitrate concentrations in groundwater during the study indicated that nitrate was entering Fish Creek from groundwater, and that the source of nitrate was commonly a septic/sewage effluent or manure source, or multiple sources, potentially including artificial nitrogen fertilizers, natural soil organic matter, and mixtures of sources. Concentrations of dissolved nitrate and orthophosphate, which are key nutrients for growth of aquatic plants, generally were low in Fish Creek and occasionally were less than reporting levels (not detected). One potential reason for the low nutrient concentrations is that nutrients were being consumed by aquatic plant life that increases during the summer growing season, as a result of the seasonal increase in temperature and larger number of daylight hours. Several aspects of Fish Creek’s hydrology contribute to higher productivity and biovolume of aquatic plants in Fish Creek than typically observed in streams of its size in Wyoming. Especially in the winter, the proportionately large, continuous gain of groundwater into Fish Creek in the perennial section keeps most of the creek free of ice. Because sunlight can still reach the streambed in Fish Creek and the water is still flowing, aquatic plants continue to photosynthesize in the winter, albeit at a lower level of productivity. Additionally, the cobble and large gravel substrate in Fish Creek provides excellent attachment points for aquatic plants, and when combined with Fish Creek’s channel stability allows rapid growth of aquatic plants once conditions allow during the spring. The aquatic plant community of Fish Creek was different than most streams in Wyoming in that it contains many different macrophytes—including macroalgae such as long streamers of Cladophora , aquatic vascular plants, and moss; most other streams in the state contain predominantly algae. From the banks of Fish Creek, the bottom of the stream sometimes appeared to be a solid green carpet. A shift was observed from higher amounts of microalgae in April/May to higher amounts macrophytes in August and October, and differences in the relative abundance of microalgae and macrophytes were statistically significant between seasons. Differences in dissolved-nitrate concentrations and in the nitrogen-to-phosphorus ratio were significantly different between seasons, as concentrations of dissolved nitrate decreased from April/May to August and October. It is likely that dissolved-nitrate concentrations in Fish Creek were lower in August and October because macrophytes were quickly utilizing the nutrient, and a negative correlation between macro-phytes and nitrate was found. Macroinvertebrates also were sampled because of their role as indicators of water quality and their documented responses to perturbation such as degradation of water quality and habitat. Statistically significant seasonal differences were noted in the macroinvertebrate community. Taxa richness and relative abundance of Ephemeroptera, Plecoptera, and Trichoptera, which tend to be intolerant of water-quality degradation, decreased from April/May to August; the same time period saw a corresponding increase in Diptera and noninsects, particularly Oligochaeta (worms) that are more tolerant. Seasonal changes in macroinvertebrate functional feeding groups were significantly different. The relative abundance of gatherer-collector and scraper feeding groups decreased from April/May to August, accompanied by an increase in filterer-collector and shredders feeding groups. Seasonal changes in feeding groups might be due to the seasonal shift in aquatic plant communities, as indicated by comparison with other streams in the area that had fewer aquatic macrophytes than Fish Creek. Statistical tests of macroinvertebrate metrics indicated few differences between years or biological sampling sites on Fish Creek, although the site farthest upstream sometimes was different not only in terms of macroinvertebrates but also in streamflow, water quality, and aquatic plants. Potential effects of contributions of additional nutrients to the Fish Creek ecosystem beyond the conditions sampled during the study period are not known. However, because virtually all of the detectable dissolved nitrate commonly was consumed by aquatic plants in August (leaving dissolved nitrate less than the reporting level in water samples), it is possible that increased nutrient contributions could cause increased growth of aquatic plants. Additional long-term monitoring of the stream, with concurrent data analysis and interpretation would be needed to determine the effects of additional nutrients on the aquatic plant community and on higher levels of the food chain.

Wyoming↗

USGS2021

Hydraulic characterization of carbonate-rock and basin-fill aquifers near Long Canyon, Goshute Valley, northeastern Nevada

Understanding groundwater flow and pumping effects near pending mining operations requires accurate subsurface hydraulic characterization. To improve conceptual models of groundwater flow and development in the complex hydrogeologic system near Long Canyon Mine, in northwestern Goshute Valley, northeastern Nevada, the U.S. Geological Survey characterized the hydraulic properties of carbonate rocks and basin-fill aquifers using an integrated analysis of steady-state and stressed aquifer conditions informed by water chemistry and aquifer-test data. Hydraulic gradients and groundwater-age data in northern Goshute Valley indicate carbonate rocks in the Pequop Mountains just west and south of the Long Canyon Mine project area constitute a more permeable and active flow system than saturated rocks in the northern Pequop Mountains, western Toano Range, and basin fill. Permeable carbonate rocks in the northern Pequop Mountains, in part, discharge to the Johnson Springs wetland complex (JSWC), where mean groundwater ages range from 500 to 2,400 years and samples all contain a small fraction of modern waters, relative to mean ages of 8,600 to more than 22,000 years for most groundwater sampled to the north and east. Recharge to the JSWC occurs from a roughly 27-square-mile area in the upgradient Pequop Mountains to the west, composed mostly of permeable carbonate rock and fractured quartzite, and bounded by low-permeability shales and marbleized and siliclastic rocks. Single-well aquifer-test analyses provided transmissivity estimates at pumped wells. Transmissivity estimates ranged from 7,000 to 400,000 feet squared per day (ft 2 /d) in carbonate rocks and from 2,000 to 80,000 ft 2 /d in basin fill near the Long Canyon Mine. Water-level drawdown from multiple-well aquifer testing and rise from unintentional leakage into the overlying basin-fill aquifer were estimated and distinguished from natural fluctuations in 93 pumping and monitoring sites using analytical water-level models. Leakage of disposed aquifer-test pumpage occurred south of the aquifer test area through an unlined irrigation ditch. Drawdown was detected at distances of as much as 3 miles (mi) from pumping wells at all but one carbonate-rock site, at basin-fill sites on the alluvial fan immediately downgradient from pumping wells, and in Big Spring and spring NS-05. Similar drawdowns in carbonate rocks within the drawdown detection area suggest all wells penetrate a highly transmissive zone (HTZ) that is bounded by low-permeability rocks. Drawdown was not detected in carbonate rocks to the west of Canyon fault, in any basin-fill sites on the valley floor east of the Hardy fault, or at volcanic sites to the north, indicating that these major fault structures and (or) permeability contrasts between hydrogeologic units impeded groundwater flow or obscured pumping signals. Alternatively, unintentional leakage might have obscured drawdown at basin-fill sites on the valley floor, where water-level rise was detected at nine sites over 3 mi. Consistent hydraulic properties were estimated by simultaneously interpreting steady-state flow during predevelopment conditions and changes in groundwater levels and springflows from the 2016 carbonate-rock aquifer test with an integrated groundwater-flow model. Hydraulic properties were distributed across carbonate rocks, basin fill, volcanic rocks, and siliciclastic rocks with a hydrogeologic framework developed from geologic mapping and hydraulic testing. Estimated transmissivity distributions spanned at least three orders of magnitude in each rock unit. In the HTZ, simulated transmissivities ranged from 10,000 to 23,000,000 ft 2 /d, with the most transmissive areas occurring around Big Spring. Comparatively low carbonate-rock transmissivities of less than 10,000 ft 2 /d were estimated in the northern Pequop Mountains and poorly defined values of less than 1,000 ft 2 /d were estimated in the western Toano Range. Transmissivities in basin fill ranged from less than 10 to 80,000 ft 2 /d and were minimally constrained by the 2016 carbonate-rock aquifer test because poorly quantified leakage affected water levels more so than pumping. The most transmissive areas were informed by single-well aquifer tests along the eastern edge of the Pequop Mountains near Long Canyon Mine and could be indicative of a hydraulic connection between basin fill and more transmissive underlying carbonate rocks. Simulated transmissivities of volcanic and low-permeability rocks mostly are less than 1,000 ft 2 /d. The estimated hydraulic-property distributions and informed interpretation of hydraulic connections among hydrogeologic units improved the characterization and representation of groundwater flow near the Long Canyon Mine.

Nevada↗

USGS2013

Naturally occurring contaminants in the Piedmont and Blue Ridge crystalline-rock aquifers and Piedmont Early Mesozoic basin siliciclastic-rock aquifers, eastern United States, 1994–2008

Groundwater quality and aquifer lithologies in the Piedmont and Blue Ridge Physiographic Provinces in the eastern United States vary widely as a result of complex geologic history. Bedrock composition (mineralogy) and geochemical conditions in the aquifer directly affect the occurrence (presence in rock and groundwater) and distribution (concentration and mobility) of potential naturally occurring contaminants, such as arsenic and radionuclides, in drinking water. To evaluate potential relations between aquifer lithology and the spatial distribution of naturally occurring contaminants, the crystalline-rock aquifers of the Piedmont and Blue Ridge Physiographic Provinces and the siliciclastic-rock aquifers of the Early Mesozoic basin of the Piedmont Physiographic Province were divided into 14 lithologic groups, each having from 1 to 16 lithochemical subgroups, based on primary rock type, mineralogy, and weathering potential. Groundwater-quality data collected by the U.S. Geological Survey (USGS) National Water-Quality Assessment (NAWQA) Program from 1994 through 2008 from 346 wells and springs in various hydrogeologic and land-use settings from Georgia through New Jersey were compiled and analyzed for this study. Analyses for most constituents were for filtered samples, and, thus, the compiled data consist largely of dissolved concentrations. Concentrations were compared to criteria for protection of human health, such as U.S. Environmental Protection Agency (USEPA) drinking water maximum contaminant levels and secondary maximum contaminant levels or health-based screening levels developed by the USGS NAWQA Program in cooperation with the USEPA, the New Jersey Department of Environmental Protection, and Oregon Health & Science University. Correlations among constituent concentrations, pH, and oxidation-reduction (redox) conditions were used to infer geochemical controls on constituent mobility within the aquifers. Of the 23 trace-element constituents evaluated, arsenic, manganese, and zinc were detected in one or more water samples at concentrations greater than established human health-based criteria. Arsenic concentrations typically were less than 1 microgram per liter (µg/L) in most groundwater samples; however, concentrations of arsenic greater than 1 µg/L frequently were detected in groundwater from clastic lacustrine sedimentary rocks of the Early Mesozoic basin aquifers and from metamorphosed clastic sedimentary rocks of the Piedmont and Blue Ridge crystalline rock aquifers. Groundwater from these rock units had elevated pH compared to other rock units evaluated in this study. Of the nine samples for which arsenic concentration was greater than 10 µg/L, six were classified as oxic and three as anoxic, and seven had pH of 7.2 or greater. Manganese concentrations typically were less than 10 µg/L in most samples; however, 8.3 percent of samples from the Piedmont and Blue Ridge crystalline-rock aquifers and 3.0 percent of samples from the Early Mesozoic basin siliciclastic rock aquifers had manganese concentrations greater than the 300-µg/L health-based screening level. The positive correlation of manganese with iron and ammonia and the negative correlation of manganese with dissolved oxygen and nitrate are consistent with the reductive dissolution of manganese oxides in the aquifer. Zinc concentrations typically were less than 10 µg/L in the groundwater samples considered in the study, but 0.4 percent and 5.5 percent of the samples had concentrations greater than the health-based screening level of 2,000 µg/L and one-tenth of the health-based screening level, respectively. The mean rank concentration of zinc in groundwater from the quartz-rich sedimentary rock lithologic group was greater than that for other lithologic groups even after eliminating samples collected from wells constructed with galvanized casing. Approximately 90 percent of 275 groundwater samples had radon-222 concentrations that were greater than the proposed alternative maximum contaminant level of 300 picocuries per liter. In contrast, only 2.0 percent of 98 samples had combined radium (radium-226 plus radium-228) concentrations greater than the maximum contaminant level of 5.0 picocuries per liter, and 0.6 percent of 310 samples had uranium concentrations greater than the maximum contaminant level of 30 µg/L. Radon concentrations were highest in the Piedmont and Blue Ridge crystalline-rock aquifers, especially in granite, and elevated median concentrations were noted in the Piedmont Early Mesozoic basin aquifers, but without the extreme maximum concentrations found in the crystalline rocks (granites). Although the siliciclastic lithologies had a greater frequency of elevated uranium concentrations, radon and radium were commonly detected in water from both siliciclastic and crystalline lithologies. Uranium concentrations in groundwater from clastic sedimentary and clastic lacustrine/evaporite sedimentary lithologic groups within the Early Mesozoic basin aquifers, which had median concentrations of 3.6 and 3.1 µg/L, respectively, generally were higher than concentrations for other siliciclastic lithologic groups, which had median concentrations less than 1 µg/L. Although 89 percent of the 260 samples from crystalline-rock aquifers had uranium concentrations less than 1 µg/L, 0.8 percent had uranium concentrations greater than the 30-µg/L maximum contaminant level, and 6.5 percent had concentrations greater than 3 µg/L.

Alabama;Delaware;Georgia;Maryl;New Jersey;North Ca↗

USGS2014

Hydrogeology and hydrology of the Punta Cabullones wetland area, Ponce, southern Puerto Rico, 2007-08

The U.S. Geological Survey, in cooperation with the Municipio Autónomo de Ponce and the Puerto Rico Department of Natural and Environmental Resources, conducted a study of the hydrogeology and hydrology of the Punta Cabullones area in Ponce, southern Puerto Rico. (Punta Cabullones is also referred to as Punta Cabullón.) The Punta Cabullones area is about 9 square miles and is an ecological system made up of a wetland, tidal flats, saltflats, mangrove forests, and a small fringing reef located a short distance offshore. The swales or depressions between successive beach ridges became development avenues for saline to hypersaline wetlands. The Punta Cabullones area was designated by the U.S. Fish and Wildlife Service as a coastal barrier in the 1980s because of its capacity to act as a buffer zone to ameliorate the impacts of natural phenomenon such as storm surges. Since 2003, Punta Cabullones has been set aside for preservation as part of the mitigation effort mandated by Federal and State laws to compensate for the potential environmental effects that might be caused by the construction of the Las Américas Transshipment Port. Total rainfall measured during 2008 within the Punta Cabullones area was 36 inches, which is slightly greater than the long-term annual average of 32 inches for the coastal plain near Ponce. Two evapotranspiration estimates, 29 and 37 inches, were obtained for the subarea of the Punta Cabullones area that is underlain by fan-delta and alluvial deposits by using two variants of the Penman semi-empirical equation. The long-term water stage and chemical character of the wetland in Punta Cabullones are highly dependent on the seasonal and annual variations of both rainfall and sea-wave activity. Also, unseasonal short-term above-normal rainfall and sea-wave events resulting from passing storms may induce substantial changes in the water stage and the chemical character of the wetland. In general, tidal fluctuations exert a minor role in modifying the water quality and stage of the wetland in Punta Cabullones. The role of the tidal fluctuations becomes important during those times when the outlets/inlets to the sea are not blocked by a sand bar and is allowed to freely flow into the wetland interior. The salinity of the wetland varies from brackish to hypersaline. The hypersaline conditions, including the occurrence of saltflats, within the Punta Cabullones wetland area result from a high evapotranspiration rate. The hypersaline conditions are further enhanced by a sand bar that blocks the inlet/outlet of the wetland’s easternmost channel, particularly during the dry season. Groundwater in Punta Cabullones mostly is present within beds of silisiclastic sand and gravel. During the study period, the depth to groundwater did not exceed 4 feet below land surface. The movement and direction of the groundwater flow in Punta Cabullones are driven by density variations that in turn result from the wide range of salinities in the groundwater. The salinity of the groundwater decreases within the first 60 to 100 feet of depth and decreases outward from a mound of hypersaline groundwater centered on piezometer nest PN2. The main groundwater types within the Punta Cabullones area vary from calcium-bicarbonate type in the northernmost part of the study area to a predominantly sodium-potassium-chloride groundwater type southward. According to stable-isotope data, groundwater within the study area is both modern meteoric water and seawater highly affected by evaporation. The chemical and stable-isotopic character of local groundwater is highly influenced by evapotranspiration because of its shallow depth. Equivalent freshwater heads indicate groundwater moves away from a mound centered on piezometer nest PN2, in a pattern similar to the spatial distribution of groundwater salinity. Vertical groundwater flow occurs in Punta Cabullones due to local differences in density. In the wetland subarea of Punta Cabullones, groundwater and surface water are hydraulically coupled. Locally, surface-hypersaline water sinks into the aquifer, providing recharge and serving as a mechanism to redistribute salt throughout the study area. The evapotranspiration in the wetland subarea is estimated at about 11 million gallons per day (Mgal/d) that is equivalent to about 12,586 acre-feet per year. The balance of evapotranspiration, in excess of the about 0.5 Mgal/d of groundwater flow within the wetland, is supplied by saline to hypersaline surface water that may include seawater and meteoric water highly affected by evaporation with dissolved salts. In one of the extreme scenarios in which no groundwater is intercepted by pumpage at the Restaurada well field, the amount of saline to hypersaline water in the wetland consumed by evapotranspiration is about 10.5 Mgal/d. In the opposite extreme in which the entire regional groundwater flow is intercepted by pumpage in the Restaurada well field, the entire evapotranpiration requirement is met by saline to hypersaline water. Hydrologic, isotopic, and chemical data indicate that all of, or a large portion of, the historical groundwater flow to Punta Cabullones is being captured by the Puerto Rico Aqueducts and Sewer Authority pumpage at the Restaurada well field at a rate of about 2 Mgal/d. As a consequence, seawater intrusion into the aquifer at the Punta Cabullones area seems to be occurring, while the current pumpage at the Restaurada well field is sustained by storage depletion of the aquifer.

Puerto Rico↗

USGS2015

Groundwater quality in Geauga County, Ohio: status, including detection frequency of methane in water wells, 2009, and changes during 1978-2009

Domestic wells that are not safeguarded by regular water-quality testing provide drinking water for 79 percent of the residents of Geauga County, in northeastern Ohio. Since 1978, the U.S. Geological Survey (USGS) has worked cooperatively with the Board of Commissioners and Geauga County Planning Commission to monitor the quality of groundwater in four commonly used aquifers in county—the glacial deposits, the Pottsville Formation, the Cuyahoga Group, and the Berea Sandstone. A 33-percent growth in population from 1980 to 2009 increased the potential for humans to influence groundwater resources by withdrawing more groundwater, disposing of more human waste near the land surface, treating an expanded network of township roads with deicing salt, and likely using more solvents, pesticides, and other chemicals on the land surface than were used in preceding decades. To describe the status of groundwater quality in 2009 and its suitability for drinking, USGS personnel collected samples of water prior to treatment from 16 wells (mostly domestic) during June 9–19. The samples were analyzed for 92 properties and constituents, 41 of which had human-health benchmarks to which analytical results could be compared to evaluate suitability for drinking. Four of these benchmarks were exceeded at the following frequencies: arsenic (2 of 16 wells, 12.5 percent), total coliform bacteria (2 of 16 wells, 12.5 percent), fecal coliform bacteria (1 of 14 wells, 7 percent), and sodium (6 of 16 wells, 38 percent). No domestic wells sampled in 2009 exceeded the health-based benchmark of 300 micrograms per liter (µg/L) for manganese, although 5 of 65 wells (8 percent) sampled since 1978 have. Analyses from domestic wells were augmented with water-quality data from seven public-supply well fields that were obtained from the Ohio Environmental Protection Agency. These public-supply data were typically collected between 2000 and 2010 and represent water samples that were collected prior to treatment or that were treated by a method that does not effectively remove the constituents of interest. Similar to the domestic-well data, these data indicated that some samples from public-supply wells have also exceeded health-based benchmarks for arsenic and sodium, along with occasional exceedances of health-based benchmarks for cadmium and lead. Concentrations of nitrate, pesticides, and volatile organic compounds in ground-water samples from domestic and (or) public-supply wells were either considerably less than the human-health benchmarks for these constituents or were not detected. Water-quality data collected in 2009 were also compared to aesthetically based benchmarks developed by the U.S. Environmental Protection Agency, called Secondary Maximum Contaminant Levels (SMCLs). Iron and manganese most frequently exceeded SMCLs (in samples from 10 of 16 domestic wells and in untreated water from 3 of 4 public-supply well fields). To evaluate the frequency of methane detection in water wells in the county, the USGS sampled 16 wells across the county and screened the samples for combustible gas within the headspace (the air above the water in a closed container). Water from three (19 percent) of the wells contained detectable combustible gas (0.10 to 0.40 percent by volume). All three detections were from wells tapping the Cuyahoga Group or the Berea Sandstone, and all detections were less than the lower explosive limit of 5 percent by volume—the concentration at which methane in air can be flammable if an ignition source is present. Analyses of dissolved gas composition in water from these three wells showed methane concentrations ranging from 0.007 to 1.8 milligrams per liter (mg/L). The primary effect of human activities on groundwater quality found during this study is the input of salinity, or chloride, near land surface. On the basis of ratios of chloride to bromide, the main sources of chloride are road salt and septic leachate rather than oil-field brines (either spilled at land surface or sprayed on roads for dust control). The correlation of chloride concentration to distance of well from road for 31 wells in the county sampled by the U.S. Geological Survey in 1999 suggests that road salt is the dominant source of chloride. The majority of constituents exceeding health-based and aesthetically based benchmarks in groundwater were those that are naturally present in aquifer rocks and sediments rather than constituents introduced by human activities. Concentrations of such natural contaminants are controlled by geochemical processes in the subsurface, particularly by oxidation-reduction (redox) reactions. The categorization of redox conditions based on the water quality of 116 samples collected from 65 wells in Geauga County during 1978 through 2009 indicates that most groundwater samples were strongly reducing (60 percent) or oxic (18 percent). Oxic waters were found only in the Pottsville Formation and Berea Sandstone and were generally associated with nitrate at concentrations of 0.38 to 6.0 mg/L. Strongly reducing waters occurred in all four commonly used aquifers and were associated with the following naturally occurring contaminants: (1) arsenic and manganese at concentrations exceeding the health-based benchmarks (10 µg/L and 300 µg/L, respectively) in some samples, (2) iron and manganese at concentrations exceeding the aesthetically based standards (300 µg/L and 50 µg/L, respectively) in most samples, and (3) total sulfides (consisting of hydrogen sulfide gas with its characteristic rotten-egg odor and [or] iron sulfide minerals that appear as finely disseminated particulates in water). Because of the association of redox conditions with specific contaminants, attempts were made to further document spatially where oxic and strongly reducing conditions occur so that contaminant occurrence can be better anticipated by planners and well owners. Within the Pottsville Formation, wells tapping strongly reducing groundwater tended to have a greater thickness of overlying low-permeability (recharge-inhibiting) material such as clay and shale than other wells tapping oxic or nitrate-reducing groundwater. In the Berea Sandstone, oxic conditions were found at well locations where either depth to groundwater was shallow (less than 45 feet [ft] below land surface) or the measured water level was within the open interval (uncased portion) of the well, whereas strongly reducing groundwater was found at well locations where depths to water were greater than 60 ft below land surface and measured water levels were 15 ft or more above the open interval of the well. To evaluate whether constituent concentrations consistently increased or decreased over time, the strength of the association between sampling year (time) and constituent concentration was statistically evaluated for 116 water-quality samples collected by the USGS in 1978, 1980, 1986, 1999, and 2009 from a total of 65 wells across the county (generally domestic wells or wells serving small businesses or churches). Results indicate that many of the constituents that have been analyzed for decades exhibited no consistent temporal trends at a statistically significant level (p-value less than 0.05); fluctuations in concentrations of these constituents represent natural variation in groundwater quality. Dissolved oxygen, calcium, and sulfate concentrations and chloride:bromide ratios increased over time in one or more aquifers, while pH and concentrations of bromide and dissolved organic carbon decreased over time. Detections of total coliform bacteria and nitrate did not become more frequent from 1986 to 2009, even though potential sources of these constituents, such as number of septic systems (linked to population) and percent developed land in the county, increased during this period.

Geauga County↗

USGS2014

Linkage of the Soil and Water Assessment Tool and the Texas Water Availability Model to simulate the effects of brush management on monthly storage of Canyon Lake, south-central Texas, 1995-2010

The U.S. Geological Survey (USGS), in cooperation with the Texas State Soil and Water Conservation Board, developed and applied an approach to create a linkage between the published upper Guadalupe River Soil Water Assessment Tool (SWAT) brush-management (ashe juniper [ Juniperus ashei ]) model and the full authorization version Guadalupe River Water Availability Model (WAM). The SWAT model was published by the USGS, and the Guadalupe River WAM is available from the Texas Commission on Environmental Quality. The upper Guadalupe River watershed is a substantial component of the Guadalupe River WAM. This report serves in part as documentation of a proof of concept on the feasibility of linking these two water-resources planning models for the purpose of simulating possible increases in water storage in Canyon Lake as a result of different brush-management scenarios. The SWAT-WAM linkage for the upper Guadalupe River is documented with a principal objective to evaluate the distributional characteristics of the monthly water storage of Canyon Lake during selected drought conditions. Focus is on the relative evaluation of select scenarios of large-scale or “extensive” brush management within the upper Guadalupe River watershed. There are six SWAT simulations for the upper Guadalupe River watershed that include a baseline (0-percent management of treatable ashe juniper, the baseline scenario from a previous study in which no percentage of ashe juniper is numerically replaced with grassland) along with five scenarios (extensions of SWAT simulations from a previous study) of 20-, 40-, 60-, 80-, and 100-percent random (numerical) replacement of treatable ashe juniper with grasslands throughout the upper Guadalupe River watershed in south-central Texas. SWAT is a process-based, semidistributed, water-balance model designed to predict the effects of landscape management decisions on water yields. A watershed is subdivided into subbasins, and each subbasin is associated with a single reach on the stream network. In general a WAM, such as the Guadalupe River WAM, provides analysis of generalized water rights in a river and reservoir framework. A WAM accommodates hydrology and water usage through several input files containing water rights, watershed parameters, and naturalized streamflow time series. A WAM is generalized for application to rivers and reservoir systems, and input datasets are uniquely developed for a river basin of concern. The extractions of SWAT output for the five extensive brush-management and baseline scenarios were offset by –21 years and, in general, the results were then mapped to the WAM input-flow file. The offset of –21 years was chosen arbitrarily for technical reasons and means that the period of monthly record 1995–2010 of the upper Guadalupe River SWAT became the synthetic period of monthly record 1974–89, hereinafter 1974–89 (synthetic) period, of the Guadalupe River WAM. The relative (between scenario to baseline) effects of extensive brush-management scenarios by using the SWAT-WAM linkage were evaluated, and two critical intermediate results were total inflow to Canyon Lake from 1995 to 2010 and the monthly storage of Canyon Lake from 1974 to 1989 (synthetic). The first quartile or lower 25th percentile of monthly storage of Canyon Lake for the baseline scenario is 381,000 acre-feet (acre-ft) for the hereinafter 1974–89 (synthetic) period. This lower quartile was chosen for analysis for two critical purposes. First, Canyon Lake is managed with a conservation pool of about 386,200 acre-ft capacity (as recognized by the WAM) and is at or near conservation capacity about 50 percent or more of the time; further, there is intrinsic data censoring that occurs for the monthly storage distribution because Canyon Lake is at or near conservation pool elevation the majority of the time. This intrinsic censoring has the effect of creating a bounded distribution with a left or low-volume tail. Statistical assessment of the brush-management scenarios beginning with the 381,000 acre-ft censoring threshold provides readily interpretable results. Second, the quantification of brush management during periods lacking abundant rainfall, which were defined in this study as months for which Canyon Lake storage was below the 25th percentile for the simulation period, are of substantial interest to water-resource managers and stakeholders in the context of water-supply enhancement. A statistical assessment of the SWAT-WAM linkage for the low-volume tail of the distribution of monthly storage of Canyon Lake is the focus of analysis and interpretation. Drought periods for the analysis are defined as the months (consecutive or not) during which Canyon Lake is below the 25th percentile of storage (381,000 acre-ft) for the baseline scenario. Such months are referred to as being within the “Drought Quartile.” The Drought Quartile is a conceptual and heuristically determined waypoint for the analysis and is not related to any administrative definition of drought by stakeholders or policy makers. The five scenarios and the baseline scenario simulated in the upper Guadalupe River SWAT were all passed through the Guadalupe River WAM by the SWAT-WAM linkage described in this report. A comparison of the mean increase per month in reservoir storage for Canyon Lake conditioned for the Drought Quartile was made. For each of the five brush-management and baseline scenarios, the months with storage below 381,000 acre-ft were extracted. The mean monthly storages during the Drought Quartile were computed for each of the five scenarios and the baseline scenario. The mean of the baseline scenario was 376,458 acre-ft and subsequently was subtracted from the mean monthly storage during the Drought Quartile for each of the five scenarios. The mean monthly offset storages of Canyon Lake during the Drought Quartile were 110 acre-ft (20 percent); 448 acre-ft (40 percent); 754 acre-ft (60 percent); 1,080 acre-ft (80 percent); and 1,090 acre-ft (100 percent). A particular mean was interpreted as follows: the value of 754 acre-ft for the 60-percent brush-management scenario implies that, on average, this scenario indicates an additional 754 acre-ft per month of storage in Canyon Lake relative to the baseline during the Drought Quartile. All of the five scenarios resulted in an increase on average to water supply relative to the baseline scenario during the Drought Quartile through the SWAT-WAM linkage.

Texas↗

USGS2014

Quality of groundwater in the Denver Basin aquifer system, Colorado, 2003-5

Groundwater resources from alluvial and bedrock aquifers of the Denver Basin are critical for municipal, domestic, and agricultural uses in Colorado along the eastern front of the Rocky Mountains. Rapid and widespread urban development, primarily along the western boundary of the Denver Basin, has approximately doubled the population since about 1970, and much of the population depends on groundwater for water supply. As part of the National Water-Quality Assessment Program, the U.S. Geological Survey conducted groundwater-quality studies during 2003–5 in the Denver Basin aquifer system to characterize water quality of shallow groundwater at the water table and of the bedrock aquifers, which are important drinking-water resources. For the Denver Basin, water-quality constituents of concern for human health or because they might otherwise limit use of water include total dissolved solids, fluoride, sulfate, nitrate, iron, manganese, selenium, radon, uranium, arsenic, pesticides, and volatile organic compounds. For the water-table studies, two monitoring-well networks were installed and sampled beneath agricultural (31 wells) and urban (29 wells) land uses at or just below the water table in either alluvial material or near-surface bedrock. For the bedrock-aquifer studies, domestic- and municipal-supply wells completed in the bedrock aquifers were sampled. The bedrock aquifers, stratigraphically from youngest (shallowest) to oldest (deepest), are the Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers. The extensive dataset collected from wells completed in the bedrock aquifers (79 samples) provides the opportunity to evaluate factors and processes affecting water quality and to establish a baseline that can be used to characterize future changes in groundwater quality. Groundwater samples were analyzed for inorganic, organic, isotopic, and age-dating constituents and tracers. This report discusses spatial and statistical distributions of chemical constituents and evaluates natural and human-related processes that affect water quality. Findings are synthesized to assess the vulnerability of the Denver Basin aquifer system to groundwater contamination. The chemistry of groundwater samples collected from the water-table wells was generally different from that of samples collected from the bedrock-aquifer wells. Samples from the water-table wells tended to have higher concentrations of total dissolved solids and most major ions. Concentrations of several constituents with potential human-health concerns, including nitrate, selenium, uranium, and arsenic, decreased with depth and were highest in samples from the water-table wells. Exceedances of drinking-water standards and water-quality benchmarks were more frequently associated with shallow groundwater samples; concentrations of total dissolved solids and sulfate exceeded water-quality benchmarks for about half or more of samples from the water-table wells. The sediments and rocks of the Denver Basin are natural sources of the trace elements selenium, uranium, and arsenic, which affect their concentrations in groundwater. Detections of organic contaminants, which are typically indicative of human sources of contamination to groundwater, were more frequent in samples from the water-table wells. Pesticide compounds and volatile organic compounds were detected in 33 and 62 percent, respectively, of water-table well samples. Detected organic contaminant concentrations were much less than the associated drinking-water standards. Samples collected from the bedrock aquifers had lower concentrations of total dissolved solids than did samples collected from the water-table wells, although within the bedrock-aquifer samples, concentrations increased from the Dawson to Denver to Arapahoe to Laramie-Fox Hills aquifers. Concentrations of total dissolved solids and many constituents varied spatially and with depth in the bedrock aquifers, likely as a result of ion-exchange and oxidation-reduction reactions, which are important processes affecting water quality. Major-ion chemistry generally evolved from a calcium-bicarbonate to calcium-sulfate composition, with some sodium-bicarbonate and sodium-sulfate facies in the deeper bedrock aquifers, likely resulting from longer residence times and more extensive water-rock interaction. Oxidation-reduction conditions generally evolved from oxic at the water table to anoxic with increasing depth in the bedrock aquifers. Most samples from the bedrock aquifers were anoxic. Exceedances of drinking-water standards and water-quality benchmarks for the bedrock aquifers occurred in 1 percent or less of samples for nitrate, selenium, or arsenic; there were no exceedances for uranium. Exceedances for total dissolved solids, sulfate, manganese, and iron were generally between about 10 and 20 percent for the bedrock-aquifer samples. Radon concentrations, which were only measured in samples collected from two of the bedrock aquifers, exceeded the lower proposed drinking-water standard for more than 90 percent of samples but exceeded the higher alternative standard for less than 5 percent of samples. Pesticide compounds and volatile organic compounds were detected in 3 and 22 percent, respectively, of bedrock-aquifer samples, all at concentrations that were that were much less than drinking-water standards. Water-quality data were synthesized to evaluate factors that affect spatial and depth variability in water quality and to assess aquifer vulnerability to contaminants from geologic materials and those of human origin. The quality of shallow groundwater in the alluvial aquifer and shallow bedrock aquifer system has been adversely affected by development of agricultural and urban areas. Land use has altered the pattern and composition of recharge. Increased recharge from irrigation water has mobilized dissolved constituents and increased concentrations in the shallow groundwater. Concentrations of most constituents associated with poor or degraded water quality in shallow groundwater decreased with depth; many of these constituents are not geochemically conservative and are affected by geochemical reactions such as oxidation-reduction reactions. Groundwater age tracers provide additional insight into aquifer vulnerability and help determine if young groundwater of potentially poor quality has migrated to deeper parts of the bedrock aquifers used for drinking-water supply. Age-tracer results were used to group samples into categories of young, mixed, and old groundwater. Groundwater ages transitioned from mostly young in the water-table wells to mostly mixed in the shallowest bedrock aquifer, the Dawson aquifer, to mostly old in the deeper bedrock aquifers. Although the bedrock aquifers are mostly old groundwater of good water quality, several lines of evidence indicate that young, contaminant-bearing recharge has reached shallow to moderate depths in some areas of the bedrock aquifers. The Dawson aquifer is the most vulnerable of the bedrock aquifers to contamination, but results indicate that the older (deeper) bedrock aquifers are also vulnerable to groundwater contamination and that mixing with young recharge has occurred in some areas. Heavy pumping has caused water-level declines in the bedrock aquifers in some parts of the Denver Basin, which has the potential to enhance the transport of contaminants from overlying units. Results of this study are consistent with the existing conceptual understanding of aquifer processes and groundwater issues in the Denver Basin and add new insight into the vulnerability of the bedrock aquifers to groundwater contamination.

Colorado↗

USGS2015

Chloride concentrations, loads, and yields in four watersheds along Interstate 95, southeastern Connecticut, 2008-11: factors that affect peak chloride concentrations during winter storms

Chloride (Cl - ) concentrations and loads and other water chemistry characteristics were assessed to evaluate potential effects of road-deicer applications on streamwater quality in four watersheds along Interstate 95 (I–95) in southeastern Connecticut from November 1, 2008, through September 30, 2011. Streamflow and water quality were studied in the Four Mile River, Oil Mill Brook, Stony Brook, and Jordan Brook watersheds, where developed land ranged from 9 to 32 percent. Water-quality samples were collected and specific conductance was measured continuously at paired water-quality monitoring sites, upstream and downstream from I–95. Specific conductance values were related to Cl - concentrations to assist in determining the effects of road-deicing operations on the levels of Cl - in the streams. Streamflow and water-quality data were compared with weather data and with the timing, amount, and composition of deicers applied to State highways. Grab samples were collected during winter stormwater-runoff events, such as winter storms or periods of rain or warm temperatures in which melting takes place. Grab samples were also collected periodically during the spring and summer and during base-flow conditions. The estimated Cl - concentrations at the eight water-quality monitoring sites during winter storms peaked as high as 270 milligrams per liter (mg/L) and were well below the U.S. Environmental Protection Agency (EPA) recommended acute chloride toxicity criterion of 860 mg/L and the chronic 4-day average toxicity criterion of 230 mg/L. Cl - concentrations in streams, particularly at sites downstream from I–95, peaked during increased streamflow in the winter and early spring as a result of deicers applied to roads and washed off by stormwater or meltwater. Cl - concentrations during most of the nonwinter seasons decreased during increases in streamflow because storm runoff was more dilute than base flow. However, peaks in specific conductance and estimated chloride concentrations at streams in more urbanized watersheds corresponded to peaks in streamflow well after winter snow or ice events; these delayed peaks in Cl - concentration likely resulted from deicer residue that remained in melting snow piles and on roadsides and (or) that were flushed from soils and shallow groundwater, then discharged downstream. Estimated peak Cl - concentrations varied with the type of winter storm event and were highest during or after winter storms of frozen precipitation and rain, in which the rain or meltwater effectively washed off the deicers. Estimated peak Cl - concentrations correlated positively with the duration of deicer application but generally not with streamflow. Estimated peak Cl - concentrations during the winter season were highest during low streamflow at most sites. Chloride concentrations varied considerably in shallow groundwater as a result of land-use differences. Cl - concentrations were very high (as high as 800 mg/L) in shallow groundwater downstream from I–95 at the Four Mile River site. Chloride/bromide mass concentration ratios and the proximity of a former landfill and sewage lagoon upstream indicate a likely source of Cl - is landfill leachate and possibly sewage leachate. Cl - loads in streams generally were highest in the winter and early spring. The estimated daily Cl - yield for the four monitoring sites downstream from I–95 ranged from 0.0004 ton per day per square mile for one of the least developed watersheds to 0.052 ton per day per square mile for the watershed with the highest percentage of urban development and impervious area. The estimated median contribution of Cl - load from atmospheric deposition was small and ranged from 0.07 percent of Cl - load at the Jordan Brook watershed to 0.57 percent at the Oil Mill Brook watershed. The Cl - loads in streams (outputs) were compared with Cl - load inputs, which include atmospheric deposition and deicer applications; Cl - load inputs were slightly larger than the Cl - load outputs at most of the sites during most years but do not account for the Cl - load in groundwater leaving the watersheds. A multiple linear regression model was developed to describe the variability of the natural log of peak specific conductance, as well as estimated Cl - concentrations. Five significant variables best explained the variability in the natural log of the peak specific conductance after deicing events: (1) snow on ground before deicing event; (2) winter precipitation with rain; (3) specific conductance in base flow; (4) State-operated road lane miles divided by watershed area; and (5) amount of Cl - from deicers applied to State-operated roads per lane mile. In this report, winter precipitation is defined as any type of precipitation, including frozen precipitation and rain, that occurs during the active deicing season, typically November through March. Frozen precipitation is defined here as snow, sleet, freezing rain, or any winter precipitation except rain. The addition of a lane mile in both directions on I–95 would result in an estimate of approximately 2 to 11 percent increase in Cl - input from deicers applied to I–95 and other roads maintained by Connecticut Department of Transportation. The largest estimated increase in Cl - load was in the watersheds with the greatest number miles of I–95 corridor relative to the total lane miles maintained by Connecticut Department of Transportation. On the basis of these estimates and the estimated peak Cl - concentrations during the study period, it is unlikely that the increased use of deicers on the additional lanes would lead to Cl - concentrations that exceed the aquatic habitat criteria.

Connecticut↗

USGS2013

Computed statistics at streamgages, and methods for estimating low-flow frequency statistics and development of regional regression equations for estimating low-flow frequency statistics at ungaged locations in Missouri

The weather and precipitation patterns in Missouri vary considerably from year to year. In 2008, the statewide average rainfall was 57.34 inches and in 2012, the statewide average rainfall was 30.64 inches. This variability in precipitation and resulting streamflow in Missouri underlies the necessity for water managers and users to have reliable streamflow statistics and a means to compute select statistics at ungaged locations for a better understanding of water availability. Knowledge of surface-water availability is dependent on the streamflow data that have been collected and analyzed by the U.S. Geological Survey for more than 100 years at approximately 350 streamgages throughout Missouri. The U.S. Geological Survey, in cooperation with the Missouri Department of Natural Resources, computed streamflow statistics at streamgages through the 2010 water year, defined periods of drought and defined methods to estimate streamflow statistics at ungaged locations, and developed regional regression equations to compute selected streamflow statistics at ungaged locations. Streamflow statistics and flow durations were computed for 532 streamgages in Missouri and in neighboring States of Missouri. For streamgages with more than 10 years of record, Kendall’s tau was computed to evaluate for trends in streamflow data. If trends were detected, the variable length method was used to define the period of no trend. Water years were removed from the dataset from the beginning of the record for a streamgage until no trend was detected. Low-flow frequency statistics were then computed for the entire period of record and for the period of no trend if 10 or more years of record were available for each analysis. Three methods are presented for computing selected streamflow statistics at ungaged locations. The first method uses power curve equations developed for 28 selected streams in Missouri and neighboring States that have multiple streamgages on the same streams. Statistical estimates on one of these streams can be calculated at an ungaged location that has a drainage area that is between 40 percent of the drainage area of the farthest upstream streamgage and within 150 percent of the drainage area of the farthest downstream streamgage along the stream of interest. The second method may be used on any stream with a streamgage that has operated for 10 years or longer and for which anthropogenic effects have not changed the low-flow characteristics at the ungaged location since collection of the streamflow data. A ratio of drainage area of the stream at the ungaged location to the drainage area of the stream at the streamgage was computed to estimate the statistic at the ungaged location. The range of applicability is between 40- and 150-percent of the drainage area of the streamgage, and the ungaged location must be located on the same stream as the streamgage. The third method uses regional regression equations to estimate selected low-flow frequency statistics for unregulated streams in Missouri. This report presents regression equations to estimate frequency statistics for the 10-year recurrence interval and for the N-day durations of 1, 2, 3, 7, 10, 30, and 60 days. Basin and climatic characteristics were computed using geographic information system software and digital geospatial data. A total of 35 characteristics were computed for use in preliminary statewide and regional regression analyses based on existing digital geospatial data and previous studies. Spatial analyses for geographical bias in the predictive accuracy of the regional regression equations defined three low-flow regions with the State representing the three major physiographic provinces in Missouri. Region 1 includes the Central Lowlands, Region 2 includes the Ozark Plateaus, and Region 3 includes the Mississippi Alluvial Plain. A total of 207 streamgages were used in the regression analyses for the regional equations. Of the 207 U.S. Geological Survey streamgages, 77 were located in Region 1, 120 were located in Region 2, and 10 were located in Region 3. Streamgages located outside of Missouri were selected to extend the range of data used for the independent variables in the regression analyses. Streamgages included in the regression analyses had 10 or more years of record and were considered to be affected minimally by anthropogenic activities or trends. Regional regression analyses identified three characteristics as statistically significant for the development of regional equations. For Region 1, drainage area, longest flow path, and streamflow-variability index were statistically significant. The range in the standard error of estimate for Region 1 is 79.6 to 94.2 percent. For Region 2, drainage area and streamflow variability index were statistically significant, and the range in the standard error of estimate is 48.2 to 72.1 percent. For Region 3, drainage area and streamflow-variability index also were statistically significant with a range in the standard error of estimate of 48.1 to 96.2 percent. Limitations on the use of estimating low-flow frequency statistics at ungaged locations are dependent on the method used. The first method outlined for use in Missouri, power curve equations, were developed to estimate the selected statistics for ungaged locations on 28 selected streams with multiple streamgages located on the same stream. A second method uses a drainage-area ratio to compute statistics at an ungaged location using data from a single streamgage on the same stream with 10 or more years of record. Ungaged locations on these streams may use the ratio of the drainage area at an ungaged location to the drainage area at a streamgage location to scale the selected statistic value from the streamgage location to the ungaged location. This method can be used if the drainage area of the ungaged location is within 40 to 150 percent of the streamgage drainage area. The third method is the use of the regional regression equations. The limits for the use of these equations are based on the ranges of the characteristics used as independent variables and that streams must be affected minimally by anthropogenic activities.

Missouri↗

USGS2015

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

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

Nebraska↗

USGS2016

Effects of water-supply reservoirs on streamflow in Massachusetts

State and local water-resource managers need modeling tools to help them manage and protect water-supply resources for both human consumption and ecological needs. The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, has developed a decision-support tool to estimate the effects of reservoirs on natural streamflow. The Massachusetts Reservoir Simulation Tool is a model that simulates the daily water balance of a reservoir. The reservoir simulation tool provides estimates of daily outflows from reservoirs and compares the frequency, duration, and magnitude of the volume of outflows from reservoirs with estimates of the unaltered streamflow that would occur if no dam were present. This tool will help environmental managers understand the complex interactions and tradeoffs between water withdrawals, reservoir operational practices, and reservoir outflows needed for aquatic habitats. A sensitivity analysis of the daily water balance equation was performed to identify physical and operational features of reservoirs that could have the greatest effect on reservoir outflows. For the purpose of this report, uncontrolled releases of water (spills or spillage) over the reservoir spillway were considered to be a proxy for reservoir outflows directly below the dam. The ratio of average withdrawals to the average inflows had the largest effect on spillage patterns, with the highest withdrawals leading to the lowest spillage. The size of the surface area relative to the drainage area of the reservoir also had an effect on spillage; reservoirs with large surface areas have high evaporation rates during the summer, which can contribute to frequent and long periods without spillage, even in the absence of water withdrawals. Other reservoir characteristics, such as variability of inflows, groundwater interactions, and seasonal demand patterns, had low to moderate effects on the frequency, duration, and magnitude of spillage. The reservoir simulation tool was used to simulate 35 single- and multiple-reservoir systems in Massachusetts over a 44-year period (water years 1961 to 2004) under two water-use scenarios. The no-pumping scenario assumes no water withdrawal pumping, and the pumping scenario incorporates average annual pumping rates from 2000 to 2004. By comparing the results of the two scenarios, the total streamflow alteration can be parsed into the portion of streamflow alteration caused by the presence of a reservoir and the additional streamflow alteration caused by the level of water use of the system. For each reservoir system, the following metrics were computed to characterize the frequency, duration, and magnitude of reservoir outflow volumes compared with unaltered streamflow conditions: (1) the median number of days per year in which the reservoir did not spill, (2) the median duration of the longest consecutive period of no-spill days per year, and (3) the lowest annual flow duration exceedance probability at which the outflows are significantly different from estimated unaltered streamflow at the 95-percent confidence level. Most reservoirs in the study do not spill during the summer months even under no-pumping conditions. The median number of days during which there was no spillage was less than 365 for all reservoirs in the study, indicating that, even under reported pumping conditions, the reservoirs refill to full volume and spill at least once during nondrought years, typically in the spring. Thirteen multiple-reservoir systems consisting of two or three hydrologically connected reservoirs were included in the study. Because operating rules used to manage multiple-reservoir systems are not available, these systems were simulated under two pumping scenarios, one in which water transfers between reservoirs are minimal and one in which reservoirs continually transferred water to intermediate or terminal reservoirs. These two scenarios provided upper and lower estimates of spillage under average pumping conditions from 2000 to 2004. For sites with insufficient data to simulate daily water balances, a proxy method to estimate the three spillage metrics was developed. A series of 4,000 Monte Carlo simulations of the reservoir water balance were run. In each simulation, streamflow, physical reservoir characteristics, and daily climate inputs were randomly varied. Tobit regression equations that quantify the relation between streamflow alteration and physical and operational characteristics of reservoirs were developed from the results of the Monte Carlo simulations and can be used to estimate each of the three spillage metrics using only the withdrawal ratio and the ratio of the surface area to the drainage area, which are available statewide for all reservoirs. A graphical user-interface for the Massachusetts Reservoir Simulation Tool was developed in a Microsoft Access environment. The simulation tool contains information for 70 reservoirs in Massachusetts and allows for simulation of additional scenarios than the ones considered in this report, including controlled releases, dam seepage and leakage, demand management plans, and alternative water withdrawal and transfer rules.

Massachusetts↗

USGS2017

Sequence stratigraphy, seismic stratigraphy, and seismic structures of the lower intermediate confining unit and most of the Floridan aquifer system, Broward County, Florida

Deep well injection and disposal of treated wastewater into the highly transmissive saline Boulder Zone in the lower part of the Floridan aquifer system began in 1971. The zone of injection is a highly transmissive hydrogeologic unit, the Boulder Zone, in the lower part of the Floridan aquifer system. Since the 1990s, however, treated wastewater injection into the Boulder Zone in southeastern Florida has been detected at three treated wastewater injection utilities in the brackish upper part of the Floridan aquifer system designated for potential use as drinking water. At a time when usage of the Boulder Zone for treated wastewater disposal is increasing and the utilization of the upper part of the Floridan aquifer system for drinking water is intensifying, there is an urgency to understand the nature of cross-formational fluid flow and identify possible fluid pathways from the lower to upper zones of the Floridan aquifer system. To better understand the hydrogeologic controls on groundwater movement through the Floridan aquifer system in southeastern Florida, the U.S. Geological Survey and the Broward County Environmental Planning and Community Resilience Division conducted a 3.5-year cooperative study from July 2012 to December 2015. The study characterizes the sequence stratigraphy, seismic stratigraphy, and seismic structures of the lower part of the intermediate confining unit aquifer and most of the Floridan aquifer system. Data obtained to meet the study objective include 80 miles of high-resolution, two-dimensional (2D), seismic-reflection profiles acquired from canals in eastern Broward County. These profiles have been used to characterize the sequence stratigraphy, seismic stratigraphy, and seismic structures in a 425-square-mile study area. Horizon mapping of the seismic-reflection profiles and additional data collection from well logs and cores or cuttings from 44 wells were focused on construction of three-dimensional (3D) visualizations of eight sequence stratigraphic cycles that compose the Eocene to Miocene Oldsmar, Avon Park, and Arcadia Formations. The mapping of these seismic-reflection and well data has produced a refined Cenozoic sequence stratigraphic, seismic stratigraphic, and hydrogeologic framework of southeastern Florida. The upward transition from the Oldsmar Formation to the Avon Park Formation and the Arcadia Formation embodies the evolution from (1) a tropical to subtropical, shallow-marine, carbonate platform, represented by the Oldsmar and Avon Park Formations, to (2) a broad, temperate, mixed carbonate-siliciclastic shallow marine shelf, represented by the lower part of the Arcadia Formation, and to (3) a temperate, distally steepened carbonate ramp represented by the upper part of the Arcadia Formation. In the study area, the depositional sequences and seismic sequences have a direct correlation with hydrogeologic units. The approximate upper boundary of four principal permeable units of the Floridan aquifer system (Upper Floridan aquifer, Avon Park permeable zone, uppermost major permeable zone of the Lower Floridan aquifer, and Boulder Zone) have sequence stratigraphic and seismic-reflection signatures that were identified on cross sections, mapped, or both, and therefore the sequence stratigraphy and seismic stratigraphy were used to guide the development of a refined spatial representation of these hydrogeologic units. In all cases, the permeability of the four permeable units is related to stratiform megaporosity generated by ancient dissolution of carbonate rock associated with subaerial exposure and unconformities at the upper surfaces of carbonate depositional cycles of several hierarchical scales ranging from high-frequency cycles to depositional sequences. Additionally, interparticle porosity also contributes substantially to the stratiform permeability in much of the Upper Floridan aquifer. Information from seismic stratigraphy allowed 3D geomodeling of hydrogeologic units—an approach never before applied to this area. Notably, the 3D geomodeling provided 3D visualizations and geocellular models of the depositional sequences, hydrostratigraphy, and structural features. The geocellular data could be used to update the hydrogeologic structure inherent to groundwater flow simulations that are designed to address the sustainability of the water resources of the Floridan aquifer system. Two kinds of pathways that could enable upward cross-formational flow of injected treated wastewater from the Boulder Zone have been identified in the 80 miles of high-resolution seismic data collected for this study: a near-vertical reverse fault and karst collapse structures. The single reverse fault, inferred to be of tectonic origin, is in extreme northeastern Broward County and has an offset of about 19 feet at the level of the Arcadia Formation. Most of the 17 karst collapse structures identified manifest as columniform, vertically stacked sagging seismic reflections that span early Eocene to Miocene age rocks equivalent to much of the Floridan aquifer system and the lower part of the overlying intermediate confining unit. In some cases, the seismic-sag structures extend upward into strata of Pliocene age. The seismic-sag structures are interpreted to have a semicircular shape in plan view on the basis of comparison to (1) other seismic-sag structures in southeastern Florida mapped with two 2D seismic cross lines or 3D data, (2) comparison to these structures located in other carbonate provinces, and (3) plausible extensional ring faults detected with multi-attribute analysis. The seismic-sag structures in the study area have heights as great as 2,500 vertical feet, though importantly, one spans about 7,800 feet. Both multi-attribute analysis and visual detection of offset of seismic reflections within the seismic-sag structures indicate faults and fractures are associated with many of the structures. Multi-attribute analysis highlighting chimney fluid pathways also indicates that the seismic-sag structures have a high probability for potential vertical cross-formational fluid flow along the faulted and fractured structures. A collapse of the seismic-sag structures within a deep burial setting evokes an origin related to hypogenic karst processes by ascending flow of subsurface fluids. In addition, paleo-epigenic karst related to major regional subaerial unconformities within the Florida Platform generated collapse structures (paleo-sinkholes) that are much smaller in scale than the cross-formational seismic-sag structures.

Florida↗

USGS2011

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

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

California↗

USGS2014

Status and understanding of groundwater quality in the Klamath Mountains study unit, 2010: California GAMA Priority Basin Project

Groundwater quality in the Klamath Mountains (KLAM) study unit was investigated as part of the Priority Basin Project of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is located in Del Norte, Humboldt, Shasta, Siskiyou, Tehama, and Trinity 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 GAMA Priority Basin Project was designed to provide a spatially unbiased, statistically robust assessment of the quality of untreated (raw) groundwater in the primary aquifer system. The assessment is based on water-quality data and explanatory factors for groundwater samples collected in 2010 by the USGS from 39 sites and on water-quality data from the California Department of Public Health (CDPH) water-quality database. The primary aquifer system was defined by the depth intervals of the wells listed in the CDPH water-quality database for the KLAM study unit. The quality of groundwater in the primary aquifer system may be different from that in the shallower or deeper water-bearing zones; shallow groundwater may be more vulnerable to surficial contamination. This study included two types of assessments: (1) a status assessment , which characterized the status of the current quality of the groundwater resource by using data from samples analyzed for volatile organic compounds, pesticides, and naturally occurring 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 KLAM study unit, not the quality of treated drinking water delivered to consumers by water purveyors. Relative-concentrations (sample concentrations divided by the health- or aesthetic-based 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 greater than (>) 1.0 indicates a concentration greater than a benchmark, and a relative-concentration less than or equal to (≤) 1.0 indicates a concentration less than or equal to a benchmark. Relative-concentrations of organic constituents were classified as “high” (relative-concentration > 1.0), “moderate” (0.1 < relative-concentration ≤ 1.0), or “low” (relative-concentration ≤ 0.1). For inorganic constituents, the boundary between low and moderate relative-concentration was set at 0.5. 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 aquifer system with a relative-concentration greater than 1.0 for a particular constituent or class of constituents; percentage is based on an areal rather than a volumetric basis. Moderate and low aquifer-scale proportions were defined as the percentages of the primary aquifer system with moderate and low relative-concentrations, respectively. The KLAM study unit includes more than 8,800 square miles (mi 2 ), but only those areas near the sampling sites, about 920 mi 2 , are included in the areal assessment of the study unit. Two statistical approaches—grid-based and spatially weighted—were used to evaluate aquifer-scale proportions for individual constituents and classes of constituents. To confirm this methodology, 90 percent confidence intervals were calculated for the grid-based high aquifer-scale proportions and were compared to the spatially weighted results, which were found to be within these confidence intervals in all cases. Grid-based results were selected for use in the status assessment unless, as was observed in a few cases, a grid-based result was zero and the spatially weighted result was not zero, in which case, the spatially weighted result was used. The status assessment showed that inorganic constituents with human-health benchmarks were detected at high relative-concentrations in 2.6 percent of the primary aquifer system and at moderate relative-concentrations in 10 percent of the system. The high aquifer-scale proportion for inorganic constituents mainly reflected the high aquifer-scale proportions of boron. Inorganic constituents with secondary maximum contaminant levels were detected at high relative-concentrations in 13 percent of the primary aquifer system and at moderate relative-concentrations in 10 percent of the system. The constituents present at high relative-concentrations included iron and manganese. Organic constituents with human-health benchmarks were not detected at high relative-concentrations, but were detected at moderate relative-concentrations in 1.9 percent of the primary aquifer system. The 1.9 percent reflected a spatially weighted moderate aquifer-scale proportion for the gasoline additive methyl tert-butyl ether. Of the 148 organic constituents analyzed, 14 constituents were detected. Only one organic constituent had a detection frequency of greater than 10 percent—the trihalomethane, chloroform. The second component of this study, the understanding assessment , identified the natural and human factors that may have affected the groundwater quality in the KLAM study unit by evaluating statistical correlations between water-quality constituents and potential explanatory factors. The potential explanatory factors evaluated were aquifer lithology, land use, hydrologic conditions, depth, groundwater age, and geochemical conditions. Results of the statistical evaluations were used to explain the occurrence and distribution of constituents in the KLAM study unit. Groundwater age distribution (modern, mixed, or pre-modern), redox class (oxic, mixed, or anoxic), and dissolved oxygen concentration were the explanatory factors that best explained occurrence patterns of the inorganic constituents. High concentrations of boron were found to be associated with groundwater classified as mixed or pre-modern with respect to groundwater age. Boron was also negatively correlated to dissolved oxygen and positively correlated to specific conductance. Iron and manganese concentrations were strongly associated with low dissolved oxygen concentrations, anoxic and mixed redox classifications, and pre-modern groundwater. Specific conductance concentrations were found to be related to pre-modern groundwater, low dissolved oxygen concentrations, and high pH. Chloroform was selected for additional evaluation in the understanding assessment because it was detected in more than 10 percent of wells sampled in the KLAM study unit. Septic tank density was the only explanatory factor that was found to relate to chloroform concentrations.

California↗

USGS2013

California State Waters Map Series: Offshore of Santa Barbara, California

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within the 3-nautical-mile limit of California’s State Waters. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data, acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. The Offshore of Santa Barbara map area lies within the central Santa Barbara Channel region of the Southern California Bight. This geologically complex region forms a major biogeographic transition zone, separating the cold-temperate Oregonian province north of Point Conception from the warm-temperate California province to the south. The map area is in the southern part of the Western Transverse Ranges geologic province, which is north of the California Continental Borderland. Significant clockwise rotation—at least 90°—since the early Miocene has been proposed for the Western Transverse Ranges province, and geodetic studies indicate that the region is presently undergoing north-south shortening. Uplift rates (as much as 2.2 mm/yr) that are based on studies of onland marine terraces provide further evidence of significant shortening. The city of Santa Barbara, the main coastal population center in the map area, is part of a contiguous urban area that extends from Carpinteria to Goleta. This urban area was developed on the coalescing alluvial surfaces, uplifted marine terraces, and low hills that lie south of the east-west-trending Santa Ynez Mountains. Several beaches line the actively utilized Santa Barbara coastal zone, including Arroyo Burro Beach Park, Leadbetter Beach, East Beach, and “Butterfly Beach.” There are ongoing coastal erosion problems associated with both development and natural processes; between 1933–1934 and 1998, cliff erosion in the map area occurred at rates of about 0.1 to 1 m/yr, the largest amount (63 m) occurring at Arroyo Burro in the western part of the map area. In addition, development of the Santa Barbara Harbor, which began in 1928, lead to shoaling west of the harbor as the initial breakwater trapped sand, as well as to coastal erosion east of the harbor. Since 1959, annual harbor dredging has mitigated at least some of the downcoast erosion problems. The Offshore of Santa Barbara map area lies in the central part of the Santa Barbara littoral cell, which is characterized by littoral drift to the east-southeast. Drift rates have been estimated to be about 400,000 tons/yr at Santa Barbara Harbor. Sediment supply to the western and central parts of the littoral cell, including the map area, is largely from relatively small transverse coastal watersheds. Within the map area, these coastal watersheds include (from east to west) San Ysidro Creek, Oak Creek, Montecito Creek, Sycamore Creek, Mission Creek, Arroyo Burro, and Atascadero Creek. The Ventura and Santa Clara Rivers, the mouths of which are about 40 to 50 km southeast of Santa Barbara, are much larger sediment sources. Still farther east, eastward-moving sediment in the littoral cell is trapped by Hueneme and Mugu Canyons and then transported to the deep-water Santa Monica Basin. The offshore part of the map area consists of a relatively flat and shallow continental shelf, which dips gently seaward (about 0.4° to 0.8°) so that water depths at the 3-nautical-mile limit of California’s State Waters are about 45 m in the east and about 75 m in the west. This part of the Santa Barbara Channel is relatively well protected from large Pacific swells from the north and northwest by Point Conception and from the south and southwest by offshore islands and banks. The shelf is underlain by variable amounts of upper Quaternary shelf, estuarine, and fluvial sediments deposited as sea level fluctuated in the late Pleistocene. Seafloor habitats in the broad Santa Barbara Channel region consist of significant amounts of soft sediment and isolated areas of rocky habitat that support kelp-forest communities nearshore and rocky-reef communities in deep water. The potential marine benthic habitat types mapped in the Offshore of Santa Barbara map area are directly related to its Quaternary geologic history, geomorphology, and active sedimentary processes. These potential habitats, which lie within the Shelf (continental shelf) megahabitat, range from soft, unconsolidated sediment to hard sedimentary bedrock. This heterogeneous seafloor provides promising habitat for rockfish, groundfish, crabs, shrimp, and other marine benthic organisms.

California↗

USGS2019

Determination of study reporting limits for pesticide constituent data for the California Groundwater Ambient Monitoring and Assessment Program Priority Basin Project, 2004–2018—Part 1: National Water Quality Schedules 2003, 2032, or 2033, and 2060

The California Groundwater Ambient Monitoring and Assessment Program Priority Basin Project (GAMA-PBP) is a long-term cooperative project designed to assess the quality of groundwater resources used for public and domestic drinking water supplies in the State of California, to monitor and evaluate changes to that quality, to investigate the human and natural factors controlling water quality, and to improve the availability of comprehensive groundwater quality data and information. Between May 18, 2004, and May 3, 2018, the GAMA-PBP collected 3001 groundwater samples for analysis of pesticide constituents by the U.S. Geological Survey (USGS) National Water Quality Laboratory (NWQL)(note that ‘pesticide constituents’ includes parent compounds and degradates). Of these samples, 2994 were analyzed for pesticide constituents on schedules 2003, 2032, or 2033 (65 to 84 constituents), and 840 were analyzed for pesticide constituents on schedule 2060 (58 constituents). The original dataset reported by the NWQL to the USGS National Water Information System (NWIS) database contained a total of 2,688 detections of 78 pesticide constituents and 253,825 non-detections. In this original dataset, 33 percent of the 3,001 samples analyzed had reported detections of one or more pesticide constituents. This report describes the GAMA-PBP data-quality objectives for pesticide data, the procedures used to establish study reporting limits, and use of those reporting limits to censor the data from the NWQL so that the final data published by the GAMA-PBP meet these data-quality objectives. The final GAMA-PBP dataset for samples collected from May 2004 to May 2018, after censoring, had a total of 1,632 detections of 37 pesticide constituents. In the final GAMA-PBP dataset, 25 percent of the 3,001 samples analyzed had detections of one or more pesticide constituents. The presence of pesticides in groundwater is commonly evaluated by calculating detection frequencies. Detection frequencies for pesticides are sensitive to detection limits and method performance for concentrations near those limits; therefore, the two primary data quality issues addressed in the GAMA-PBP data-quality objectives for pesticides are (1) establishing criteria for classifying data from the laboratory as detections or non-detections for the purpose of data reporting by the project and (2) accounting for changes in analytical methods or method performance over time. The GAMA-PBP addresses these issues by developing study reporting limits that are used as the boundary between detections and non-detections for the reporting of GAMA-PBP results. These reporting limits are defined from method detection limits (MDLs) provided by the NWQL, unless examination of results from laboratory set blanks (LSBs) and GAMA-PBP field blanks indicates that a higher concentration censoring limit is warranted. The GAMA-PBP selected the MDL as the primary choice for defining study reporting limits for consistency with U.S. Environmental Protection Agency (EPA) guidelines for reporting detections of pesticides and other organic constituents. A five-step procedure is used to develop study reporting limits and censor the GAMA-PBP dataset accordingly. The effect of the censoring at each step is described to provide information about the relative effect of each step on the overall censoring of the dataset. Steps 1 and 2 can be implemented at the time the data are received, whereas steps 3−5 require information accumulated over an extended period. Step 1: Reject results that were most likely the result of specific contamination instances attributable to unusual field or laboratory conditions during sample collection or processing. Two such instances were identified, leading to rejection of 25 detections, which were assigned a data-quality indicator code of “Q” for “reviewed and rejected” in the NWIS database. Step 2: Use the NWQL MDLs in effect at the time each sample was analyzed as the reporting limit. A total of 506 detections were censored on this basis. Step 3: Use the maximum MDL established by the NWQL during July 2004–August 2018 (MDLmax) as the reporting limit. The rationale for using the MDLmax as the reporting limit is based primarily on the observation that the concentrations of MDLs generally increased over time. A total of 438 detections were censored on this basis. Step 4: Use the LSBs to identify periods of greater potential laboratory contamination bias and define raised reporting limits to be used during those periods. These periods were defined by using a moving average detection frequency approach. For consistency with the NWQL procedures for defining raised reporting limits on the basis of detections in LSBs, the raised reporting limits were defined as equal to three times the highest concentration measured in an LSB during the period. A total of 25 detections in groundwater samples analyzed during periods of increased laboratory contamination bias were censored. Step 5: Use the LSBs and field blanks to identify potential contamination bias from field or laboratory processes outside of the time periods identified in step 4. The NWQL protocols were used to define the MDLs from blanks analyzed outside of the periods identified in step 4. If an MDL defined from blanks was greater than the MDLmax, the MDL defined from blanks was used to censor the data. One constituent had a study reporting limit defined on this basis, and a total of 62 detections in groundwater samples were censored. As of 2019, the USGS NWIS database does not have the capability to store both the original value reported by the NWQL and the final value published by the GAMA-PBP that reflects application of the quality-control censoring described in this report. In the interim, while this capability is developed, the 1,031 results censored in steps 2−5 are blocked from public release in NWIS, and the GAMA-PBP has published the original and final values in a USGS data release accompanying this report. The entire GAMA-PBP final dataset for pesticide constituents on schedules 2003, 2032, or 2033, or on schedule 2060 is publicly available in that USGS data release, through the USGS GAMA-PBP public web portal, and through the California State Water Resources Control Board GAMA public groundwater information system.

California↗