USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Reports”

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 919 records · Page 51Linked to original sources

Sources and preparation of data for assessing trends in concentrations of pesticides in streams of the United States, 1992-2006

This report provides a water-quality data set of 44 commonly used pesticides and 8 pesticide degradates suitable for a national assessment of trends in pesticide concentrations in streams of the United States. Water-quality samples collected from January 1992 through August 2006 at stream-water sites of the U.S. Geological Survey National Water-Quality Assessment Program and the National Stream Quality Accounting Network Program were compiled, reviewed, selected, and prepared for trend analysis as described in this report. Samples analyzed at the U.S. Geological Survey National Water Quality Laboratory by a gas chromatography/mass spectrometry analytical method were the most extensive in time and space and were selected for national trend analysis. The selection criteria described in the report produced a trend data set of 16,869 pesticide samples at 201 stream and river sites.

Scientific Investigations Report↗

A deposit model for magmatic iron-titanium-oxide deposits related to Proterozoic massif anorthosite plutonic suite

This descriptive model for magmatic iron-titanium-oxide (Fe-Ti-oxide) deposits hosted by Proterozoic age massif-type anorthosite and related rock types presents their geological, mineralogical, geochemical, and geoenvironmental attributes. Although these Proterozoic rocks are found worldwide, the majority of known deposits are found within exposed rocks of the Grenville Province, stretching from southwestern United States through eastern Canada; its extension into Norway is termed the Rogaland Anorthosite Province. This type of Fe-Ti-oxide deposit dominated by ilmenite rarely contains more than 300 million tons of ore, with between 10- to 45-percent titanium dioxide (TiO 2 ), 32- to 45-percent iron oxide (FeO), and less than 0.2-percent vanadium (V). The origin of these typically discordant ore deposits remains as enigmatic as the magmatic evolution of their host rocks. The deposits clearly have a magmatic origin, hosted by an age-constrained unique suite of rocks that likely are the consequence of a particular combination of tectonic circumstances, rather than any a priori temporal control. Principal ore minerals are ilmenite and hemo-ilmenite (ilmenite with extensive hematite exsolution lamellae); occurrences of titanomagnetite, magnetite, and apatite that are related to this deposit type are currently of less economic importance. Ore-mineral paragenesis is somewhat obscured by complicated solid solution and oxidation behavior within the Fe-Ti-oxide system. Anorthosite suites hosting these deposits require an extensive history of voluminous plagioclase crystallization to develop plagioclase-melt diapirs with entrained Fe-Ti-rich melt rising from the base of the lithosphere to mid- and upper-crustal levels. Timing and style of oxide mineralization are related to magmatic and dynamic evolution of these diapiric systems and to development and movement of oxide cumulates and related melts. Active mines have developed large open pits with extensive waste-rock piles, but because of the nature of the ore and waste rock, the major environmental impacts documented at the mine sites are reported to be waste disposal issues and somewhat degraded water quality.

Scientific Investigations Report↗

Estimation of natural historical flows for the Manitowish River near Manitowish Waters, Wisconsin

The Wisconsin Department of Natural Resources is charged with oversight of dam operations throughout Wisconsin and is considering modifications to the operating orders for the Rest Lake Dam in Vilas County, Wisconsin. State law requires that the operation orders be tied to natural low flows at the dam. Because the presence of the dam confounds measurement of natural flows, the U.S. Geological Survey, in cooperation with the Wisconsin Department of Natural Resources, installed streamflow-gaging stations and developed two statistical methods to improve estimates of natural flows at the Rest Lake Dam. Two independent methods were used to estimate daily natural flow for the Manitowish River approximately 1 mile downstream of the Rest Lake Dam. The first method was an adjusted drainage-area ratio method, which used a regression analysis that related measured water yield (flow divided by watershed area) from short-term (2009–11) gaging stations upstream of the Manitowish Chain of Lakes to the water yield from two nearby long-term gaging stations in order to extend the flow record (1991–2011). In this approach, the computed flows into the Chain of Lakes at the upstream gaging stations were multiplied by a coefficient to account for the monthly hydrologic contributions (precipitation, evaporation, groundwater, and runoff) associated with the additional watershed area between the upstream gaging stations and the dam at the outlet of the Chain of Lakes (Rest Lake Dam). The second method used to estimate daily natural flow at the Rest Lake Dam was a water-budget approach, which used lake stage and dam outflow data provided by the dam operator. A water-budget model was constructed and then calibrated with an automated parameter-estimation program by matching simulated flow-duration statistics with measured flow-duration statistics at the upstream gaging stations. After calibration of the water-budget model, the model was used to compute natural flow at the dam from 1973 to 2011. Daily natural flows at the dam, as computed by the adjusted drainage-area ratio method and the water-budget method, were used to compute monthly flow-duration values for the period of historical data available for each method. Monthly flow-durations provide a means for evaluating the frequency and range in flows that have been observed for each month over the course of many years. Both methods described the pattern and timing of measured high-flow and low-flow events at the upstream gaging stations. The adjusted drainage-area ratio method generally had smaller residual errors across the full range of observed flows and had smaller monthly biases than the water-budget method. Although it is not possible to evaluate which method may be more "correct" for estimating monthly natural flows at the dam, comparisons between the results of each method indicate that the adjusted drainage-area ratio method may be susceptible to biases at high flows due to isolated storms outside of the Manitowish River watershed. Conversely, it appears that the water-budget method may be susceptible to biases at low flows because of its sensitivity to the accuracy of reported lake stage and outflows, as well as effects of upstream diversions that could not be fully compensated for with this method. Results from both methods are useful for understanding the natural flow patterns at the dam. Flows for both methods have similar patterns, with high median flows in spring and low median flows in late summer. Similarly, the range from monthly high-flow durations to low-flow durations increases during spring, decreases during summer, and increases again during fall. These seasonal patterns illustrate a challenge with interpreting a single value of natural low flow. That is, a natural low flow computed for September is not representative of a natural low flow in April. Moreover, alteration of natural flows caused by storing water in the Chain of Lakes during spring and releasing it in fall causes a change in the timing of high and low flows compared with natural conditions. That is, the lowest reported dam outflows occurred in spring and highest reported outflows occurred in fall, which is opposite the natural patterns.

