Search USGS⌕ Search

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.

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

At least 1,531 records · Page 85Linked to original sources
USGS2004

Application of health-based screening levels to ground-water quality data in a state-scale pilot effort

A state-scale pilot effort was conducted to evaluate a Health-Based Screening Level (HBSL) approach developed for communicating findings from the U.S. Geological Survey (USGS) National Water-Quality Assessment Program in a human-health context. Many aquifers sampled by USGS are used as drinking-water sources, and water-quality conditions historically have been assessed by comparing measured contaminant concentrations to established drinking-water standards and guidelines. Because drinking-water standards and guidelines do not exist for many analyzed contaminants, HBSL values were developed collaboratively by the USGS, U.S. Environmental Protection Agency (USEPA), New Jersey Department of Environmental Protection, and Oregon Health & Science University, using USEPA toxicity values and USEPA Office of Water methodologies. The main objective of this report is to demonstrate the use of HBSL approach as a tool for communicating water-quality data in a human-health context by conducting a retrospective analysis of ground-water quality data from New Jersey. Another important objective is to provide guidance on the use and interpretation of HBSL values and other human-health benchmarks in the analyses of water-quality data in a human-health context. Ground-water samples collected during 1996-98 from 30 public-supply, 82 domestic, and 108 monitoring wells were analyzed for 97 pesticides and 85 volatile organic compounds (VOCs). The occurrence of individual pesticides and VOCs was evaluated in a human-health context by calculating Benchmark Quotients (BQs), defined as ratios of measured concentrations of regulated compounds (that is, compounds with Federal or state drinking-water standards) to Maximum Contaminant Level (MCL) values and ratios of measured concentrations of unregulated compounds to HBSL values. Contaminants were identified as being of potential human-health concern if maximum detected concentrations were within a factor of 10 of the associated MCL or HBSL (that is, maximum BQ value (BQmax) greater than or equal to 0.1) in any well type (public supply, domestic, monitoring). Most (57 of 77) pesticides and VOCs with human-health benchmarks were detected at concentrations well below these levels (BQmax less than 0.1) for all three well types; however, BQmax values ranged from 0.1 to 3,000 for 6 pesticides and 14 VOCs. Of these 20 contaminants, one pesticide (dieldrin) and three VOCs (1,2-dibromoethane, tetrachloroethylene, and trichloroethylene) both (1) were measured at concentrations that met or exceeded MCL or HBSL values, and (2) were detected in more than 10 percent of samples collected from raw ground water used as sources of drinking water (public-supply and (or) domestic wells) and, therefore, are particularly relevant to human health. The occurrence of multiple pesticides and VOCs in individual wells also was evaluated in a human-health context because at least 53 different contaminants were detected in each of the three well types. To assess the relative human-health importance of the occurrence of multiple contaminants in different wells, the BQ values for all contaminants in a given well were summed. The median ratio of the maximum BQ to the sum of all BQ values for each well ranged from 0.83 to 0.93 for all well types, indicating that the maximum BQ makes up the majority of the sum for most wells. Maximum and summed BQ values were statistically greater for individual public-supply wells than for individual domestic and monitoring wells. The HBSL approach is an effective tool for placing water-quality data in a human-health context. For 79 of the 182 compounds analyzed in this study, no USEPA drinking-water standards or guidelines exist, but new HBSL values were calculated for 39 of these 79 compounds. The new HBSL values increased the number of detected pesticides and VOCs with human-health benchmarks from 65 to 77 (of 97 detected compounds), thereby expanding the basis for interpreting contaminant-occu

Scientific Investigations Report↗
USGS2021

Hydrogeology and model-simulated groundwater availability in the Salt Fork Red River aquifer, southwestern Oklahoma, 1980–2015

