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

Water budgets and groundwater volumes for abandoned underground mines in the Western Middle Anthracite Coalfield, Schuylkill, Columbia, and Northumberland Counties, Pennsylvania: Preliminary estimates with identification of data needs

This report, prepared in cooperation with the Pennsylvania Department of Environmental Protection (PaDEP), the Eastern Pennsylvania Coalition for Abandoned Mine Reclamation, and the Dauphin County Conservation District, provides estimates of water budgets and groundwater volumes stored in abandoned underground mines in the Western Middle Anthracite Coalfield, which encompasses an area of 120 square miles in eastern Pennsylvania. The estimates are based on preliminary simulations using a groundwater-flow model and an associated geographic information system that integrates data on the mining features, hydrogeology, and streamflow in the study area. The Mahanoy and Shamokin Creek Basins were the focus of the study because these basins exhibit extensive hydrologic effects and water-quality degradation from the abandoned mines in their headwaters in the Western Middle Anthracite Coalfield. Proposed groundwater withdrawals from the flooded parts of the mines and stream-channel modifications in selected areas have the potential for altering the distribution of groundwater and the interaction between the groundwater and streams in the area. Preliminary three-dimensional, steady-state simulations of groundwater flow by the use of MODFLOW are presented to summarize information on the exchange of groundwater among adjacent mines and to help guide the management of ongoing data collection, reclamation activities, and water-use planning. The conceptual model includes high-permeability mine voids that are connected vertically and horizontally within multicolliery units (MCUs). MCUs were identified on the basis of mine maps, locations of mine discharges, and groundwater levels in the mines measured by PaDEP. The locations and integrity of mine barriers were determined from mine maps and groundwater levels. The permeability of intact barriers is low, reflecting the hydraulic characteristics of unmined host rock and coal. A steady-state model was calibrated to measured groundwater levels and stream base flow, the latter at many locations composed primarily of discharge from mines. Automatic parameter estimation used MODFLOW-2000 with manual adjustments to constrain parameter values to realistic ranges. The calibrated model supports the conceptual model of high-permeability MCUs separated by low-permeability barriers and streamflow losses and gains associated with mine infiltration and discharge. The simulated groundwater levels illustrate low groundwater gradients within an MCU and abrupt changes in water levels between MCUs. The preliminary model results indicate that the primary result of increased pumping from the mine would be reduced discharge from the mine to streams near the pumping wells. The intact barriers limit the spatial extent of mine dewatering. Considering the simulated groundwater levels, depth of mining, and assumed bulk porosity of 11 or 40 percent for the mined seams, the water volume in storage in the mines of the Western Middle Anthracite Coalfield was estimated to range from 60 to 220 billion gallons, respectively. Details of the groundwater-level distribution and the rates of some mine discharges are not simulated well using the preliminary model. Use of the model results should be limited to evaluation of the conceptual model and its simulation using porous-media flow methods, overall water budgets for the Western Middle Anthracite Coalfield, and approximate storage volumes. Model results should not be considered accurate for detailed simulation of flow within a single MCU or individual flooded mine. Although improvements in the model calibration were possible by introducing spatial variability in permeability parameters and adjusting barrier properties, more detailed parameterizations have increased uncertainty because of the limited data set. The preliminary identification of data needs includes continuous streamflow, mine discharge rate, and groundwater levels in the mines and adjacent areas. Data collected whe

Pennsylvania↗

Hydrogeologic framework and occurrence, movement, and chemical characterization of groundwater in Dixie Valley, west-central Nevada

Dixie Valley, a primarily undeveloped basin in west-central Nevada, is being considered for groundwater exportation. Proposed pumping would occur from the basin-fill aquifer. In response to proposed exportation, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation and Churchill County, conducted a study to improve the understanding of groundwater resources in Dixie Valley. The objective of this report is to characterize the hydrogeologic framework, the occurrence and movement of groundwater, the general water quality of the basin-fill aquifer, and the potential mixing between basin-fill and geothermal aquifers in Dixie Valley. Various types of geologic, hydrologic, and geochemical data were compiled from previous studies and collected in support of this study. Hydrogeologic units in Dixie Valley were defined to characterize rocks and sediments with similar lithologies and hydraulic properties influencing groundwater flow. Hydraulic properties of the basin-fill deposits were characterized by transmissivity estimated from aquifer tests and specific-capacity tests. Groundwater-level measurements and hydrogeologic-unit data were combined to create a potentiometric surface map and to characterize groundwater occurrence and movement. Subsurface inflow from adjacent valleys into Dixie Valley through the basin-fill aquifer was evaluated using hydraulic gradients and Darcy flux computations. The chemical signature and groundwater quality of the Dixie Valley basin-fill aquifer, and potential mixing between basin-fill and geothermal aquifers, were evaluated using chemical data collected from wells and springs during the current study and from previous investigations. Dixie Valley is the terminus of the Dixie Valley flow system, which includes Pleasant, Jersey, Fairview, Stingaree, Cowkick, and Eastgate Valleys. The freshwater aquifer in the study area is composed of unconsolidated basin-fill deposits of Quaternary age. The basin-fill hydrogeologic unit can be several orders of magnitude more transmissive than surrounding and underlying consolidated rocks and Dixie Valley playa deposits. Transmissivity estimates in the basin fill throughout Dixie Valley ranged from 30 to 45,500 feet squared per day; however, a single transmissivity value of 0.1 foot squared per day was estimated for playa deposits. Groundwater generally flows from the mountain range uplands toward the central valley lowlands and eventually discharges near the playa edge. Potentiometric contours east and west of the playa indicate that groundwater is moving eastward from the Stillwater Range and westward from the Clan Alpine Mountains toward the playa. Similarly, groundwater flows from the southern and northern basin boundaries toward the basin center. Subsurface groundwater flow likely enters Dixie Valley from Fairview and Stingaree Valleys in the south and from Jersey and Pleasant Valleys in the north, but groundwater connections through basin-fill deposits were present only across the Fairview and Jersey Valley divides. Annual subsurface inflow from Fairview and Jersey Valleys ranges from 700 to 1,300 acre-feet per year and from 1,800 to 2,300 acre-feet per year, respectively. Groundwater flow between Dixie, Stingaree, and Pleasant Valleys could occur through less transmissive consolidated rocks, but only flow through basin fill was estimated in this study. Groundwater in the playa is distinct from the freshwater, basin-fill aquifer. Groundwater mixing between basin-fill and playa groundwater systems is physically limited by transmissivity contrasts of about four orders of magnitude. Total dissolved solids in playa deposit groundwater are nearly 440 times greater than total dissolved solids in the basin-fill groundwater. These distinctive physical and chemical flow restrictions indicate that groundwater interaction between the basin fill and playa sediments was minimal during this study period (water years 2009–11). Groundwater in Dixie Valley generally can be characterized as a sodium bicarbonate type, with greater proportions of chloride north of the Dixie Valley playa, and greater proportions of sulfate south of the playa. Analysis of major ion water chemistry data sampled during the study period indicates that groundwater north and south of Township 22N differ chemically. Dixie Valley groundwater quality is marginal when compared with national primary and secondary drinking-water standards. Arsenic and fluoride concentrations exceed primary drinking water standards, and total dissolved solids and manganese concentrations exceed secondary drinking water standards in samples collected during this study. High concentrations of boron and tungsten also were observed. Chemical comparisons between basin-fill and geothermal aquifer water indicate that most basin-fill groundwater sampled could contain 10–20 percent geothermal water. Geothermal indicators such as high temperature, lithium, boron, chloride, and silica suggest that mixing occurs in many wells that tap the basin-fill aquifer, particularly on the north, south, and west sides of the basin. Magnesium-lithium geothermometers indicate that some basin-fill aquifer water sampled for the current study likely originates from water that was heated above background mountain-block recharge temperatures (between 3 and 15 degrees Celsius), highlighting the influence of mixing with warm water that was possibly derived from geothermal sources.

Nevada↗

Comparison of Hydrologic and Water-Quality Characteristics of Two Native Tallgrass Prairie Streams with Agricultural Streams in Missouri and Kansas