Wisconsin↗

A deposit model for magmatic iron-titanium-oxide deposits related to Proterozoic massif anorthosite plutonic suites

This descriptive model for magmatic iron-titanium-oxide (Fe-Ti-oxide) deposits hosted by Proterozoic age massif-type anorthosite and related rock types presents their geological, mineralogical, geochemical, and geoenvironmental attributes. Although these Proterozoic rocks are found worldwide, the majority of known deposits are found within exposed rocks of the Grenville Province, stretching from southwestern United States through eastern Canada; its extension into Norway is termed the Rogaland Anorthosite Province. This type of Fe-Ti-oxide deposit dominated by ilmenite rarely contains more than 300 million tons of ore, with between 10- to 45-percent titanium dioxide (TiO 2 ), 32- to 45-percent iron oxide (FeO), and less than 0.2-percent vanadium (V). The origin of these typically discordant ore deposits remains as enigmatic as the magmatic evolution of their host rocks. The deposits clearly have a magmatic origin, hosted by an age-constrained unique suite of rocks that likely are the consequence of a particular combination of tectonic circumstances, rather than any a priori temporal control. Principal ore minerals are ilmenite and hemo-ilmenite (ilmenite with extensive hematite exsolution lamellae); occurrences of titanomagnetite, magnetite, and apatite that are related to this deposit type are currently of less economic importance. Ore-mineral paragenesis is somewhat obscured by complicated solid solution and oxidation behavior within the Fe-Ti-oxide system. Anorthosite suites hosting these deposits require an extensive history of voluminous plagioclase crystallization to develop plagioclase-melt diapirs with entrained Fe-Ti-rich melt rising from the base of the lithosphere to mid- and upper-crustal levels. Timing and style of oxide mineralization are related to magmatic and dynamic evolution of these diapiric systems and to development and movement of oxide cumulates and related melts. Active mines have developed large open pits with extensive waste-rock piles, but because of the nature of the ore and waste rock, the major environmental impacts documented at the mine sites are reported to be waste disposal issues and somewhat degraded water quality.

Scientific Investigations Report↗

Methods for estimating regional skewness of annual peak flows in parts of the Great Lakes and Ohio River Basins, based on data through water year 2013

Bulletin 17C (B17C) recommends fitting the log-Pearson Type III (LP−III) distribution to a series of annual peak flows at a streamgage by using the method of moments. The third moment, the skewness coefficient (or skew), is important because the magnitudes of annual exceedance probability (AEP) flows estimated by using the LP−III distribution are affected by the skew; interest is focused on the right-hand tail of the distribution, which represents the larger annual peak flows that correspond to small AEPs. For streamgages having modest record lengths, the skew is sensitive to extreme events like large floods, which cause a sample to be highly asymmetrical or “skewed.” For this reason, B17C recommends using a weighted-average skew computed from the station skew for a given streamgage and a regional skew. This report generates an estimate of regional skew for a study area encompassing most of the Great Lakes Basin (hydrologic unit 04) and part of the Ohio River Basin (hydrologic unit 05). A total of 551 candidate streamgages that were unaffected by extensive regulation, diversion, urbanization, or channelization were considered for use in the skew analysis; after screening for redundancy and pseudo record length greater than 36 years, 368 streamgages were selected for use in the study. Flood frequencies for candidate streamgages were analyzed by employing the Expected Moments Algorithm, which extends the method of moments so that it can accommodate interval, censored, and historic/paleo flow data, as well as the Multiple Grubbs-Beck test to identify potentially influential low floods in the data series. Bayesian weighted least squares/Bayesian generalized least squares regression was used to develop a regional skew model for the study area that would incorporate possible variables (basin characteristics) to explain the variation in skew in the study area. Twelve basin characteristics were considered as possible explanatory variables; however, none produced a pseudo coefficient of determination greater than 5 percent; as a result, these characteristics did not help to explain the variation in skew in the study area. Therefore, a constant model having a regional skew coefficient of 0.086 and an average variance of prediction ( AVP new ) (which corresponds to the mean square error [MSE]) of 0.13 at a new streamgage was selected. The AVP new corresponds to an effective record length of 54 years, a marked improvement over the Bulletin 17B national skew map, whose reported MSE of 0.302 indicated a corresponding effective record length of only 17 years.

