SEARCH · Search USGS

Results for “Scientific Investigations Report”

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,423 records · Page 79Linked to original sources

USGS2011

Use of instantaneous streamflow measurements to improve regression estimates of index flow for the summer month of lowest streamflow in Michigan

In Michigan, index flow Q50 is a streamflow characteristic defined as the minimum of median flows for July, August, and September. The state of Michigan uses index flow estimates to help regulate large (greater than 100,000 gallons per day) water withdrawals to prevent adverse effects on characteristic fish populations. At sites where long-term streamgages are located, index flows are computed directly from continuous streamflow records as GageQ50. In an earlier study, a multiple-regression equation was developed to estimate index flows IndxQ50 at ungaged sites. The index equation explains about 94 percent of the variability of index flows at 147 (index) streamgages by use of six explanatory variables describing soil type, aquifer transmissivity, land cover, and precipitation characteristics. This report extends the results of the previous study, by use of Monte Carlo simulations, to evaluate alternative flow estimators, DiscQ50, IntgQ50, SiteQ50, and AugmQ50. The Monte Carlo simulations treated each of the available index streamgages, in turn, as a miscellaneous site where streamflow conditions are described by one or more instantaneous measurements of flow. In the simulations, instantaneous flows were approximated by daily mean flows at the corresponding site. All estimators use information that can be obtained from instantaneous flow measurements and contemporaneous daily mean flow data from nearby long-term streamgages. The efficacy of these estimators was evaluated over a set of measurement intensities in which the number of simulated instantaneous flow measurements ranged from 1 to 100 at a site. The discrete measurement estimator DiscQ50 is based on a simple linear regression developed between information on daily mean flows at five or more streamgages near the miscellaneous site and their corresponding GageQ50 index flows. The regression relation then was used to compute a DiscQ50 estimate at the miscellaneous site by use of the simulated instantaneous flow measurement. This process was repeated to develop a set of DiscQ50 estimates for all simulated instantaneous measurements, a weighted DiscQ50 estimate was formed from this set. Results indicated that the expected value of this weighted estimate was more precise than the IndxQ50 estimate for all measurement intensities evaluated. The integrated index-flow estimator, IntgQ50, was formed by computing a weighted average of the index estimate IndxQ50 and the DiscQ50 estimate. Results indicated that the IntgQ50 estimator was more precise than the DiscQ50 estimator at low measurement intensities of one to two measurements. At greater measurement intensities, the precision of the IntgQ50 estimator converges to the DiscQ50 estimator. Neither the DiscQ50 nor the IntgQ50 estimators provided site-specific estimates. In particular, although expected values of DiscQ50 and IntgQ50 estimates converge with increasing measurement intensity, they do not necessarily converge to the site-specific value of Q50. The site estimator of flow, SiteQ50, was developed to facilitate this convergence at higher measurement intensities. This is accomplished by use of the median of simulated instantaneous flow values for each measurement intensity level. A weighted estimate of the median and information associated with the IntgQ50 estimate was used to form the SiteQ50 estimate. Initial simulations indicate that the SiteQ50 estimator generally has greater precision than the IntgQ50 estimator at measurement intensities greater than 3, however, additional analysis is needed to identify streamflow conditions under which instantaneous measurements will produce estimates that generally converge to the index flows. A preliminary augmented index regression equation was developed, which contains the index regression estimate and two additional variables associated with base-flow recession characteristics. When these recession variables were estimated as the medians of recession parameters compute

Scientific Investigations Report↗

USGS2014

Conceptual model of the uppermost principal aquifer systems in the Williston and Powder River structural basins, United States and Canada