This report presents the results of a study by the U.S. Geological Survey, in cooperation with the Missouri Department of Natural Resources, to analyze and compare hydrologic and water-quality characteristics of tallgrass prairie and agricultural basins located within the historical distribution of tallgrass prairie in Missouri and Kansas. Streamflow and water-quality data from two remnant, tallgrass prairie basins (East Drywood Creek at Prairie State Park, Missouri, and Kings Creek near Manhattan, Kansas) were compared to similar data from agricultural basins in Missouri and Kansas. Prairie streams, especially Kings Creek in eastern Kansas, received a higher percentage of base flow and a lower percentage of direct runoff than similar-sized agricultural streams in the region. A larger contribution of direct runoff from the agricultural streams made them much flashier than prairie streams. During 22 years of record, the Kings Creek base-flow component averaged 66 percent of total flow, but base flow was only 16 to 26 percent of flows at agricultural sites of various record periods. The large base-flow component likely is the result of greater infiltration of precipitation in prairie soils and the resulting greater contribution of groundwater to streamflow. The 1- and 3-day annual maximum flows were significantly greater at three agricultural sites than at Kings Creek. The effects of flashier agricultural streams on native aquatic biota are unknown, but may be an important factor in the sustainability of some native aquatic species. There were no significant differences in the distribution of dissolved-oxygen concentrations at prairie and agricultural sites, and some samples from most sites fell below the 5 milligrams per liter Missouri and Kansas standard for the protection of aquatic life. More than 10 percent of samples from the East Drywood Creek prairie stream were less than this standard. These data indicate low dissolved-oxygen concentrations during summer low-flow periods may be a natural phenomenon for small prairie streams in the Osage Plains. Nutrient concentrations including total nitrogen, ammonia, nitrate, and total phosphorus were significantly less in base-flow and runoff samples from prairie streams than from agricultural streams. The total nitrogen concentration at all sites other than one of two prairie sampling sites were, on occasion, above the U.S. Environmental Protection Agency recommended criterion for total nitrogen for the prevention of nutrient enrichment, and typically were above this recommended criterion in runoff samples at all sites. Nitrate and total phosphorus concentrations in samples from the prairie streams generally were below the U.S. Environmental Protection Agency recommended nutrient criteria in base-flow and runoff samples, whereas samples from agricultural sites generally were below the criteria in base-flow samples and generally above in runoff samples. The lower concentrations of nutrient species in prairie streams is likely because prairies are not fertilized like agricultural basins and prairie basins are able to retain nutrients better than agricultural basins. This retention is enhanced by increased infiltration of precipitation into the prairie soils, decreased surface runoff, and likely less erosion than in agricultural basins. Streamflow in the small native prairie streams had more days of zero flow and lower streamflow yields than similar-sized agricultural streams. The prairie streams were at zero flow about 50 percent of the time, and the agricultural streams were at zero flow 25 to 35 percent of the time. Characteristics of the prairie basins that could account for the greater periods of zero flow and lower yields when compared to agricultural streams include greater infiltration, greater interception and evapotranspiration, shallower soils, and possible greater seepage losses in the prairie basins. Another difference between the prairie and agricultural strea

Scientific Investigations Report↗

Assessing potential effects of highway runoff on receiving-water quality at selected sites in Oregon with the Stochastic Empirical Loading and Dilution Model (SELDM)

In 2012, the U.S. Geological Survey and the Oregon Department of Transportation began a cooperative study to demonstrate use of the Stochastic Empirical Loading and Dilution Model (SELDM) for runoff-quality analyses in Oregon. SELDM can be used to estimate stormflows, constituent concentrations, and loads from the area upstream of a stormflow discharge site, from the site of interest and in the receiving waters downstream of the discharge. SELDM also can be used to assess the potential effectiveness of best management practices (BMP) for mitigating potential effects of runoff in receiving waters. Nominally, SELDM is a highway-runoff model, but it is well suited for analysis of runoff from other land uses as well. This report provides case studies and examples to demonstrate stochastic-runoff modeling concepts and to demonstrate application of the model. Basin characteristics from six Oregon highway study sites were used to demonstrate various applications of the model. The highway catchment and upstream basin drainage areas of these study sites ranged from 3.85 to 11.83 acres and from 0.16 to 6.56 square miles, respectively. The upstream basins of two sites are urbanized, and the remaining four sites are less than 5 percent impervious. SELDM facilitates analysis by providing precipitation, pre-storm streamflow, and other variables by region or from hydrologically similar sites. In Oregon, there can be large variations in precipitation and streamflow among nearby sites. Therefore, spatially interpolated geographic information system data layers containing storm-event precipitation and pre-storm streamflow statistics specific to Oregon were created for the study using Kriging techniques. Concentrations and loads of cadmium, chloride, chromium, copper, iron, lead, nickel, phosphorus, and zinc were simulated at the six Oregon highway study sites by using statistics from sites in other areas of the country. Water‑quality datasets measured at hydrologically similar basins in the vicinity of the study sites in Oregon were selected and compiled to estimate stormflow-quality statistics for the upstream basins. The quality of highway runoff and some upstream stormflow constituents were simulated by using statistical moments (average, standard deviation, and skew) of the logarithms of data. Some upstream stormflow constituents were simulated by using transport curves, which are relations between stormflow and constituent concentrations. Stochastic analyses were done by using SELDM to demonstrate use of the model and to illustrate the types of information that stochastic analyses may provide: 1. An analysis was done to demonstrate use of dilution factors as an initial reconnaissance tool for comparing relative risk among sites. 2. An analysis of hardness-dependent, water-quality criteria was done to illustrate the effects of variations in hardness and flow on the application and interpretation of such criteria. This analysis shows that hardness-dependent criteria can vary by an order of magnitude among storm events because hardness is diluted by stormflows. 3. An analysis of uncertainties in input and output values was done to demonstrate that properly selected robust datasets are needed to represent conditions at a site of interest. This analysis shows that the rate of water-quality exceedances that are measured or simulated may depend on sample size and the luck of the draw. 4. An analysis was done to demonstrate that SELDM and other Monte Carlo models may generate extreme values from input statistics, which may or may not be feasible based on physicochemical or hydrological limits. 5. An analysis of BMP modeling methods was done to demonstrate use of the model for estimating treatment requirements for meeting water-quality objectives. 6. An analysis of the use of grab sampling and nonstochastic upstream modeling methods was done to evaluate the potential effects on modeling outcomes. Additional analyses using surrogate water-quality datasets for the upstream basin and highway catchment were provided for six Oregon study sites to illustrate the risk-based information that SELDM will produce. These analyses show that the potential effects of highway runoff on receiving-water quality downstream of the outfall depends on the ratio of drainage areas (dilution), the quality of the receiving water upstream of the highway, and the concentration of the criteria of the constituent of interest. These analyses also show that the probability of exceeding a water-quality criterion may depend on the input statistics used, thus careful selection of representative values is important.

Oregon↗

Simulation of groundwater flow and chloride transport in the “1,500-foot” sand, “2,400-foot” sand, and “2,800-foot” sand of the Baton Rouge area, Louisiana