Great Lakes and Ohio River Basins↗

Methods for estimating regional skewness of annual peak flows in parts of eastern New York and Pennsylvania, based on data through water year 2013

Bulletin 17C (B17C) recommends fitting the log-Pearson Type III (LP−III) distribution to a series of annual peak flows at a streamgage by using the method of moments. The third moment, the skewness coefficient (or skew), is important because the magnitudes of annual exceedance probability (AEP) flows estimated by using the LP–III distribution are affected by the skew; interest is focused on the right-hand tail of the distribution, which represents the larger annual peak flows that correspond to small AEPs. For streamgages having modest record lengths, the skew is sensitive to extreme events like large floods, which cause a sample to be highly asymmetrical or “skewed.” For this reason, B17C recommends using a weighted-average skew computed from the skew of the annual peak flows for a given streamgage and a regional skew. This report presents an estimate of regional skew for a study area encompassing parts of eastern New York and Pennsylvania. A total of 232 candidate U.S. Geological Survey streamgages that were unaffected by extensive regulation, diversion, urbanization, or channelization were considered for use in the skew analysis; after screening for redundancy and pseudo record length ( P RL ) of at least 36 years, 183 streamgages were selected for use in the study. Flood frequencies for candidate streamgages were analyzed by employing the expected moments algorithm, which extends the method of moments so that it can accommodate interval, censored, and historical/paleo flow data, as well as the multiple Grubbs-Beck test to identify potentially influential low floods in the data series. Bayesian weighted least squares/Bayesian generalized least squares regression was used to develop a regional skew model for the study area that would incorporate possible variables (basin characteristics) to explain the variation in skew in the study area. Ten basin characteristics were considered as possible explanatory variables; however, none produced a pseudo coefficient of determination greater than 1 percent; as a result, these characteristics did not help to explain the variation in skew in the study area. Therefore, a constant model that had a regional skew coefficient of 0.32 and an average variance of prediction at a new streamgage ( AVP new , which corresponds to the mean square error [MSE] of 0.11) was selected. The AVP new corresponds to an effective record length of 68 years, a marked improvement over the Bulletin 17B national skew map, whose reported MSE of 0.302 indicated a corresponding effective record length of only 17 years.

New York, Pennsylvania↗

Creek and quarry water quality at Pipestone National Monument and pilot study of pathogen detection methods in waterfall mist at Winnewissa Falls, Pipestone, Minnesota, 2018–19

Pipestone National Monument is a 301-acre site sacred to many Native American Tribes, providing cultural exhibits and walking trails to Pipestone Creek, Winnewissa Falls, and historical pipestone quarries for numerous visitors each year. However, the Minnesota Pollution Control Agency has determined turbidity and fecal coliform bacteria occur in Pipestone Creek in high enough numbers to be a potential health hazard. Concerns also were raised about exposure risk from waterfall mist to visitors and staff. The U.S. Geological Survey and the National Park Service collaborated on a study to collect 21 water-quality samples from 8 creek sites and 3 quarries in 2018 and analyzed them for over 250 water-quality parameters and contaminants. Additional samples were collected in August 2019 to assess the waterfall mists from Winnewissa Falls. Nutrient concentrations in the creek and quarries were elevated in 2018, indicating they are affected by agricultural inputs. All sample concentrations for nitrate and total nitrogen in Pipestone Creek exceeded Minnesota standards and U.S. Environmental Protection Agency nutrient criteria. Minnesota standards and U.S. Environmental Protection Agency nutrient criteria for total phosphorus also were exceeded in some of the quarry samples. Twenty of 210 micropollutants had measurable concentrations: 13 pesticides, 5 pharmaceuticals, and 2 other types of micropollutants. Atrazine, deethylatrazine, and metolachlor ethanesulfonic acid were detected in all 21 samples collected during the study. The five pharmaceuticals detected were acetaminophen, gabapentin, gemfibrozil, metformin, and oxycodone. Gabapentin (10 of 21 samples) and metformin (8 of 21 samples) were the most commonly detected pharmaceuticals. None of the detected micropollutant concentrations exceeded any Minnesota standards or U.S. Environmental Protection Agency aquatic life benchmarks, except the acute toxicity benchmark for nonvascular plants for atrazine. Two cyanotoxins, anatoxin-a and microcystin, were detected, but concentrations were below U.S. Environmental Protection Agency guidelines for swimming or recreation. Notably, total coliform, fecal coliform, and Escherichia coli were detected in all creek samples, and concentrations generally decreased downstream, suggesting contamination potentially occurred upstream from the monument. Mycobacterium avium ssp. paratuberculosis was not detected in any creek sediment samples but was detected in three water samples from the creek. Three organisms were detected in the 2019 water and mist sampling from Winnewissa Falls. Two of these organisms can cause illness in humans ( Cryptosporidium and Legionella ), and a third (ruminant Bacteroides ) is an indicator of manure contamination. Despite few samples, pathogen-positive water samples and air sampling demonstrated the feasibility and utility of the mist sampling approach outlined in this report.