The three uppermost principal aquifer systems of the Northern Great Plains—the glacial, lower Tertiary, and Upper Cretaceous aquifer systems—are described in this report and provide water for irrigation, mining, public and domestic supply, livestock, and industrial uses. These aquifer systems primarily are present in two nationally important fossil-fuelproducing areas: the Williston and Powder River structural basins in the United States and Canada. The glacial aquifer system is contained within glacial deposits that overlie the lower Tertiary and Upper Cretaceous aquifer systems in the northeastern part of the Williston structural basin. Productive sand and gravel aquifers exist within this aquifer system. The Upper Cretaceous aquifer system is contained within bedrock lithostratigraphic units as deep as 2,850 and 8,500 feet below land surface in the Williston and Powder River structural basins, respectively. Petroleum extraction from much deeper formations, such as the Bakken Formation, is rapidly increasing because of recently improved hydraulic fracturing methods that require large volumes of relatively freshwater from shallow aquifers or surface water. Extraction of coalbed natural gas from within the lower Tertiary aquifer system requires removal of large volumes of groundwater to allow degasification. Recognizing the importance of understanding water resources in these energy-rich basins, the U.S. Geological Survey (USGS) Groundwater Resources Program ( http://water.usgs.gov/ogw/gwrp/ ) began a groundwater study of the Williston and Powder River structural basins in 2011 to quantify this groundwater resource, the results of which are described in this report. The overall objective of this study was to characterize, quantify, and provide an improved conceptual understanding of the three uppermost and principal aquifer systems in energy-resource areas of the Northern Great Plains to assist in groundwater-resource management for multiple uses. The study area includes parts of Montana, North Dakota, South Dakota, and Wyoming in the United States and Manitoba and Saskatchewan in Canada. The glacial aquifer system is contained within glacial drift consisting primarily of till, with smaller amounts of glacial outwash sand and gravel deposits. The lower Tertiary and Upper Cretaceous aquifer systems are contained within several formations of the Tertiary and Cretaceous geologic systems, which are hydraulically separated from underlying aquifers by a basal confining unit. The lower Tertiary and Upper Cretaceous aquifer systems each were divided into three hydrogeologic units that correspond to one or more lithostratigraphic units. The period prior to 1960 is defined as the predevelopment period when little groundwater was extracted. From 1960 through 1990, numerous flowing wells were installed near the Yellowstone, Little Missouri and Knife Rivers, resulting in local groundwater declines. Recently developed technologies for the extraction of petroleum resources, which largely have been applied in the study area since about 2005, require millions of gallons of water for construction of each well, with additional water needed for long-term operation; therefore, the potential for an increase in groundwater extraction is high. In this study, groundwater recharge and discharge components were estimated for the period 1981–2005. Groundwater recharge primarily occurs from infiltration of rainfall and snowmelt (precipitation recharge) and infiltration of streams into the ground (stream infiltration). Total estimated recharge to the Williston and Powder River control volumes is 4,560 and 1,500 cubic feet per second, respectively. Estimated precipitation recharge is 26 and 15 percent of total recharge for the Williston and Powder River control volumes, respectively. Estimated stream infiltration is 71 and 80 percent of total recharge for the Williston and Powder River control volumes, respectively. Groundwater discharge primarily is to streams and springs and is estimated to be about 97 and 92 percent of total discharge for the Williston and Powder River control volumes, respectively. Most of the remaining discharge results from pumped and flowing wells. Groundwater flow in the Williston structural basin generally is from the west and southwest toward the east, where discharge to streams occurs. Locally, in the uppermost hydrogeologic units, groundwater generally is unconfined and flows from topographically high to low areas, where discharge to streams occurs. Groundwater flow in the Powder River structural basin generally is toward the north, with local variations, particularly in the upper Fort Union aquifer, where flow is toward streams.

Manitoba;Montana;North Dakota;Saskatchewan;South D↗

USGS2014

Seismic instrumentation plan for the Hawaiian Volcano Observatory

The seismic network operated by the U.S. Geological Survey’s Hawaiian Volcano Observatory (HVO) is the main source of authoritative data for reporting earthquakes in the State of Hawaii, including those that occur on the State’s six active volcanoes (Kīlauea, Mauna Loa, Hualālai, Mauna Kea, Haleakalā, Lō‘ihi). Of these volcanoes, Kīlauea and Mauna Loa are considered “very high threat” in a report on the rationale for a National Volcanic Early Warning System (NVEWS) (Ewert and others, 2005). This seismic instrumentation plan assesses the current state of HVO’s seismic network with respect to the State’s active volcanoes and calculates the number of stations that are needed to upgrade the current network to provide a seismic early warning capability for forecasting volcanic activity. Further, the report provides proposed priorities for upgrading the seismic network and a cost assessment for both the installation costs and maintenance costs of the improved network that are required to fully realize the potential of the early warning system. HVO has operated seismometers on the Island of Hawai‘i since 1912. Currently, the seismic network includes more than 70 stations from four different organizations. Generally, the Island of Hawai‘i has most of the seismic stations in the network (and most of the activity), with the density of seismic stations increasing from the northern part of the island to the south-southeast. The strength of the current network, based on theoretical detection and location capabilities, is at the summit of Kīlauea Volcano and Kīlauea’s upper East Rift Zone and Pu‘u ‘Ō‘ō—where few, if any, upgrades need to be made to the seismic network. The network in the region between Kīlauea and Mauna Loa is slightly weaker, as is the summit of Mauna Loa. In general, the rift zones of each volcano are more poorly monitored seismically than the summits and thus require a greater number of stations to achieve a volcanic early warning capability for monitoring seismicity. Priorities for new seismic installations on the volcanoes depend on several factors, including current activity, historical activity, population exposure, and current network quality. On Kīlauea, new installations on the middle East Rift Zone, lower East Rift Zone, and lower Southwest Rift Zone appear to be the highest priorities. On Mauna Loa, improvements to the summit seismic network should be prioritized based on the analysis of the data, followed by the installation of a sparse network on both rift zones. Once installed, the next priority would be to create denser seismic networks on the rift zones, particularly where eruptions could quickly threaten populated areas (middle Northeast Rift Zone, lower Southwest Rift Zone). On Hualālai, analysis of the data indicates that the Northwest Rift Zone is the most important priority, particularly where it runs through the population center of Kalaoa. Hualālai’s South Rift Zone appears to be the lowest priority for additional seismic instrumentation of any rift zone on Kīlauea, Mauna Loa, or Hualālai because of its low historical activity and lack of population exposure. Mauna Kea and Haleakalā have less active historical eruptive activity and thus have more modest proposed upgrades to seismic instrumentation. The installation of new seismic stations is only the first part of building a volcanic early warning capability for seismicity in the State of Hawaii. Additional personnel will likely be required to study the volcanic processes at work under each volcano, analyze the current seismic activity at a level sufficient for early warning, build new tools for monitoring, maintain seismic computing resources, and maintain the new seismic stations.

Hawaii↗

USGS2011

Use of diverse geochemical data sets to determine sources and sinks of nitrate and methane in groundwater, Garfield County, Colorado, 2009

Previous water-quality assessments reported elevated concentrations of nitrate and methane in water from domestic wells screened in shallow zones of the Wasatch Formation, Garfield County, Colorado. In 2009, the U.S. Geological Survey, in cooperation with the Colorado Department of Public Health and Environment, analyzed samples collected from 26 domestic wells for a diverse set of geochemical tracers for the purpose of determining sources and sinks of nitrate and methane in groundwater from the Wasatch Formation. Nitrate concentrations ranged from less than 0.04 to 6.74 milligrams per liter as nitrogen (mg/L as N) and were significantly lower in water samples with dissolved-oxygen concentrations less than 0.5 mg/L than in samples with dissolved-oxygen concentrations greater than or equal to 0.5 mg/L. Chloride/bromide mass ratios and tracers of groundwater age (tritium, chlorofluorocarbons, and sulfur hexafluoride) indicate that septic-system effluent or animal waste was a source of nitrate in some young groundwater (less than 50 years), although other sources such as fertilizer also may have contributed nitrate to the groundwater. Nitrate and nitrogen gas (N2) concentrations indicate that denitrification was the primary sink for nitrate in anoxic groundwater, removing 99 percent of the original nitrate content in some samples that had nitrate concentrations greater than 10 mg/L as N at the time of recharge. Methane concentrations ranged from less than 0.0005 to 32.5 mg/L and were significantly higher in water samples with dissolved-oxygen concentrations less than 0.5 mg/L than in samples with dissolved-oxygen concentrations greater than or equal to 0.5 mg/L. High methane concentrations (greater than 1 mg/L) in some samples were biogenic in origin and appeared to be derived from a relatively deep source on the basis of helium concentrations and isotopic data. One such sample had water-isotopic and major-ion compositions similar to that of produced water from the underlying Mesaverde Group, which was the primary natural-gas producing interval in the study area. Methane in the Mesaverde Group was largely thermogenic in origin so biogenic methane in the sample probably was derived from deeper zones in the Wasatch Formation. The primary methane sink in the aquifer appeared to be methane oxidation on the basis of dissolved-oxygen and methane concentrations and methane isotopic data. The diverse data sets used in this study enhance previous water-quality assessments by providing new and more complete insights into the sources and sinks of nitrate and methane in groundwater. Field measurements of dissolved oxygen in groundwater were useful indicators of the Wasatch Formation's vulnerability to nitrate and methane contamination or enrichment. Results from this study also provide new evidence for the movement of water, ions, and gases into the shallow Wasatch Formation from sources such as the Mesaverde Group and deeper Wasatch Formation.

Scientific Investigations Report↗

USGS2012

Evaluation of volatile organic compound (VOC) blank data and application of study reporting levels to groundwater data collected for the California GAMA Priority Basin Project, May 2004 through September 2010

Volatile organic compounds (VOCs) were analyzed in quality-control samples collected for the California Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project. From May 2004 through September 2010, a total of 2,026 groundwater samples, 211 field blanks, and 109 source-solution blanks were collected and analyzed for concentrations of 85 VOCs. Results from analyses of these field and source-solution blanks and of 2,411 laboratory instrument blanks during the same time period were used to assess the quality of data for the 2,026 groundwater samples. Eighteen VOCs were detected in field blanks or source-solution blanks: acetone, benzene, bromodichloromethane, 2-butanone, carbon disulfide, chloroform, 1,1-dichloroethene, dichloromethane, ethylbenzene, tetrachloroethene, styrene, tetrahydrofuran, toluene, trichloroethene, trichlorofluoromethane, 1,2,4-trimethylbenzene, m - and p -xylenes, and o -xylene. The objective of the evaluation of the VOC-blank data was to determine if study reporting levels (SRLs) were needed for any of the VOCs detected in blanks to ensure the quality of the data from groundwater samples. An SRL is equivalent to a raised reporting level that is used in place of the reporting level used by the analyzing laboratory [long‑term method detection level (LT-MDL) or laboratory reporting level (LRL)] to reduce the probability of reporting false-positive detections. Evaluation of VOC-blank data was done in three stages: (1) identification of a set of representative quality‑control field blanks (QCFBs) to be used for calculation of SRLs and identification of VOCs amenable to the SRL approach, (2) evaluation of potential sources of contamination to blanks and groundwater samples by VOCs detected in field blanks, and (3) selection of appropriate SRLs from among four potential SRLs for VOCs detected in field blanks and application of those SRLs to the groundwater data. An important conclusion from this study is that to ensure the quality of the data from groundwater samples, it was necessary to apply different methods of determining SRLs from field blank data to different VOCs, rather than use the same method for all VOCs. Four potential SRL values were defined by using three approaches: two values were defined by using a binomial probability method based on one-sided, nonparametric upper confidence limits, one was defined as equal to the maximum concentration detected in the field blanks, and one was defined as equal to the maximum laboratory method detection level used during the period when samples were collected for the project. The differences in detection frequencies and concentrations among different types of blanks (laboratory instrument blanks, source-solution blanks, and field blanks collected with three different sampling equipment configurations) and groundwater samples were used to infer the sources and mechanisms of contamination for each VOC detection in field blanks. Other chemical data for the groundwater samples (oxidation-reduction state, co-occurrence of VOCs, groundwater age) and ancillary information about the well sites (land use, presence of known sources of contamination) were used to evaluate whether the patterns of detections of VOCs in groundwater samples before and after application of potential SRLs were plausible. On this basis, the appropriate SRL was selected for each VOC that was determined to require an SRL. The SRLs for ethylbenzene [0.06 microgram per liter (μg/L)], m - and p -xylenes (0.33 μg/L), o -xylene (0.12 μg/L), toluene (0.69 μg/L), and 1,2,4-trimethylbenzene (0.56 μg/L) corresponded to the highest concentrations detected in the QCFBs and were selected because they resulted in the most censoring of groundwater data. Comparisons of hydrocarbon ratios in groundwater samples and blanks and comparisons between detection frequencies of the five hydrocarbons in groundwater samples and different types of blanks suggested three dominant sources of contamination that affected groundwater samples and blanks: (1) ethylbenzene, m - and p -xylenes, o -xylene, and toluene from fuel or exhaust components sorbed onto sampling lines, (2) toluene from vials and the source blank water, and (3) 1,2,4-trimethylbenzene from materials used for collection of samples for radon-222 analysis.

California↗

USGS2015

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

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

California↗

USGS2024

Status of water-level altitudes and long-term and short-term water-level changes in the Chicot and Evangeline (undifferentiated) and Jasper aquifers, greater Houston area, Texas, 2023

Since the early 1900s, groundwater withdrawn from the primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers—has been an important source of water in the greater Houston area, Texas. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of water-level altitudes and water-level changes in these aquifers in the greater Houston area. In this report, the Chicot and Evangeline aquifers are treated as a single aquifer for the purposes of providing annual assessments of regional-scale water-level altitudes and water-level changes over time. In 2023, shaded depictions of water-level altitudes for the Chicot and Evangeline aquifers (undifferentiated) ranged from about 286 feet (ft) below the North American Vertical Datum of 1988 (NAVD 88) to about 169 ft above NAVD 88. The largest decline in water-level altitudes indicated by the 1977–2023 long-term water-level-change map was in south-central Montgomery County southeast of The Woodlands. In comparison, the 1990–2023 long-term water-level-change map depicts the largest declines in water-level altitudes in localized areas at or near certain wells in parts of northwestern Harris County and south-central Montgomery County. The largest rise in water-level altitudes for 1977–2023 is depicted in a relatively large area in southeastern Harris County, whereas the largest rise in water-level altitudes for 1990–2023 is depicted in a relatively large area in central Harris County. The 5-year short-term water-level-change map depicts the largest declines at three wells in northern Fort Bend County, one well in western Harris County, and three wells in south-central Montgomery County and the largest rise at one well in central Harris County. The 1-year short-term water-level-change map depicts the largest declines at one well in northern Fort Bend County and two wells in southwestern Harris County and the largest rises at one well in northern Brazoria County and one well in south-central Montgomery County. In 2023, shaded depictions of water-level altitudes for the Jasper aquifer ranged from about 242 ft below NAVD 88 to about 218 ft above NAVD 88. The 2000–23 long-term water-level-change map depicts water-level declines throughout the study area where water-level-measurement data from the aquifer were collected, with the largest declines in north-central Harris County and south-central Montgomery County south of The Woodlands. The 5-year short-term water-level-change map depicts the largest declines at two wells in central Montgomery County near Conroe and two wells in south-central Montgomery County southeast of The Woodlands and the largest rise at one well in western Montgomery County. The 1-year short-term water-level-change map depicts the largest declines at four wells in south-central Montgomery County southeast of The Woodlands and one well in central Montgomery County near Conroe and the largest rises at two wells in western Montgomery County.

Texas↗

USGS2010

Estimated Withdrawals and Use of Water in Colorado, 2005

The future health and economic welfare of the people and environment of Colorado depend on a continuous supply of fresh water. Detailed, comprehensive information on the use of water from Colorado's diverse surface-water and groundwater resources is important to water managers and planners by providing information they need to quantify current stresses and estimate and plan for future water needs. As part of the U.S. Geological Survey's (USGS) National Water Use Information Program (NWUIP), Statewide water withdrawal and water-use data have been collected or estimated and summarized in this report by county and by four-digit hydrologic unit code for the following seven water-use categories: irrigation (crop and golf course), public supply, self-supplied domestic, self-supplied industrial, livestock, mining, and thermoelectric power generation. A summary for instream water use for hydroelectric power generation also is included. This report is published in cooperation with the Colorado Water Conservation Board. In 2005, an estimated 13,581.22 million gallons per day (Mgal/d) was withdrawn from groundwater and surface-water sources in Colorado for the seven water-use categories. Withdrawals from surface water represented about 11,035 Mgal/d, or 81.3 percent of the total, whereas withdrawals from groundwater sources represented an estimated 2,546 Mgal/d or 18.7 percent of the total. Irrigation (combined crop and golf course) totaled 12,362.49 Mgal/d or 91 percent of the total water withdrawals in the State of Colorado. Crop irrigation accounted for 99.7 percent (12,321.85 Mgal/d) of the irrigation, whereas the 243 turf golf courses in Colorado accounted for 0.3 percent (40.64 Mgal/d) of the total irrigation water withdrawals. Total withdrawals for the other water-use categories were public supply, 864.17 Mgal/d; self-supplied domestic, 34.43 Mgal/d; self-supplied industrial, 142.44 Mgal/d; livestock, 33.06 Mgal/d; mining, 21.42 Mgal/d (includes both fresh and saline water); and thermoelectric, 123.21 Mgal/d. The counties with the largest total withdrawals (greater than 500 Mgal/d) were Mesa, Weld, Rio Grande, Montrose, Gunnison, and Saguache. Counties with the smallest total withdrawals (less than 5 Mgal/d) were Clear Creek, Gilpin, and San Juan. Four-digit hydrologic unit codes with the greatest withdrawals were 1019 (South Platte River Basin), 1301 (Rio Grande Basin), and 1102 (Arkansas River Basin); the high withdrawal rates were driven by crop irrigation withdrawals. Total instream water use for hydroelectric power generation was 5,253.60 Mgal/d. Groundwater withdrawals were estimated for 2004 for the bedrock and overlying alluvial aquifers in the Denver Basin for irrigation, public supply, commercial/industrial, household use only, and domestic/livestock water-use categories. Withdrawals were estimated for input into the USGS Denver Basin model by using the equations in the Senate Bill 96-074 groundwater model. The greatest withdrawals were for public supply. The smallest withdrawals were for household-use-only wells. Douglas County had the greatest groundwater withdrawals (183.98 Mgal/d), whereas Broomfield County had the smallest (3.09 Mgal/d). Of the seven Denver Basin aquifers, the Lower Arapahoe aquifer had the greatest total estimated withdrawals (287.11 Mgal/d), with Douglas County having the greatest public-supply withdrawal of any county (95.29 Mgal/d) from this aquifer. The Upper Dawson aquifer was the least used of the Denver Basin aquifers, based on estimated withdrawals of 17.64 Mgal/d. As part of the Colorado Statewide Water Supply Initiative (SWSI), forecasts of future water demand were made based on information such as population, climate, and then-current (2000) water-use information and did not include the effects of future water conservation. Categories compared between estimates in the SWSI baseline forecasted water demand and the USGS water-use compilation were limited to county population and w

Scientific Investigations Report↗

USGS2011

Crims Island-Restoration and monitoring of juvenile salmon rearing habitat in the Columbia River Estuary, Oregon, 2004-10

Under the 2004 Biological Opinion for operation of the Federal Columbia River Power System released by the National Marine Fisheries Service, the U.S. Army Corps of Engineers (USACE), the Bonneville Power Administration (BPA), and the Bureau of Reclamation (Reclamation) were directed to restore more than 4,047 hectares (10,000 acres) of tidal marsh in the Columbia River estuary by 2010. Restoration of Crims Island near Longview, Washington, restored 38.1 hectares of marsh and swamp in the tidal freshwater portion of the lower Columbia River. The goal of the restoration was to improve habitat for juveniles of Endangered Species Act (ESA)-listed salmon stocks and ESA-listed Columbian white-tailed deer. The U.S. Geological Survey (USGS) monitored and evaluated the fisheries and aquatic resources at Crims Island in 2004 prior to restoration (pre-restoration), which began in August 2004, and then post-restoration from 2006 to 2009. This report summarizes pre- and post-restoration monitoring data used by the USGS to evaluate project success. We evaluated project success by examining the interaction between juvenile salmon and a suite of broader ecological measures including sediments, plants, and invertebrates and their response to large-scale habitat alteration. The restoration action at Crims Island from August 2004 to September 2005 was to excavate a 0.6-meter layer of soil and dig channels in the interior of the island to remove reed canary grass and increase habitat area and tidal exchange. The excavation created 34.4 hectares of tidal emergent marsh where none previously existed and 3.7 hectares of intertidal and subtidal channels. Cattle that had grazed the island for more than 50 years were relocated. Soil excavated from the site was deposited in upland areas next to the tidal marsh to establish an upland forest. Excavation deepened and widened an existing T-shaped channel to increase tidal flow to the interior of the island. The western arm of the existing 'T-channel' was extended westward and connected to Bradbury Slough to create a second outlet to the main river. New intertidal channels were constructed from the existing 'T-channel' and tidal mudflats became inundated at high tide to increase rearing habitat for juvenile salmonids. The restoration action resulted in a 95-percent increase in available juvenile salmon rearing habitat. We collected juvenile salmon and other fishes at Crims Island and a nearby reference site using beach seines and fyke nets annually from March through August during all years. Benthic invertebrates were collected with sediment corers and drift invertebrates were collected with neuston nets. Juvenile salmon stomach contents were sampled using lavage. Vegetation and sediments characteristics were surveyed and we conducted a topographic/bathymetric survey using a RTK (real time kinematic) GPS (global positioning system). The fish assemblage at Crims Island, composed primarily of threespine stickleback (Gasterosteus aculeatus), non-native banded killifish (Fundulus diaphanus), peamouth chub (Mylocheilus caurinus), subyearling Chinook salmon (Oncorhynchus tshawytscha) (hereinafter referred to as subyearlings), and small numbers of juvenile chum salmon (Oncorhynchus keta), did not differ appreciably pre- and post-restoration. Subyearlings were the primary salmonid collected and were seasonally abundant from April through May during all years. The abundance of juvenile salmon declined seasonally as water temperature exceeded 20 degrees C in the Reference site by mid-June; however, subyearlings persisted at the Mainstem site and in subtidal channels of the Restoration site through the summer in water temperatures exceeding 22 degrees C. Residence times of subyearlings in Crims Island backwaters generally were short consisting of one or two tidal cycles. Median residence time was longer in the Restoration site than in the Reference site pre- and post-restoration. Small (mean = 55.7 millimeters) subyea

Scientific Investigations Report↗

USGS2014

Numerical simulation of groundwater flow in the Columbia Plateau Regional Aquifer System, Idaho, Oregon, and Washington

A three-dimensional numerical model of groundwater flow was constructed for the Columbia Plateau Regional Aquifer System (CPRAS), Idaho, Oregon, and Washington, to evaluate and test the conceptual model of the system and to evaluate groundwater availability. The model described in this report can be used as a tool by water-resource managers and other stakeholders to quantitatively evaluate proposed alternative management strategies and assess the long‑term availability of groundwater. The numerical simulation of groundwater flow in the CPRAS was completed with support from the Groundwater Resources Program of the U.S. Geological Survey Office of Groundwater. The model was constructed using the U.S. Geological Survey modular three-dimensional finite-difference groundwater-flow model, MODFLOW-NWT. The model uses 3-kilometer (9,842.5 feet) grid cells that subdivide the model domain by 126 rows and 131 columns. Vertically, the model domain was subdivided into six geologic model units. From youngest to oldest, the units are the Overburden, the Saddle Mountains Basalt, the Mabton Interbed, the Wanapum Basalt, the Vantage Interbed, and the Grande Ronde Basalt. Natural recharge was estimated using gridded historical estimates of annual precipitation for the period 1895–2007. Pre-development recharge was estimated to be the average natural recharge for this period. Irrigation recharge and irrigation pumping were estimated using a remote-sensing based soil-water balance model for the period 1985–2007. Pre-1985 irrigation recharge and pumping were estimated using previously published compilation maps and the history of large-scale irrigation projects. Pumping estimates for municipal, industrial, rural, residential, and all other uses were estimated using reported values and census data. Pumping was assumed to be negligible prior to 1920. Two models were constructed to simulate groundwater flow in the CPRAS: a steady-state predevelopment model representing conditions before large-scale pumping and irrigation altered the system, and a transient model representing the period 1900–2007. Automated parameter-estimation techniques (steady-state predevelopment model) and traditional trial-and-error (transient model) methods were used for calibration. To calibrate the steady-state and transient models, 10,525 and 46,460 water level measurements, respectively, and 50 base-flow estimates were used. The steady-state model simulated the shape, slope, and trends of a potentiometric surface that was generally consistent with mapped water levels. For the transient model, the mean and median difference between simulated and measured hydraulic heads is -10 and 4 ft, respectively, with a standard deviation of 164 ft over a 5,648 ft range of measured heads. The residuals for the simulation period show that 52 percent of the simulated heads exceeded measured heads with a median residual value of 43 ft, and 48 percent were less than measured heads with a median residual value of -76 ft. The CPRAS model was constructed to derive components of the groundwater budget and help understand the interactions of stresses, such as recharge, groundwater pumping, and commingling wells on the groundwater and surface-water system. Through these applications, the model can be used to identify trends in groundwater storage and use, and quantify groundwater availability. The annual groundwater budgets showed several patterns of change over the simulation period. Groundwater pumping was negligible until the 1950s and began to increase significantly during the 1970s and 1980s. Recharge was highly variable due to the interannual variability of precipitation, but began to increase in the late 1940s due to the increase in surface-water irrigation projects. Groundwater contributions to streamflow (base flow) followed recharge closely. However, in areas of significant groundwater-level decline, base flow is reduced. Groundwater pumping had the greatest effect on water levels, followed by irrigation enhanced recharge. Commingling was a larger factor in structurally complex upland areas where hydraulic-head gradients are naturally high. Groundwater pumping has increased substantially over the past 40–50 years; this increase resulted in declining water levels at depth and decreased base flows over much of the study area. The effects of pumping are mitigated somewhat by the increase of surface-water irrigation, especially in the shallow Overburden unit, and commingling wells in some areas. During dry to average years, groundwater pumping causes a net loss of groundwater in storage and current condition (2000–2007) groundwater pumping exceeds recharge in all but the wettest of years.

Idaho, Oregon, Washington↗

USGS2008

Simulation of Flow, Sediment Transport, and Sediment Mobility of the Lower Coeur d'Alene River, Idaho

A one-dimensional sediment-transport model and a multi-dimensional hydraulic and bed shear stress model were developed to investigate the hydraulic, sediment transport, and sediment mobility characteristics of the lower Coeur d?Alene River in northern Idaho. This report documents the development and calibration of those models, as well as the results of model simulations. The one-dimensional sediment-transport model (HEC-6) was developed, calibrated, and used to simulate flow hydraulics and erosion, deposition, and transport of sediment in the lower Coeur d?Alene River. The HEC-6 modeled reach, comprised of 234 cross sections, extends from Enaville, Idaho, on the North Fork of the Coeur d?Alene River and near Pinehurst, Idaho, on the South Fork of the river to near Harrison, Idaho, on the main stem of the river. Bed-sediment samples collected by previous investigators and samples collected for this study in 2005 were used in the model. Sediment discharge curves from a previous study were updated using suspended-sediment samples collected at three sites since April 2000. The HEC-6 was calibrated using river discharge and water-surface elevations measured at five U.S. Geological Survey gaging stations. The calibrated HEC-6 model allowed simulation of management alternatives to assess erosion and deposition from proposed dredging of contaminated streambed sediments in the Dudley reach. Four management alternatives were simulated with HEC-6. Before the start of simulation for these alternatives, seven cross sections in the reach near Dudley, Idaho, were deepened 20 feet?removing about 296,000 cubic yards of sediments?to simulate dredging. Management alternative 1 simulated stage-discharge conditions from 2000, and alternative 2 simulated conditions from 1997. Results from alternatives 1 and 2 indicated that about 6,500 and 12,300 cubic yards, respectively, were deposited in the dredged reach. These figures represent 2 and 4 percent, respectively, of the total volume of dredged sediments removed before the start of simulation. In alternatives 3 and 4, the incoming total sediment discharges from the South Fork of the river were decreased by one-half. Management alternative 3 simulated stage-discharge conditions from 2000, and alternative 4 simulated conditions from 1997. Reducing incoming sediment discharge from the South Fork did not affect the streambed and deposition in the Dudley and downstream reaches, probably because the distance between the South Fork and the Dudley reach is long enough for sediment supply, transport capacity, and channel geometry to be balanced before reaching the Dudley and downstream reaches. Development and calibration of a multi-dimensional hydraulic and bed shear stress model (FASTMECH) allowed simulation of water-surface elevation, depth, velocity, bed shear stress, and sediment mobility in the Dudley reach (5.3 miles). The computational grid incorporated bathymetric and Light Detection and Ranging (LIDAR) data, with a node spacing of about 2.5 meters. With the exception of the fourth FASTMECH calibration simulation, results from the FASTMECH calibration simulations indicated that flow depths, flow velocities, and bed shear stresses increased as river discharge increased. Water-surface elevations in the fourth calibration simulation were about 2 feet higher than those in the other simulations because high lake levels in Coeur d?Alene Lake caused backwater conditions. Average simulated velocities along the thalweg ranged from about 3 to 5.3 feet per second, and maximum simulated velocities ranged from 3.9 to 7 feet per second. In the dredged reach, average simulated velocity along the thalweg ranged from 3.5 to 6 feet per second. The model also simulated several back-eddies (flow reversal); the largest eddy encompassed about one-third of the river width. Average bed shear stresses increased more than 200 percent from the first to the last simulation. Simulated sediment mobility, asses

Scientific Investigations Report↗

USGS2010

Potentiometric Surfaces and Water-Level Trends in the Cockfield (Upper Claiborne) and Wilcox (Lower Wilcox) Aquifers of Southern and Northeastern Arkansas, 2009

Eocene-age sand beds near the base of the Cockfield Formation of Claiborne Group constitute the aquifer known locally as the Cockfield aquifer. Upper-Paleocene age sand beds within the lower parts of the Wilcox Group constitute the aquifer known locally as the Wilcox aquifer. In 2005, reported water withdrawals from the Cockfield aquifer in Arkansas totaled 16.1 million gallons per day, while reported water withdrawals from the Wilcox aquifer in Arkansas totaled 27.0 million gallons per day. Major withdrawals from these units were for industrial and public water supplies with lesser but locally important withdrawals for commercial, domestic, and agricultural uses. During February 2009, 56 water-level measurements were made in wells completed in the Cockfield aquifer and 57 water-level measurements were made in wells completed in the Wilcox aquifer. The results from the 2009 water-level measurements are presented in potentiometric-surface maps and in combination with previous water-level measurements. Trends in water-level change over time within the two aquifers are investigated using water-level difference maps and well hydrographs. Water-level difference maps were constructed for each aquifer using the difference between depth to water measurements made in 2003 to 2009. Well hydrographs for each aquifer were constructed for wells with 20 or more years of historical water-level data. The hydrographs were evaluated individually using linear regression to calculate the annual rise or decline in water levels, and by aggregating the regression results by county and statistically summarizing for the range, mean, and median water-level change in each county. The 2009 potentiometric surface of the Cockfield aquifer map indicates the regional direction of groundwater flow generally towards the east and southeast, except in two areas of intense groundwater withdrawals that have developed into cones of depression. The lowest water-level altitude measured was 43 feet and the highest water-level altitude measured was 351 feet. A water-level difference map was constructed from 54 wells completed in the Cockfield aquifer within Arkansas. The largest rise in water level was 14.9 feet and the largest decline was 27.4 feet. Seven wells had a rise in water level, and the remaining 47 wells had a decline in water level. Hydrographs for 33 wells completed in the Cockfield aquifer were developed. Hydrographs indicate water-level changes in individual wells ranged from rises of 0.33 feet per year to declines of 1.21 feet per year over the 20-year period (1990-2009). County summaries of the linear regression analysis indicate Cleveland and Columbia Counties have mean annual rises. Arkansas, Ashley, Bradley, Calhoun, Chicot, Desha, Drew, Lincoln, and Union Counties have mean annual declines. The potentiometric surface for the Wilcox aquifer is presented using two maps, one for a southern area and another for a northeastern area, because of the absence of water-level data in the central part of the State. The direction of groundwater flow in the southern area is generally the east, except around two cones of depression and around two mounds of elevated water levels. Water-level altitudes in the southern area range from 147 feet to 400 feet. The direction of groundwater flow in the northeastern area is generally to the south and southeast except in an area of intense groundwater withdrawals that has altered the flow to a westerly direction. Two water-level difference maps were constructed using water-level altitudes measured in 2003 to 2009 from 53 wells completed in the Wilcox aquifer within southern and northeastern Arkansas. In the southern area the largest rise in water level was 16.0 feet and the largest decline was 17.7 feet. Eight wells in the southern area had rising water levels and the remaining five wells had declining water levels. In the northeastern area, the largest rise in water level was 1.3 feet and the larg

Scientific Investigations Report↗

USGS2022

Hydrogeology and groundwater quality in the San Agustin Basin, New Mexico, 1975–2019

This report describes the findings of a U.S. Geological Survey study, completed in cooperation with the Bureau of Land Management, focused on better understanding the present-day (1975–2019) hydrogeology and groundwater quality of the San Agustin Basin in west-central New Mexico to support sustainable groundwater resource management. The basin hosts a relatively undeveloped basin-fill and alluvium aquifer system and is topographically divided into east and west subbasins by the McClure Hills. Groundwater chemistry and groundwater elevation data were compiled, collected, and interpreted in the context of groundwater flow and quality. The analyses presented in this report consider groundwater chemistry data collected within the last decade (2010–19) and groundwater elevation data collected from 1975 through 2019 to provide insight into present-day conditions. Groundwater elevations show that groundwater typically moves from the highlands to the lowlands, with a prominent east to west regional trend. Groundwater elevations were lowest in the southwestern portion of the west subbasin, where estimated flow directions suggest underflow through the local highlands into the northern East Fork Gila River watershed, which is further supported by historical groundwater elevation data from the northern East Fork Gila River watershed. Gradual groundwater elevation gradients (about 2 feet per mile) near the east and west subbasin divide suggest that groundwater slowly flows from the east subbasin to the west subbasin. Quantitative analyses of groundwater chemistry data show that groundwater in both subbasins has similar chemical characteristics. A systematic east to west groundwater evolution in water chemistry was not observed despite evidenced subbasin connectivity. The absence of this pattern suggests that groundwater mixing is regionally prevalent, sediment reactivity is low and variable, and (or) recharge conditions are comparable in both subbasins. Groundwater chemistry was generally independent of aquifer type, suggesting that the aquifers are hydrologically well connected. Corrected carbon-14 groundwater age estimates in the basin ranged from 232 to 13,916 years before present with a median of 5,409 years. A wide range of groundwater ages is therefore present in the basin, with waters commonly being thousands of years old, thereby supporting generally slow regional groundwater movement. A component of relatively young groundwater, for which estimated ages could not be accurately computed, is also present in the basin, and it may commonly mix with older waters. The spatial distribution of categorical and quantitative groundwater ages indicates that most recharge likely occurs in the highlands through mountain-block recharge and as focused recharge within arroyos, although evidence of modern (1953 and after) groundwater was minimal at sampled sites. Median annual gradients (groundwater elevation change over time) indicate that most groundwater elevations in the lowlands changed little (−0.2 to 0.2 foot per year) from 1975 through 2019. Groundwater elevations in the highlands varied more annually, which is likely due to recharge from precipitation events. These more variable groundwater elevations in the highlands compared with the lowlands, along with groundwater ages, provide further evidence that most groundwater recharge takes place in the highlands, with minimal recharge in the lowlands. Median groundwater elevation change for all sites was −0.05 foot per year. Temporal consistency of lowland groundwater elevations suggests that regional groundwater dynamics have been more or less stable through time under current climate and development conditions, although median annual gradients indicate that groundwater elevations may have slightly declined on average between 1975 and 2019.

New Mexico↗

USGS2007

Estimated water use and availability in the Pawtuxet and Quinebaug River basins, Rhode Island, 1995-99

Water availability became a concern in Rhode Island during a drought in 1999, and an investigation was needed to assess demands on the hydrologic system from withdrawals during periods of little to no precipitation. The low water levels during the drought prompted the U.S. Geological Survey and the Rhode Island Water Resources Board to begin a series of studies on water use and availability in each drainage area in Rhode Island for 1995–99. The study area for this report, which includes the Pawtuxet River Basin in central Rhode Island (231.6 square miles) and the Quinebaug River Basin in western Rhode Island (60.97 square miles), was delineated as the surface-water drainage areas of these basins. During the study period from 1995 through 1999, two major water suppliers withdrew an average of 71.86 million gallons per day (Mgal/d) from the Pawtuxet River Basin; of this amount, about 35.98 Mgal/d of potable water were exported to other basins in Rhode Island. The estimated water withdrawals from minor water suppliers were 0.026 Mgal/d in the Pawtuxet River Basin and 0.003 Mgal/d in the Quinebaug River Basin. Total self-supply withdrawals were 2.173 Mgal/d in the Pawtuxet River Basin and 0.360 Mgal/d in the Quinebaug River Basin, which has no public water supply. Total water use averaged 18.07 Mgal/d in the Pawtuxet River Basin and 0.363 Mgal/d in the Quinebaug River Basin. Total return flow in the Pawtuxet River Basin was 30.64 Mgal/d, which included about 12.28 Mgal/d that were imported from other basins in Rhode Island. Total return flow was 0.283 Mgal/d in the Quinebaug River Basin. During times of little to no recharge in the form of precipitation, the surface- and ground-water flows are from storage primarily in the stratified sand and gravel deposits; water also flows through the till deposits, but at a slower rate. The ground water discharging to the streams during times of little to no recharge from precipitation is referred to as base flow. The PART program, a computerized hydrograph-separation application, was used to analyze the data collected at two selected index stream-gaging stations to determine water availability on the basis of the 75th, 50th, and 25th percentiles of the total base flow; the base flow for the 7-day, 10-year low-flow scenario; and the base flow for the Aquatic Base Flow scenario for both stations. The index stream-gaging stations used in the analysis were the Branch River at Forestdale, Rhode Island (period of record 1957–1999) and the Nooseneck River at Nooseneck, Rhode Island (period of record 1964–1980). A regression equation was used to estimate unknown base-flow contributions from sand and gravel deposits at the two stations. The base-flow contributions from sand and gravel deposits and till deposits at the index stations were computed for June, July, August, and September within the periods of record, and divided by the area of each type of surficial deposit at each index station. These months were selected because they define a period when there is usually an increased demand for water and little to no precipitation. The base flows at the stream-gaging station Branch River at Forestdale, Rhode Island were lowest in August at the 75th, 50th, and 25th percentiles (29.67, 21.48, and 13.30 Mgal/d, respectively). The base flows at the stream-gaging station Nooseneck River at Nooseneck, Rhode Island were lowest in September at the 75th percentile (3.551 Mgal/d) and lowest in August at the 50th and 25th percentiles (2.554 and 1.811 Mgal/d). The base flows per unit area for the index stations were multiplied by the areas of sand and gravel and till in the studyarea subbasins to determine the amount of available water for each scenario. The water availability in the Pawtuxet River Basin at the 50th percentile ranged from 126.5 Mgal/d in August to 204.7 Mgal/d in June, and the total gross water availability for the 7-day, 10-year low-flow scenario at the 50th percentile ranged from 112.2 Mgal/d in August to 190.4 Mgal/d in June. The Scituate Reservoir safe yield was 83 Mgal/d in all scenarios. Water availability in the Quinebaug River Basin ranged from 13.94 Mgal/d in August to 30.53 Mgal/d in June at the 50th percentile. The total gross water availability for the 7-day, 10-year low-flow scenario at the 50th percentile ranged from 14.26 Mgal/d in August to 42.69 Mgal/d in June. Because water withdrawals and use are greater during the summer than other times of the year, water availability in June, July, August, and September was compared to water withdrawals in the basin and subbasins. The ratios of water withdrawn to water available were calculated for the 75th, 50th, and 25th percentiles for the subbasins; the closer the ratio is to 1, the closer the withdrawals are to the estimated water available, and the less net water is available. Withdrawals in July were higher than in the other summer months in both basins. In the Pawtuxet River Basin, the ratios were close to 1 in July for the estimated gross yield (from sand and gravel and from till and from the Scituate Reservoir safe yield), 7-day, 10-year low-flow scenario, and Aquatic Base Flow scenario at the 75th percentile and in August for all three scenarios at the 50th and 25th percentiles. In the Quinebaug River Basin, the ratios were close to 1 in August for the estimated gross yield; 7-day, 10-year low-flow scenario; and Aquatic Base Flow scenario. A long-term water budget was calculated for 1941 through 1999 to identify and assess the basin and subbasin inflow and outflows for the Pawtuxet and Quinebaug River Basins. The water withdrawals and return flows used in the budget were from 1995 through 1999. Inflow was assumed to be equal to outflow; total inflows and outflows were 574.9 Mgal/d in the Pawtuxet River Basin and 148.4 Mgal/d in the Quinebaug River Basin. Precipitation and return flow were 95 and 5 percent of the estimated inflows to the Pawtuxet River Basin, respectively. Precipitation was 100 percent of the estimated inflow to the Quinebaug River Basin; return flow was less than 1 percent of the inflow. Evapotranspiration, streamflow, and water withdrawals were 46, 41, and 13 percent, respectively, of the estimated outflows in the Pawtuxet River Basin. Evapotranspiration and streamflow were 49 and 51 percent, respectively, of the estimated outflows in the Quinebaug River Basin. Water withdrawals were less than 1 percent of outflows in the Quinebaug River Basin.

Rhode Island↗

USGS2024

Peak streamflow trends in South Dakota and their relation to changes in climate, water years 1921–2020

Peak-flow (flood) frequency analysis is essential to water-resources management applications, including the design of critical infrastructure such as bridges and culverts, and floodplain mapping. Federal guidelines for performing peak-flow flood frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C, which documents the guidelines for determining annual peak streamflow frequency, is that, for basins without major hydrologic alterations (for example, regulation, diversion, and urbanization), statistical properties of the distribution of annual peak streamflows are stationary; that is, the mean, variance, and skew are constant through time. Nonstationarity is a statistical property of a peak-flow series such that the long-term (on the order of decades) distributional properties change one or more times either gradually or abruptly through time. Individual nonstationarities may be attributed to one source such as flow regulation, land-use change, or climate but are often the result of a combination of sources, making detection and attribution of nonstationarities challenging. In response to a growing concern regarding nonstationarity in peak streamflows in the region, the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, assessed the potential nonstationarity in peak streamflows in the north-central United States. This chapter characterizes the effects of natural hydroclimatic shifts and potential climate change on annual peak streamflows in the State of South Dakota. Annual peak and daily streamflow as well as model-simulated gridded climatic data were examined for temporal monotonic trends, change points, and other statistical properties indicative of changing climatic and environmental conditions. Changes in annual peak and daily flows were evaluated among 13, 35, and 81 qualifying U.S. Geological Survey streamgages for the 75-, 50-, and 30-year trend periods through water year 2020 (the period from October 1, 2019, to September 30, 2020) in South Dakota, respectively. No qualifying streamgages were in the 100-year trend period in the State. Statistical tests for autocorrelation (independent and identically distributed assumption), monotonic trends, and change points in the median and scale are analyzed to evaluate potential stationarity violations (nonstationarity) for performing at-site peak-flow flood-frequency analysis. The trends are reported using a likelihood approach as an alternative to simply reporting significant trends with an arbitrary p -value cutoff point. A distinct east-west spatial pattern of likely upward and downward monotonic trends and change points, respectively, was detected in 75- and 50-year trend periods, but an inconsistent spatial pattern was detected in the 30-year trend period. Additionally, change points in the median annual peak streamflows were detected in the late 1970s and early 1980s in the western part of the State, but in the east, the change point was more commonly detected in 1992–93. A similar east-west spatial pattern of likely upward and downward trends was detected in the annual peak-flow timing, the day of the year of the annal peak streamflow. In the western part of the State, the annual peak streamflows are arriving earlier, but in the east, the annual peak streamflows are arriving later. A peaks-over-threshold (POT) analysis where, on average, there are two events per year (POT2) and four events per year (POT4) was also used to evaluate changes in the frequency (count) of daily streamflows exceeding the threshold. Similar to detected changes in the annual peak streamflow, an east-west likely upward or downward change corresponding to an increase or decrease, respectively, in the frequency of daily streamflow greater than a POT2 and POT4 threshold was detected. A monthly water-balance model was used to evaluate hydroclimatic variation in annual and seasonal precipitation, snowfall, potential evapotranspiration, and soil moisture storage for all qualifying streamgages in the 75-, 50-, and 30-year trend periods. Detected trends in the annual hydroclimatic metrics for the 75- and 50-year trend periods indicate a spatially consistent statewide increase in precipitation, decrease in snowfall, increase in potential evapotranspiration, and increase in soil moisture storage. Furthermore, detected trends in seasonal precipitation in the 75- and 50-year trend periods highlight a pronounced change in precipitation in winter and later into the summer season, especially in the 50-year trend period in the eastern part of the State. Statewide increases in seasonal soil moisture storage were also detected, highlighting year-round increasing flood magnitudes, particularly in the eastern part of the State. Based on the results of these stationarity tests for the qualifying streamgages in South Dakota among the 75-, 50-, and 30-year trend periods, consistent temporal and spatial patterns of nonstationarity were detected among the 75- and 50-year trend periods. Furthermore, when nonstationarity is detected in daily streamflow, increased streamflow and volume (increasing frequency in POT), as well as potentially bridge scour, may have implications on culvert and highway design in the eastern part of South Dakota. Thus, when performing at-site peak-flow flood-frequency analyses in South Dakota, potential nonstationarities and alternative approaches are important considerations.

South Dakota↗

USGS2012

Occurrence of anthropogenic organic compounds and nutrients in source and finished water in the Sioux Falls area, South Dakota, 2009-10

Anthropogenic organic compounds (AOCs) in drinking-water sources commonly are derived from municipal, agricultural, and industrial wastewater sources, and are a concern for water-supply managers. A cooperative study between the city of Sioux Falls, S. Dak., and the U.S. Geological Survey was initiated in 2009 to (1) characterize the occurrence of anthropogenic organic compounds in the source waters (groundwater and surface water) to water supplies in the Sioux Falls area, (2) determine if the compounds detected in the source waters also are present in the finished water, and (3) identify probable sources of nitrate in the Big Sioux River Basin and determine if sources change seasonally or under different hydrologic conditions. This report presents analytical results of water-quality samples collected from source waters and finished waters in the Sioux Falls area. The study approach included the collection of water samples from source and finished waters in the Sioux Falls area for the analyses of AOCs, nutrients, and nitrogen and oxygen isotopes in nitrate. Water-quality constituents monitored in this study were chosen to represent a variety of the contaminants known or suspected to occur within the Big Sioux River Basin, including pesticides, pharmaceuticals, sterols, household and industrial products, polycyclic aromatic hydrocarbons, antibiotics, and hormones. A total of 184 AOCs were monitored, of which 40 AOCs had relevant human-health benchmarks. During 11 sampling visits, 45 AOCs (24 percent) were detected in at least one sample of source or finished water, and 13 AOCs were detected in at least 20 percent of all samples. Concentrations of detected AOCs were all less than 1 microgram per liter, except for two AOCs in multiple samples from the Big Sioux River, and one AOC in finished-water samples. Concentrations of AOCs were less than 0.1 microgram per liter in more than 75 percent of the detections. Nutrient concentrations varied seasonally in source-water samples from surface water and groundwater. In the Big Sioux River, nitrite plus nitrate concentrations were typically less than 1 milligram per liter as nitrogen, and reached a maximum of 4.06 milligrams per liter as nitrogen following a June 2010 storm. Nitrite plus nitrate concentrations in groundwater ranged from less than 0.1 to 0.701 milligram per liter as nitrogen. Eight of the AOCs detected have a human-health benchmark that could be used to evaluate the concentrations in a human-health context. Four AOCs had maximum concentrations within an order of magnitude of the benchmark, indicating that additional monitoring of the compound may be warranted. Three herbicides (atrazine, metolachlor, and prometon) and one degradate (deethylatrazine) were detected in finished-water samples as frequently as in source-water samples. The concentrations of herbicides in source water varied by an order of magnitude from the period of peak use (early summer) to the winter months. Groundwater and finished-water concentrations of atrazine were similar for the six sampling dates when groundwater was the only source water used. Upstream wastewater discharges contributed a fairly small percentage of the flow to the Big Sioux River near Sioux Falls, but several AOCs associated with wastewater were frequently detected. The interpretation of all potential sources of nitrogen cannot be accomplished by use of nitrogen and oxygen isotopes in nitrate alone, but provides a qualitative indication that very little nitrate originates from excess fertilizer runoff, and most nitrate originates from municipal wastewater effluent, manure runoff (either from field application or feeding operations), or fertilizers mineralized by processes in the soil.

South Dakota↗

USGS2012

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

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

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

USGS2014

A precipitation-runoff model for simulating natural streamflow conditions in the Smith River watershed, Montana, water years 1996-2008

This report documents the construction of a precipitation-runoff model for simulating natural streamflow in the Smith River watershed, Montana. This Precipitation-Runoff Modeling System model, constructed in cooperation with the Meagher County Conservation District, can be used to examine the general hydrologic framework of the Smith River watershed, including quantification of precipitation, evapotranspiration, and streamflow; partitioning of streamflow between surface runoff and subsurface flow; and quantifying contributions to streamflow from several parts of the watershed. The model was constructed by using spatial datasets describing watershed topography, the streams, and the hydrologic characteristics of the basin soils and vegetation. Time-series data (daily total precipitation, and daily minimum and maximum temperature) were input to the model to simulate daily streamflow. The model was calibrated for water years 2002–2007 and evaluated for water years 1996–2001. Though water year 2008 was included in the study period to evaluate water-budget components, calibration and evaluation data were unavailable for that year. During the calibration and evaluation periods, simulated-natural flow values were compared to reconstructed-natural streamflow data. These reconstructed-natural streamflow data were calculated by adding Bureau of Reclamation’s depletions data to the observed streamflows. Reconstructed-natural streamflows represent estimates of streamflows for water years 1996–2007 assuming there was no agricultural water-resources development in the watershed. Additional calibration targets were basin mean monthly solar radiation and potential evapotranspiration. The model estimated the hydrologic processes in the Smith River watershed during the calibration and evaluation periods. Simulated-natural mean annual and mean monthly flows generally were the same or higher than the reconstructed-natural streamflow values during the calibration period, whereas they were lower during the evaluation period. The shape of the annual hydrographs for the simulated-natural daily streamflow values matched the shape of the hydrographs for the reconstructed-natural values for most of the calibration period, but daily streamflow values were underestimated during the evaluation period for water years 1996–1998. The model enabled a detailed evaluation of the components of the water budget within the Smith River watershed during the water year 1996–2008 study period. During this study period, simulated mean annual precipitation across the Smith River watershed was 16 inches, out of which 14 inches evaporated or transpired and 2 inches left the basin as streamflow. Per the precipitation-runoff model simulations, during most of the year, surface runoff rarely (less than 2 percent of the time during water years 2002–2008) makes up more than 10 percent of the total streamflow. Subsurface flow (the combination of interflow and groundwater flow) makes up most of the total streamflow (99 or more percent of total streamflow for 71 percent of the time during water years 2002–2008).

Montana↗