Groundwater withdrawals since the 1940s have lowered water levels, altered groundwater-flow directions, and caused saltwater to intrude within some freshwater-containing sands of the fluvial-deltaic Southern Hills regional aquifer system beneath Baton Rouge, Louisiana. New interpretations of stratigraphic correlations amongst geophysical well logs were utilized to revise a hydrogeologic framework that delineates the depth and thickness variations of aquifers and confining units in the Southern Hills regional aquifer system. A groundwater-flow and chloride-transport model incorporating the revised framework was constructed to assess the effects of groundwater withdrawals on the rate and pathways of saltwater migration in the “1,500-foot” sand, “2,400-foot” sand, and the “2,800-foot” sand. Groundwater withdrawals reported since 1940 were compiled to specify annual average withdrawal rates through 2016 for 722 wells. Regional groundwater flow throughout the Southern Hills regional aquifer system was first simulated with MODFLOW, and flow-model parameters were calibrated to 8,810 water levels observed through 2016 by using the parameter-estimation code PEST++. Saltwater transport was subsequently simulated for the “1,500-foot” sand, “2,400-foot” sand, and the “2,800-foot” sand by using the variable-density code, SEAWAT. Chloride-concentration measurements were used as a proxy for saltwater to formulate the concentration initial conditions and calibrate the transport-model parameters. Three groundwater-management scenarios were simulated to evaluate the effects of different groundwater withdrawals on future groundwater levels and saltwater concentrations in the “1,500-foot” sand, “2,400-foot” sand, and “2,800-foot” sand. All three scenarios simulated the period from 2017 through 2112 (96 years), and the water levels and concentrations simulated for 2047 and 2112 were compared among the scenarios. The first scenario simulated a continuation of groundwater withdrawals at 2016 rates and represents the “status quo” of groundwater withdrawals. The second scenario simulated the effects of discontinuing 10,620 gallons per minute (gal/min) of withdrawals from the “2,800-foot” sand, and the third scenario simulated reallocating 2,000 gal/min of withdrawals from the “1,500-foot” sand to the “2,800-foot” sand. Continuation of the “status quo” withdrawals results in lower water levels by 2047 around groundwater-withdrawal locations in the “1,500-foot” sand, “2,400-foot” sand, and “2,800-foot” sand. By 2112, water levels recover to higher levels as flow in the aquifer approaches equilibrium. Saltwater within the “1,500-foot” sand would continue migrating toward public-supply wells located 2.4 miles (mi) north of the Baton Rouge Fault, but a “scavenger well” that removes relatively concentrated water from the base of the “1,500-foot” sand attenuates chloride concentrations at the public-supply wells. Saltwater within the “2,400-foot” sand would continue to encroach on a well with large withdrawals and farther east within an area about 1 mi north of the Baton Rouge Fault. Saltwater within the “2,800-foot” sand would migrate northward toward withdrawal wells located about 3 mi north of the industrial district. Cessation of 10,620 gal/min of industrial withdrawals from the “2,800-foot” sand about 12 mi northwest of the industrial district (scenario 2) would cause a substantial water-level recovery in the “2,800-foot” sand in the area of discontinued withdrawals. Groundwater levels 3 mi north of the industrial district would be 25–30 feet higher in 2047 than predicted for the “status quo” withdrawals. Saltwater encroachment toward wells north of the industrial district would be slowed because of the decreased hydraulic gradient. Reallocating 2,000 gal/min of withdrawals from the “1,500-foot” sand to the “2,800-foot” sand 12 mi northwest of the industrial district (scenario 3) would have a negligible effect on water levels and chloride concentrations in the “1,500-foot” sand 15 mi to the south-southeast where saltwater is encroaching toward wells in the “1,500-foot” sand. Within the “2,800-foot” sand, the area of saltwater encroachment is only 3 mi from increased withdrawals in the “2,800-foot” sand, and water levels would be about 5 feet lower in 2047 than for the “status quo” scenario. A larger hydraulic gradient would cause slightly faster saltwater transport and higher chloride concentrations within this area of the “2,800-foot” sand.

Louisiana↗

Net Acid Production, Acid Neutralizing Capacity, and Associated Mineralogical and Geochemical Characteristics of Animas River Watershed Igneous Rocks Near Silverton, Colorado

This report presents results from laboratory and field studies involving the net acid production (NAP), acid neutralizing capacity (ANC), and magnetic mineralogy of 27 samples collected in altered volcanic terrain in the upper Animas River watershed near Silverton, Colo., during the summer of 2005. Sampling focused mainly on the volumetrically important, Tertiary-age volcanic and plutonic rocks that host base- and precious-metal mineralization in the study area. These rocks were analyzed to determine their potential for neutralization of acid-rock drainage. Rocks in the study area have been subjected to a regional propylitic alteration event, which introduced calcite, chlorite (clinochlore), and epidote that have varying amounts and rates of acid neutralizing capacity (ANC). Locally, hydrothermal alteration has consumed any ANC and introduced minerals, mainly pyrite, that have a high net acid production (NAP). Laboratory studies included hydrogen pyroxide (H2O2) acid digestion and subsequent sodium hydroxide (NaOH) titration to determine NAP, and sulfuric acid (H2SO4) acid titration experiments to determine ANC. In addition to these environmental rock-property determinations, mineralogical, chemical, and petrographic characteristics of each sample were determined through semiquantitative X-ray diffractometry (Rietveld method), optical mineralogy, wavelength dispersive X-ray fluorescence, total carbon-carbonate, and inductively coupled plasma?mass spectrometric analysis. An ANC ranking was assigned to rock samples based on calculated ANC quantity in kilograms/ton (kg/t) calcium carbonate equivalent and ratios of ANC to NAP. Results show that talus near the southeast Silverton caldera margin, composed of andesite clasts of the Burns Member of the Silverton Volcanics, has the highest ANC (>100 kg/t calcium carbonate equivalent) with little to no NAP. The other units found to have moderate to high ANC include (a) andesite lavas and volcaniclastic rocks of the San Juan Formation, west and northwest of the Silverton caldera, and (b) the Picayune Megabreccia Member of Sapinero Mesa Tuff along the western San Juan caldera margin. Sultan Mountain stock, composed of granitoid intrusive rocks, was shown to have low ANC and moderate NAP. Sequential leachate analyses on a suite of whole-rock samples from the current and a previous study indicate that host rock composition and mineralogy control leachate compositions. The most mafic volcanic samples had high leachate concentrations for Mg, Fe, and Ca, whereas silicic volcanic samples had lower ferromagnesiun compositions. Samples with high chlorite abundance also had high leachable Mg concentrations. Trace-element substitution, such as Sr for Ca in plagioclase, controls high Sr concentrations in those samples with high plagioclase abundance. High Ti abundance in leachate was observed in those samples with high magnetite concentrations. This is likely due to samples containing intergrown magnetite-ilmenite. Whole rocks having high trace-element concentrations have relatively high leachate trace-element abundances. Some lavas of the San Juan Formation and Burns Member of the Silverton Volcanics had elevated Zn-, Cd-, and Pb-leachate concentrations. Manganese was also elevated in one San Juan Formation sample. Other San Juan Formation and Burns Member lavas had low to moderate trace-element abundances. One sample of the pyroxene andesite member of the Silverton Volcanics had elevated concentrations for As and Mo. Most other pyroxene andesite member samples had low leachate trace-element abundances. Mine-waste-leachate analyses indicated that one mine-waste sample had elevated concentrations of Cu (1.5 orders of magnitude), Zn (1 order of magnitude), As (1 order of magnitude), Mo (1.5 to 2 orders of magnitude), Cd (1 to 2 orders of magnitude), and Pb (2 to 3 orders of magnitude) compared to whole rocks. These data indicate the importance of whole-rock geochemistry or leachate analys

Scientific Investigations Report↗

Effects of removing Good Hope Mill Dam on selected physical, chemical, and biological characteristics of Conodoguinet Creek, Cumberland County, Pennsylvania

The implications of dam removal on channel characteris-tics, water quality, benthic invertebrates, and fish are not well understood because of the small number of removals that have been studied. Comprehensive studies that document the effects of dam removal are just beginning to be published, but most research has focused on larger dams or on the response of a sin-gle variable (such as benthic invertebrates). This report, pre-pared in cooperation with the Conodoguinet Creek Watershed Association, provides an evaluation of how channel morphol-ogy, bed-particle-size distribution, water quality, benthic inver-tebrates, fish, and aquatic habitat responded after removal of Good Hope Mill Dam (a small 'run of the river' dam) from Conodoguinet Creek in Cumberland County, Pa. Good Hope Mill Dam was a 6-foot high, 220-foot wide concrete structure demolished and removed over a 3-day period beginning with the initial breach on November 2, 2001, at 10:00 a.m. eastern standard time. To isolate the effects of dam removal, data were collected before and after dam removal at five monitoring stations and over selected reaches upstream, within, and downstream of the impoundment. Stations 1, 2, and 5 were at free-flowing control locations 4.9 miles upstream, 2.5 miles upstream, and 5 miles downstream of the dam, respec-tively. Stations 3 and 4 were located where the largest responses were anticipated, 115 feet upstream and 126 feet downstream of the dam, respectively Good Hope Mill Dam was not an effective barrier to sedi-ment transport. Less than 3 inches of sediment in the silt/clay-size range (less than 0.062 millimeters) coated bedrock within the 7,160-foot (1.4-mile) impoundment. The bedrock within the impoundment was not incised during or after dam removal, and the limited sediment supply resulted in no measurable change in the thalweg elevation downstream of the dam. The cross-sec-tional areas at stations 3 and 4, measured 17 days and 23 months after dam removal, were within 3 percent of the area measured before removal. Some of the impounded silt/clay at station 3 and other sed-iment in the work area downstream of the dam were initially entrained over the 3-day removal period and deposited on sub-strate at station 4. Remaining silt/clay at station 3 and deposits at station 4 were transported downstream by the flows mea-sured over the 23 months after removal (daily mean flow ranged from 38 to 5,180 cubic feet per second). The median bed-parti-cle size at station 3 increased by approximately 32 millimeters in the 23-month period after removal. Bed-particle-size distri-bution at station 4 became finer when silt/clay was initially deposited but coarsened as high flows flushed it downstream; median bed-particle size was 77.7 millimeters before removal compared to 31.3 millimeters 17 days after removal and 99 mil-limeters 23 months after removal. Good Hope Mill Dam had either no effect on water-quality characteristics or the effect was so small it was masked by sea-sonal and periodic variability. Measurements of daily mean temperature, dissolved-oxygen concentration, pH, and specific conductance on a short time scale (every 15 minutes) indicate the daily range of temperature was suppressed under impounded conditions and daily extremes of temperature, dis-solved-oxygen concentration, pH, and specific conductance at station 2 were out of phase by approximately 12 hours with station 3. Once the dam was removed, the pattern at station 3 shifted and converged with the pattern at station 2. The offset before removal may be related to a lag time resulting from a decrease in velocity through the impoundment. Total nitrogen and suspended-sediment concentrations increased upon the initial dam breach but were within the range of concentrations measured from March 2001 through April 2002 over varying flow conditions at station 1. Total nitrogen concentration at station 4 was 4.66 milligrams per liter upon the initial breach of the dam,