Minnesota↗

2023 Earthquake Ground-Motion Workshop for the Central and Eastern United States, with a focus on the Gulf and Atlantic Coastal Plains—Agenda and abstracts

The U.S. Geological Survey held a virtual workshop December 7–8, 2023, to share research and ideas about earthquake ground motions in the Central and Eastern United States, with a focus on the Atlantic and Gulf Coastal Plains. The workshop was organized to learn about potential regionalization of ground-motion characteristics (source, path, and site), consider new explanatory variables for site response, and hear and discuss updates on ground-motion research on the Atlantic and Gulf Coastal Plains. The workshop was organized into a series of contributed presentations and three panel discussions held during 2 days. This report documents the agenda, contributed abstracts, and panel summaries.

Scientific Investigations Report↗

Mercury and methylmercury in water and bottom sediments of wetlands at Lostwood National Wildlife Refuge, North Dakota, 2003-04

Certain ecosystem types, particularly wetlands, have environmental characteristics that can make them particularly sensitive to mercury inputs and that can result in large mercury concentrations in fish or other aquatic biota. To provide information needed to make effective management decisions to decrease human and wildlife exposure to methylmercury in northern prairie pothole wetlands, the U.S. Geological Survey, in cooperation with the North Dakota Department of Health, conducted a study to assess mercury and methylmercury concentrations in wetlands at the Lostwood National Wildlife Refuge (the Refuge) in northwest North Dakota. In April 2003 and 2004, water and bottom-sediment samples were collected from 44 individual wetlands that were classified as one of four wetland types. Many factors that may affect methylmercury production were considered in the study. The prairie pothole wetlands at the Refuge had large ranges in major environmental characteristics. Hydrologic differences, most notably semiannual wetting and drying cycles, that are intrinsic to prairie pothole wetlands affected methylmercury concentrations. This likely resulted from the stimulation of anaerobic microbial activity following reflooding of soils, particularly soils containing substantial organic carbon. Among the four wetland types considered for this study, seasonal and semipermanent wetlands generally had the largest methylmercury concentrations. Regardless of wetland type, however, methylmercury concentrations at the Refuge are large in relation to reported concentrations for natural aquatic systems.

North Dakota↗

Preliminary Volcano-Hazard Assessment for Gareloi Volcano, Gareloi Island, Alaska

Gareloi Volcano (178.794 degrees W and 51.790 degrees N) is located on Gareloi Island in the Delarof Islands group of the Aleutian Islands, about 2,000 kilometers west-southwest of Anchorage and about 150 kilometers west of Adak, the westernmost community in Alaska. This small (about 8x10 kilometer) volcano has been one of the most active in the Aleutians since its discovery by the Bering expedition in the 1740s, though because of its remote location, observations have been scant and many smaller eruptions may have gone unrecorded. Eruptions of Gareloi commonly produce ash clouds and lava flows. Scars on the flanks of the volcano and debris-avalanche deposits on the adjacent seafloor indicate that the volcano has produced large landslides in the past, possibly causing tsunamis. Such events are infrequent, occurring at most every few thousand years. The primary hazard from Gareloi is airborne clouds of ash that could affect aircraft. In this report, we summarize and describe the major volcanic hazards associated with Gareloi.

Scientific Investigations Report↗

Summary of the Georgia Agricultural Water Conservation and Metering Program and evaluation of methods used to collect and analyze irrigation data in the middle and lower Chattahoochee and Flint River basins, 2004-2010