The 1973 Oklahoma Water Law (82 OK Stat § 82-1020.5) requires that the Oklahoma Water Resources Board (OWRB) conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin (hereinafter referred to as an “aquifer”). The maximum annual yield allocated per acre of land is known as the equal-proportionate-share (EPS) pumping rate. At present (2021), the OWRB has not yet established a maximum annual yield and EPS pumping rate for the Salt Fork Red River aquifer. To provide updated information to the OWRB that could support evaluation and determination of an appropriate maximum annual yield, the U.S. Geological Survey (USGS), in cooperation with the OWRB, conducted a hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater availability in the Salt Fork Red River aquifer. The Salt Fork Red River aquifer in Greer, Harmon, and Jackson Counties of southwestern Oklahoma is composed of about 274.5 square miles of alluvium and terrace deposits associated with the Salt Fork Red River. The mean annual recharge rate to the Salt Fork Red River aquifer for the period 1980–2015 was estimated to be about 2.94 inches per year, or 10.0 percent of the mean annual precipitation for the same period (29.4 inches per year). This 1980–2015 mean annual recharge rate is equivalent to a mean annual recharge rate of about 38,000 acre-feet per year (acre-ft/yr) for the Salt Fork Red River aquifer excluding about 19,764 acres comprising the Mulberry Creek and Horse Creek terraces. The mean annual recharge rates upgradient and downgradient from USGS streamgage 07300500 Salt Fork Red River at Mangum, Okla. (hereinafter referred to as the “Mangum gage”), apportioned by aquifer area (41.5 and 58.5 percent, respectively), were about 16,000 and 22,000 acre-ft/yr, respectively. Mean annual groundwater use for the study period (1980–2015) was 3,532.7 acre-ft/yr; about 77 percent of that groundwater use was for irrigation, and about 23 percent was for public supply. Most groundwater use for irrigation was associated with wells in the Martha terrace. A hydrogeologic framework was developed for the Salt Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of the numerical groundwater-flow model of the Salt Fork Red River aquifer described in this report. A conceptual model for the Salt Fork Red River aquifer that reasonably represents the groundwater-flow system was developed to constrain the construction and calibration of the numerical model. The conceptual-model water budget estimated mean annual inflows to, and outflows from, the Salt Fork Red River aquifer for the period 1980–2015 and included a subaccounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from the Mangum gage. The numerical groundwater-flow model of the Salt Fork Red River aquifer was constructed by using MODFLOW-2005 with the Newton formulation solver. The model of the Salt Fork Red River aquifer was spatially discretized into 1,050 rows, 1,125 columns, about 170,000 active cells measuring 200 by 200 feet (ft), and a single convertible layer. The model was temporally discretized into 432 monthly transient stress periods (each with two time steps to improve model stability). An initial steady-state stress period represented mean annual inflows to, and outflows from, the aquifer and produced a solution that was used as the initial condition for subsequent transient stress periods as well as some groundwater-availability scenarios. The model was calibrated to water-table-altitude observations at selected wells and base-flow observations at selected streamgages. The simulated saturated thickness of the Salt Fork Red River aquifer was determined by subtracting the altitude of the aquifer base from the simulated water-table altitude at the end of the numerical-model period (2015). The simulated saturated thickness was more than 75 ft in a paleochannel in the Dodson terrace near the Texas border. The mean aquifer thickness (sum of saturated and unsaturated) was 49.62 ft, and the mean saturated thickness was 28.55 ft. A simulated mean transmissivity of 1,024 feet squared per day was computed from the calibrated hydraulic conductivity and saturated thickness of each cell. The simulated available water in storage at the end of the numerical-model period (2015) was 526,117 acre-feet (acre-ft); about 42 percent of that total was available upgradient from the Mangum gage, and about 58 percent of that total was available downgradient from the Mangum gage (including the Mangum terrace). Three types of groundwater-availability scenarios were run using the calibrated numerical model. These scenarios were used to (1) estimate the EPS pumping rate that ensures a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought on base flow and groundwater storage. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were about 0.51, 0.48, and 0.48 acre-foot per acre per year, respectively. Given the 155,929-acre modeled aquifer area, these rates correspond to annual yields of about 78,800, 74,900, and 74,700 acre-ft/yr, respectively. For the 20-year EPS scenario, decreasing and increasing recharge by 10 percent resulted in a 6-percent change in the EPS pumping rate in both cases; for the 40- and 50-year EPS scenarios, decreasing and increasing recharge by 10 percent resulted in a 7-percent change in the EPS pumping rate in both cases. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Salt Fork Red River aquifer and base flows in the Salt Fork Red River. The effects of well withdrawals were evaluated by quantifying differences in groundwater storage and base flow in four 50-year scenarios, which applied (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2015), (3) 2015 pumping rates, and (4) increasing demand pumping rates at simulated wells. The increasing demand pumping rates assumed a cumulative 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwestern Oklahoma. Groundwater storage after 50 years with no pumping was 535,000 acre-ft, or 8,900 acre-ft (1.7 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-table-altitude increase of 0.48 ft. Groundwater storage after 50 years of pumping at the mean rate for the study period (1980–2015) was 519,900 acre-ft, or 6,200 acre-ft (1.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.34 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 513,100 acre-ft, or 13,000 acre-ft (2.5 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.71 ft. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was 509,700 acre-ft, or 16,500 acre-ft (3.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.89 ft. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage. The period January 1983–December 1992 was chosen as the simulated drought period. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought) at the end of the simulated drought period (1992). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–92). Upstream inflows from the Salt Fork Red River, Turkey Creek, and Bitter Creek were reduced by 75 percent. Groundwater storage at the end of the drought period (1992) was 479,200 acre-ft, or 53,200 acre-ft (10.0 percent) less than the groundwater storage of the calibrated numerical model at the end of the drought period. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 2.9 ft. At the end of the 10-year hypothetical drought period, simulated base flows at the Mangum gage and USGS streamgage 07301110 Salt Fork Red River near Elmer, Okla., had decreased by about 80 and 70 percent, respectively.

Oklahoma↗
USGS2006

Water-quality trend analysis and sampling design for the Devils Lake Basin, North Dakota, January 1965 through September 2003

This report presents the results of a study conducted by the U.S. Geological Survey, in cooperation with the North Dakota State Water Commission, the Devils Lake Basin Joint Water Resource Board, and the Red River Joint Water Resource District, to analyze historical water-quality trends in three dissolved major ions, three nutrients, and one dissolved trace element for eight stations in the Devils Lake Basin in North Dakota and to develop an efficient sampling design to monitor the future trends. A multiple-regression model was used to detect and remove streamflow-related variability in constituent concentrations. To separate the natural variability in concentration as a result of variability in streamflow from the variability in concentration as a result of other factors, the base-10 logarithm of daily streamflow was divided into four components-a 5-year streamflow anomaly, an annual streamflow anomaly, a seasonal streamflow anomaly, and a daily streamflow anomaly. The constituent concentrations then were adjusted for streamflow-related variability by removing the 5-year, annual, seasonal, and daily variability. Constituents used for the water-quality trend analysis were evaluated for a step trend to examine the effect of Channel A on water quality in the basin and a linear trend to detect gradual changes with time from January 1980 through September 2003. The fitted upward linear trends for dissolved calcium concentrations during 1980-2003 for two stations were significant. The fitted step trends for dissolved sulfate concentrations for three stations were positive and similar in magnitude. Of the three upward trends, one was significant. The fitted step trends for dissolved chloride concentrations were positive but insignificant. The fitted linear trends for the upstream stations were small and insignificant, but three of the downward trends that occurred during 1980-2003 for the remaining stations were significant. The fitted upward linear trends for dissolved nitrite plus nitrate as nitrogen concentrations during 1987-2003 for two stations were significant. However, concentrations during recent years appear to be lower than those for the 1970s and early 1980s but higher than those for the late 1980s and early 1990s. The fitted downward linear trend for dissolved ammonia concentrations for one station was significant. The fitted linear trends for total phosphorus concentrations for two stations were significant. Upward trends for total phosphorus concentrations occurred from the late 1980s to 2003 for most stations, but a small and insignificant downward trend occurred for one station. Continued monitoring will be needed to determine if the recent trend toward higher dissolved nitrite plus nitrate as nitrogen and total phosphorus concentrations continues in the future. For continued monitoring of water-quality trends in the upper Devils Lake Basin, an efficient sampling design consists of five major-ion, nutrient, and trace-element samples per year at three existing stream stations and at three existing lake stations. This sampling design requires the collection of 15 stream samples and 15 lake samples per year rather than 16 stream samples and 20 lake samples per year as in the 1992-2003 program. Thus, the design would result in a program that is less costly and more efficient than the 1992-2003 program but that still would provide the data needed to monitor water-quality trends in the Devils Lake Basin.

North Dakota↗
USGS2010

Occurrence of organic wastewater-indicator compounds in urban streams of the Atlanta area, Georgia, 2003-2006

Between March 2003 and January 2006, 863 water samples were collected from streams in seven urban watersheds with varying land uses within or near the City of Atlanta, Georgia. Sixty-four sampling sites representing three site types were established in those watersheds. The first type consisted of sites within three watersheds not affected by combined sewer overflows; these were designated as the control basins. The second and third site types were established in four watersheds and were designated as sites upstream or downstream from combined sewer outfalls. Stream samples collected during the study were analyzed for major ions, nutrients, trace metals, and 60 organic compounds commonly found in wastewater (organic wastewater-indicator compounds, OWICs). Inorganic constituents were analyzed to discern possible relations between OWICs and urban runoff, sewage effluent, or combined sewer overflows (CSOs). The OWICs were grouped into nine compound classes based either on an already accepted class of compounds (such as pesticide) or the manner in which the compounds are used (such as automotive uses). The compounds benzo(a)pyrene, 4-cumylphenol, 3-tert-butyl-4-hydroxyanisole (BHA), isophorone, isoquinoline, metolachlor, metalaxyl, and 4-octylphenol were not detected in any sample collected during the study. As many as 33 OWICs were detected above study reporting levels in water samples collected from streams in the Intrenchment Creek, Peachtree Creek, Proctor Creek, and South River watersheds (basins with CSOs), a number markedly higher than the number detected in water samples from the control basins (watersheds without CSOs). Several compounds known to disrupt the endocrine systems of aquatic biota were among the compounds detected. The median numbers of OWICs detected in base-flow samples from the control basins ranged from 3 to 4 and 7 to 9 in stormflow samples, while the median in base-flow samples from the four CSO-affected watersheds ranged from 4 to 16 and 11 to 19 in stormflow samples. The detection frequencies and concentrations of OWICs in water samples varied depending on flow conditions during sample collection; however, regardless of flow condition, the total OWICs concentrations were strongly related to the numbers of OWICs detected in these samples. In addition, the median number of detectable OWICs and total OWIC concentrations increased linearly with increasing impervious area and stream flashiness (flashiness indicates the rapidity with which streamflow responds to high rainfall amounts). Four compounds-tris(2-butoxyethyl) phosphate (TBEP), tris(2-chloroethyl) phosphate (TCEP), bromacil, and cholesterol-were detected at concentrations greater than study reporting levels in at least 45 percent of all samples collected during the study. On a broad scale, the seasonal distributions of OWICs detected in samples collected during the study period were consistent with use patterns or urban activity, but were markedly different and more variable within individual watersheds. Seven of the nine OWIC classes were detected with greater frequency in base-flow samples from sites downstream from CSOs than from those upstream from CSOs or from control-basin sites; these OWICs also were detected in a greater percentage of base-flow samples from upstream than control-basin sites. Polycyclic aromatic hydrocarbon (PAH) and automotive-use compounds were detected with similar frequency in base-flow samples from all sites. The pesticides and industrial-use compounds were detected with similar frequency in base-flow samples from sites upstream and downstream from CSOs. The compounds 7-acetyl-1,1,3,4,4,6-hexamethyl tetrahydronaphthalene (AHTN), 1,3,4,6,7,8-hexahydro-4,6,6,7,8,8-hexamethyl cyclopenta-g-2-benzopyran (HHCB), and triclosan were detected with greater frequency in base-flow than in stormflow samples from all sites. The detection frequencies of these three compounds were particularly high in base-flow samples from downstream sites, especially those from the Intrenchment Creek watershed. In stormflow samples, only the industrial-use compounds were detected with similar frequency among samples from control basins, and upstream and downstream sites. Caffeine, camphor, and menthol were detected in a greater percentage of stormflow than base-flow samples from all sites. The disinfectant byproduct bromoform was detected with the highest frequency in base-flow samples from the upstream sites, particularly those from the South River watershed. Typically, compounds in the pesticide class were detected with similar frequency in base-flow and stormflow samples from upstream and downstream sites, although bromacil and carbaryl were the exceptions. Bromacil was detected with greatest frequency in base-flow samples from all sites, but was detected in a larger percentage of base-flow samples from upstream and downstream sites, especially upstream and downstream sites in the Proctor Creek and South River watersheds. Carbaryl, however, was detected in a greater percentage of stormflow samples from all sites. Although collectively the industrial-use compounds were detected in more stormflow than base-flow samples, tetrachloroethene (PCE) was detected in more base-flow than stormflow samples at all sites. More specifically, PCE was detected with the highest frequency in samples from the upstream sites, particularly those from the Proctor Creek watershed. The similarity in the pattern and distribution of OWICs in samples at sites upstream and downstream from known CSO outfalls indicates that CSOs were not the dominant source of OWICs during the study period. Other sources may include non-sewage discharges-both permitted, permitted but out of compliance, and non-permitted, contaminated groundwater from leaking sewer lines or septic systems, sanitary-sewer overflows, or dry-weather runoff from outdoor water use. These OWICs may be better suited for identifying sewage-contaminated groundwater than sewage-contaminated surface water because groundwater is not typically affected by the OWICs that are more common in urban runoff.

Georgia↗
USGS2005

Ground-water hydrology of the Willamette basin, Oregon

The Willamette Basin encompasses a drainage of 12,000 square miles and is home to approximately 70 percent of Oregon's population. Agriculture and population are concentrated in the lowland, a broad, relatively flat area between the Coast and Cascade Ranges. Annual rainfall is high, with about 80 percent of precipitation falling from October through March and less than 5 percent falling in July and August, the peak growing season. Population growth and an increase in cultivation of crops needing irrigation have produced a growing seasonal demand for water. Because many streams are administratively closed to new appropriations in summer, ground water is the most likely source for meeting future water demand. This report describes the current understanding of the regional ground-water flow system, and addresses the effects of ground-water development. This study defines seven regional hydrogeologic units in the Willamette Basin. The highly permeable High Cascade unit consists of young volcanic material found at the surface along the crest of the Cascade Range. Four sedimentary hydrogeologic units fill the lowland between the Cascade and Coast Ranges. Young, highly permeable coarse-grained sediments of the upper sedimentary unit have a limited extent in the floodplains of the major streams and in part of the Portland Basin. Extending over much of the lowland where the upper sedimentary unit does not occur, silts and clays of the Willamette silt unit act as a confining unit. The middle sedimentary unit, consisting of permeable coarse-grained material, occurs beneath the Willamette silt and upper sedimentary units and at the surface as terraces in the lowland. Beneath these units is the lower sedimentary unit, which consists of predominantly fine-grained sediments. In the northern part of the basin, lavas of the Columbia River basalt unit occur at the surface in uplands and beneath the basin-fill sedimentary units. The Columbia River basalt unit contains multiple productive water-bearing zones. A basement confining unit of older marine and volcanic rocks of low permeability underlies the basin and occurs at land surface in the Coast Range and western part of the Cascade Range. Most recharge in the basin is from infiltration of precipitation, and the spatial distribution of recharge mimics the distribution of precipitation, which increases with elevation. Basinwide annual mean recharge is estimated to be 22 inches. Rain and snowmelt easily recharge into the permeable High Cascade unit and discharge within the High Cascade area. Most recharge in the Coast Range and western part of the Cascade Range follows short flowpaths through the upper part of the low permeability material and discharges to streams within the mountains. Consequently, recharge in the Coast and Ranges is not available as lateral ground-water flow into the lowland, where most ground-water use occurs. Within the lowland, annual mean recharge is 16 inches and most recharge occurs from November to April, when rainfall is large and evapotranspiration is small. From May to October recharge is negligible because precipitation is small and evapotranspiration is large. Discharge of ground water is mainly to streams. Ground-water discharge is a relatively large component of flow in streams that drain the High Cascade unit and parts of the Portland Basin where permeable units are at the surface. In streams that do not head in the High Cascade area, streamflow is generally dominated by runoff of precipitation. Ground-water in the permeable units in the lowland discharges to the major streams where there is a good hydraulic connection between aquifers and streams. Ground-water discharge to smaller streams, which flow on the less permeable Willamette silt unit, is small and mostly from the Willamette silt unit. Most ground-water withdrawals occur within the lowland. Irrigation is the largest use of ground water, accounting for 240,000 acre feet of withdrawals, or 81 p

Scientific Investigations Report↗
USGS2024

Analysis of water use associated with hydraulic fracturing and determination of baseline water quality in watersheds within the shale play of eastern Ohio, 2021–23

The U.S. Geological Survey, in cooperation with the Ohio Department of Natural Resources, performed a two-part study to (1) assess water use and temporal trends and changes in streamflow, and to (2) characterize 2021–23 baseline water quality, as they relate to oil and gas extraction activities in selected eastern Ohio watersheds. Between calendar years 2010 and 2019, hydraulic fracturing water withdrawals totaling about 27,168 million gallons were reported at 643 locations in Ohio. In 2021, wells developed with hydraulic fracturing were the source of most of the oil and gas produced in Ohio. Daily streamflow time-series data from seven study gages and two reference gages were used to assess temporal trends and changes in streamflow. The study gages were in basins with reported water withdrawals for hydraulic fracturing. The reference gages, which have long periods of record and were subject to minimal streamflow regulation, were in nearby basins with no hydraulic fracturing water withdrawals. Trend slopes for the period of record annual minimum and median daily streamflows and for annual daily streamflow nonexceedance probabilities less than 0.9 were all uniformly positive at the study and reference gages. This trend indicates a consistently increasing pattern over the periods of record, except for high flows. In addition, analyses of annual streamflow statistics showed no general indication that low flows or extreme low flows at the reference or study gages have lowered, become more frequent, or lengthened in duration since 2010, when records for hydraulic fracturing water withdrawals began in Ohio. In fact, in almost all cases, the opposite was indicated. Nonexceedance percentiles of daily streamflows were compared between the full and pre-2012 periods of record for reference and study gages. The streamflows associated with nonexceedance percentiles in the lower quartile of daily streamflows determined for the full period of record were larger than or equal to those determined for the pre-2012 period of record for all study and reference gages. This indicates that low flows did not decrease during the post-2011 period of record when water was withdrawn for hydraulic fracturing. Water-quality data were collected eight times at each of eight sampling sites (six of which were colocated with the study gages). Sampling was done during a variety of flow conditions to assess baseline water quality. In 2021, the 8 sampling sites had drainage basins that were wholly or partially within 7 of the 10 most active counties in Ohio for oil and gas development. As part of the record of baseline conditions, water-quality data were used to assess (1) water types based on major-ion chemistry; (2) sources of salinity to streams; (3) exceedances of aquatic life use criteria; and (4) the correlations between water chemistry and drainage-basin characteristics, such as density of oil and gas wells, density of wastewater treatment plants, or the percentage of different types of land cover (agriculture, developed, forest). Seven of the water-quality sampling sites were designated as coal-mine impacted based on criteria developed for assessing mine-drainage impacts in Ohio. Mine drainage from historical coal mining in the region likely affected the quality of these streams and complicated the use of some constituents typically used as indicators of oil and gas influence. Based on major-ion chemistry, three main types of water were in the study area―sulfate (three sites), calcium-bicarbonate (one site), and mixed bicarbonate-chloride (four sites) type waters. One site had samples with a higher proportion of sodium and chloride ions than other stream samples, indicating potential contamination with oil-field brine or road salt. Binary mixing curves revealed that 11 samples from 4 of the sampling sites likely contained a component of brine. The results of the baseline assessment of surface-water quality in the study area showed no exceedances of Ohio Environmental Protection Agency aquatic life use criteria. Spearman’s rank correlation coefficients indicated no significant positive correlations with the density of vertical or horizontal oil and gas wells.

Ohio↗
USGS2009

Occurrence of volatile organic compounds in selected urban streams in the United States, 1995-2003

As part of the U.S. Geological Survey's (USGS) National Water-Quality Assessment (NAWQA) Program, urban indicator sites were monitored to (1) characterize the stream quality from drainage basins with predominantly residential and commercial land use, and (2) determine which selected natural and anthropogenic factors affect stream quality. A total of 869 water samples were collected from 37 urban streams during 1995-2003 and were analyzed for 87 volatile organic compounds (VOCs). The occurrence of VOCs in urban streams is described in this report for (1) all samples as a single dataset, (2) all samples grouped by streamflow pentiles, and (3) all samples grouped by warmer (April through September) and cooler (October through March) months by the detection frequency and (or) concentration of (a) any VOC, (b) VOC groups, and (c) individual compounds. An assessment level of 0.02 microgram per liter (ug/L) was used to compute the detection frequencies and concentrations of VOCs. Concentrations of VOCs were compared to (1) U.S. Environmental Protection Agency's (USEPA) drinking-water Maximum Contaminant Levels (MCLs) or Drinking Water Advisories, (2) Health-Based Screening Levels (HBSLs) developed by the USGS in collaboration with the USEPA and other agencies, and (3) USEPA and Canadian aquatic-life criteria. One or more VOCs were detected in 97.1 percent of 869 samples, and one or more VOCs were detected frequently (greater than 80 percent) at all sites. The median total VOC concentration for all samples was 0.57 ug/L, and total VOC concentrations in a single sample ranged from not detected to 698 ug/L. About 85 percent of the samples contained two or more VOCs, and about one-half contained five or more VOCs. The gasoline hydrocarbons were the most frequently occurring VOC group followed by solvents, trihalomethanes (THMs), gasoline oxygenates, organic synthesis compounds, fumigants, and refrigerants. Concentration ranges for most VOC groups were distributed over at least two orders of magnitude. Fifty-seven of the 87 VOCs analyzed were detected in at least one sample at an assessment level of 0.02 ug/L. More than one-half of the 30 VOCs not detected in samples were organic synthesis compounds. Fifteen compounds had detection frequencies greater than or equal to 10 percent. With the exception of toluene and chloroform, the median concentration of each VOC for all samples was less than the assessment level. Furthermore, the median concentrations of detections for the 15 most frequently occurring VOCs ranged from 0.03 to 3.9 ug/L, and typically were less than or equal to 0.10 ug/L. The 869 samples from the 37 sites were stratified into five streamflow pentiles (less than 20, 20-less than 40, 40-less than 60, 60-less than 80, and greater than or equal to 80 percent of estimated long-term streamflow statistics) for comparison of the occurrence of VOCs. The detection frequency of one or more VOCs by streamflow pentile varied only slightly from 96.7 to 97.7 percent. The median total VOC concentrations in samples for the five streamflow pentiles ranged from 0.39 to 1.0 ug/L. Two or more VOCs were present in more than 80 percent of samples in each of the five pentiles. The gasoline hydrocarbons, solvents, THMs, and gasoline oxygenates occurred frequently (greater than 30 percent) in all streamflow pentiles, in contrast to the organic synthesis compounds, fumigants, and refrigerants that occurred less frequently in urban streams under all streamflow conditions. The median total VOC concentrations for gasoline hydrocarbons, solvents, gasoline oxygenates, and organic synthesis compounds generally increased as streamflow increased. In contrast, the median total VOC concentrations for THMs and fumigants generally decreased as streamflow increased. The median total VOC concentrations for refrigerants showed no pattern as streamflow increased. Because differences between VOC occurrence and streamflow pentiles were small for most compariso

Scientific Investigations Report↗
USGS2010

Processes of Terrace Formation on the Piedmont of the Santa Cruz River Valley During Quaternary Time, Green Valley-Tubac Area, Southeastern Arizona

In this report we describe a series of stepped Quaternary terraces on some piedmont tributaries of the Santa Cruz River valley in southeastern Arizona. These terraces began to form in early Pleistocene time, after major basin-and-range faulting ceased, with lateral planation of basin fill and deposition of thin fans of alluvium. At the end of this cycle of erosion and deposition, tributaries of the Santa Cruz River began the process of dissection and terrace formation that continues to the present. Vertical cutting alternated with periods of equilibrium, during which streams cut laterally and left thin deposits of channel fill. The distribution of terraces was mapped and compiled with adjacent mapping to produce a regional picture of piedmont stream history in the middle part of the Santa Cruz River valley. For selected tributaries, the thickness of terrace fill was measured, particle size and lithology of gravel were determined, and sedimentary features were photographed and described. Mapping of terrace stratigraphy revealed that on two tributaries, Madera Canyon Wash and Montosa Canyon Wash, stream piracy has played an important role in piedmont landscape development. On two other tributaries, Cottonwood Canyon Wash and Josephine Canyon Wash, rapid downcutting preempted piracy. Two types of terraces are recognized: erosional and depositional. Gravel in thin erosional terraces has Trask sorting coefficients and sedimentary structures typical of streamflood deposits, replete with bar-and-swale surface topography on young terraces. Erosional-terrace fill represents the channel fill of the stream that cuts the terrace; the thickness of the fill indicates the depth of channel scour. In contrast to erosional terraces, depositional terraces show evidence of repeated deposition and net aggradation, as indicated by their thickness (as much as 20+ m) and weakly bedded structure. Depositional terraces are common below mountain-front canyon mouths where streams drop their load in response to abrupt flattening of gradients and expansion of channel banks, and they extend down the piedmont along Josephine Canyon Wash. Gravel in depositional terraces also has sorting coefficients typical of streamflood deposits. Sedimentary features in both types of terraces are consistent with deposition by flash floods in ephemeral streams, suggesting the climate was arid. Bedding and clast armor are weakly developed, clast clusters and imbrication are common, and crossbedding is generally absent. Debris-flow deposits, even near the mountain front, are surprisingly rare. On the tectonically stable piedmont of southeastern Arizona, stream piracy and climate change are the most likely agents of terrace formation. Both piracy and climate change can cause rapid changes in discharge and sediment supply, which initiate cycles of incision, lateral cutting, and aggradation. Increased stream discharge initiates downcutting, but increased sediment supply interrupts downcutting and causes streams to cut laterally and aggrade. At times, on Madera Canyon Wash and Montosa Canyon Wash, stream piracy affected stream discharge and sediment supply, but on Cottonwood Canyon Wash and Josephine Canyon Wash, only climate change could have initiated terrace cutting. Terraces probably formed during extended arid intervals when sparse vegetation and flashy stream discharge combined to increase sediment supply. In most cases, sediment supply was sufficient to promote lateral cutting but not long-term aggradation. Thus, most streams formed erosional terraces. The middle Pleistocene Josephine Canyon Wash formed a depositional terrace because it had a source of abundant unconsolidated sediment.

Scientific Investigations Report↗
USGS2015

Estimated water use in Arkansas, 2010

The Arkansas Natural Resources Commission (ANRC) conducts an annual inventory of reported groundwater and surface-water withdrawals in Arkansas in cooperation with the U.S. Geological Survey (USGS). This report describes withdrawals from groundwater and surface-water resources in Arkansas for 2010. The report compiles withdrawals by county for 10 categories of water use—public supply, domestic (self-supplied), commercial (self-supplied), industrial (self-supplied), mining, livestock, aquaculture, irrigation, duck (hunting) clubs, and thermoelectric power generation. Water-use trends in Arkansas from 1965 to 2010 and sources of groundwater withdrawals also are described. During 2010, total withdrawals from groundwater and surface-water sources in Arkansas were 11,300 million gallons per day (Mgal/d). Of the total withdrawn, about 69 percent (7,790 Mgal/d) was from groundwater. Public-supply systems served about 94 percent of Arkansas’ population. Public-supply total withdrawals were 429 Mgal/d, with about 69 percent from surface-water sources. The statewide average of per capita residential use from public-supply systems was about 155 gallons per day (gal/d). The domestic (self-supplied) water use was 12.8 Mgal/d. Total commercial (self-supplied) water use was 11.7 Mgal/d. Total industrial (self-supplied) water use was 276 Mgal/d. Total mining water use was 44.3 Mgal/d. Total livestock water use was 39.0 Mgal/d. Total aquaculture water use was 268 Mgal/d. Irrigation water use totaled 8,720 Mgal/d. Total duck (hunting) club water use was 216 Mgal/d. Total thermoelectric power water use was 1,540 Mgal/d. The three water-use categories with the largest withdrawals and their effects on total water use were examined. Total water use in Arkansas has increased about 428 percent between 1965 and 2010. Total groundwater use increased about 533 percent and total surface-water use increased about 289 percent. Since about 2000, total water use in Arkansas has plateaued, peaking in 2005. The plateauing of total water use is the result of decreasing surface-water use in the irrigation, thermoelectric power, and public-supply categories offsetting the continuing increases in groundwater use in the irrigation category. An examination of total water use and irrigation water-use changes over time demonstrates the increasing dominance of groundwater withdrawals for irrigation on Arkansas’ total water use. Total irrigation water use in Arkansas between 1965 and 2010 has increased about 652 percent. Withdrawals for thermoelectric power water use in Arkansas have continuously accounted for about half of the State’s total surface-water use for the period 1965–2010. Thermoelectric power water use in Arkansas between 1965 and 2010 increased 264 percent. The percent of Arkansas’ population served by public water suppliers has continued to increase while the percentage of the State’s total water use withdrawn by public water suppliers has remained relatively constant. Public-supply water use in Arkansas between 1965 and 2010 increased about 238 percent. Regardless of continuing increases in population, since about 2000, total public-supply water use in Arkansas has plateaued at about 425 Mgal/d. Groundwater withdrawals comprised about 69 percent of the total amount of water used in Arkansas in 2010. Four aquifers in Arkansas account for more than 99 percent of the total groundwater withdrawals. The aquifers in deposits of Quaternary age supplied about 97 percent of all groundwater withdrawals. The Sparta-Memphis aquifer supplied about 2.5 percent of all groundwater withdrawals, the Wilcox aquifer supplied about 0.5 percent of all groundwater withdrawals, and the Paleozoic aquifer supplied about 0.3 percent of all groundwater withdrawals.

Arkansas↗
USGS2009

Ground-Water Budgets for the Wood River Valley Aquifer System, South-Central Idaho, 1995-2004

The Wood River Valley contains most of the population of Blaine County and the cities of Sun Valley, Ketchum, Haley, and Bellevue. This mountain valley is underlain by the alluvial Wood River Valley aquifer system which consists of a single unconfined aquifer that underlies the entire valley, an underlying confined aquifer that is present only in the southernmost valley, and the confining unit that separates them. The entire population of the area depends on ground water for domestic supply, either from domestic or municipal-supply wells, and rapid population growth since the 1970s has caused concern about the long-term sustainability of the ground-water resource. To help address these concerns this report describes a ground-water budget developed for the Wood River Valley aquifer system for three selected time periods: average conditions for the 10-year period 1995-2004, and the single years of 1995 and 2001. The 10-year period 1995-2004 represents a range of conditions in the recent past for which measured data exist. Water years 1995 and 2001 represent the wettest and driest years, respectively, within the 10-year period based on precipitation at the Ketchum Ranger Station. Recharge or inflow to the Wood River Valley aquifer system occurs through seven main sources (from largest to smallest): infiltration from tributary canyons, streamflow loss from the Big Wood River, areal recharge from precipitation and applied irrigation water, seepage from canals and recharge pits, leakage from municipal pipes, percolation from septic systems, and subsurface inflow beneath the Big Wood River in the northern end of the valley. Total estimated mean annual inflow or recharge to the aquifer system for 1995-2004 is 270,000 acre-ft/yr (370 ft3/s). Total recharge for the wet year 1995 and the dry year 2001 is estimated to be 270,000 acre-ft/yr (370 ft3/s) and 220,000 acre-ft/yr (300 ft3/s), respectively. Discharge or outflow from the Wood River Valley aquifer system occurs through five main sources (from largest to smallest): Silver Creek streamflow gain, ground-water pumpage, Big Wood River streamflow gain, direct evapotranspiration from riparian vegetation, and subsurface outflow (treated separately). Total estimated mean 1995-2004 annual outflow or discharge from the aquifer system is 250,000 acre-ft/yr (350 ft3/s). Estimated total discharge is 240,000 acre-ft/yr (330 ft3/s) for both the wet year 1995 and the dry year 2001. The budget residual is the difference between estimated ground-water inflow and outflow and encompasses subsurface outflow, ground-water storage change, and budget error. For 1995-2004, mean annual inflow exceeded outflow by 20,000 acre-ft/yr (28 ft3/s); for the wet year 1995, mean annual inflow exceeded outflow by 30,000 acre-ft/yr (41 ft3/s); for the dry year 2001, mean annual outflow exceeded inflow by 20,000 acre-ft/yr (28 ft3/s). These values represent 8, 13, and 8 percent, respectively, of total outflows for the same periods. It is difficult to differentiate the relative contributions of the three residual components, although the estimated fluctuations between the wet and dry year budgets likely are primarily caused by changes in ground-water storage. The individual components in the wet and dry year ground-water budgets responded in a consistent manner to changes in precipitation and temperature. Although the ground-water budgets for the three periods indicated that ground-water storage is replenished in wet years, statistical analyses by Skinner and others (2007) suggest that such replenishment is not complete and over the long term more water is removed from storage than is replaced. In other words, despite restoration of water to ground-water storage in wet years, changes have occurred in either recharge and (or) discharge to cause ground-water storage to decline over time. Such changes may include, but are not limited to: lining or abandoning canals and ditches, conversion of surface-water irriga

Scientific Investigations Report↗
USGS2009

Numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units in the Rapid City area, South Dakota

The city of Rapid City and other water users in the Rapid City area obtain water supplies from the Minnelusa and Madison aquifers, which are contained in the Minnelusa and Madison hydrogeologic units. A numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units in the Rapid City area was developed to synthesize estimates of water-budget components and hydraulic properties, and to provide a tool to analyze the effect of additional stress on water-level altitudes within the aquifers and on discharge to springs. This report, prepared in cooperation with the city of Rapid City, documents a numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units for the 1,000-square-mile study area that includes Rapid City and the surrounding area. Water-table conditions generally exist in outcrop areas of the Minnelusa and Madison hydrogeologic units, which form generally concentric rings that surround the Precambrian core of the uplifted Black Hills. Confined conditions exist east of the water-table areas in the study area. The Minnelusa hydrogeologic unit is 375 to 800 feet (ft) thick in the study area with the more permeable upper part containing predominantly sandstone and the less permeable lower part containing more shale and limestone than the upper part. Shale units in the lower part generally impede flow between the Minnelusa hydrogeologic unit and the underlying Madison hydrogeologic unit; however, fracturing and weathering may result in hydraulic connections in some areas. The Madison hydrogeologic unit is composed of limestone and dolomite that is about 250 to 610 ft thick in the study area, and the upper part contains substantial secondary permeability from solution openings and fractures. Recharge to the Minnelusa and Madison hydrogeologic units is from streamflow loss where streams cross the outcrop and from infiltration of precipitation on the outcrops (areal recharge). MODFLOW-2000, a finite-difference groundwater-flow model, was used to simulate flow in the Minnelusa and Madison hydrogeologic units with five layers. Layer 1 represented the fractured sandstone layers in the upper 250 ft of the Minnelusa hydrogeologic unit, and layer 2 represented the lower part of the Minnelusa hydrogeologic unit. Layer 3 represented the upper 150 ft of the Madison hydrogeologic unit, and layer 4 represented the less permeable lower part. Layer 5 represented an approximation of the underlying Deadwood aquifer to simulate upward flow to the Madison hydrogeologic unit. The finite-difference grid, oriented 23 degrees counterclockwise, included 221 rows and 169 columns with a square cell size of 492.1 ft in the detailed study area that surrounded Rapid City. The northern and southern boundaries for layers 1-4 were represented as no-flow boundaries, and the boundary on the east was represented with head-dependent flow cells. Streamflow recharge was represented with specified-flow cells, and areal recharge to layers 1-4 was represented with a specified-flux boundary. Calibration of the model was accomplished by two simulations: (1) steady-state simulation of average conditions for water years 1988-97 and (2) transient simulations of water years 1988-97 divided into twenty 6-month stress periods. Flow-system components represented in the model include recharge, discharge, and hydraulic properties. The steady-state streamflow recharge rate was 42.2 cubic feet per second (ft3/s), and transient streamflow recharge rates ranged from 14.1 to 102.2 ft3/s. The steady-state areal recharge rate was 20.9 ft3/s, and transient areal recharge rates ranged from 1.1 to 98.4 ft3/s. The upward flow rate from the Deadwood aquifer to the Madison hydrogeologic unit was 6.3 ft3/s. Discharge included springflow, water use, flow to overlying units, and regional outflow. The estimated steady-state springflow of 32.8 ft3/s from seven springs was similar to the simulated springflow of 31.6 ft3/s, which included 20.5 ft3

South Dakota↗
USGS2015

Comparisons of estimates of annual exceedance-probability discharges for small drainage basins in Iowa, based on data through water year 2013

Traditionally, the Iowa Department of Transportation has used the Iowa Runoff Chart and single-variable regional-regression equations (RREs) from a U.S. Geological Survey report (published in 1987) as the primary methods to estimate annual exceedance-probability discharge (AEPD) for small (20 square miles or less) drainage basins in Iowa. With the publication of new multi- and single-variable RREs by the U.S. Geological Survey (published in 2013), the Iowa Department of Transportation needs to determine which methods of AEPD estimation provide the best accuracy and the least bias for small drainage basins in Iowa. Twenty five streamgages with drainage areas less than 2 square miles (mi 2 ) and 55 streamgages with drainage areas between 2 and 20 mi 2 were selected for the comparisons that used two evaluation metrics. Estimates of AEPDs calculated for the streamgages using the expected moments algorithm/multiple Grubbs-Beck test analysis method were compared to estimates of AEPDs calculated from the 2013 multivariable RREs; the 2013 single-variable RREs; the 1987 single-variable RREs; the TR-55 rainfall-runoff model; and the Iowa Runoff Chart. For the 25 streamgages with drainage areas less than 2 mi 2 , results of the comparisons seem to indicate the best overall accuracy and the least bias may be achieved by using the TR-55 method for flood regions 1 and 3 (published in 2013) and by using the 1987 single-variable RREs for flood region 2 (published in 2013). For drainage basins with areas between 2 and 20 mi 2 , results of the comparisons seem to indicate the best overall accuracy and the least bias may be achieved by using the 1987 single-variable RREs for the Southern Iowa Drift Plain landform region and for flood region 3 (published in 2013), by using the 2013 multivariable RREs for the Iowan Surface landform region, and by using the 2013 or 1987 single-variable RREs for flood region 2 (published in 2013). For all other landform or flood regions in Iowa, use of the 2013 single-variable RREs may provide the best overall accuracy and the least bias. An examination was conducted to understand why the 1987 single-variable RREs seem to provide better accuracy and less bias than either of the 2013 multi- or single-variable RREs. A comparison of 1-percent annual exceedance-probability regression lines for hydrologic regions 1-4 from the 1987 single-variable RREs and for flood regions 1-3 from the 2013 single-variable RREs indicates that the 1987 single-variable regional-regression lines generally have steeper slopes and lower discharges when compared to 2013 single-variable regional-regression lines for corresponding areas of Iowa. The combination of the definition of hydrologic regions, the lower discharges, and the steeper slopes of regression lines associated with the 1987 single-variable RREs seem to provide better accuracy and less bias when compared to the 2013 multi- or single-variable RREs; better accuracy and less bias was determined particularly for drainage areas less than 2 mi 2 , and also for some drainage areas between 2 and 20 mi 2 . The 2013 multi- and single-variable RREs are considered to provide better accuracy and less bias for larger drainage areas. Results of this study indicate that additional research is needed to address the curvilinear relation between drainage area and AEPDs for areas of Iowa.

Iowa↗
USGS2009

Geography of Alaska lake districts: Identification, description, and analysis of lake-rich regions of a diverse and dynamic state

Lakes are abundant landforms and important ecosystems in Alaska, but are unevenly distributed on the landscape with expansive lake-poor regions and several lake-rich regions. Such lake-rich areas are termed lake districts and have landscape characteristics that can be considered distinctive in similar respects to mountain ranges. In this report, we explore the nature of lake-rich areas by quantitatively identifying Alaska’s lake districts, describing and comparing their physical characteristics, and analyzing how Alaska lake districts are naturally organized and correspond to climatic and geophysical characteristics, as well as studied and managed by people. We use a digital dataset (National Hydrography Dataset) of lakes greater than 1 hectare, which includes 409,040 individual lakes and represents 3.3 percent of the land-surface area of Alaska. The selection criteria we used to identify lake districts were (1) a lake area (termed limnetic ratio, in percent) greater than the mean for the State, and (2) a lake density (number of lakes per unit area) greater than the mean for the State using a pixel size scaled to the area of interest and number of lakes in the census. Pixels meeting these criteria were grouped and delineated and all groups greater than 1,000 square kilometers were identified as Alaska’s lake districts. These lake districts were described according to lake size-frequency metrics, elevation distributions, geology, climate, and ecoregions to better understand their similarities and differences. We also looked at where lake research and relevant ecological monitoring has occurred in Alaska relative to lake districts and how lake district lands and waters are currently managed. We identified and delineated 20 lake districts in Alaska representing 16 percent of the State, but including 65 percent of lakes and 75 percent of lake area. The largest lake districts identified are the Yukon-Kuskokwim Delta, Arctic Coastal Plain, and Iliamna lake districts with high limnetic ratios of 19, 17, and 21 percent, respectively. The three smallest districts we considered were Tetlin in the eastern interior, Menhiskof on the Alaska Peninsula, and Matanuska–Susitna at the head of Cook Inlet with limnetic ratios of 14, 9, and 9 percent, respectively. Lake density and limnetic ratio were poorly related among lake districts, such that some districts had a few large lakes like Iliamna with Lakes Iliamna and Becharof—the two largest in the State, compared to other districts with many very small lakes like Yukon-Kuskokwim Delta with 111,130 lakes and 63 percent of these less than 10 hectares. Most lake districts are in regions with relatively low precipitation, but temperature regimes varied widely among lake districts. Approximately one-half of lake districts were glaciated during the Pleistocene and similar numbers occur in regions classified as having continuous, discontinuous, and sporadic permafrost, or perennially unfrozen soils. Most districts are at low elevations (less than 250 meters) with two important exceptions being Tetlin with a mean elevation of 530 meters and Ahtna with a mean elevation of 760 meters. These higher elevation districts, particularly Ahtna, had distinct characteristics from other lake districts such as continuous permafrost and Pleistocene glaciation. Several lake districts share similar boundaries to defined ecoregions with lake districts occurring in less than one-half of these 32 ecoregions of Alaska. Most lake districts are lands fully or partly managed by the U.S. Fish and Wildlife Service and the National Park Service, with other land management by the Bureau of Land Management and State and borough government. Much of the U.S. Geological Survey’s lake water-quality sampling efforts has been done in the Arctic Coastal Plain, Matanuska-Susitna, and Iliamna districts but no recorded collections in nine lake districts. Similarly, most lake limnological studies in Alaska were site-specific and represent only a small portion of Alaska’s lake districts. This identification, characterization, and analysis of lake-rich regions may help provide a template to guide future limnological and other scientific research for Alaska.

Alaska↗
USGS2015

Water-budgets and recharge-area simulations for the Spring Creek and Nittany Creek Basins and parts of the Spruce Creek Basin, Centre and Huntingdon Counties, Pennsylvania, Water Years 2000–06

This report describes the results of a study by the U.S. Geological Survey in cooperation with ClearWater Conservancy and the Pennsylvania Department of Environmental Protection to develop a hydrologic model to simulate a water budget and identify areas of greater than average recharge for the Spring Creek Basin in central Pennsylvania. The model was developed to help policy makers, natural resource managers, and the public better understand and manage the water resources in the region. The Groundwater and Surface-water FLOW model (GSFLOW), which is an integration of the Precipitation-Runoff Modeling System (PRMS) and the Modular Groundwater Flow Model (MODFLOW-NWT), was used to simulate surface water and groundwater in the Spring Creek Basin for water years 2000–06. Because the groundwater and surface-water divides for the Spring Creek Basin do not coincide, the study area includes the Nittany Creek Basin and headwaters of the Spruce Creek Basin. The hydrologic model was developed by the use of a stepwise process: (1) develop and calibrate a PRMS model and steady-state MODFLOW-NWT model; (2) re-calibrate the steady-state MODFLOW-NWT model using potential recharge estimates simulated from the PRMS model, and (3) integrate the PRMS and MODFLOW-NWT models into GSFLOW. The individually calibrated PRMS and MODFLOW-NWT models were used as a starting point for the calibration of the fully coupled GSFLOW model. The GSFLOW model calibration was done by comparing observations and corresponding simulated values of streamflow from 11 streamgages and groundwater levels from 16 wells. The cumulative water budget and individual water budgets for water years 2000–06 were simulated by using GSFLOW. The largest source and sink terms are represented by precipitation and evapotranspiration, respectively. For the period simulated, a net surplus in the water budget was computed where inflows exceeded outflows by about 1.7 billion cubic feet (0.47 inches per year over the basin area); storage increased by about the same amount to balance the budget. The rate and distribution of recharge throughout the Spring Creek, Nittany Creek, and Spruce Creek Basins is variable as a result of the high degree of hydrogeologic heterogeneity and karst features. The greatest amount of recharge was simulated in the carbonate-bedrock valley, near the toe slopes of Nittany and Tussey Mountains, in the Scotia Barrens, and along the area coinciding with the Gatesburg Formation. Runoff extremes were observed for water years 2001 (dry year) and 2004 (wet year). Simulated average recharge rates (water reaching the saturated zone as defined in GSFLOW) for 2001 and 2004 were 5.4 in/yr and 22.0 in/yr, respectively. Areas where simulations show large variations in annual recharge between wet and dry years are the same areas where simulated recharge was large. Those areas where rates of groundwater recharge are much higher than average, and are capable of accepting substantially greater quantities of recharge during wet years, might be considered critical for maintaining the flow of springs, stream base flow, or the source of water to supply wells. The slopes of the Bald Eagle, Tussey, and Nittany Mountains are relatively insensitive to variations in recharge, primarily because of reduced infiltration rates and steep slopes.

Pennsylvania↗
USGS2015

Data Collection and Simulation of Ecological Habitat and Recreational Habitat in the Shenandoah River, Virginia

This report presents updates to methods, describes additional data collected, documents modeling results, and discusses implications from an updated habitat-flow model that can be used to predict ecological habitat for fish and recreational habitat for canoeing on the main stem Shenandoah River in Virginia. Given a 76-percent increase in population predictions for 2040 over 1995 records, increased water-withdrawal scenarios were evaluated to determine the effects on habitat and recreation in the Shenandoah River. Projected water demands for 2040 vary by watershed: the North Fork Shenandoah River shows a 55.9-percent increase, the South Fork Shenandoah River shows a 46.5-percent increase, and the main stem Shenandoah River shows a 52-percent increase; most localities are projected to approach the total permitted surface-water and groundwater withdrawals values by 2040, and a few localities are projected to exceed these values. The habitat model used for this study evaluates the suitability of ecological habitat, represented by fish, and recreational habitat, represented by canoeing, based on depth, velocity, and substrate conditions, which are weighted for the physical habitat types (riffles, runs, or pools) present within a stretch of river. Weighted usable-habitat area in the Lockes Mill reach was maximized for adult smallmouth bass and sub-adult smallmouth bass ( Micropterus dolomieu ) and river chub ( Nocomis micropogon ) when streamflows were equal to median flow (900 cubic feet per second) for summer months. Ecological maximum weighted usable-habitat areas for smaller fish, such as spotfin or satinfin shiner ( Cyprinella spp.), margined madtom ( Noturus insignis ), and juvenile redbreast sunfish ( Lepomis auritus ) occurred with 10 th percentile flows (482 cubic feet per second) and lower. Recreational weighted usable-habitat areas for canoeing were maximized when streamflows were above the 75 th percentile (1,410 cubic feet per second). During historic droughts, streamflows were less than the 10 th percentile, and adult smallmouth bass and sub-adult smallmouth bass habitat was below normal for the majority of days during at least 2 months of the summer. When streamflows were less than the lowest 7-day average in a 10-year period, or 7Q10 flow (357 cubic feet per second), margined madtom, river chub, and sub-adult redbreast sunfish habitat areas were below normal as well. Streamflows that limit most fish species habit availability range from 300 to 500 cubic feet per second. For the drought years simulated, flows that were equal to or less than the 10 th percentile for summer months did not provide adequate depth for canoe passage through riffle habitats. A modeling limitation for higher flows than those studied during development of the habitat-suitability criteria is that modeled habitat availability will decrease as flows increase. Time-series analyses were used to investigate changes in habitat availability with increased water withdrawals of 10, 20, and almost 50 percent (48.6 percent) up to the 2040 amounts projected by local water supply plans. Adult and sub-adult smallmouth bass frequently had habitat availability outside the normal range for habitat conditions during drought years, yet 10- or 20-percent increases in withdrawals did not contribute to a large reduction in habitat. When withdrawals were increased by 50 percent, there was an additional decrease in habitat. During 2002 drought scenarios, reduced habitat availability for sub-adult redbreast sunfish or river chub was only slightly evident with 50-percent increased withdrawal scenarios. Recreational habitat represented by canoeing decreased lower than normal during the 2002 drought. For a recent normal year, like 2012, increased water-withdrawal scenarios did not affect habitat availability for fish such as adult and sub-adult smallmouth bass, sub-adult redbreast sunfish, or river chub. Canoeing habitat availability was within the normal range most of 2012, and increased water-withdrawal scenarios showed almost no affect. For both ecological fish habitat and recreational canoeing habitat, the antecedent conditions (habitat within normal range of habitat or below normal) appear to govern whether additional water withdrawals will affect habitat availability. As human populations and water demands increase, many of the ecological or recreational stresses may be lessened by managing the timing of water withdrawals from the system.

Virginia↗
USGS2005

Historical perspective of statewide streamflows during the 2002 and 1977 droughts in Colorado

Since 1890, Colorado has experienced a number of widespread drought periods; the most recent statewide drought began during 1999 and includes 2002, a year characterized by precipitation, snowpack accumulation, and streamflows that were much lower than normal. Because the drought of 2002 had a substantial effect on streamflows in Colorado, the U.S. Geological Survey, in cooperation with the Colorado Water Conservation Board, began a study in 2004 to analyze statewide streamflows during 2002 and develop a historical perspective of those streamflows. The purpose of this report is to describe an analysis of streamflows recorded throughout Colorado during the drought of 2002, as well as other drought years such as 1977, and to provide some historical perspective of drought-diminished streamflows in Colorado. Because most streamflows in Colorado are derived from melting of mountain snowpacks during April through July, streamflows primarily were analyzed for the snowmelt (high-flow) period, but streamflows also were analyzed for the winter (low-flow) period. The snowmelt period is defined as April 1 through September 30 and the winter period is defined as October 1 through March 31. Historical daily average streamflows were analyzed on the basis of 7, 30, 90, and 180 consecutive-day periods (N-day) for 154 selected stations in Colorado. Methods used for analysis of the N-day snowmelt and winter streamflows include evaluation of trends in the historical streamflow records, computation of the rank of each annual N-day streamflow value for each station, analysis for years other than 2002 and 1977 with drought-diminished streamflows, and frequency analysis (on the basis of nonexceedance probability) of the 180-day streamflows. Ranking analyses for the N-day snowmelt streamflows indicated that streamflows during 2002 were ranked as the lowest or second lowest historical values at 114-123 stations, or about 74-80 percent of the stations; by comparison, the N-day snowmelt streamflows during 1977 were ranked as the lowest or second lowest historical values at 69-87 stations, or about 47-59 percent of the stations. Many of the stations in the mountainous headwaters where snowmelt streamflows were ranked lowest during 2002 were ranked second lowest during 1977. These results indicate that snowmelt streamflows during 2002 were considerably more diminished than those during 1977. The 180-day snowmelt streamflows were ranked among the five lowest historical values at about 90 percent of the stations during 2002 and were ranked among the five lowest historical values at about 77 percent of the stations during 1977. Other years during which the 180-day snowmelt streamflows were ranked among the five lowest values at a substantial percentage of stations include 1934, 1954, 1963, and 1981, but the percentages of stations with 180-day snowmelt streamflows ranked among the five lowest values were smaller during those years than during 2002 and 1977. Frequency analysis of snowmelt streamflows indicated that recurrence intervals for the 180-day snowmelt streamflows during 2002 were greater than 50 years for about 57 percent of the stations and were more than 100 years for about 14 percent of the stations. By comparison, recurrence intervals for the 180-day snowmelt streamflows during 1977 were greater than 50 years only for about 15 percent of the stations and were more than 100 years only for about 1 percent of the stations. Generally, snowmelt streamflows during 2002 were more diminished and have higher recurrence intervals than snowmelt streamflows during 1977. The N-day winter streamflows during 2002 and 1977 were not ranked among the five lowest historical values at about 86-103 stations, or about 58-70 percent of the stations, compared to about 10-27 percent of the stations for the N-day snowmelt streamflows. These results indicate that winter streamflows during the 2002 and 1977 droughts were diminished to a lesser extent than t

Scientific Investigations Report↗
USGS2013

Characteristics of sediment transport at selected sites along the Missouri River during the high-flow conditions of 2011

During 2011, many tributaries in the Missouri River Basin experienced near record peak streamflow and caused flood damage to many communities along much of the Missouri River from Montana to the confluence with the Mississippi River. The large runoff event in 2011 provided an opportunity to examine characteristics of sediment transport in the Missouri River at high-magnitude streamflow and for a long duration. The purpose of this report is to describe sediment characteristics during the 2011 high-flow conditions at six selected sites on the Missouri River, two in the middle region of the basin between Lake Sakakawea and Lake Oahe in North Dakota, and four downstream from Gavins Point Dam along the Nebraska-South Dakota and Nebraska-Iowa borders. A wider range in suspended-sediment concentration was observed in the middle segment of the Missouri River compared to sites in the lower segment. In the middle segment of the Missouri River, suspended-sediment concentrations increased and peaked as flows increased and started to plateau; however, while flows were still high and steady, suspended-sediment concentrations decreased and suspended-sediment grain sizes coarsened, indicating the decrease possibly was related to fine-sediment supply limitations. Measured bedload transport rates in the lower segment of the Missouri River (sites 3 to 6) were consistently higher than those in the middle segment (sites 1 and 2) during the high-flow conditions in 2011. The median bedload transport rate measured at site 1 was 517 tons per day and at site 2 was 1,500 tons per day. Measured bedload transport rates were highest at site 3 then decreased downstream to site 5, then increased at site 6. The median bedload transport rates were 22,100 tons per day at site 3; 5,640 tons per day at site 4; 3,930 tons per day at site 5; and 8,450 tons per day at site 6. At the two sites in the middle segment of the Missouri River, the greatest bedload was measured during the recession of the streamflow hydrograph. A similar pattern was observed at sites 3–5 in the lower segment of the Missouri River, where the greatest bedload was measured later in the event on the recession of the streamflow hydrograph, although the change in bedload was not as dramatic as observed at the sites in the middle segment of the Missouri River. With the exception of site 3, the total-sediment load on the Missouri River was highest at the beginning of the high-flow event and decreased as streamflow decreased. In the middle segment of the Missouri River, measured total-sediment load ranged from 2,320 to 182,000 tons per day at site 1 and from 3,190 to 279,000 tons per day at site 2. In the lower segment of the Missouri River, measured total-sediment load ranged from 50,600 to 223,000 tons per day at site 4; from 23,500 to 403,000 tons per day at site 5; and from 52,700 to 273,000 tons per day at site 6. The total-sediment load was dominated by suspended sediment at all of the sites measured on the Missouri River in 2011. In general, the percentage of total-sediment load that was bedload increased as the streamflow decreased, although this pattern was more prevalent at sites in the middle segment than those in the lower segment. The suspended-sediment load comprised an average of 93 percent of the total load, with the exception of site 3, where the suspended-sediment load comprised only 72 percent of the total-sediment load.

Missouri River↗
USGS2021

Occurrence and distribution of mercury in streams and reservoirs in the Triangle Area of North Carolina, July 2007–June 2009

During the time period 2001–2006, the U.S. Geological Survey reported mercury-concentration measurements that exceeded the North Carolina water-quality criterion (NCWQC) of 0.012 microgram per liter for total recoverable mercury in streams and reservoirs across the Triangle Area of North Carolina. Mercury data were sparse, however, generally consisting of only one or two water samples per year. Additional monitoring and data analysis were needed to better determine the occurrence and distribution of mercury in the Triangle Area for all seasons and waterbody types as well as associations between mercury concentrations and water-quality and land-use parameters. Water at fifteen reservoir and 14 stream sites across the Triangle Area was sampled at various times between August 2007 and June 2009, with water samples collected from both the surfaces and bottoms of the water columns in reservoirs and from the surfaces of streams. A bed sediment sample was also collected at all reservoir sites and at all but one stream site. A total of 301 water samples was collected at reservoir sites. Filtered and total recoverable mercury were detected in at least one water sample collected from each reservoir site. A total of 77 water samples was collected from stream sites with filtered mercury detected in samples from one-half of these sites, and total recoverable mercury detected in at least one water sample from all but two sites. Total recoverable and filtered mercury concentrations exceeded the NCWQC for mercury more frequently in reservoir than in stream samples. Differences in sampling frequencies among seasons and between streams and reservoirs, however, may have negatively biased overall estimates of mercury concentrations in streams relative to reservoirs. Filtered mercury concentrations in surface-water samples from reservoirs and total recoverable mercury concentrations in bottom samples from reservoirs were highest in the fall, whereas no seasonal trends in filtered or total recoverable mercury were detected from stream samples. Total mercury concentrations were calculated for the bulk sample on the basis of the percentage of the grains in the bulk sample whose diameters that were smaller than 0.0625 millimeters. Total mercury concentrations in bed sediment were generally higher for samples from reservoir sites compared to streams sites, although the highest total mercury concentration in bed sediment was from a stream site. Concentrations of total recoverable mercury in water samples from stream sites all fell within the general range for streams and lakes without on-site significant anthropogenic sources (for example, mercury mines or industrial pollution), whereas samples collected from eight reservoir sites had total mercury concentrations in a range characteristic of sites affected by mercury mines or industrial pollution. Results suggested that litterfall may be a source of mercury in streams, whereas atmospheric deposition is likely a dominant source for reservoirs; however, high concentrations of filtered and total recoverable mercury concentrations in the fall season in some reservoir-water samples may warrant further analysis of potential hydrologic factors. Mercury concentrations in all water and bed sediment samples were below levels expected to cause adverse effects to humans and aquatic biota, indicating that mercury levels at the study sites in the Triangle Area were unlikely to cause an immediate health risk to humans or aquatic organisms. The high variability among several sample replicates for total recoverable mercury, however, indicated that inferences from total recoverable mercury concentrations can be tenuous.

North Carolina↗