Pennsylvania↗

Groundwater quality and simulation of sources of water to wells in the Marsh Creek valley at the U.S. Geological Survey Northern Appalachian Research Laboratory, Tioga County, Pennsylvania

This report provides a November 2010 snapshot of groundwater quality and an analysis of the sources of water to wells at the U.S. Geological Survey (USGS) Northern Appalachian Research Laboratory (NARL) near Wellsboro, Pennsylvania. The laboratory, which conducts fisheries research, currently (2011) withdraws 1,000 gallons per minute of high-quality groundwater from three wells completed in the glacial sand and gravel aquifer beneath the Marsh Creek valley; a fourth well that taps the same aquifer provides the potable supply for the facility. The study was conducted to document the source areas and quality of the water supply for this Department of Interior facility, which is surrounded by the ongoing development of natural gas from the Marcellus Shale. Groundwater samples were collected from the four wells used by the NARL and from two nearby domestic-supply wells. The domestic-supply wells withdraw groundwater from bedrock of the Catskill Formation. Samples were analyzed for major ions, nutrients, trace metals, radiochemicals, dissolved gases, and stable isotopes of oxygen and hydrogen in water and carbon in dissolved carbonate to document groundwater quality. Organic constituents (other than hydrocarbon gases) associated with hydraulic fracturing and other human activities were not analyzed as part of this assessment. Results show low concentrations of all constituents. Only radon, which ranged from 980 to 1,310 picocuries per liter, was somewhat elevated. These findings are consistent with the pristine nature of the aquifer in the Marsh Creek valley, which is the reason the laboratory was sited at this location. The sources of water and areas contributing recharge to wells were identified by the use of a previously documented MODFLOW groundwater-flow model for the following conditions: (1) withdrawals of 1,000 to 3,000 gallons per minute from the NARL wells, (2) average or dry hydrologic conditions, and (3) withdrawals of 1,000 gallons per minute from a new well 3,500 feet to the southwest that was drilled to provide water for Marcellus gas-well operations. Results of simulations indicate that during average hydrologic conditions, infiltration from Straight Run, a tributary to Marsh Creek, provides nearly all the water to the NARL wells. During dry conditions, the areas contributing recharge expand such that Asaph Run contributes about half of the water to the NARL wells when withdrawals are 1,000 or 2,000 gallons per minute. The addition of a simulated withdrawal of 1,000 gallons per minute from the nearby new well does not substantially affect the sources of water captured by the NARL wells. These results are subject to some limitations. The water-quality samples represent a snapshot of groundwater chemistry for only one hydrologic condition; the concentrations of some constituents may change temporally. In addition, samples were collected and analyzed for hydrocarbon gases, but not organic constituents associated with hydraulic fracturing; additional sampling for these constituents would provide a more complete water-quality baseline. The sources contributing water to the NARL wells and the new well were simulated by use of a simplified one-layer model of the glacial sand and gravel aquifer for steady-state conditions that in reality are never achieved. Steady-state simulations of dry hydrologic conditions show that it is possible for the NARL wells to capture water from Asaph Run; however, maps of simulated groundwater time-of-travel indicate that a dry period of unusually long duration would be required. A better analysis could be done by recalibrating the groundwater-flow model with a finite-difference grid having multiple layers, cells smaller than the 200-foot by 200-foot cells used in this study, and transient stress periods.

Pennsylvania↗

Selected approaches to estimate water-budget components of the High Plains, 1940 through 1949 and 2000 through 2009

The High Plains aquifer, underlying almost 112 million acres in the central United States, is one of the largest aquifers in the Nation. It is the primary water supply for drinking water, irrigation, animal production, and industry in the region. Expansion of irrigated agriculture throughout the past 60 years has helped make the High Plains one of the most productive agricultural regions in the Nation. Extensive withdrawals of groundwater for irrigation have caused water-level declines in many parts of the aquifer and increased concerns about the long-term sustainability of the aquifer. Quantification of water-budget components is a prerequisite for effective water-resources management. Components analyzed as part of this study were precipitation, evapotranspiration, recharge, surface runoff, groundwater discharge to streams, groundwater fluxes to and from adjacent geologic units, irrigation, and groundwater in storage. These components were assessed for 1940 through 1949 (representing conditions prior to substantial groundwater development and referred to as "pregroundwater development" throughout this report) and 2000 through 2009. Because no single method can perfectly quantify the magnitude of any part of a water budget at a regional scale, results from several methods and previously published work were compiled and compared for this study when feasible. Results varied among the several methods applied, as indicated by the range of average annual volumes given for each component listed in the following paragraphs. Precipitation was derived from three sources: the Parameter-Elevation Regressions on Independent Slopes Model, data developed using Next Generation Weather Radar and measured precipitation from weather stations by the Office of Hydrologic Development at the National Weather Service for the Sacramento-Soil Moisture Accounting model, and precipitation measured at weather stations and spatially distributed using an inverse-distance-weighted interpolation method. Precipitation estimates using these sources, as a 10-year average annual total volume for the High Plains, ranged from 192 to 199 million acre-feet (acre-ft) for 1940 through 1949 and from 185 to 199 million acre-ft for 2000 through 2009. Evapotranspiration was obtained from three sources: the National Weather Service Sacramento-Soil Moisture Accounting model, the Simplified-Surface-Energy-Balance model using remotely sensed data, and the Soil-Water-Balance model. Average annual total evapotranspiration estimated using these sources was 148 million acre-ft for 1940 through 1949 and ranged from 154 to 193 million acre-ft for 2000 through 2009. The maximum amount of shallow groundwater lost to evapotranspiration was approximated for areas where the water table was within 5 feet of land surface. The average annual total volume of evapotranspiration from shallow groundwater was 9.0 million acre-ft for 1940 through 1949 and ranged from 9.6 to 12.6 million acre-ft for 2000 through 2009. Recharge was estimated using two soil-water-balance models as well as previously published studies for various locations across the High Plains region. Average annual total recharge ranged from 8.3 to 13.2 million acre-ft for 1940 through 1949 and from 15.9 to 35.0 million acre-ft for 2000 through 2009. Surface runoff and groundwater discharge to streams were determined using discharge records from streamflow-gaging stations near the edges of the High Plains and the Base-Flow Index program. For 1940 through 1949, the average annual net surface runoff leaving the High Plains was 1.9 million acre-ft, and the net loss from the High Plains aquifer by groundwater discharge to streams was 3.1 million acre-ft. For 2000 through 2009, the average annual net surface runoff leaving the High Plains region was 1.3 million acre-ft and the net loss by groundwater discharge to streams was 3.9 million acre-ft. For 2000 through 2009, the average annual total estimated groundwater pumpage volume from two soil-water-balance models ranged from 8.7 to 16.2 million acre-ft. Average annual irrigation application rates for the High Plains ranged from 8.4 to 16.2 inches per year. The USGS Water-Use Program published estimated total annual pumpage from the High Plains aquifer for 2000 and 2005. Those volumes were greater than those estimated from the two soil-water-balance models. Total groundwater in storage in the High Plains aquifer was estimated as 3,173 million acre-ft prior to groundwater development and 2,907 million acre-ft in 2007. The average annual decrease of groundwater in storage between 2000 and 2007 was 10 million acre-ft per year.