Since receiving jurisdiction from the State Legislature in June 2003 to implement the Georgia Agricultural Water Conservation and Metering Program, the Georgia Soil and Water Conservation Commission (Commission) by year-end 2010 installed more than 10,000 annually read water meters and nearly 200 daily reporting, satellite-transmitted, telemetry sites on irrigation systems located primarily in southern Georgia. More than 3,000 annually reported meters and 50 telemetry sites were installed during 2010 alone. The Commission monitored rates and volumes of agricultural irrigation supplied by groundwater, surface-water, and well-to-pond sources to inform water managers on the patterns and amounts of such water use and to determine effective and efficient resource utilization. Summary analyses of 4 complete years of irrigation data collected from annually read water meters in the middle and lower Chattahoochee and Flint River basins during 2007-2010 indicated that groundwater-supplied fields received slightly more irrigation depth per acre than surface-water-supplied fields. Year 2007 yielded the largest disparity between irrigation depth supplied by groundwater and surface-water sources as farmers responded to severe-to-exceptional drought conditions with increased irrigation. Groundwater sources (wells and well-to-pond systems) outnumbered surface-water sources by a factor of five; each groundwater source applied a third more irrigation volume than surface water; and, total irrigation volume from groundwater exceeded that of surface water by a factor of 6.7. Metered irrigation volume indicated a pattern of low-to-high water use from northwest to southeast that could point to relations between agricultural water use, water-resource potential and availability, soil type, and crop patterns. Normalizing metered irrigation-volume data by factoring out irrigated acres allowed irrigation water use to be expressed as an irrigation depth and nearly eliminated the disparity between volumes of applied irrigation derived from groundwater and surface water. Analysis of per-acre irrigation depths provided a commonality for comparing irrigation practices across the entire range of field sizes in southern Georgia and indicated underreporting of irrigated acres for some systems. Well-to-pond systems supplied irrigation at depths similar to groundwater and can be combined with groundwater irrigation data for subsequent analyses. Average irrigation depths during 2010 indicated an increase from average irrigation depths during 2008 and 2009, most likely the result of relatively dry conditions during 2010 compared to conditions in 2008 and 2009. Geostatistical models facilitated estimation of irrigation water use for unmetered systems and demonstrated usefulness in redesigning the telemetry network. Geospatial analysis evaluated the ability of the telemetry network to represent annually reported water-meter data and presented an objective, unbiased method for revising the network.

Georgia↗

Indian National Gas Hydrate Program Expedition 01 report

Gas hydrate is a naturally occurring “ice-like” combination of natural gas and water that has the potential to serve as an immense resource of natural gas from the world’s oceans and polar regions. However, gas-hydrate recovery is both a scientific and a technical challenge and much remains to be learned about the geologic, engineering, and economic factors controlling the ultimate energy resource potential of gas hydrate. The amount of natural gas contained in the world’s gas-hydrate accumulations is enormous, but these estimates are speculative and range over three orders of magnitude from about 2,800 to 8,000,000 trillion cubic meters of gas. By comparison, conventional natural gas accumulations (reserves and undiscovered, technically recoverable resources) for the world are estimated at approximately 440 trillion cubic meters. Gas recovery from gas hydrate is hindered because the gas is in a solid form and because gas hydrate commonly occurs in remote Arctic and deep marine environments. Proposed methods of gas recovery from gas hydrate generally deal with disassociating or “melting” in situ gas hydrate by heating the reservoir beyond the temperature of gas-hydrate formation, or decreasing the reservoir pressure below hydrate equilibrium. The pace of energy-related gas hydrate assessment projects has accelerated over the past several years. The Indian National Gas Hydrate Program Expedition 01 was designed to study the gas-hydrate occurrences off the Indian Peninsula and along the Andaman convergent margin with special emphasis on understanding the geologic and geochemical controls on the occurrence of gas hydrate in these two diverse settings. During Indian National Gas Hydrate Program Expedition 01, dedicated gas-hydrate coring, drilling, and downhole logging operations were conducted from 28 April 2006 to 19 August 2006.

Krishna-Godivari Basin↗

Effects of water-supply reservoirs on streamflow in Massachusetts

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

Massachusetts↗

Comparison of U.S. Geological Survey and Bureau of Reclamation water-use reporting in the Colorado River Basin

The use of water in the United States is arguably one of the most important factors determining water availability at any specific place and time. Numerous local, State, and Federal entities develop, compile, and report water-use data, which can lead to confusing or conflicting information. This report was authored jointly by the U.S. Geological Survey (USGS) and Bureau of Reclamation (Reclamation) to compare and contrast the two agencies’ water-use information programs in the Colorado River Basin. The report also describes the legal drivers for each program, clarifies confusing terminology, compares the methods used, and contrasts the information reported by each agency. This detailed comparison demonstrates that these two Federal agencies have different missions, different programmatic drivers, and different user communities, all of which lead to different approaches to water-use data collection, analysis, and reporting. This report highlights those differences and explains why the USGS and Reclamation programs exist and how the data serve different user communities. Even though the two water-use programs are different by design and purpose, the program comparison presented in this report has identified opportunities for closer coordination and sharing of information between the USGS and Reclamation, as well as program components where agency collaboration can improve water-use estimate methodologies. This comparison effort emphasizes that it is incumbent upon each agency to clearly define the meaning of the terms used and the appropriate application of the reported information to avoid confusion or the accidental misuse of the information. An additional benefit of this comparison effort is the formation of a joint USGS/Reclamation water-use team that will continue to investigate opportunities to expand and coordinate future water-use data compilation and reporting.

Arizona, Colorado, New Mexico, Utah, Wyoming↗

Sediment lithology and borehole erosion testing, American and Sacramento Rivers, California

Executive Summary A field investigation intended to measure the potential for erosion of sediments beside the American and Sacramento Rivers near Sacramento, California, is described. The study featured two primary components: (1) drilling and soil sampling to reveal lithology, down to depths matching the local river thalweg, where possible, and (2) borehole erosion tests (BETs) as described by Briaud and others (2017) at many of the same locations. The latter test involved drilling a vertical hole, measuring its diameter profile, inserting a hollow drilling rod to almost the bottom of the hole, and pumping fluid through the drilling rod at a known discharge for a chosen time interval. The hole was then resurveyed to establish an erosion rate (change in borehole radius divided by duration of flow event) as a function of depth, and the test was repeated. This test was performed with water as the erosive fluid at 12 locations, with 1 test repeated with drilling mud. Lithology holes were drilled at these same locations and an additional five locations. Drilling operations took place on river left and river right on the American River and river left (left bank, when looking downstream) on the Sacramento River. The drilling to acquire sediment samples and reveal lithology involved the use of a mobile drilling rig equipped with a 6-inch (in.) auger, a 140-pound pneumatic hammer to drive split spoon and Calmod samplers, and a piston to push Shelby tube samplers to obtain samples of clayey material. Blow count (hammer blows per 6-in. sampler advance) was recorded while sampling, and the process was logged using standard U.S. Army Corps of Engineers (USACE), Sacramento District procedures. Sediment samples were identified and described in the field per ASTM D2488 and then delivered to a USACE laboratory and to Texas A&M University for additional laboratory analysis. The BETs were performed with the same drilling rig that performed the drilling for definition of lithology. In most instances, tests were limited to regions above the water table, to avoid slumping of the borehole and heaving sands pushing into the hole. Most of the tests featured sediments that were primarily silty sand or sandy silt. The testing procedure involved comparing borehole profiles before and after passing an assumed constant discharge through a drilling rod to the bottom of the drilled hole. Discharge and water losses were logged during the testing procedure, and water losses into the walls of the drilled hole were typically less than 5 percent of the introduced volume. For the tests performed with water, the coefficient of variation of the discharge ranged from 4.5 to 28 percent, with a mean of 13 percent, but the mean discharge appeared to be reasonably steady over the typical test duration of 10–30 minutes. It was thus assumed that discharge was constant and water losses during the tests were neglected. Coefficients of variation of the discharge for the three tests performed with drilling mud were much higher (20–50 percent), but erosion rates were much smaller. Resolution of the borehole caliper-reported diameter was 0.1 in. and several of the tests lasted for 10 minutes. With boreholes measured twice, before and after each test, and averaged, these numbers correspond to an apparent erosion rate (radius change divided by test duration) of 0.3 inches per hour (in/hr), which is a theoretical lower bound on what could be measured with this approach and equipment. In practice, 0.5 in/hr appears to be a more realistic lower bound on the detectable erosion rate, based on inspection of computed changes and erosion rates. Three flow speeds (5, 8, and 12 feet per second; ft/s) were targeted for the tests. Because of equipment limitations, it was not possible in the field to reach an average of 12 ft/s throughout any given borehole, although much higher flow speeds were reached locally in some cases. Most tests featured at least two different flow rates, and the borehole was typically surveyed at least twice for each condition, to allow averaging to reduce the influence of random diameter measurement errors. Errors arising from out-of-round boreholes appeared to be uncommon. Briaud and others (2017) recommend stepped increases in the flow rate during a borehole test. This approach was taken during initial testing but proved to be problematic. The drilled hole would be enlarged by the first (smaller) discharge, and then it would be difficult to reach the desired higher flow speed because of the larger annulus between the drilled hole and the drilling rod that supplied the water for testing. This was largely solved by starting with a high discharge and, in many cases, maintaining it for subsequent tests with the average flow speed decreasing as the hole enlarged. Several different measures of erosion rate were computed and investigated by comparison to lithological profiles. The vertically averaged erosion rate for each hole was computed, but this result does not reveal vertical variability of erodibility; and the mean flow speed within the hole is not a good representation of the speed when attempting to determine a relationship between erosion rate and flow speed. Instead, for each 6-inch layer within the hole, vertically averaged erosion rates and local flow speeds were computed and plotted. Where possible, the soil type for each layer was identified. For later laboratory analysis, project protocol dictated collection of Shelby tube samples whenever clay was encountered. Plots of erosion rate versus flow speed displayed scatter that indicate that several other factors influence the erosion potential of the soil. Blow count was not a good predictor variable; it is better correlated with soil type than erodibility. Soils were classified as sand, silt, or clay, depending on which soil type dominated within a sample. In general, those classified as sand and silt did not reveal clear patterns allowing erosion rate to be computed directly from flow speed, but the test results define the range and bounds on the erosion rate. Results for clay were slightly clearer with the erosion rate increasing with flow speed, once a threshold had been reached. In this case, the erosion rate appeared to change near a speed of 7 ft/s; above this threshold, erosion rates jumped from less than 2 in/hr to greater than 3 in/hr. Even for soils with similar classifications, large differences in erodibility were observed between sites and in different layers within an individual hole. One potential means of dealing with this problem would be to perform more tests at each site to allow establishment of relationships between flow speed and erodibility for individual layers within a borehole. The maximum number of tests performed at a site in this study was four, but in some cases, results are available for only one or two flow events. Comparison of data to a set of Erosion Function Apparatus tests that provide better resolution of the vertical variation in the erosion rate versus flow speed relationship would allow further investigation of this idea. It was hypothesized that drilling mud could expand the utility of the test in soft sands by reducing the likelihood of slumping that would be interpreted as erosion. The one test that was performed with drilling mud indicated that it greatly reduced the erosion rate of the soils encountered. It yielded very different results from the test performed at the same site with water. Erosion rate is often expressed as a function of shear stress applied to a soil. In order to compute shear stress on the walls of the drilled hole, one must assume a form for the relationship between flow speed and shear stress and select a friction factor that is often estimated empirically from head loss, observed water-surface profiles, surface roughness, or other data not available in this report. One methodology for computing shear stress from flow speed is discussed in this report, but the test results have been presented in terms of erosion rate versus flow speed to avoid assuming values that are not verifiable via the field data collected in this study. Erosion rate was computed from directly measured values (sequential borehole profiles) and flow speed was computed directly from measured quantities (discharge and borehole geometry). The BET has seen limited application, primarily in clayey soils, whereas most of the soils encountered in this study were primarily sand or silt. The objective of the BET is to determine the erodibility of in situ soil below the ground or riverbed surface. The BET is simple in principle and has the advantage of revealing erodibility of in situ sediments below the ground or riverbed surface; it appears to be very useful in clayey soils, based on previously published work, but is more difficult to apply in sandy soils where slumping and water losses within the hole during testing are more likely to occur. The BET did reveal a large variation in the results both laterally and vertically, even for the same soil-type classification. It is thus recommended that the results be applied considering these spatial variations rather than attempting to universally assign an erosion-rate relationship to a particular soil type. Results have been provided showing the results by site and by sediment classification (sand, silt, and clay), to allow either approach. Where possible, it is important to rely on site-specific results because the erosion-rate relationship for a given soil type varied by site. Data collected during this project have been made publicly available online via the U.S. Geological Survey (USGS) Sciencebase database. The measured borehole profiles, discharge, lithology log sheets, and photos are available in the data release that accompanies this report (see Work and Livsey (2019) in the “Selected References” section for the appropriate link).