High Plains Aquifer↗

Alkalic-type epithermal gold deposit model

This report summarizes the primary characteristics of alkalic-type epithermal gold (Au) deposits and provides an updated descriptive model. These deposits, primarily of Mesozoic to Neogene age, are among the largest epithermal gold deposits in the world. Considered a subset of low-sulfidation epithermal deposits, they are spatially and genetically linked to small stocks or clusters of intrusions containing high alkali-element contents. Deposits occur as disseminations, breccia-fillings, and veins and may be spatially and genetically related to skarns and low-grade porphyry copper (Cu) or molybdenum (Mo) systems. Gold commonly occurs as native gold, precious metal tellurides, and as sub-micron gold in arsenian pyrite. Quartz, carbonate, fluorite, adularia, and vanadian muscovite/roscoelite are the most common gangue minerals. Alkalic-type gold deposits form in a variety of geological settings including continent-arc collision zones and back-arc or post-subduction rifts that are invariably characterized by a transition from convergent to extensional or transpressive tectonics. The geochemical compositions of alkaline igneous rocks spatially linked with these deposits span the alkaline-subalkaline transition. Their alkali enrichment may be masked by potassic alteration, but the unaltered or least altered rocks (1) have chondrite normalized patterns that are commonly light rare earth element (LREE) enriched, (2) are heavy rare earth element (HREE) depleted, and (3) have high large ion lithophile contents and variable enrichment of high-field strength elements. Radiogenic isotopes suggest a mantle derivation for the alkalic magmas but allow crustal contamination. Oxygen and hydrogen isotope compositions show that the fluids responsible for deposit formation are dominantly magmatic, although meteoric or other external fluids (seawater, evolved groundwater) also contributed to the ore-forming fluids responsible for these deposits. Carbon and sulfur isotope compositions in vein-hosted carbonates and sulfide gangue minerals, respectively, coincide with magmatic values, although a sedimentary source of carbon and sulfur is evident in several deposits. Deep-seated structures are critical for the upwelling of hydrous alkalic magmas and for focusing magmatic-hydrothermal fluids to the site of precious metal deposition. The source of gold, silver (Ag), tellurium (Te), vanadium (V), and fluorine (F) was probably the alkalic igneous rocks themselves, and the coexistence of native gold, gold tellurides, and roscoelite in several deposits is primarily a function of similar physicochemical conditions during deposition (for example, overlapping pH and oxygen fugacity ( f O2). Potential environmental impacts related to the mining and processing of alkalic-type epithermal gold deposits include acid mine drainage with high levels of metals, especially zinc (Zn), copper, lead (Pb), and arsenic. However, because alkalic-type gold deposits typically contain carbonates, which contribute calcium and magnesium ions that increase water hardness, aquatic life may be afforded some protection. Impacts vary widely as a function of host rocks, climate, topography, and mining methods. Geologic mapping to (1) highlight the distribution of potassic alteration; (2) define fault density and orientation of structures; (3) determine the distribution of alkaline rocks and hydrothermal breccias; and (4) identify uniquely colored gangue minerals, such as fluorite and roscoelite, will be critical to exploration and future discoveries. Geophysical techniques that identify potassium (K) anomalies (for example, radiometric and spectroscopic surveys), as well as magnetic, resistivity, aeromagnetic, and gravity surveys, may help locate zones of high-permeability that control advecting hydrothermal fluids. Geochemical surveys that include analyses for Au, Ag, barium, Te, K, F, V, Mo, and mercury, which are key elements in these deposits, should be undertaken along with the measurement of other pathfinder elements such as arsenic, bismuth, Cu, iron, nickel, Pb, antimony, selenium, and Zn.

Scientific Investigations Report↗

Water resources of Monroe County, New York, water years 2003-08: Streamflow, constituent loads, and trends in water quality

This report, the sixth in a series published since 1994, presents analyses of hydrologic data in Monroe County for the period October 2002 through September 2008. Streamflows and water quality were monitored at nine sites by the Monroe County Department of Health and the U.S. Geological Survey. Streamflow yields (flow per unit area) were highest in Northrup Creek, which had sustained flows from year-round inflow from the village of Spencerport wastewater-treatment plant and seasonal releases from the New York State Erie (Barge) Canal. Genesee River streamflow yields also were high, at least in part, as a result of higher rainfall and lower evapotranspiration rates in the upper part of the Genesee River Basin than in the other study basins. The lowest streamflow yields were measured in Honeoye Creek, which reflected a decrease in flows due to the withdrawals from Hemlock and Canadice Lakes for the city of Rochester water supply. Water samples collected at nine monitoring sites were analyzed for nutrients, chloride, sulfate, and total suspended solids. The loads of constituents, which were computed from the concentration data and the daily flows recorded at each of the monitoring sites, are estimates of the mass of the constituents that was transported in the streamflow. Annual yields (loads per unit area) also were computed to assess differences in constituent transport among the study basins. All urban sites - Allen Creek and the two downstream sites on Irondequoit Creek - had seasonally high concentrations and annual yields of chloride. Chloride loads are attributed to the application of road-deicing salts to the county's roadways and are related to population and road densities. The less-urbanized sites in the study - Genesee River, Honeoye Creek, and Oatka Creek - had relatively low concentrations and yields of chloride. The highest concentrations and yields of sulfate were measured in Black Creek, Oatka Creek, and Irondequoit Creek at Railroad Mills and are attributable to dissolution of sulfate from gypsum (calcium sulfate) deposits in Silurian shale bedrock that crops out upstream from these monitoring sites. Northrup Creek had the highest concentrations of phosphorus, orthophosphate, and nitrogen, and high yields of nitrate plus nitrite nitrogen and ammonia plus organic nitrogen. These results are attributed to discharges from the Spencerport wastewater-treatment plant (which ceased operation in June 2008), diversions from the New York State Erie (Barge) Canal, and manure and fertilizers applied to agricultural fields. Concentrations and yields of nitrate plus nitrite nitrogen also were high in Oatka Creek and Black Creek; basins with substantial agricultural land uses. Allen Creek had the second highest yield of ammonia plus organic nitrogen. Honeoye Creek, which drains a relatively undeveloped basin, had the lowest yields of nitrogen constituents. The second highest median concentrations and highest sample concentrations of phosphorus and orthophosphate, as well as the highest phosphorus yields, were measured in the Genesee River. A comparison of the yields computed for the two downstream sites on Irondequoit Creek - above Blossom Road and at Empire Boulevard - permitted an assessment of the mitigative effects of the Ellison Park wetland on constituent loads, which would otherwise be transported to Irondequoit Bay. These effects also include those provided by a flow-control structure (installed mid-way through the wetland during February 1997), which was designed to increase the dispersal and short-term detention of stormflows in the wetland. The wetland decreased yields of particulate constituents - phosphorus and ammonia plus organic nitrogen - but had little effect on the yields of dissolved constituents - chloride, sulfate, and nitrate plus nitrite nitrogen. Trends in flow-adjusted concentrations were identified at all sites for most of the nutrient constituents that were evaluated. All of the linear time tren

Scientific Investigations Report↗

Characterization of stormwater runoff from bridge decks in eastern Massachusetts, 2014–16

The quality of stormwater runoff from bridge decks (hereafter referred to as “bridge-deck runoff”) was characterized in a field study from August 2014 through August 2016 in which concentrations of suspended sediment (SS) and total nutrients were monitored. These new data were collected to supplement existing highway-runoff data collected in Massachusetts which were deficient in bridge-deck runoff concentration data. Monitoring stations were installed at three bridges maintained by the Massachusetts Department of Transportation in eastern Massachusetts (State Route 2A in the city of Boston, Interstate 90 in the town of Weston, and State Route 20 near Quinsigamond Village in the city of Worcester). The bridges had annual average daily traffic volumes from 21,200 to 124,000 vehicles per day; the land use surrounding the monitoring stations was 25 to 67 percent impervious. Automatic-monitoring techniques were used to collect more than 160 flow-proportional composite samples of bridge-deck runoff. Samples were analyzed for concentrations of SS, loss on ignition of suspended solids (LOI), particulate carbon (PC), total phosphorus (TP), total dissolved nitrogen (DN), and particulate nitrogen (PN). The distribution of particle size of SS also was determined for composite samples. Samples of bridge-deck runoff were collected year round during rain, mixed precipitation, and snowmelt runoff and with different dry antecedent periods throughout the 2-year sampling period. At the three bridge-deck-monitoring stations, median concentrations of SS in composite samples of bridge-deck runoff ranged from 1,490 to 2,020 milligrams per liter (mg/L); however, the range of SS in individual composites was vast at 44 to 142,000 mg/L. Median concentrations of SS were similar in composite samples collected from the State Route 2A and Interstate 90 bridge (2,010 and 2,020 mg/L, respectively), and lowest at the State Route 20 bridge (1,490 mg/L). Concentrations of coarse sediment (greater than 0.25 millimeters in diameter) dominated the SS matrix by more than an order of magnitude. Concentrations of LOI and PC in composite samples ranged from 15 to 1,740 mg/L and 6.68 to 1,360 mg/L, respectively, and generally represented less than 10 and 3 percent of the median mass of SS, respectively. Concentrations of TP in composite samples ranged from 0.09 to 7.02 mg/L; median concentrations of TP ranged from 0.505 to 0.69 mg/L and were highest on the bridge on State Route 2A in Boston. Concentrations of total nitrogen (TN) (sum DN and PN) in composite samples were variable (0.36 to 29 mg/L). Median DN (0.64 to 0.90 mg/L) concentrations generally represented about 40 percent of the TN concentration at each bridge and were similar to annual volume-weighted mean concentrations of nitrogen in precipitation in Massachusetts. Nonparametric statistical methods were used to test for differences between sample constituent concentrations among the three bridges. These results indicated that there are no statistically significant differences for concentrations of SS, LOI, PC, and TP among the three bridges (one-way analysis of variance test on rank-transformed data, 95-percent confidence level). Test results for concentrations of TN in composite samples indicated that concentrations of TN collected on State Route 20 near Quinsigamond Village were significantly higher than those concentrations collected on State Route 2A in Boston and Interstate 90 near Weston. Median concentrations of TN were about 93 and 55 percent lower at State Route 2A and at Interstate 90, respectively, compared to the median concentrations of TN at State Route 20. Samples of sediment were collected from five fixed locations on each bridge on three occasions during dry weather to calculate semiquantitative distributions of sediment yields on the bridge surface relative to the monitoring location. Mean yields of bridge-deck sediment during this study for State Route 2A in Boston, Interstate 90 near Weston, and State Route 20 near Quinsigamond Village were 1,500, 250, and 5,700 pounds per curb-mile, respectively. Sediment yields at each sampling location varied widely (26 to 25,000 pounds per curb-mile) but were similar to yields reported elsewhere in Massachusetts and the United States. Yields calculated for each sampling location indicated that the sediment was not evenly distributed across each bridge in this study for plausible reasons such as bridge slope, vehicular tracking, and bridge deterioration. Bridge-deck sediment quality was largely affected by the distribution of sediment particle size. Concentrations of TP in the fine sediment-size fraction (less than 0.0625 millimeter in diameter) of samples of bridge-deck sediment were about 6 times greater than in the coarse size fraction. Concentrations for many total-recoverable metals were 2 to 17 times greater in the fine size fraction compared to concentrations in the coarse size fraction (greater than or equal to 0.25 millimeter in diameter), and concentrations of total-recoverable copper and lead in the fine size fraction were 2 to 65 times higher compared to concentrations in the intermediate (greater than or equal to 0.0625 to 0.25 millimeter in diameter) or the coarse size fraction. However, the proportion of sediment particles less than 0.0625 millimeter in diameter in composite samples of bridge-deck runoff was small (median values range from 4 to 8 percent at each bridge) compared to the larger sediment particle-size mass. As a result, more than 50 percent of the sediment-associated TP, aluminum, chromium, manganese, and nickel was estimated to be associated with the coarse size fraction of the SS load. In contrast, about 95 percent of the estimated sediment-associated copper concentration was associated with the fine size fraction of the SS load. Version 1.0.2 of the Stochastic Empirical Loading and Dilution Model was used to simulate long-term (29–30-year) concentrations and annual yields of SS, TP, and TN in bridge-deck runoff and in discharges from a hypothetical stormwater treatment best-management practice structure. Three methods (traditional statistics, robust statistics, and L-moments) were used to calculate statistics for stochastic simulations because the high variability in measured concentration values during the field study resulted in extreme simulated concentrations. Statistics of each dataset, including the average, standard deviation, and skew of the common (base 10) logarithms, for each of the three bridges, and for a lumped dataset, were calculated and used for simulations; statistics representing the median of statistics calculated for the three bridges also were used for simulations. These median statistics were selected for the interpretive simulations so that the simulations could be used to estimate concentrations and yields from other, unmonitored bridges in Massachusetts. Comparisons of the standard and robust statistics indicated that simulation results with either method would be similar, which indicated that the large variability in simulated results was not caused by a few outliers. Comparison to statistics calculated by the L-moments methods indicated that L-moments do not produce extreme concentrations; however, they also do not produce results that represent the bulk of concentration data. The runoff-quality risk analysis indicated that bridge-deck runoff would exceed discharge standards commonly used for large, advanced wastewater treatment plants, but that commonly used stormwater best-management practices may reduce the percentage of exceedances by one-half. Results of simulations indicated that long-term average yields of TN, TP, and SS may be about 21.4, 6.44, and 40,600 pounds per acre per year, respectively. These yields are about 1.3, 3.4, and 16 times simulated ultra-urban highway yields in Massachusetts; however, simulations indicated that use of a best-management practice structure to treat bridge-deck runoff may reduce discharge yields to about 10, 2.8, and 4,300, pounds per acre per year, respectively.

Massachusetts↗

Evaluation of a mass-balance approach to determine consumptive water use in northeastern Illinois

A principal component of evaluating and managing water use is consumptive use. This is the portion of water withdrawn for a particular use, such as residential, which is evaporated, transpired, incorporated into products or crops, consumed by humans or livestock, or otherwise removed from the immediate water environment. The amount of consumptive use may be estimated by a water (mass)-balance approach; however, because of the difficulty of obtaining necessary data, its application typically is restricted to the facility scale. The general governing mass-balance equation is: Consumptive use = Water supplied - Return flows . This study explored a mass-balance field-based computation of consumptive use in a residential setting at the scale of a sanitary sewer service area (sewershed). In addition, the feasibility (cost and difficulty) and relative uncertainties (accuracies) associated with applying the approach at this scale were evaluated. The study was conducted during 2011–13 within a 3.5-square mile (mi 2 ) sewershed confined to a predominantly residential area of Elk Grove Village, Illinois. Following background evaluation of the geohydrologic setting, sewershed infrastructure, and possible components of supplied and returned water, the identified primary components were 1. public water deliveries by the Elk Grove Village Department of Public Works, 2. self-served groundwater withdrawals in an included unincorporated neighborhood with public sanitary sewer service, 3. return flows to the sanitary sewer system, and 4. direct return of water discharged from swimming pools to Salt Creek. Water volumes principally were reported for deliveries, measured for sanitary sewer returns by using an acoustic Doppler current-velocity meter, and estimated for domestic withdrawals and swimming pool discharges to storm sewers. All water volumes required some degree of estimation. Observation wells were installed adjacent to sewer pipelines (lines) to determine the depth of the water table relative to that of the sewer lines and to collect water samples for detection of optical brighteners, as they are routinely discharged as clotheswashing waste to sanitary sewers. These data provided qualitative information on gains (inflow and infiltration) and losses (exfiltration) of sewer flow by pipe leakage, which might otherwise not be considered in the sewer flow return measurements. Hydrographs of sewer flow also were evaluated to identify and estimate storm-associated inputs to sewer flow. The volume of sanitary sewer return flow (778 million gallons per year [Mgal/yr]) was determined to substantially exceed the volume of supplied water (566 Mgal/yr), thus, for this study setting, voiding the utility of the applied mass-balance approach for estimating consumptive water use. Mass-balance components, including sanitary sewer flow and supplied-water use, were estimated within reasonable limits of uncertainty. Evidence of a water table that is typically shallower than the area’s sewer lines, yet is sometimes depressed near more deeply buried sewer lines, suggests groundwater infiltration into the sewers contributes to the excess volume of return flow. Technical obstacles and project resources precluded accurate quantification of infiltration volumes and other gains and losses to sanitary sewer flow. As estimated from various simplified methods, a minimum of 26 percent of return flow measured in the sanitary sewer represented groundwater infiltration and stormwater inflow; separately, about 2 percent of return flow was estimated as inflow. On the basis of the alternative winter base-rate method, consumptive use in the sewershed was estimated as 13 percent, which compares favorably with that used by the State of Illinois for Lake Michigan allocation accounting (10 percent) and other States and Canadian Provinces in the Great Lakes region (generally 10-15 percent). The study also provided other findings considered useful to studies of water use and to performance evaluation of sanitary sewer infrastructure. In urban residential settings, the comparatively small volumes of nonpublic sources of water (self-supplied) and direct (nonsanitary) return flow potentially can be ignored in the estimation of consumptive use. An acoustic Doppler current-velocity meter can be used in sanitary sewers to accurately measure discharge and reasonably estimate storm-associated inflows. Hourly to daily patterns of water use can be readily identified and quantified in the return flow record for the sanitary sewers. Relative volumes of infiltration gains (and exfiltration losses) can be substantial, even in sewer systems of communities making significant investments in system upgrades to limit sewer line leakage. Monitoring of optical brighteners in groundwater (and potentially in sanitary sewer flow) can provide a useful means of identifying probable leakage from (and to) sewer lines. Accurate quantification of gains and losses to sanitary sewer flow at the sewershed scale will require additional research effort and technical advances. Under ideal conditions, accurate quantification of consumptive use at the sewershed scale by the described mass-balance approach might be possible. Under most prevailing conditions, quantification likely would be more costly and time consuming than that of the present study, given the freely contributed technical support of the host community and relatively appropriate conditions of the study area. Essentials to quantification of consumptive use are a fully cooperative community, storm and sanitary sewers that are separate, and newer sewer infrastructure and (or) a robust program for limiting infiltration, exfiltration, and inflow.

Illinois↗

Groundwater discharge characteristics for selected streams within the Loup River Basin, Nebraska, 2014–16

Streams in the Loup River Basin are sensitive to groundwater withdrawals because of the close hydrologic connection between groundwater and surface water. Groundwater discharge is the primary component of streamflow in the Loup River Basin and constitutes more than 90 percent of streamflow in the central part of the Sand Hills. To improve the understanding of geologic controls and various climatic and land-use changes on groundwater discharge, the U.S. Geological Survey (USGS), in cooperation with the Upper Loup Natural Resources District (NRD), the Lower Loup NRD, and the Nebraska Environmental Trust, studied the spatial and temporal characteristics of groundwater discharge within the Loup River Basin. This report documents the methods of data collection and analysis, which include the collection of approximately 350 river miles of aerial thermal infrared imagery and continuous groundwater-level and temperature data from six streamflow-gaging stations within the Loup River Basin. The results from the stream reconnaissance and examination of aerial thermal infrared imagery demonstrated the influence of the surficial and subsurface geology on the spatial characteristics of groundwater discharge to streams in the Loup River Basin. At the headwaters of the South Loup River, streamflow is sustained and increased from focused groundwater discharge emanating from Quaternary deposits at many small (less than 0.1 cubic foot per second) focused points. The volume of water produced from this dense network of focused groundwater discharge points along the North Fork South Loup River is sufficient to provide approximately 40 percent of the flow measured at the South Loup River at Arnold, Nebraska streamflow-gaging station (USGS station 06781600) during the irrigation season. Approximately 5 miles downstream from the South Loup River at Arnold, Nebr., streamflow-gaging station, the river incises into Pliocene-age sand and gravel deposits, which provide additional groundwater discharge to the stream. The streamflow of the South Loup River increases by a factor of 5 across a 62-mile reach of the middle South Loup River. Increases in streamflow along the upper Dismal River result from a dense network of focused groundwater discharge points within semiconsolidated Pliocene-age deposits. Below the Dismal River near Thedford, Nebr., streamflow-gaging station (USGS station 06775900), the Dismal River incises into the Ogallala Formation over a short reach before flowing over coarser, more permeable Quaternary-age alluvial deposits. Diffuse groundwater discharge sustains and increases the streamflow of the lower Dismal River in this reach. Groundwater sapping was evident on some stream reaches and has increased the size and flow of focused groundwater discharge points. Previous researchers have documented streambed incision and groundwater sapping on the upper Dismal River that have created and enlarged focused groundwater discharge points capturing additional groundwater. Similar processes appear to have played a role in the formation of larger focused groundwater discharge points, which sustain the flow of the middle South Loup River. The constant flow of groundwater into the South Loup River has removed finer-grained Quaternary sediments and further exposed Pliocene-age gravel deposits. Headward erosion is evident where some of the large focused groundwater discharge points have incised their own draws and terminate in bowl-like depressions away from the stream. Within the Loup River NRDs, the percentage of groundwater-irrigated land in a stream basin is one factor that affects groundwater discharge to streams. A striking example was at the South Loup River at Saint Michael, Nebr., groundwater and streamflow-gaging station (USGS station 06784000) where the shallow groundwater levels declined below the level of the stream during the middle to late part of the growing season (July to September) when consumptive groundwater use was at its peak. The South Loup River Basin above the South Loup River at Saint Michael, Nebr., streamflow-gaging station has the highest percentage of groundwater-irrigated row crops of all the basins examined in this study. Continuous groundwater and surface-water levels measured at the North Loup River at the Taylor, Nebr., streamflow-gaging station (USGS station 06786000) indicate that the stream is receiving groundwater throughout the year; however, when consumptive groundwater use peaks during the middle to late part of the growing season (July to September), the difference in elevation between the groundwater level and the stream elevation decreases, which indicates a reduction in the amount of groundwater discharge received.

Nebraska↗

Geology, Water, and Wind in the Lower Helmand Basin, Southern Afghanistan

This report presents an overview of the geology, hydrology, and climate of the lower Helmand Basin, a large, closed, arid basin in southern Afghanistan. The basin is drained by the Helmand River, the only perennial desert stream between the Indus and Tigris-Euphrates Rivers. The Helmand River is the lifeblood of southern Afghanistan and has supported desert civilizations in the Sistan depression for over 6,000 years. The Helmand Basin is a structurally closed basin that began to form during the middle Tertiary as a consequence of the collision of several Gondwanaland fragments. Aeromagnetic studies indicate the basin is 3-5 kilometers deep over basement rocks. Continued subsidence along basin-bounding faults in Iran and Pakistan throughout the Neogene has formed the Sistan depression in the southwest corner of the basin. Lacustrine, eolian, and fluvial deposits are commonly exposed in the basin and were intruded by latest Miocene-middle Quaternary volcanoes, which indicates that depositional environments in the lower Helmand Basin have not substantially changed for nearly 10 million years. Lakes expanded in the Sistan depression during the Quaternary; however, the size and extent of these pluvial lakes are unknown. Climate conditions in the lower Helmand Basin likely mirrored climate changes in the Rajasthan Desert to the east and in Middle Eastern deserts to the west: greater aridity during global episodes of colder temperatures and increased available moisture during episodes of warmer temperatures. Eolian processes are unusually dominant in shaping the landscape in the basin. A strong wind blows for 120 days each summer, scouring dry lakebeds and creating dune fields from annual flood deposits. Nearly one-third of the basin is mantled with active or stabilized dunes. Blowing winds combined with summer temperatures over 50? Celsius and voluminous insect populations hatched from the deltaic wetlands create an environment referred to as the 'most odious place on earth' by 19th century visitors. During dry years, large plumes of dust originating from Sistan are recorded by weather satellites. The Helmand River drains about 40 percent of Afghanistan and receives most of its moisture from melting snow and spring storms. Similar to many desert streams, the Helmand and its main tributary, the Arghandab River, are characterized by large fluctuations in monthly and annual discharges. Water from the Helmand accumulates in several hamuns (shallow lakes) in the Sistan depression. The wetlands surrounding these hamuns are the largest in western Asia and are directly affected by droughts and floods on the Helmand. Average annual discharge on the Helmand is about 6.12 million megaliters (million cubic meters), and the annual discharge varies by a factor of five. In 2005, the region was just beginning to recover from the longest drought (1998-2005) of record back to 1830. Annual peak discharges range from less than 80 cubic meters per second in 1971 to nearly 19,000 cubic meters per second in 1885. Large floods fill each hamun to overflowing to create one large lake that overflows into the normally dry Gaud-i Zirreh basin. The interaction of flooding, active subsidence, and wind erosion causes frequent channel changes on the Helmand delta. A major development effort on the Helmand River was initiated after World War II with substantial aid from the United States. Two dams and several major canals were completed in the 1950s; however, poor drainage conditions on the newly prepared agricultural fields caused extensive waterlogging and salinization. New drains were installed and improved agricultural methods were implemented in the 1970s, and some lands became more productive. Since 1980, Afghanistan has endured almost constant war and civil and political strife. In 2005, the country was on a path to rebuild much of its technical infrastructure. Revitalization of agricultural lands in the lower Helmand Basin and improved managem

Scientific Investigations Report↗

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota↗

Temporal Differences in Flow Depth and Velocity Distributions and Hydraulic Microhabitats Near Bridges of the Lower Platte River, Nebraska, 1934-2006

As part of a collaborative study of the cumulative impacts on stream and riparian ecology of water and channel management practices in the lower Platte River, Nebraska, this report describes a study by the U.S. Geological Survey in cooperation with the Lower Platte South Natural Resources District that summarizes: (1) temporal differences in distribution of streamflow depth, velocity, and microhabitats among five discrete 11-water-year periods 1934-44, 1951-61, 1966-76, 1985-95, and 1996-2006, and (2) the effects of bridge proximity on distribution of streamflow depth, velocity, and microhabitat of the Platte River when cross sections were measured at a similar discharge. The scope of the study included the four presently (2008) active streamflow-gaging stations located near bridges over the lower Platte River at North Bend, near Leshara, near Ashland, and at Louisville, Nebraska, and the most downstream streamflow-gaging station within the central Platte River segment near Duncan, Nebraska. Generally, in cases where temporal differences in streamflow depth and velocity were evident, at least one of the water-year periods from 1934 through 1995 had deeper streamflow than the recent water-year period (1996-2006). Temporal differences in distributions of streamflow depth were not strongly associated with differences in either climatic conditions or the maximum peak flow that occurred prior to the latest discharge measurement during each period. The relative cross-sectional area of most hydraulic niches did not differ among the water-year periods. Part of this apparent uniformity likely was an artifact of the broad microhabitat classification used for this study. In cases where temporal differences in relative cross-sectional area of hydraulic niches were evidenced, the differences occurred during high- and low-flow conditions, not during median flow conditions. The temporal differences in relative cross-sectional area were found more frequently for hydraulic niches defined by moderate and fast velocities than for hydraulic niches defined by slow velocities. Generally, any significant increase or decrease in the relative cross-sectional areas of hydraulic niches during the water-year periods from 1934 through 1995 had disappeared during the most recent water-year period, 1996-2006. Deep-Swift niche was the predominant hydraulic niche for all near-bridge sites on the lower Platte River for high- and median-flow conditions. The Deep-Swift niche also was the predominant niche for the near-bridge sites near Ashland and at Louisville for low-flow conditions; for the near-bridge sites at North Bend and near Leshara, streamflow cross-sectional areas during low-flow conditions were shared among the Shallow-Moderate, Intermediate-Moderate, Intermediate-Swift, and Deep-Swift hydraulic niches. For the near-bridge site near Duncan, the site farthest downstream in the central Platte River system, the Deep-Swift hydraulic niche was predominant only during high-flow conditions; during median- and low-flow conditions the relative cross-sectional area was shared among the Shallow-Slow, Shallow-Moderate, Intermediate-Moderate, and Intermediate-Swift hydraulic niches. Significant temporal differences in the relative cross-sectional area of the Deep-Swift hydraulic niche were found for sites near the two farthest downstream bridges near Ashland and at Louisville, but only for low-flow conditions. The Deep-Swift microhabitat was of special interest because it is the preferred hydraulic habitat during the adult life of the endangered pallid sturgeon (Scaphirhynchus albus). Temporal differences in relative cross-sectional areas of the Glide low-flow geomorphic microhabitat that contained the Deep-Swift hydraulic niche also indicated that relative cross-sectional areas of the Glide during the 1951-61 and 1996-2006 water-year periods were lower than during the 1966-76 period. The temporal differences indicated that any significant temporal chang

Scientific Investigations Report↗

Application of a stream-aquifer model to Monument Creek for development of a method to estimate transit losses for reusable water, El Paso County, Colorado

The U.S. Geological Survey, in cooperation with Colorado Springs Utilities, the Colorado Water Conservation Board, and the El Paso County Water Authority, began a study in 2004 to (1) apply a stream-aquifer model to Monument Creek, (2) use the results of the modeling to develop a transit-loss accounting program for Monument Creek, (3) revise the existing transit-loss accounting program for Fountain Creek to incorporate new water-management strategies and allow for incorporation of future changes in water-management strategies, and (4) integrate the two accounting programs into a single program with a Web-based user interface. The purpose of this report is to present the results of applying a stream-aquifer model to the Monument Creek study reach. Transit losses were estimated for reusable-water flows in Monument Creek that ranged from 1 to 200 cubic feet per second (ft 3 /s) and for native streamflows that ranged from 0 to 1,000 ft 3 /s. Transit losses were estimated for bank-storage, channel-storage, and evaporative losses. The same stream-aquifer model used in the previously completed (1988) Fountain Creek study was used in the Monument Creek study. Sixteen model nodes were established for the Monument Creek study reach, defining 15 subreaches. Channel length, aquifer length, and aquifer width for the subreaches were estimated from available topographic and geologic maps. Thickness of alluvial deposits and saturated thickness were estimated using lithologic and water-level data from about 100 wells and test holes in or near the Monument Creek study reach. Estimated average transmissivities for the subreaches ranged from 2,000 to 12,000 feet squared per day, and a uniform value of 0.20 was used for storage coefficient. Qualitative comparison of recorded and simulated streamflow at the downstream node for the calibration and verification simulations indicated that the two streamflows compared reasonably well. No adjustments were made to the model parameters. Differences between recorded and simulated streamflow volumes for all calibration and verification simulations ranged from about –8.8 to 7.5 percent; the total error for all simulations was about –0.7 percent. The model was used to estimate bank-storage losses for 10 to 15 native streamflows for each reusable-water flow of 1, 3, 5, 7, 10, 15, 20, 30, 40, 50, 100, and 200 ft 3 /s. Then the 10 to 15 bank-storage loss values were used in least-squares linear regression to estimate a relation between bank-storage loss and native streamflow for each of the 12 reusable-water flow rates. The 12 regression relations then were used to develop “look-up” tables of bank-storage loss for reusable-water flows ranging from 1 to 200 ft 3 /s (in 1-ft 3 /s increments). Additional model simulations indicated that (1) when the ratio of downstream native streamflow to upstream native streamflow was less than 1, bank-storage loss generally increased and (2) when the ratio of downstream native streamflow to upstream native streamflow was larger than 1, bank-storage loss generally decreased. These results were used to develop a bank-storage loss adjustment factor based on the ratio of native streamflow at the downstream node to native streamflow at the upstream node. The model also was used to estimate a recovery period, which is the length of time needed for the bank-storage loss to return to the stream. The recovery period was 1 day for six subreaches; 2 days for four subreaches; between 3 and 12 days for four subreaches; and 28 days for one subreach. Channel-storage losses are about 10 percent of the reusable-water flow for most of the subreaches, except for two subreaches, where the channel-storage losses are about 20 percent, and one subreach, where the losses are about 30 percent, owing to the greater channel lengths. Evaporative losses were estimated by the use of monthly pan-evaporation data and the incremental increase in stream width resulting from any reusable-water flows. Monthly pan-evaporation data were converted to a daily rate. The daily rate, when multiplied by the stream-width increase (in feet) that results from reusable-water flow and by the subreach length (in miles) gives the daily evaporative loss in cubic feet per second.

Colorado↗