California↗

Streamflow characterization and effects of Interstate 26 highway construction on water-quality and channel morphology conditions of the French Broad River in Buncombe and Henderson Counties, North Carolina, 2019–22

The North Carolina Department of Transportation began a series of highway construction projects to upgrade Interstate 26 (I–26) in Buncombe and Henderson Counties, North Carolina, starting in early 2020. The North Carolina Department of Transportation and the U.S. Geological Survey partnered in 2019 to monitor and assess water-quality and geomorphic conditions throughout the French Broad River corridor in areas near Interstate 26 highway construction activities. This report summarizes findings from multi-year water-quality and channel morphology monitoring from 2019 through 2022 in the French Broad River corridor near ongoing highway construction activities. Dissolved oxygen, pH, specific conductance, turbidity, and water temperature were monitored at six gaging locations, three of which were near highway sections undergoing construction during the study period from June 1, 2019, through December 31, 2022. Water quality was compared across three monitoring locations before and during construction using Brunner-Dette-Munk statistical tests, a nonparametric form of analysis of variance used for unbalanced datasets. Most Brunner-Dette-Munk test results detected insignificant effects of highway construction on water-quality conditions when controlling for monitoring location, except for the test of daily maximum specific conductance values. However, the effects of highway construction could not be isolated from other external factors that may have affected specific conductance, including an extended period of low streamflow and possible road salt application after winter storms. Additional Brunner-Dette-Munk tests, controlling for flow magnitude (streamflow quantiles) or hydrograph position (“rising limb,” “falling limb,” or “base flow”), did not suggest that highway construction affected water-quality conditions in the study area. Geomorphic surveys of the streambed and streambanks were performed in a 1-kilometer reach from 2019 through 2022 to assess geomorphic changes during the reconstruction of a bridge spanning the French Broad River. Bathymetric surveys of the streambed and light detection and ranging (lidar) scans of the streambanks were collected from a canoe using motion-compensated, real-time positioning solutions. Initial geomorphic surveys were used as a baseline to which subsequent repeat surveys were compared using geomorphic change detection analysis. Detectable elevation changes were determined using statistically robust thresholds, ensuring only changes exceeding measurement uncertainty were considered for interpretive analyses. Change detection results identified that 16.4 percent of the reach experienced detectable geomorphic change during the study period, with inter-survey changes ranging from 3.5 to 12.4 percent of the total reach. Most detectable change is attributed to the installation and reconfiguration of temporary construction pads; other changes are attributed to enlargement of a scour hole and upstream deposition near the bridge. Although anthropogenic activities and changes in streamflow regime contributed to these observations, the approach used in this study provides confidence that reported values reflect real geomorphic changes and not noise.

North Carolina↗

Preliminary volcano-hazard assessment for the Tanaga volcanic cluster, Tanaga Island, Alaska

Summary of Volcano Hazards at Tanaga Volcanic Cluster The Tanaga volcanic cluster lies on the northwest part of Tanaga Island, about 100 kilometers west of Adak, Alaska, and 2,025 kilometers southwest of Anchorage, Alaska. The cluster consists of three volcanoes-from west to east, they are Sajaka, Tanaga, and Takawangha. All three volcanoes have erupted in the last 1,000 years, producing lava flows and tephra (ash) deposits. A much less frequent, but potentially more hazardous phenomenon, is volcanic edifice collapse into the sea, which likely happens only on a timescale of every few thousands of years, at most. Parts of the volcanic bedrock near Takawangha have been altered by hydrothermal activity and are prone to slope failure, but such events only present a local hazard. Given the volcanic cluster's remote location, the primary hazard from the Tanaga volcanoes is airborne ash that could affect aircraft. In this report, we summarize the major volcanic hazards associated with the Tanaga volcanic cluster.

Alaska↗

Modeling uncertainty in coal resource assessments, with an application to a central area of the Gillette coal field, Wyoming

Standards for the public disclosure of mineral resources and reserves do not require the use of any specific methodology when it comes to estimating the reliability of the resources. Unbeknownst to most intended recipients of resource appraisals, such freedom commonly results in subjective opinions or estimations based on suboptimal approaches, such as use of distance methods. This report presents the results of a study of the third of three coal deposits in which drilling density has been increased one order of magnitude in three stages. Applying geostatistical simulation, the densest dataset was used to check the results obtained by modeling the sparser drillings. We have come up with two summary displays of results based on the same simulations, which individually and combined provide a better assessment of uncertainty than traditional qualitative resource classifications: (a) a display of cell 90 percent confidence interval versus cumulative cell tonnage, and (b) a histogram of total resources. The first graph allows classification of data into any number of bins with dividers to be decided by the assessor on the basis of a discriminating variable that is statistically accepted as a measure of uncertainty, thereby improving the quality and flexibility of the modeling. The second display expands the scope of the modeling by providing a quantitative measure of uncertainty for total tonnage, which is a fundamental concern for stockholders, geologists, and decision makers. Our approach allows us to correctly model uncertainty issues not possible to predict with distance methods, such as (a) different levels of uncertainty for individual beds with the same pattern and density of drill holes, (b) different local degrees of reduction of uncertainty with drilling densification reflecting fluctuation in the complexity of the geology, (c) average reduction in uncertainty at a disproportionately lesser rate than the reduction in area per drill hole, (d) the proportional effect of higher uncertainty in areas of higher tonnages, despite a regular drilling pattern, (e) the possibility of a local increase in uncertainty despite drilling densification to reflect a more complex geology as the deposit is known in more detail, and (f) for exactly the same drilling pattern, tonnage per individual beds with different uncertainty than the aggregated tonnage. These results should be considered realistic improvements over distance methods used for quantitative classification of uncertainty in coal resource, such as U.S. Geological Survey Circular 891. 1 The approach should be a welcome addition to the toolkit of Competent Persons preparing public disclosures according to international mineral codes such as those promoted by the Combined Reserves International Reporting Standards Committee 2 and the Joint Ore Reserve Committee. 3 1 Wood, G.H., Jr., Kehn, T.M., Carter, M.D., and Culbertson, W.C., 1983, Coal resources classification system of the U.S. Geological Survey: U.S. Geological Survey Circular 891, 65 p. 2 CRIRSCO (Combined Reserves International Reporting Standards Committee), 2013, International reporting template for the public reporting of exploration results, mineral resources and mineral reserves: Accessed February 2014 at http://www.crirsco.com/crirsco_template_may2013.pdf. 3 JORC (Joint Ore Reserves Committee), 2012, Australasian code for reporting of exploration results, mineral resources and ore reserves: Accessed September 2014 at http://www.jorc.org/docs/jorc_code2012.pdf.

Wyoming↗