Search USGS⌕ Search

SEARCH · Search USGS

Results for “Scientific Monthly”

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 703 records · Page 39Linked to original sources
USGS2011

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

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

Scientific Investigations Report↗
USGS2024

Reservoir evolution, downstream sediment transport, downstream channel change, and synthesis of geomorphic responses of Fall Creek and Middle Fork Willamette River to water years 2012–18 streambed drawdowns at Fall Creek Lake, Oregon

Executive Summary Chapter A. Introduction Fall Creek Dam impounds Fall Creek Lake, a 10-kilometer-long reservoir in western Oregon and is operated by the U.S. Army Corps of Engineers (USACE) primarily for flood-risk management (or flood control) in late autumn through early spring months, as well as for water quality, irrigation, recreation, and habitat in late spring through early autumn. Since 2011 (water year [WY] 2012), Fall Creek Lake has been temporarily drawn down each year to facilitate downstream passage of juvenile spring Chinook salmon ( Oncorhynchus tshawytscha ) through the 55-meter (m) high dam. This annual dam operation is temporary, typically lasting about 1–2 weeks from WY 2012 through 2020 (drawdown operations in WY 2022–24 have increased to more than 6 weeks). Drawdown of the reservoir results in lake levels being lowered to the elevation near the historical, pre-dam streambed. The annual streambed drawdowns of WY 2012–18 have improved fish passage and led the USACE to formally adopt streambed drawdowns as part of annual operations at Fall Creek Dam. However, temporarily lowering the lake to streambed creates free-flowing conditions in the reservoir that result in the erosion and episodic export of predominantly sand and finer-grained sediments (less than 2 millimeters [mm]) to the lower gravel-bed reaches of Fall Creek and the Middle Fork Willamette River. The introduction of large volumes of sand and finer-grain sediment into the dam-regulated reaches downstream from Fall Creek Dam prompted questions about the geomorphic responses to annual streambed drawdowns within Fall Creek Lake and downstream reaches along Fall Creek and the Middle Fork Willamette River. The U.S. Geological Survey (USGS) in partnership with USACE initiated a comprehensive geomorphic and sediment transport investigation to assess the coupled processes of reservoir erosion, sediment evacuation from Fall Creek Lake, and patterns of sediment transport and deposition in reaches downstream from the Fall Creek Dam that have resulted from annual streambed drawdowns. The purpose of this report is to systematically describe the processes of sediment erosion, transport, and deposition at Fall Creek Lake and geomorphic interactions between reaches upstream and downstream from Fall Creek Dam that relate to dam operations. Specifically, this report focuses on evaluating geomorphic responses to streambed drawdowns from WY 2012 through 2018 and placing drawdown-induced geomorphic responses within the broader context of physiographic and historical conditions and dam operations of Fall Creek and Middle Fork Willamette Rivers. Key objectives for this study were to characterize changes in reservoir morphology and substrate at Fall Creek Lake, describe the character and temporal pattern of sediment transport downstream from Fall Creek Dam, characterize geomorphic changes in channel reaches downstream from the Fall Creek Dam, and relate these data to the annual streambed drawdowns of WY 2012–18. This study uses multiple independent monitoring and measurement approaches to assess site, reach, and river-scale geomorphic responses to drawdowns to inform dam and reservoir management. Patterns and processes of reservoir evolution were assessed with geomorphic mapping and volumetric analyses of topography through comparison of multiple digital surface models (DSMs). Just downstream from Fall Creek Dam, analyses of sediment export from the reservoir focused on suspended sediment but also incorporated bedload analyses to assess sediment sizes. Geomorphic assessments downstream from the dam used reach-scale and site-scale approaches to document changes in channel morphology and substrate, including site measurements of sand and finer-grained sediment deposition and in-channel bed-material, volumetric change analyses from comparison of digital elevation models (DEMs), and repeat geomorphic mapping. Findings from this study inform river management and dam operations by providing an understanding of (1) coupled upstream-downstream geomorphic responses to the Fall Creek Lake streambed drawdowns, (2) geomorphic responses of Fall Creek Lake streambed drawdowns in comparison to drawdowns at other large dams, (3) controls on reservoir erosion and downstream geomorphic responses, and (4) implications for future hydrogeomorphic changes that may result from continued drawdowns and monitoring activities to assess those changes. Chapter B. Reservoir Morphology and Evolution Related to Dam Operations at Fall Creek Lake To understand the volume and distribution of sediment accumulation in Fall Creek Lake since dam closure in 1965, decadal-scale sedimentation patterns (spanning approximately 1965–2016) are evaluated using a combination of storage curve analyses and geomorphic mapping. Short-term (drawdown event-scale) patterns of erosion, sedimentation, and sediment export downstream are evaluated using a combination of geomorphic mapping and change detection analyses that quantify the distribution and total volume of sediment erosion and deposition within Fall Creek Lake. Geomorphic mapping of reservoir topography and analyses of historical datasets reveals four categories of landforms and sediment processes within Fall Creek Lake related to lake level operations: lacustrine sedimentation expressed in the reservoir floor, fluvial erosion and deposition within historical stream channels during streambed drawdowns, channel-like features created by erosion within the reservoir floor during streambed drawdowns, and erosion on reservoir hillslopes. Where the reservoir floor is mapped for this study as pelagic (deep water), deposition up to 3 meters (m) thick by lacustrine processes and burial of pre-dam topography with deposits thinning toward the edges of the valley floor and upstream areas of reservoir are observed. Despite over 50 years of sediment accumulation since dam construction, the main stream channels of Fall and Winberry Creeks (or reservoir thalwegs) through the reservoir are well defined, though their distinct morphology is likely influenced by a long history of recurring historical drawdowns to or near streambed since dam construction. Unregulated streamflow and sediment transport through the reservoir primarily are confined to these channels during the streambed drawdown periods. Erosional channel-like features created by drawdowns are carved through underlying, unconsolidated reservoir floor sediments and are most prominent in the lower reservoir below minimum conservation pool (the low pool elevation during winter flood season); sediment generated from the formation of these drawdown channels is more likely to be transported through and out of the reservoir than sediment deposits along the reservoir hillslopes at the valley margins that are separated from main channels by areas of low-gradient reservoir floor. Morphologic changes in the lower reservoir topography between January 2012 and November 2016 indicate overall net erosion of about 129,500 cubic meters (m 3 ). The most prominent geomorphic changes occurred along the main channels of Fall and Winberry Creeks near the Fall Creek Dam where incision, lateral migration, and slumping banks resulted in vertical and lateral adjustments to channel position, whereas most changes fell below the detectable limit on higher-elevation reservoir floor surfaces except where erosion occurred along features mapped as drawdown channels. Chapter C. Sediment Delivery from Fall Creek Lake and Transport through Downstream Reaches USGS implemented a sediment monitoring program in WY 2013–18 to evaluate the quantity and character of reservoir sediment exported from Fall Creek Lake during streambed drawdowns. Turbidity and suspended sediments were monitored annually autumn through spring to span the WY 2013–18 streambed drawdowns; however, unequal monitoring timeframes each year reduced the ability to compare results and factors affecting sediment export from the reservoir difficult between years. These data were originally measured to develop regressions and compute suspended-sediment loads (SSL). Bedload sediment monitoring from a cableway at the Fall Creek streamgage was completed in the autumn-winter of WY 2013 and 2017. The limited number of samples and presumed variability in sediment supply from the reservoir precluded construction of streamflow and bedload discharge relations to compute more than instantaneous bedload. Sand and finer-grained silts and clays were transported from the reservoir in suspension, though some coarser grains (up to 32 mm) were also mobilized and transported downstream from the dam as bedload. Observations of increased sediment transport downstream from Fall Creek Dam coincided with lake levels approaching about 3 m (10 feet [ft] or elevation 690 ft) above the streambed regulating outlets. Suspended-sediment loads computed for the full monitoring periods WY 2013–18 at the Fall Creek streamgage, located 1.4 kilometers (km) downstream from Fall Creek Dam, range from 54,700 metric tons (t) in WY 2013 to 13,900 t in WY 2018. Although the total annual SSL varied from year to year, the overall seasonal patterns of suspended sediment transport throughout each year were similar during monitoring in WY 2013-18. Suspended-sediment loads were low prior to the drawdown, then increased rapidly as lake levels lowered and approached the streambed. In the weeks following the drawdown period, as pool levels were increased, SSL remained slightly elevated above pre-drawdown levels but generally declined through the following winter and spring except during streamflow-driven pulses of suspended-sediment transport. WY 2013 had the greatest total computed SSL for each streambed drawdown and partial-year monitoring period. SSL computed for the partial-year period have generally decreased since WY 2013 and have varied by about 6,800 t with the exception of WY 2014. WY 2014 SSL reflects anomalously low sediment export due to low streamflows and freezing conditions that stabilized reservoir floor deposits. Bedload measurements in the short 1.4-km reach between Fall Creek Dam and the Fall Creek streamgage showed an inverse correlation between bedload transport rates and discharge, which probably reflects diminishing supply of coarse-sized sediment. Sand was more abundant (60–100 percent) than gravel in bedload samples confirming sand and finer-grained sediment dominated sediment evacuated from the reservoir during streambed drawdowns at Fall Creek Lake. Chapter D. Geomorphic Responses to Fall Creek Lake Streambed Drawdowns Downstream from Fall Creek Dam In the days, weeks, and months following streambed drawdown operations at Fall Creek Dam through WY 2018, sites downstream from the dam displayed a variety of geomorphic responses to reservoir sediment delivery within the main channel and overbank areas. Evaluation of streambed elevations at two streamgages located 1.4 km downstream from the dam on Fall Creek and 16.3 km downstream from the dam on the Middle Fork Willamette River indicated the effects of drawdown sediment on bed elevations were modest and transient. Repeat particle size measurements (October 2015 and September 2016) at five sites along Fall Creek and the Middle Fork Willamette River showed similar grain-sized distributions that do not reveal substantial deposition of fine-grained sediment related to the WY 2016 streambed drawdown. Altogether, these findings indicate that transport capacity in the main, low-flow channels of Fall Creek and Middle Fork Willamette River during WY 2012–18 was sufficient to mobilize and evacuate reservoir sediments from streambed drawdowns or other bank material and tributary sources. However, other monitoring for this study indicate low-velocity zones in off-channel areas are prime locations for sand and finer-grain sediment deposition. Patterns of overbank sediment accumulation indicate that the magnitude and timing of overbank deposition on bars and low-elevation floodplain varies with proximity to the dam, geomorphic setting, streamflows, and other factors. Sand and finer-grained reservoir sediments carried as suspended-sediment load in the reaches downstream from Fall Creek Dam were deposited in overbank areas as observed with clay-horizon markers during WY 2016–17. Overbank deposition quantified with Geomorphic Change Detection (GCD) software evaluated landform-scale patterns of erosion and deposition using repeat light detection and ranging (lidar) surveys at two sites in the Upper Fall Creek reach and one site in the Jasper reach for 3 years (2012–15) and one site in the Clearwater reach for 6 years (2009–15). Deposition thickness and spatial patterns from the GCD analysis were variable; some sites had dispersed but measurable deposition while at others, deposition was highly localized and exceeded 1 m in depth. Patterns of overbank deposition illustrate interactions among bar morphology, local hydraulics, and suspended-sediment transport dynamics that can create patches of highly localized deposition. The measured deposition at the two Fall Creek GCD sites likely resulted from reservoir sediments released from Fall Creek Lake during streambed drawdowns in WY 2016 and 2017 because the limited sediment inputs from bank material (geomorphically laterally stable reach) or tributaries (no significant tributaries) provided few other sediment sources. On the Middle Fork Willamette River, observed patterns of overbank deposition could reflect sediment sourced from upstream tributaries, bank erosion, or Fall Creek Lake streambed drawdown operations. Despite the introduction of several thousand tons of reservoir sediment delivered from the Fall Creek Lake streambed drawdowns to below-dam river corridors, reach-scale mapping of channel features downstream from Fall Creek Dam shows minimal evidence of changes in channel planform or landforms that can be attributed to a drawdowns in WY 2012–16. On Upper Fall Creek reach, widespread increases in gravel bars or other in-channel sediment did not result from the five streambed drawdowns. The main changes attributable to sediment releases from Fall Creek Lake were localized increases in vegetated bar area, particularly on channel margin areas where sand and finer-grain sediment was deposited and rapidly colonized by vegetation. The area of mapped secondary water features decreased between 2005 and 2016, but that may be due to lower discharges depicted in the 2016 aerial photographs and less mapped area of inundation. Primary changes along the Lower Fall Creek reach include a 6.4 percent decrease in area of secondary water features between 2011 and 2016, and a nearly twofold increase in the area of unvegetated bars. Immediately downstream from the Fall Creek confluence, there were negligible changes in the location and areas of vegetated bars and the main wetted channel between 2005 and 2016, and local increases in bar area cannot be attributed solely to deposition of reservoir sediments from Fall Creek Lake because (1) areas along the Middle Fork Willamette River just upstream from the Fall Creek confluence display similar type and magnitude of changes and (2) some of the increases at the confluence area pre-date the drawdowns. The cumulative effect of sediment releases from Fall Creek Lake streambed drawdowns from WY 2012 to 2016 on downstream channel planform and landforms are modest compared to the river-scale transformations and planform changes that occurred in the decades following dam construction. Chapter E. Discussion of Geomorphic Responses of Fall Creek and Middle Fork Willamette River to Streambed Drawdowns at Fall Creek Lake Multiple aspects of Fall Creek Dam infrastructure and operations exert first-order controls on the magnitudes of reservoir erosion that occur during the streambed drawdowns and ultimately determine the sediment delivery to downstream reaches. Key aspects of the dam and its operations that are most relevant to assessing geomorphic responses to streambed drawdowns include the (1) dam infrastructure, including configuration and size of regulating outlets and their proximity to the streambed which dictates the capacity and competence of the river to deliver sediment to downstream reaches and mode of sediment transport as suspended-sediment load or bedload; (2) frequency of historical drawdowns and long-term, year-round dam operations and lake level management, which partly dictate reservoir morphology and locations and magnitudes of readily erodible materials; (3) dam operations and hydroclimatic conditions during the streambed drawdown (including length of the drawdown and streamflows entering the reservoir), which directly control the timing, duration and magnitude of reservoir erosion and sediment evacuation; and (4) dam operations following the streambed drawdown operation that regulate streamflows (and thereby sediment transport conditions) downstream of Fall Creek Dam which primarily reflect interactions between hydroclimatic conditions and flood control operations. Patterns of sediment erosion and evacuation observed in this study at Fall Creek Lake from WY 2012–18 suggest that reservoir erosion during annual streambed drawdowns may remain similar or decrease in future years assuming (1) annual streambed drawdown operations are implemented in similar manner as the WY 2012–18 drawdowns (in terms of duration, late autumn or early winter implementation, rate of pool-level lowering to reach streambed, and other factors), (2) streambed drawdowns coincide with similar conditions as were observed WY 2012–18 (similar sediment yield into reservoir, low reservoir inflows, limited precipitation, moderate air temperature), and (3) no major geomorphic changes in the main reservoir channels of Fall and Winberry Creeks occur (for example, channel avulsion). Under such conditions, it is hypothesized that the stream channel within the reservoir would achieve a quasi-equilibrium state with respect to annual influx and export of sediment and aided by the substantial amount of in-channel bedrock, will remain laterally stable without erosion across reservoir deposits. Patterns of sediment transport measured at the Fall Creek streamgage downstream from Fall Creek Dam provide insight into the potential effects of future streambed drawdowns at Fall Creek Lake. Analyses of suspended sediment measured in WY 2013–18 show a major reduction in suspended-sediment loads between WY 2013 and later years, indicating streamflows transporting sediment through the reservoir to downstream reaches during streambed drawdowns have become supply limited. The 6-year suspended-sediment monitoring and sampling program is insufficient to make predictions about future sediment transport conditions because of uneven monitoring periods and varying controls on reservoir sediment erosion. It is likely that future suspended-sediment loads will be variable but similar to those observed in WY 2015–18 if operational, climatic, and geomorphological factors remain similar to those monitored WY 2015–18. Suspended-sediment loads downstream from Fall Creek Lake will likely remain highest when regulating outlets are fully open and Fall Creek is free flowing with the reservoir fully drained with little to no residual pool. Over time, it is possible that the suspended-sediment loads would reflect mobilization of reservoir sediment deposited in the previous year rather than erosion of sediment deposited years or decades earlier. Bedload is likely to remain a small fraction of the total sediment load evacuated from the reservoir and is relatively modest compared with pre-dam bedload transport rates because most coarse sediment remains trapped by the dam. If sediment releases from Fall Creek Lake and ensuing streamflow conditions follow a similar pattern in the future as was assessed in this study spanning WY 2012–18, near-term geomorphic adjustments downstream of the dam are expected to be modest. Barring major operational, climatic, and geomorphological changes, local site-scale deposition on bars, overbank areas, or off-channel features that persists several months after the streambed drawdown will likely continue to be highly variable, ranging from negligible to several centimeters of deposition. At the landform-scale, low velocity areas nearest to Fall Creek Dam will likely continue to undergo rapid deposition immediately during and after a streambed drawdown event, similar to patterns observed for WY 2012–18. Some of the sediment entering these off-channel features and margin areas may be temporarily stored, then later remobilized and dispersed farther downstream. But if newly deposited sediment persists through the following spring, there is a greater likelihood that local vegetation will establish, reinforce deposited material, and trap sediment during later drawdowns. The reach-scale geomorphic changes may become more apparent if (1) streambed drawdowns continued for several decades, and geomorphic changes were measured at decadal scales or (2) the amount of sediment introduced to downstream reaches substantially increased and (or) sediment transport capacity decreased. The continued streamflow regulation of Fall Creek Dam after sediment releases provides an opportunity to strategically manage streamflows during and after the streambed drawdowns to minimize downstream sediment impacts and ensure other operational thresholds are satisfied. This study provides a comprehensive foundation of datasets and geomorphic analyses to inform dam operations at Fall Creek Lake, monitor sediment transport downstream, and consider operational schemes for future drawdowns. The datasets from this study also provide baselines of sediment transport and geomorphic conditions to assess future changes in reservoir and downstream environments. Future monitoring could be tailored to address specific questions regarding the long-term geomorphic effects of streambed drawdowns on fluvial habitats, flood hazards, cultural resources, or downstream water quality. Future monitoring activities could focus on the relevant geomorphic processes and spatial domains within the three categories used for this study: (1) reservoir erosion and net sediment evacuation, (2) sediment delivery to downstream reaches, including magnitude and temporal pattern of sediment transport, and (3) geomorphic responses of downstream reaches to sediment delivery. Specifically, high priority future monitoring activities could include: Repeat topographic or photographic surveys in the reservoir to characterize changes occurring within individual drawdowns, to quantify sediment export, to determine temporal changes in reservoir storage, and to identify locations of erosion and deposition. Continuous, year-round turbidity monitoring supplemented with suspended-sediment measurements at a streamflow-gaging station immediately downstream from the dam to quantify sediment export. Repeat geomorphic monitoring, mapping, or modeling in downstream reaches to track changes in channel and over bank features using a combination of site- and reach-scale monitoring approaches. This could support assessments of sediment deposition and ensuing vegetation encroachment on flood hazards and habitats and examine how sediment transport and depositional processes may be affected by different sediment supply, streamflow, or dam management scenarios.

Oregon↗
USGS2010

Determination of time-of-travel, dispersion characteristics, and oxygen reaeration coefficients during low streamflows--Lower Tacony/Frankford Creek, Philadelphia, Pennsylvania

Time-of-travel, dispersion characteristics, and oxygen reaeration coefficients were determined by use of dye and gas tracing for a 2-mile reach of Tacony/Frankford Creek in Philadelphia, southeastern Pennsylvania. The reach frequently has concentrations of dissolved oxygen (DO) below the water-quality standard of 4 milligrams per liter during warm months. Several large combined sewer overflows (CSOs), including one of the largest in Philadelphia (former Wingohocking Creek), discharge to the study reach in this urbanized watershed, affecting water quality and the timing and magnitude of storm peaks. In addition, a dam that commonly results in backwater conditions and reduced natural reaeration is present a few hundred feet from the end of the study reach. Time-of-travel and reaeration data were collected under base-flow conditions in August and September 2009 for three sub-reaches from Roosevelt Boulevard (U.S. Route 1) to Castor Avenue. Determination of traveltimes to the centroid of the dye cloud were needed for calculation of the reaeration coefficients. Results of the dye study in Tacony/Frankford Creek indicate that traveltimes were affected by the presence of man-made structures, such as the large scour hole and pool developed at the outfall of the T14 CSO and the dam, both of which reduce stream velocities. Mean stream velocities during the dye-tracer tests ranged from a maximum of 0.44 to 0.04 foot per second through a large pool. The dispersion efficiency of the stream was determined from relations between normalized unit concentrations to time to peak for use in water-quality modeling. Oxygen reaeration coefficients determined by a constant rate-injection method using propane as the tracer gas were as low as 0.04 unit per hour in a long pool affected by backwater conditions behind a dam. The highest reaeration coefficient was 2.29 units per hour for a steep-gradient reach with multiple winding channels through gravel deposits, just downstream of a large scour pool developed at the outlet of the T14 CSO. Reaeration coefficients determined from the field tracer-gas method were compared to values calculated by two other methods, one that is based on theoretical equations using physical properties of the stream as variables and the other that is based on equations using the timing of measured daily maximum DO concentrations in the stream. Reaeration coefficients from the two alternate methods were most similar to values determined from the field tracer-gas method for the upstream portion of the study reach, characterized by free-flowing riffle and pools. Values of reaeration coefficients determined by the tracer-gas method were 2 to 10 times higher than values determined by 2 alternate methods for most subreaches hydraulically affected by man-made structures. In addition to the tracer gas, propane, the gas analysis also included methane, ethane, and ethene, of which only methane was measured in concentrations above a few micrograms per liter. Methane, thought to occur naturally or because of ongoing processes in the stream, was measured in concentrations ranging from 6.6 to 78 micrograms per liter; the concentrations were greatest in sub-reaches dominated by pools.

Pennsylvania↗
USGS2022

Tracking heat in the Willamette River system, Oregon

The Willamette River Basin in northwestern Oregon is home to several cold-water fish species whose habitat has been altered by the Willamette Valley Project, a system of 13 dams and reservoirs operated by the U.S. Army Corps of Engineers. Water-resource managers use a variety of flow- and temperature-management strategies to ameliorate the effects of upstream Willamette Valley Project dams on the habitat and viability of these anadromous and native fish. In this study, new capabilities were added to the CE-QUAL-W2 two-dimensional flow and water-quality model to inform those flow- and temperature-management strategies by tracking the quantities and ages of water and heat from individual upstream sources to downstream locations in the Willamette River system. Specifically, the fraction of water and heat attributable to upstream dam releases or other water inputs, and the fraction of heat sourced from environmental heat fluxes across the water and sediment surfaces, were tracked and quantified in the river at all locations and times simulated by the model. Applying the updated CE-QUAL-W2 models to the Willamette River system for the months of March through October in the years 2011 (cool and wet), 2015 (hot and dry), and 2016 (warm and somewhat dry) demonstrated that the influence of upstream dam releases on downstream water temperature diminished within a few days as water moved downstream. At sites that are roughly two or more days of travel from upstream dams (Albany and downstream), the July–August fraction of riverine heat content that could be tracked back to upstream dam releases typically diminished to less than 20 percent, despite the fact that roughly 50 percent of July–August streamflow could be attributed to upstream dam releases at the same sites. In contrast, the fraction of riverine heat content that could be attributed to environmental energy fluxes continued to increase with downstream distance, from about 59 to 67 percent at Albany during July–August to 62 to 73 percent at Keizer and 68 to 79 percent at Newberg. At locations sufficiently far downstream, upstream dam releases affect water temperature mainly through a decrease in travel time (less time for environmental heat fluxes to warm the river during summer) and an increase in thermal mass (more water to dilute and buffer incoming heat fluxes) rather than through the simple transport of heat content (water temperature) released from the dams. This concept was explored not only for the baseline conditions that occurred in March–October of 2011, 2015, and 2016, but also for a hypothetical high-flow release during August 2016 and an actual high-flow release during August 2017. In these high-flow releases, an extra 2,500 cubic feet per second (roughly) was released from Dexter Dam on the Middle Fork Willamette River, and downstream effects were measured (2017, actual) and simulated (2016, hypothetical). Results of the simulations were consistent with insights gained from the baseline conditions, such that temperature changes caused by flow augmentation were substantial in upstream reaches (measured cooling of about 1.5 °C near Harrisburg [43 miles downstream] and Albany [84 miles downstream] in 2017, and cooling of about 0.5 °C near Albany in 2016) and diminished farther downstream, but still measurable (more than a few tenths of a degree Celsius) even at Newberg, which is about 154 miles downstream. The direct downstream effects of dam releases on the river heat content attributable to those releases were increased by the hypothetical flow augmentation, with increases of 20 percent at Harrisburg and 12 percent at Keizer. Even with a decreased influence of environmental energy fluxes on river heat content, however, the fraction of heat content attributable to such fluxes was still more than 50 percent at and downstream of Albany and more than 70 percent at Newberg, where the river temperature was less affected by upstream dam-release temperatures and instead was affected primarily by a decreased travel time and increased thermal mass.

Oregon↗
USGS2020

Hydrogeology and groundwater flow in alluvial deposits, north Summerset, South Dakota

The city of Summerset is a growing community in west South Dakota. The Sun Valley Estates subdivision in the north part of the city was developed on unconsolidated deposits surrounded by steep terrain. During years with greater than normal precipitation, particularly in 2019, groundwater levels increased in the unconsolidated deposits and caused damage to stormwater systems, sewer infrastructure, and houses with basements. The U.S. Geological Survey, in cooperation with the City of Summerset, completed a study of the hydrogeology and groundwater flow in the alluvial aquifer part of the unconsolidated deposits in north Summerset to understand the groundwater system in the area and to provide hydrogeologic information in support of future development planning. The study area included most of the Sun Valley Estates subdivision in the north part of the city of Summerset in the east Black Hills of west South Dakota. About 0.7 square mile of water-bearing alluvial deposits is included in the study area. Precipitation in the study area from 2017 to 2019 was compared to the monthly normal values at a nearby climate site. The largest departure from normal was in May 2019 with precipitation exceeding the monthly normal by about 5 inches (in.). All months in 2019, except March, exceeded the monthly normal precipitation. Cumulative departure from normal precipitation in 2019 increased from about 4 in. greater than normal in January to about 18 in. greater than normal in December. The geologic setting of the study area is characterized by the surrounding Black Hills. Unconsolidated Quaternary-age deposits overlie consolidated to partially consolidated Mesozoic-age and Paleozoic-age shales, sandstones, and limestones. Surficial deposits of alluvium and other unconsolidated deposits are the primary surficial geologic units in the study area and form the components of the alluvium hydrogeologic unit of the study area. Results from previous studies of alluvium along nearby Rapid Creek estimated hydraulic conductivity to range from 89 to 2,292 feet per day (ft/d), transmissivity to range from 1,001 to 32,083 feet squared per day, and storage coefficients to range from 0.0002 to 0.16. Hydraulic conductivity and transmissivity generally decreased downstream along Rapid Creek (west to east). Slug tests were completed August 16, 2019, at two observation wells completed in the alluvial aquifer in the Sun Valley Estates subdivision to determine hydraulic conductivity. Hydraulic conductivity estimated from AQTESOLV curve-fitting analysis using the Bouwer-Rice method for all slug-in and slug-out trials from two observation wells in the study ranged from 0.20 to 0.26 ft/d for well 441318103220001 (SunValley1 well) and from 0.54 to 14 ft/d for well 441319103215701 (SunValley2 well). The mean, median, and standard deviation of all trials at both wells were 4.3 ft/d, 0.8 ft/d, and 5.6 ft/d, respectively. The extent of the alluvial aquifer was determined by geologic maps and lithologic logs. Alluvial deposits in the study area extend to about 1 mile in the north–south direction and about 1.5 miles in the southeast–northwest direction. The direction of groundwater flow was estimated using water-level records and topographic maps. The resulting potentiometric map indicated that groundwater in the alluvial aquifer under the Sun Valley Estates subdivision originates from higher elevations of the west part of the area of interest and from streams in the southeast part. Recharge and evapotranspiration estimates were results from a Soil-Water Balance model that calculated a matrix of recharge for 2019 with values ranging from 0 to 11.4 in. and an annual mean value of 5.1 in. across the study area. Soil-Water Balance-estimated potential evapotranspiration for 2019 ranged from 28.90 to 28.75 in. and the estimated annual mean was 28.86 in. across the study area. Estimated groundwater budget components for the alluvial aquifer in the area of interest included inflows and outflows. Total estimated groundwater budget components for inflows for 2019 were about 66 percent from recharge, 33 percent from streamflow, and 1 percent from inflow from adjacent aquifers. Total estimated outflows were about 99-percent evapotranspiration and less than 1-percent outflow to adjacent aquifers.

South Dakota↗
USGS2014

Turbidity and suspended sediment in the upper Esopus Creek watershed, Ulster County, New York

Suspended-sediment concentrations (SSCs) and turbidity were measured for 2 to 3 years at 14 monitoring sites throughout the upper Esopus Creek watershed in the Catskill Mountains of New York State. The upper Esopus Creek watershed is part of the New York City water-supply system that supplies water to more than 9 million people every day. Turbidity, caused primarily by high concentrations of inorganic suspended particles, is a potential water-quality concern because it colors the water and can reduce the effectiveness of drinking-water disinfection. The purposes of this study were to quantify concentrations of suspended sediment and turbidity levels, to estimate suspended-sediment loads within the upper Esopus Creek watershed, and to investigate the relations between SSC and turbidity. Samples were collected at four locations along the main channel of Esopus Creek and at all of the principal tributaries. Samples were collected monthly and during storms and were analyzed for SSC and turbidity in the laboratory. Turbidity was also measured every 15 minutes at six of the sampling stations with in situ turbidity probes. The largest tributary, Stony Clove Creek, consistently produced higher SSCs and turbidity than any of the other Esopus Creek tributaries. The rest of the tributaries fell into two groups: those that produced moderate SSCs and turbidity and those that produced low SSCs and turbidity. Within those two groups the tributary that produced the highest SSCs and turbidity varied from year to year depending on the hydrologic conditions within each subwatershed. During the 3-year study, Stony Clove Creek accounted for an average of 40 percent of the annual suspended-sediment load measured at the upper Esopus Creek watershed outlet at Coldbrook, more than all of the other measured tributaries combined. The other tributaries to the upper Esopus Creek, taken together, accounted for an average of about 20 percent of the load at Coldbrook during 2010 and 2011, when most of the tributaries were sampled. Woodland Creek, the third largest tributary in the watershed, also accounted for a substantial amount of the load at Coldbrook, an average of 10 percent during the 3 years. Stony Clove Creek appeared to be a persistent source of sediment to Esopus Creek; it had the highest sediment yield (load per unit area) of all monitoring sites, including the outlet at Coldbrook. Discharge, SSC, and turbidity were strongly related at the Coldbrook site but not at every monitoring site. In general, relations between discharge and SSC and turbidity were strongest at sites with high SSCs, with the exception of Stony Clove Creek. Stony Clove Creek had high SSCs and turbidity regardless of discharge, and although concentrations and turbidity values generally increased with increasing discharge, the relation was not strong. Five of the six sites used to investigate the relations between SSC and laboratory turbidity had a coefficient of determination (r 2 ) greater than 0.7. Relations were not as strong between SSC and the turbidity measured by in situ probes because the period of record was shorter and therefore the sample sizes were smaller. Data from in situ turbidity probes were strongly related to turbidity data measured in the laboratory for all but one of the monitoring sites where the relation was strongly leveraged by one sample. Although the in situ turbidity probes appeared to provide a good surrogate for SSC and could allow more accurate calculations of suspended-sediment load than discrete suspended-sediment samples alone, more data would be required to define the regression models throughout the range in discharge, SSCs, and turbidity levels that occur at each monitoring site. Nonetheless, the in situ probes provided much greater detail about the relation between discharge and turbidity than did the grab samples and storm samples measured in the laboratory.

New York↗
USGS2012

Description of 2005-10 domestic water use for selected U.S. cities and guidance for estimating domestic water use

Domestic water-use and related socioeconomic and climatic data for 2005-10 were used in an analysis of 21 selected U.S. cities to describe recent domestic per capita water use, investigate variables that potentially affect domestic water use, and provide guidance for estimating domestic water use. Domestic water use may be affected by a combination of several factors. Domestic per capita water use for the selected cities ranged from a median annual average of 43 to 177 gallons per capita per day (gpcd). In terms of year-to-year variability in domestic per capita water use for the selected cities, the difference from the median ranged from ± 7 to ± 26 percent with an overall median variability of ± 14 percent. As a percentage of total annual water use, median annual domestic water use for the selected cities ranged from 33 to 71 percent with an overall median of 57 percent. Monthly production and water sales data were used to calculate daily per capita water use rates for the lowest 3 consecutive months (low-3) and the highest 3 consecutive months (high-3) of usage. Median low-3 domestic per capita water use for 16 selected cities ranged from 40 to 100 gpcd. Median high-3 domestic per capita water use for 16 selected cities ranged from 53 to 316 gpcd. In general, the median domestic water use as a percentage of the median total water use for 16 selected cities was similar for the low-3 and high-3 periods. Statistical analyses of combined data for the selected cities indicated that none of the socioeconomic variables, including cost of water, were potentially useful as determinants of domestic water use at the national level. However, specific socioeconomic variables may be useful for the estimation of domestic water use at the State or local level. Different socioeconomic variables may be useful in different States. Statistical analyses indicated that specific climatic variables may be useful for the estimation of domestic water use for some, but not all, of the selected cities. National average public supply per capita water use declined from 185 gpcd in 1990 to 171 gpcd in 2005. National average domestic delivery per capita water use declined from 105 gpcd in 1990 to 99 gpcd in 2005. Average State domestic delivery per capita water use ranged from 51 to 189 gpcd in 2005. The average annual total per capita water use in 19 selected cities that provided data for each year declined from 167 gpcd in 2006 to 145 gpcd in 2010. During this time period, average per capita water use measured during the low-3 period each year declined from 115 to 102 gpcd, and average per capita use measured during the high-3 period declined from 250 to 211 gpcd. Continued collection of data on water deliveries to domestic populations, as well as updated estimates of the population served by these deliveries, is recommended for determination of regional and temporal trends in domestic per capita water use. Declines in various measures of per capita water use have been observed in recent years for several States with municipal water use data-collection programs. Domestic self-supplied water use historically has not been metered. Estimates of self-supplied domestic water use are made using estimates of the population that is not served by public water suppliers and per capita coefficients. For 2005, the average State domestic self-supplied per capita use in the United States ranged from 50 to 206 gpcd. The median domestic self-supplied per capita use was 76 gpcd for States in which standard coefficients were used, and 98 gpcd for States in which coefficients were based on domestic deliveries from public supply. In specific areas with scarce resources or increasing numbers of households with private wells, an assessment of domestic water use may require metering of households or development of more specific per capita coefficients to estimate water demand.

California;Kansas;Mississippi;Montana;Texas;Wiscon↗
USGS2021

Assessment of water availability in the Osage Nation using an integrated hydrologic-flow model

The Osage Nation of northeastern Oklahoma, conterminous with Osage County, covers about 2,900 square miles. The area is primarily rural with 62 percent of the land being native prairie grass, and much of the area is used for cattle ranching and extraction of petroleum and natural gas. Protection of water rights are important to the Osage Nation because of its reliance on cattle ranching and the potential for impairment of water quality by petroleum extraction. Additionally, the potential for future population increases, demands for water from neighboring areas such as the Tulsa metropolitan area, and expansion of petroleum and natural-gas extraction on water resources of this area further the need for the Osage Nation to better understand its water availability. Therefore, the U.S. Geological Survey, in cooperation with the Osage Nation, completed a hydrologic investigation to assess the status and availability of surface-water and groundwater resources in the Osage Nation. A transient integrated hydrologic-flow model was constructed using the U.S. Geological Survey fully integrated hydrologic-flow model called the MODFLOW One-Water Hydrologic Model. The integrated hydrologic-flow model, called the Osage Nation Integrated Hydrologic Model (ONIHM), was constructed and uses an orthogonal grid of 276 rows and 289 columns, and each grid cell measures 1,312.34 feet (ft; 400 meters) per side, with eight variably thick vertical layers that represented the alluvial and bedrock aquifers within the study area, including the alluvial aquifer, the Vamoosa-Ada aquifer, and the minor Pennsylvanian bedrock aquifers, and the confining units. Landscape and groundwater-flow processes were simulated for two periods: (1) the 1950–2014 period from January 1950 through September 2014 and (2) the forecast period from October 2014 through December 2099. The 1950–2014 period ONIHM simulated past conditions using measured or estimated inputs, and the forecast-period ONIHM simulated three separate potential forecast conditions under constant dry, average, or wet climate conditions using calibrated input values from the 1950–2014 period ONIHM. The 1950–2014 period ONIHM was calibrated by linking the Parameter Estimation software (PEST) with the MODFLOW One-Water Hydrologic Model. PEST uses statistical parameter estimation techniques to identify the best set of parameter values to minimize the difference between measured or estimated calibration targets and their simulated equivalent values (residuals). Tikhonov regularization and singular-value decomposition-assist features of PEST were used during the calibration process. The 1950–2014 period ONIHM was calibrated to 713 measured groundwater levels at 195 wells; 95,636 estimated monthly mean groundwater levels at 124 wells; 5,307 measured streamflows at 13 streamgages; and 8,679 simulated mean monthly streamflows at 10 streamgages extracted from a surface-water model by adjusting 231 parameters. The estimated groundwater-level observations and streamflows were included as observations to improve the spatial and temporal density of observation targets during calibration. The best set of parameter values obtained during the calibration process of the 1950–2014 model was then used as the input parameter values for the forecast model simulations. A comparison of the calibration targets to their corresponding simulated values indicated that the model adequately reproduced streamflows and groundwater levels for some streamgages and wells and underestimated streamflows and groundwater levels at other locations. Measured and simulated streamflows correlated adequately with a coefficient of determination of 0.938, as did water levels with a coefficient of determination of 0.795. The 1950–2014 period ONIHM underestimated certain groundwater levels and streamflows, but generally measured or estimated calibration targets correlated well with simulated equivalents, which indicated that the model can adequately simulate the response of the hydrologic system to stresses in the 1950–2014 and forecast periods. In the 1950–2014 period ONIHM, the calibrated mean horizontal hydraulic conductivity for layer 1 alluvial aquifer was 30.7 feet per day, and the seven lower layers had a calibrated mean horizontal hydraulic conductivity of less than 3.3 feet per day. The mean calibrated groundwater-level residual was 16.6 ft, and the mean calibrated streamflow residual of the Arkansas River at Ralston, Oklahoma, streamgage (U.S. Geological Survey station 07152500) was within 6 percent (373 cubic feet per second) of mean measured streamflow for the 1950–2014 period ONIHM. The ONIHM simulated landscape fluxes of precipitation; groundwater applied by irrigation wells; evapotranspiration from precipitation, groundwater, and irrigation; runoff from precipitation; and deep percolation from precipitation. The largest loss of water from the landscape was evapotranspiration from precipitation with a calibrated mean annual outflow of 32 inches (in.): mean annual precipitation was about 36 in. Calibrated mean annual runoff and deep percolation (recharge to the water table) rates were 4.7 inches per year (in/yr) and 0.70 in/yr, respectively, for the 1950–2014 period ONIHM. The calibrated 1950–2014 period ONIHM groundwater fluxes included net farm net recharge (calculated as the difference between the inflow of recharge to the water table and the outflow of evapotranspiration from the water table such that negative values indicate that evapotranspiration from the water table was greater than deep percolation [recharge to the water table] and vice versa). Net farm net recharge was the largest flux from the groundwater system with a mean annual net outflow of 153.4 cubic feet per second. Stream leakage was the largest flux to the groundwater system with a mean annual net inflow of 152.5 cubic feet per second, indicating that, on average, the groundwater/surface-water interaction was a “losing” system where stream water leaked into the subsurface and recharged the water table. Simulated monthly trends demonstrated that net stream leakage was the largest inflow to the groundwater-flow system for 10 of the 12 months; for the other 2 months (January and March), farm net recharge (January) and net storage (March) were the largest inflow to the groundwater-flow system. A saline groundwater interface map was created for the study and compared to the water levels from the final stress period of the 1950–2014 model to identify the presence of fresh/marginal groundwater throughout the study area. Fresh/marginal groundwater was characterized as groundwater with less than 1,500 milligrams per liter of total dissolved solids. Fresh/marginal groundwater thickness ranged from 0 to 438.2 ft within the study area. The thickest regions of fresh/marginal groundwater were in the eastern part of the study area near Sand Creek, Bird Creek, and Hominy Creek and in the Arkansas River alluvial aquifer in the region downstream from the Arkansas River at Ralston, Okla. Like the 1950–2014 model, forecast model results for the landscape indicated that transpiration from precipitation was the largest flux out of the landscape for all three forecasts, constituting 77, 73, and 58 percent of precipitation for the dry, average, and wet forecasts, respectively. The dry and average forecast landscape fluxes demonstrated similar trends and magnitudes, whereas the wet forecast landscape fluxes indicated the largest changes compared to the average forecast fluxes. Most notably, runoff increased from a mean of 1.1 and 1.6 in/yr for the dry and average forecasts, respectively, to 10 in/yr for the wet forecast. Similar changes occurred for the other wet forecast landscape fluxes. The calibrated 1950–2014 period ONIHM simulated three forecasts to assess the effects of potential climatic changes on the hydrologic system from October 2014 to December 2099. The three forecasts simulated theoretical dry, average, and wet conditions using precipitation and potential evapotranspiration datasets from selected years in the calibrated 1950–2014 period ONIHM. Annual precipitation amounts were 26.89, 35.47, and 50.73 in. for the dry, average, and wet forecasts, respectively. Groundwater-flow component forecast results indicated that stream leakage is always a net inflow to the groundwater-flow system for dry, average, and wet conditions, meaning the study area stream network is always predominantly a “losing” regime where stream water infiltrates into the underlying aquifer. Storage was only a net outflow from the groundwater-flow system and indicated a replenishment to groundwater storage that resulted in an increase in groundwater levels only during the wet forecast. Further, these gains in groundwater storage for the wet forecast occurred only during February through June. Mean fresh/marginal groundwater saturated thicknesses were 125 and 126 ft for the dry and average forecast conditions, respectively, and wet forecast average thickness was 145 ft and ranged from 0 to 443 ft. The spatial extents of fresh/marginal groundwater at the end of the dry, average, and wet forecast model periods (December 2099) did not change substantially from the end of the 1950–2014 model period (September 2014).

Kansas, Oklahoma↗
USGS2013

Streamflow monitoring and statistics for development of water rights claims for Wild and Scenic Rivers, Owyhee Canyonlands Wilderness, Idaho, 2012

The U.S. Geological Survey, in cooperation with the Bureau of Land Management (BLM), collected streamflow data in 2012 and estimated streamflow statistics for stream segments designated "Wild," "Scenic," or "Recreational" under the National Wild and Scenic Rivers System in the Owyhee Canyonlands Wilderness in southwestern Idaho. The streamflow statistics were used by BLM to develop and file a draft, federal reserved water right claim in autumn 2012 to protect federally designated "outstanding remarkable values" in the stream segments. BLM determined that the daily mean streamflow equaled or exceeded 20 and 80 percent of the time during bimonthly periods (two periods per month) and the bankfull streamflow are important streamflow thresholds for maintaining outstanding remarkable values. Prior to this study, streamflow statistics estimated using available datasets and tools for the Owyhee Canyonlands Wilderness were inaccurate for use in the water rights claim. Streamflow measurements were made at varying intervals during February–September 2012 at 14 monitoring sites; 2 of the monitoring sites were equipped with telemetered streamgaging equipment. Synthetic streamflow records were created for 11 of the 14 monitoring sites using a partial‑record method or a drainage-area-ratio method. Streamflow records were obtained directly from an operating, long-term streamgage at one monitoring site, and from discontinued streamgages at two monitoring sites. For 10 sites analyzed using the partial-record method, discrete measurements were related to daily mean streamflow at a nearby, telemetered “index” streamgage. Resulting regression equations were used to estimate daily mean and annual peak streamflow at the monitoring sites during the full period of record for the index sites. A synthetic streamflow record for Sheep Creek was developed using a drainage-area-ratio method, because measured streamflows did not relate well to any index site to allow use of the partial-record method. The synthetic and actual daily mean streamflow records were used to estimate daily mean streamflow that was exceeded 80, 50, and 20 percent of the time (80-, 50-, and 20-percent exceedances) for bimonthly and annual periods. Bankfull streamflow statistics were calculated by fitting the synthetic and actual annual peak streamflow records to a log Pearson Type III distribution using Bulletin 17B guidelines in the U.S. Geological Survey PeakFQ program. The coefficients of determination (R 2 ) for the regressions between the monitoring and index sites ranged from 0.74 for Wickahoney Creek to 0.98 for the West Fork Bruneau River and Deep Creek. Confidence in computed streamflow statistics is highest among other sites for the East Fork Owyhee River and the West Fork Bruneau River on the basis of regression statistics, visual fit of the related data, and the range and number of streamflow measurements. Streamflow statistics for sites with the greatest uncertainty included Big Jacks, Little Jacks, Cottonwood, Wickahoney, and Sheep Creeks. The uncertainty in computed streamflow statistics was due to a number of factors which included the distance of index sites relative to monitoring sites, relatively low streamflow conditions that occurred during the study, and the limited number and range of streamflow measurements. However, the computed streamflow statistics are considered the best possible estimates given available datasets in the remote study area. Streamflow measurements over a wider range of hydrologic and climatic conditions would improve the relations between streamflow characteristics at monitoring and index sites. Additionally, field surveys are needed to verify if the streamflows selected for the water rights claims are sufficient for maintaining outstanding remarkable values in the Wild and Scenic rivers included in the study.

Idaho;Nevada;Oregon↗
USGS2014

Regional regression equations for the estimation of selected monthly low-flow duration and frequency statistics at ungaged sites on streams in New Jersey

Regional regression equations were developed for estimating monthly flow-duration and monthly low-flow frequency statistics for ungaged streams in Coastal Plain and non-coastal regions of New Jersey for baseline and current land- and water-use conditions. The equations were developed to estimate 87 different streamflow statistics, which include the monthly 99-, 90-, 85-, 75-, 50-, and 25-percentile flow-durations of the minimum 1-day daily flow; the August–September 99-, 90-, and 75-percentile minimum 1-day daily flow; and the monthly 7-day, 10-year (M7D10Y) low-flow frequency. These 87 streamflow statistics were computed for 41 continuous-record streamflow-gaging stations (streamgages) with 20 or more years of record and 167 low-flow partial-record stations in New Jersey with 10 or more streamflow measurements. The regression analyses used to develop equations to estimate selected streamflow statistics were performed by testing the relation between flow-duration statistics and low-flow frequency statistics for 32 basin characteristics (physical characteristics, land use, surficial geology, and climate) at the 41 streamgages and 167 low-flow partial-record stations. The regression analyses determined drainage area, soil permeability, average April precipitation, average June precipitation, and percent storage (water bodies and wetlands) were the significant explanatory variables for estimating the selected flow-duration and low-flow frequency statistics. Streamflow estimates were computed for two land- and water-use conditions in New Jersey—land- and water-use during the baseline period of record (defined as the years a streamgage had little to no change in development and water use) and current land- and water-use conditions (1989–2008)—for each selected station using data collected through water year 2008. The baseline period of record is representative of a period when the basin was unaffected by change in development. The current period is representative of the increased development of the last 20 years (1989–2008). The two different land- and water-use conditions were used as surrogates for development to determine whether there have been changes in low-flow statistics as a result of changes in development over time. The State was divided into two low-flow regression regions, the Coastal Plain and the non-coastal region, in order to improve the accuracy of the regression equations. The left-censored parametric survival regression method was used for the analyses to account for streamgages and partial-record stations that had zero flow values for some of the statistics. The average standard error of estimate for the 348 regression equations ranged from 16 to 340 percent. These regression equations and basin characteristics are presented in the U.S. Geological Survey (USGS) StreamStats Web-based geographic information system application. This tool allows users to click on an ungaged site on a stream in New Jersey and get the estimated flow-duration and low-flow frequency statistics. Additionally, the user can click on a streamgage or partial-record station and get the “at-site” streamflow statistics. The low-flow characteristics of a stream ultimately affect the use of the stream by humans. Specific information on the low-flow characteristics of streams is essential to water managers who deal with problems related to municipal and industrial water supply, fish and wildlife conservation, and dilution of wastewater.

New Jersey↗
USGS2023

Critical minerals in subduction-related magmatic-hydrothermal systems of the United States

During the World War and Cold War eras (1910s–1990s), domestic consumption of numerous mineral commodities relied increasingly on imported supplies. Consumption reliance has since expanded to include 50 “critical minerals” (elements and mineral commodities) that are mostly to entirely imported and subject to curtailment by suppliers or supply chain disruption. New domestic supplies of critical minerals are being pursued by mining companies and by several federal departments and agencies. Information on domestic deposits and resources of critical minerals is being compiled by the U.S. Geological Survey Mineral Resources Program, which has organized investigations by mineral system, deposit type, and commodity. Production, reserves, resources, and inventories of 21 critical minerals in domestic magmatic-hydrothermal deposits related to subduction-generated magmatism, and in tailings, slag, slimes, and electrolyte from copper concentrators, smelters, and refineries that processed some deposits, are largely restricted to Western States and Alaska. The critical mineral commodities Al, Sb, As, Bi, Co, fluorite, Ga, Ge, In, Mn, Ni, Nb, Pd, Pt, potash, Re, Ta, Te, Sn, W, and V are variably concentrated in porphyry/skarn copper-(molybdenum), skarn-replacement-vein (S-R-V) tungsten, polymetallic sulfide S-R-V intermediate sulfidation (IS), high-sulfidation gold-silver, low-sulfidation gold-silver, and lithocap alunite deposit types. These deposit types occur in porphyry copper-molybdenum-gold, alkalic porphyry, porphyry tin (granite related), and reduced intrusion-related mineral systems. Production, reserves, and resources of Co, Ni, Nb, Pd, Pt, Ta, Sn, and V in subduction-related deposits in Western States are insignificant to small, mostly equivalent to months to a few years of recent annual domestic consumption (2016–2020). Significant inventories, equivalent to 2 or more years of consumption of aluminum, antimony, potash, and tungsten in unmined S-R-V tungsten, polymetallic sulfide S-R-V-IS, and lithocap alunite deposits vary from approximately 2 to 8 years. Several decades of consumption of arsenic, bismuth, fluorite, gallium, germanium, and indium exist in some polymetallic sulfide S-R-V-IS and lithocap alunite deposit types. Based on concentrations of critical minerals in reserves, resources, drill holes, and deposit domains (ore types), and in captive refinery records, the largest domestic inventories of Sb, As, Bi, Re, and Te, and possibly Ga, Ge, In, Sn, and W, are in porphyry copper-molybdenum (Cu-Mo) deposits in Alaska, Idaho, Utah, and Arizona, and in interim products of processing porphyry Cu-Mo deposit ores for recovery of copper and molybdenum. Concentrations of critical minerals in archival specimens and sample collections, although somewhat biased by collection and conservation decisions and categorization, are broadly proportionate to those in reserves, resources, and drill holes. These concentrations imply significant inventories of some critical minerals in deposits for which production, resources, and refinery records are unavailable or incomplete. Because of the large masses of ores mined and processed annually (hundreds of millions of metric tons) and in reserves and resources (hundreds of millions of metric tons to billions of metric tons), calculated inventories of critical minerals in porphyry Cu-Mo deposits are equivalent to decades and centuries of recent consumption. However, these inventories should not be considered consumable supplies without reserve definition and development of economically viable mining plans and recovery techniques. An expeditious strategy for elimination or reduction of import reliance is recovery, and improved recovery efficiency, of Sb, As, Bi, Re, and Te, and possibly Ga, Ge, In, Ni, Sn, Ti, and W; during concentration and refining of copper and molybdenum minerals in ores of operating porphyry Cu-Mo mines; and in unmined porphyry Cu-Mo resources. These chalcophile, siderophile, and lithophile critical minerals, often undetectable in ore, are concentrated (hundreds of parts per million to percents) in slimes and electrolyte during copper electrorefining or could be recovered, in part, during sulfide concentration and smelting. Other than rhenium (recovered during molybdenum refining) and tellurium, all have been routinely discarded. Subsidization (for example, commodity price guarantees, tax credits, recovery technology development), political initiative, and (or) sustained market favorability could support new production of critical mineral commodities from subduction-related magmatic-hydrothermal deposits in Western States. In addition, insufficient domestic refining capacity could relegate the large inventories of critical minerals in porphyry Cu-Mo reserves and resources (for example, Pebble, Alaska; Resolution and Copper World [Rosemont], Arizona) to exportation in concentrates and importation insecurity, fortifying their present status.

Scientific Investigations Report↗
USGS2010

Effects of groundwater levels and headwater wetlands on streamflow in the Charlie Creek basin, Peace River watershed, west-central Florida

The Charlie Creek basin was studied from April 2004 to December 2005 to better understand how groundwater levels in the underlying aquifers and storage and overflow of water from headwater wetlands preserve the streamflows exiting this least-developed tributary basin of the Peace River watershed. The hydrogeologic framework, physical characteristics, and streamflow were described and quantified for five subbasins of the 330-square mile Charlie Creek basin, allowing the contribution of its headwaters area and tributary subbasins to be separately quantified. A MIKE SHE model simulation of the integrated surface-water and groundwater flow processes in the basin was used to simulate daily streamflow observed over 21 months in 2004 and 2005 at five streamflow stations, and to quantify the monthly and annual water budgets for the five subbasins including the changing amount of water stored in wetlands. Groundwater heads were mapped in Zone 2 of the intermediate aquifer system and in the Upper Floridan aquifer, and were used to interpret the location of artesian head conditions in the Charlie Creek basin and its relation to streamflow. Artesian conditions in the intermediate aquifer system induce upward groundwater flow into the surficial aquifer and help sustain base flow which supplies about two-thirds of the streamflow from the Charlie Creek basin. Seepage measurements confirmed seepage inflow to Charlie Creek during the study period. The upper half of the basin, comprised largely of the Upper Charlie Creek subbasin, has lower runoff potential than the lower basin, more storage of runoff in wetlands, and periodically generates no streamflow. Artesian head conditions in the intermediate aquifer system were widespread in the upper half of the Charlie Creek basin, preventing downward leakage from expansive areas of wetlands and enabling them to act as headwaters to Charlie Creek once their storage requirements were met. Currently, the dynamic balance between wetland storage, rainfall-runoff processes, and groundwater-level differences in the upper basin allow it to generate approximately half of the streamflow from the Charlie Creek basin. Therefore, future development in the upper basin that would alter the hydraulic connectivity of wetlands during high flow conditions or expand recharging groundwater conditions could substantially affect streamflow in Charlie Creek. LIDAR (Light detection and ranging) based topographic maps and integrated modeling results were used to quantify the water stored in wetlands and other topographic depressions, and to describe the network of shallow stream channels connecting wetlands to Charlie Creek and its tributaries over distances of several thousand feet. Peak flows at all but one streamflow station were underpredicted in MIKE SHE simulations, possibly because the hydraulics of surface channels connecting wetlands to stream channels were not explicitly simulated in the model. Explicitly simulating the smaller channels connecting wetlands and stream channels should improve the ability of future watershed models to simulate peak flows in streams with headwater wetlands. The runoff potential was greater in the lower half of the Charlie Creek basin than in the upper half, and the streambed of Charlie Creek had greater potential to both directly gain streamflow from groundwater and lose streamflow to groundwater. Charlie Creek is more incised into the surficial aquifer in the lower basin than in the upper basin, and the streambed intersects the top of the intermediate aquifer system at two known locations. Groundwater levels in the intermediate aquifer system varied widely in the lower half of the basin from artesian conditions inducing upward flow toward the surficial aquifer and streams, to recharging conditions allowing downward flow and stream leakage. Recharge areas were greatest in May 2004 when rainfall was at a seasonal low and irrigation pumping was at a seasonal high. Recharge conditions

Scientific Investigations Report↗
USGS2014

Equations for estimating selected streamflow statistics in Rhode Island

Regional regression equations were developed for estimating selected natural—unaffected by alteration—streamflows of specific flow durations and low-flow frequency statistics for ungaged stream sites in Rhode Island. Selected at-site streamflow statistics are provided for 41 long-term streamgages, 21 short-term streamgages, and 135 partial-record stations in Rhode Island, eastern Connecticut, and southeastern and south-central Massachusetts. The regression equations for estimating selected streamflow statistics and the at-site statistics estimated for each of the 197 sites may be used by Federal, State, and local water managers in addressing water issues in and near Rhode Island. Multiple and simple linear regression equations were developed to estimate the 99-, 98-, 95-, 90-, 85-, 80-, 75-, 70-, 60-, 50-, 40-, 30-, 25-, 20-, 15-, 10-, 5-, 2-, and 1-percent flow durations and the 7Q2 (7-day, 2-year) and 7Q10 (7-day, 10-year) low-flow-frequency statistics. An additional 49 selected statistics, for which regression equations were not developed, also were estimated for the long- and short-term streamgages and partial-record stations for flow durations between the 99.99 and 0.01 percent and for the mean annual, mean monthly, and median monthly streamflows. A total of 70 selected streamflow statistics were estimated for 41 long-term streamgages, 21 short-term streamgages, and 135 partial-record stations in and near Rhode Island. Estimates of the long-term streamflow statistics for the 21 short-term streamgages and 135 partial-record stations were developed by the Maintenance of Variance Extension, type 1 (MOVE.1), record-extension technique. The equations used to estimate selected streamflow statistics were developed by relating the 19 flow-duration and 2 low-flow-frequency statistics to 31 different basin characteristics (physical, land-cover, and climatic) at the 41 long-term and 19 of 21 short-term streamgages (a total of 60 streamgages) in and near Rhode Island. The 135 partial-record stations were not used in the regression analyses. The regression analyses were done by using a user-weighted least-squares technique in the weighted-multiple-linear regression program for the 90- to 1-percent flow-duration statistics. For the 99-, 98-, and 95-percent flow durations and the 7Q2 and 7Q10 statistics, left-censored regression analyses were used to account for zero flows at a few streamgages. The regression analyses determined that two basin characteristics—drainage area and stream density—were the only significant explanatory variables for 16 of the 19 flow-duration and the 2 low-flow regression equations. For the 10-, 15-, and 20-percent flow-duration regression equations, drainage area was the only significant explanatory variable. The standard error of the estimate for the 21 regression equations ranged from 17.58 to 141.83 percent. The 99- to 85-percent flow durations and the low-flow statistics 7Q2 and 7Q10 had the highest standard errors of the estimate, ranging from 48.68 to 141.83 percent. The standard error of the estimate for the medium- to high-flow statistics—the 80- to 1-percent flow durations—ranged from 17.58 to 37.65 percent, with the standard errors for the 60- to 1-percent flow durations all being less than about 21 percent. Data also are provided to allow the user to calculate the 90-percent prediction intervals for the 21 streamflow statistics. The equations, which are based on data from streams with little to no flow alterations, will provide an estimate of the natural flows for a selected site. They will not estimate flows for altered sites with dams, surface-water withdrawals, groundwater withdrawals (pumping wells), diversions, and wastewater discharges. If the equations are used to estimate streamflow statistics for altered sites, the user should adjust the flow estimates for the alterations. The regression equations should be used only for ungaged sites with drainage areas between 0.52 and 294 square miles and stream densities between 0.94 and 3.49 miles per square mile; these are the ranges of the explanatory variables in the equations.

Rhode Island↗
USGS2007

Investigation of Ground-Water Contamination at Solid Waste Management Unit 12, Naval Weapons Station Charleston, North Charleston, South Carolina

The U.S. Geological Survey and the Naval Facilities Engineering Command Southeast investigated natural and engineered remediation of chlorinated volatile organic compound ground-water contamination at Solid Waste Management Unit 12 at the Naval Weapons Station Charleston, North Charleston, South Carolina. The primary contaminants of interest are tetrachloroethene, 1,1,1-trichloroethane, trichloroethene, cis-1,2-dichloroethene, vinyl chloride, 1,1-dichloroethane, and 1,1-dichloroethene. In general, the hydrogeology of Solid Waste Management Unit 12 consists of a surficial aquifer, composed of sand to clayey sand, overlain by dense clay that extends from about land surface to a depth of about 8 to 10 feet and substantially limits local recharge. During some months in the summer, evapotranspiration and limited local recharge result in ground-water level depressions in the forested area near wells 12MW-12S and 12MW-17S, seasonally reflecting the effects of evapotranspiration. Changes in surface-water levels following Hurricane Gaston in 2004 resulted in a substantial change in the ground-water levels at the site that, in turn, may have caused lateral shifting of the contaminant plume. Hydraulic conductivity, determined by slug tests, is higher along the axis of the plume in the downgradient part of the forests than adjacent to the plume, implying that there is some degree of lithologic control on the plume location. Hydraulic conductivity, hydraulic gradient, sulfur-hexafluoride measurements, and historical data indicate that ground-water flow rates are substantially slower in the forested area relative to upgradient areas. The ground-water contamination, consisting of chlorinated volatile organic compounds, extends eastward in the surficial aquifer from the probable source area near a former underground storage tank. Engineered remediation approaches include a permeable reactive barrier and phytoremediation. The central part of the permeable reactive barrier along the main axis of the contaminant plume appears to be actively removing contamination; however, ground-water contamination is moving around the southern end of the permeable reactive barrier. Changes in the contaminant concentrations along the path of ground-water transport reflect a complex variety of influences. Potential influences include dechlorination, sorption and desorption, transpirative removal by trees, lateral shifting of the plume, and the presence of zones of differing concentrations possibly reflecting one or more pulse releases of contamination from the source area. Near the source area at well 12MW-10S, volatile organic compound concentrations of cis-1,2-dichlorothene, vinyl chloride, 1,1-dichloroethane, and 1,1,1-trichloroethane continued an irregular decline, while tetrachloroethene and 1,1-dichloroethene showed marked fluctuations in concentration during 2005 and 2006. Volatile organic compound concentrations at well 12MW-03S continued to show decreasing concentrations with the June 2006 concentrations being the lowest yet recorded at that well for several volatile organic compounds. Concentration and delta carbon 13 data indicate that in the upgradient part of the plume, tetrachloroethene is being degraded to trichloroethene, which is being degraded to cis-1,2-dichloroethene, and cis-1,2-dichloroethene is accumulating faster than it is being depleted. Ground-water volatile organic compound concentrations also changed in some wells in the forested area in the midpart of the plume. Increasing tetrachloroethene and decreasing trichloroethene and 1,1-dichloroethene concentrations were observed at wells 12MW-05S and 12MW-29S, possibly reflecting a lateral shift in the axis of the contamination plume or an advancing contamination pulse. Substantial decreases in contamination occur in the forested area downgradient from well 12MW-05S. Probable major loss mechanisms in this area include evapotranspiration and sorption.

North Carolina, South Carolina↗
USGS2018

Preliminary synthesis and assessment of environmental flows in the middle Verde River watershed, Arizona

A 3-year study was undertaken to evaluate the suitability of the available modeling tools for characterizing environmental flows in the middle Verde River watershed of central Arizona, describe riparian vegetation throughout the watershed, and estimate sediment mobilization in the river. Existing data on fish and macroinvertebrates were analyzed in relation to basin characteristics, flow regimes, and microhabitat, and a pilot study was conducted that sampled fish and macroinvertebrates and the microhabitats in which they were found. The sampling for the pilot study took place at five different locations in the middle Verde River watershed. This report presents the results of this 3-year study. The Northern Arizona Groundwater Flow Model (NARGFM) was found to be capable of predicting long-term changes caused by alteration of regional recharge (such as may result from climate variability) and groundwater pumping in gaining, losing, and dry reaches of the major streams in the middle Verde River watershed. Over the period 1910 to 2006, the model simulated an increase in dry reaches, a small increase in reaches losing discharge to the groundwater aquifer, and a concurrent decrease in reaches gaining discharge from groundwater. Although evaluations of the suitability of using the NARGFM and Basin Characteristic Model to characterize various streamflow intervals showed that smallerscale basin monthly runoff could be estimated adequately at locations of interest, monthly stream-flow estimates were found unsatisfactory for determining environmental flows. Orthoimagery and Moderate Resolution Imaging Spectroradiometer data were used to quantify stream and riparian vegetation properties related to biotic habitat. The relative abundance of riparian vegetation varied along the main channel of the Verde River. As would be expected, more upland plant species and fewer lowland species were found in the upper-middle section compared to the lower-middle section, and vice-versa. Vegetation changes within the upper-middle and lower-middle reaches are related to differences in climate and hydrology. In general, the riparian vegetation of the middle Verde River watershed is that of a healthy ecosystem’s mixed age, mixed patch structure, likely a result of the mostly unaltered disturbance regime. The frequency of in-river hydrogeomorphic features (pool, riffle, run) varied along the middle Verde River channel. There was a greater abundance of riffle habitat in the upper-middle reach; the lower-middle reach included more pool habitat. The Oak Creek tributary was more homogenous in geomorphic stream habitat composition than West Clear Creek, where runs dominated the upper reaches and pools dominated many of the lower reaches. On the basis of the period of record and discharges recorded at 15-minute intervals, five flows were found to reach the gravel-transport threshold. Sediment mobilization computed with flows averaged over daily time steps yielded just three flows that reached the gravel-transport threshold, and monthly averaged flows yielded none. In the middle Verde River watershed, 15-minute data should be used when possible to evaluate sediment transport in the river system. Data from more than 300 fish surveys conducted from 1992 to 2011 were analyzed using two schemes, one that divided the river into five reaches based on basin characteristics, and a second that divided the river into five reaches based on degree of flow alteration (specifically, diversions). Fish community metrics and assemblage data were used to analyze patterns of species composition and abundance in the two approaches. Overall, native and non-native species were regularly interacting and probably competing for similar resources. Fish abundances were also analyzed in response to floods and other flow metrics. Although the data are limited, native fish abundances increased more rapidly than non-native fish abundances in response to large floods. The basin-characteristic reach analysis showed native fish in greater abundance in the upper-middle reaches of the Verde River watershed and generally decreasing with downstream distance. The median relative abundance of native fish decreased by 50 percent from reach 1 to reach 5. Using the reach scheme based on degree of flow alteration, nondiverted reaches were found to have a greater abundance of native fish than diverted reaches. In heavily diverted reaches, non-native species outnumbered native species. Fish metrics and stream-flow metrics for the 30, 90, and 365-day periods before collection were computed and the results analyzed statistically. Only abundance of all fish species was associated with the 30-day flow metrics. The 90-day flow metrics were generally positively associated with fish metrics, whereas the 365-day flow metrics had more negative correlations. In particular, significant relations were found between fish metrics and the magnitude and frequency of high flows, including maximum monthly flow, median annual number of high-flow events, and median annual maximum streamflow. Native sucker (Catostomidae) populations tended to decrease in periods of extended base flow, and fish in the non-native sunfish family (Centrarchidae) decreased in periods of flashy, high magnitude flows. A pilot study surveyed fish at five locations in the upper part of the middle Verde River watershed as a means to measure microhabitat availability and quantify native and non-native fish use of that available microhabitat. Results indicated that native and non-native species exhibit some clear differences in microhabitat use. Although at least some native and non-native fish were found in each velocity, depth, and substrate category, preferential microhabitat use was common. On a percentage basis, non-native species had a strong preference for slow-moving and deeper water with silt and sand substrate, with a secondary preference for faster moving and very shallow water and a coarse gravel substrate. Native species showed a general preference for somewhat faster, moderate depth water over coarse gravel and had no clear secondary preference. Macroinvertebrate-variables index period, high-flow year, and collection location (upper-middle Verde River, lowermiddle Verde River, or Verde River tributaries) were found to be important explanatory variables in differentiating among community metrics. Overall richness (number of unique taxa), Shannon’s diversity index, and the percent of the most dominant taxa were all highly correlated, but their response to each macroinvertebrate variable was different. The percentage of mayfly (order Ephemeroptera) taxa was significantly higher in Oak Creek and the upper-middle and lower-middle Verde River reaches, locations which have higher flows and more urbanization than other reaches. When community metrics were related to hydrologic metrics, caddisfly (order Trichoptera) populations appeared to increase and mayfly populations to decrease in response to less flashy and more stable streamflows. Conversely, caddisfly populations appeared to decrease and mayfly populations to increase in response to greater flow variability. Six locations along the Verde River were sampled for macroinvertebrates as part of a pilot study associated with this report—(1) below Granite Creek, (2) near Campbell Ranch, (3) at the U.S. Geological Survey Paulden gage, (4) at the Perkinsville Bridge, (5) at the USGS Clarkdale gage, and (6) near the Reitz Ranch property. A nonmetric multidimensional scaling ordination of macroinvertebrate assemblages showed that the Verde River below Granite Creek site was different from the five other sites and that the Perkinsville Bridge and near Reitz Ranch samples had similar community structure. The near Campbell Ranch and Paulden gage locations had similar microhabitat characteristics, with the exception of riparian cover, yet the assemblage structure was very different. The different community composition at Verde River below Granite Creek was likely due to it having the smallest substrate sizes, lowest velocities, shallowest depths, and most riparian cover of the six sites.

Arizona↗
USGS2008

A comparison of runoff quantity and quality from two small basins undergoing implementation of conventional- and low-impact-development (LID) strategies: Cross Plains, Wisconsin, water years 1999–2005

Environmental managers are often faced with the task of designing strategies to accommodate development while minimizing adverse environmental impacts. Low-impact development (LID) is one such strategy that attempts to mitigate environmental degradation commonly associated with impervious surfaces. The U.S. Geological Survey, in cooperation with the Wisconsin Department of Natural Resources, studied two residential basins in Cross Plains, Wis., during water years 1999–2005. A paired-basin study design was used to compare runoff quantity and quality from the two basins, one of which was developed in a conventional way and the other was developed with LID. The conventional-developed basin (herein called “conventional basin”) consisted of curb and gutter, 40-foot street widths, and a fully connected stormwater-conveyance system. The LID basin consisted of grassed swales, reduced impervious area (32-foot street widths), street inlets draining to grass swales, a detention pond, and an infiltration basin. Data collected in the LID basin represented predevelopment through near-complete build-out conditions. Smaller, more frequent precipitation events that produced stormwater discharge from the conventional basin were retained in the LID basin. Only six events with precipitation depths less than or equal to 0.4 inch produced measurable discharge from the LID basin. Of these six events, five occurred during winter months when underlying soils are commonly frozen, and one was likely a result of saturated soil from a preceding storm. In the conventional basin, the number of discharge events, using the same threshold of precipitation depth, was 180, with nearly one-half of those resulting from precipitation depths less than 0.2 inch. Precipitation events capable of producing appreciable discharge in the LID basin were typically those of high intensity or precipitation depth or those that occurred after soils were already saturated. Total annual discharge volume measured from the conventional basin ranged from 1.3 to 9.2 times that from the LID basin. Development of the LID basin did not appreciably alter the hydrologic response to precipitation characterized during predevelopment conditions. Ninety-five percent or more of precipitation in the LID basin was retained during each year of construction from predevelopment through near-complete build-out, surpassing the 90-percent benchmark established for new development by the Wisconsin Department of Natural Resources. The amount of precipitation retained in the conventional basin did not exceed 94 percent and fell below the 90-percent standard 2 of the 6 years monitored. Much of the runoff in the LID basin was retained by an infiltration basin, the largest control structure used to mitigate storm-runoff quantity and quality. The infiltration basin also was the last best-management practice (BMP) used to treat runoff before it left the LID basin as discharge. From May 25, 2002, to September 30, 2005, only 24 of 155 precipitation events exceeded the retention/ infiltrative capacity of the infiltration basin. The overall reduction in runoff volume from these few events was 51 percent. The effectiveness of the infiltration basin decreased as precipitation intensities exceeded 0.5 inch per hour. Annual loads were estimated to characterize the overall effectiveness of low-impact design practices for mitigating delivery of total solids, total suspended solids, and total phosphorus. Annual loads of these three constituents were greater in the LID basin than in the conventional basin in 2000 and 2004. Seventy percent or more of all constituent annual loads were associated with two discharge events in 2000, and a single discharge event produced 50 percent or more of constituent annual loads in 2004. Each of these discharge events was associated with considerable precipitation depths and (or) intensities, ranging from 4.89 to 6.21 inches and from 1.13 to 1.2 inches per hour, respectively. These same storms did not contribute as much of the annual load in the conventional basin. With large storms and saturated soils, the ability of low-impact design techniques to reduce runoff, and thus constituent loads, can be greatly diminished. For both the LID and conventional basins, the temperature of runoff was largely affected by ambient air temperatures. However, the temperature of discharge from the LID basin increased upon runoff cessation. This increase is likely due to solar heating of water that is temporarily stored in the detention pond and infiltration basin.

Wisconsin↗
USGS2024

Geomorphic change, hydrology, and hydraulics of Caulks Creek, Wildwood, Missouri

Caulks Creek is a small stream that flows through the city of Wildwood in western St. Louis County, Missouri. The U.S. Geological Survey, in cooperation with the city of Wildwood, has documented historical and recent geomorphic change along Caulks Creek, simulated the hydrologic and hydraulic response of Caulks Creek to a variety of design storm scenarios, and simulated bank retreat resulting from fluvial erosion and mass failure processes. Six study reaches were selected for monitoring geomorphic change based on known locations of erosion issues documented by the city of Wildwood. Recent short-term rates and patterns of geomorphic change in the study reaches, with a focus on bank retreat, were determined from repeat terrestrial light detection and ranging surveys and field observations of the six study reaches. Historical aerial photographs of the study reaches were analyzed to determine long-term rates of bank retreat and channel widening. In general, rapid bank retreat and widening was observed at the outer banks of meander bends and where banks are unforested. Short-term bank retreat varied substantially within individual study reaches, across the study area, and during the study period from no change to as much as 16 feet of retreat between consecutive surveys (5 to 8 months). The field surveys and visual observations indicated that bank retreat occurs episodically owing to a combination of fluvial erosion and mass failure processes, as well as freezing and thawing cycles. Long-term rates of bank retreat ranged from 0.6 to 4.4 feet per year. Hydrologic and hydraulic models of Caulks Creek were used to quantify the peak, volume, and timing of the flow response and the spatial distribution of hydraulic drivers of erosion (velocity and shear stress) along Caulks Creek for design storm scenarios that represent current (as of this publication) and projected future climate. The projected climate conditions resulted in higher peak flows compared to current conditions, including 6 to 21 percent for the year 2050 and 10 to 42 percent for the year 2099 at the downstream end of the study area. Additionally, for a given design storm, projected climate change is predicted to result in faster flows with greater shear stress, as well as more within-stream variability in velocity and shear stress. Many factors affect the velocity and shear stress at a given location, but in general, somewhat fast velocities and high shear stresses tended to occur where the channel is relatively narrow and straight. The velocity and shear stress in the study reaches (known areas of widening and bank retreat) were not particularly high, at least in part owing to the relatively large widths and high sinuosity of the present-day channel in these reaches. The potential mitigating effect of adding runoff storage to the basin also was examined for a selection of design storm scenarios. Additional runoff storage was more effective at mitigating peak flows and total runoff volumes for higher-frequency, lower-intensity storms than for lower-frequency, higher-intensity storms. The additional storage also resulted in an overall reduction in velocity (by as much as 28 percent) and shear stress (by as much as 40 percent) in the study area. However, the effect of the additional storage on peak flows, total runoff volumes, velocity, and shear stress decreased with distance downstream through the study area. For the simulated scenarios, added runoff storage was effective at mitigating the increases in peak flows, total runoff volumes, velocity, and shear stress caused by projected climate change. Lastly, the bank stability and toe erosion model (BSTEM) was used to predict bank erosion and potential bank failure surfaces at five locations along Caulks Creek for a selection of design storm scenarios. The lower-frequency, higher-magnitude design storms resulted in more bank retreat than the higher-frequency, lower-magnitude design storms, though the magnitude of the difference was site dependent. Although scenarios with additional storage were not directly simulated in BSTEM, it is likely that the additional storage would result in reduced bank retreat compared to the same design storm with existing storage, based on the hydraulic modeling results for scenarios with added runoff storage.

Missouri↗
USGS2013

Effects of surface applications of biosolids on groundwater quality and trace-element concentrations in crops near Deer Trail, Colorado, 2004-2010

The U.S. Geological Survey (USGS), in cooperation with Metro Wastewater Reclamation District (Metro District), studied biosolids composition and the effects of biosolids applications on groundwater quality and trace-element concentrations in crops of the Metro District properties near Deer Trail, Colorado, during 2004 through 2010. Priority parameters for each monitoring component included the nine trace elements regulated by Colorado for biosolids (arsenic, cadmium, copper, lead, mercury, molybdenum, nickel, selenium, and zinc); other constituents also were analyzed. All concentrations for the priority parameters in monthly biosolids samples were less than Colorado regulatory limits, and the concentrations were relatively consistent. Biosolids likely were the largest source of nitrogen and phosphorus on the Metro District properties. Plutonium isotopes were not detected in the biosolids, but many organic wastewater compounds (organic wastewater compounds: wastewater indicators, pharmaceuticals, and hormones) were detected in substantial concentrations relative to minimum reporting levels and various surface-water concentrations. Bismuth, copper, mercury, nitrogen, phosphorus, silver, biogenic sterols, detergent degradates, disinfectants, fire retardants, fragrances, pharmaceuticals, and plasticizers would be the most likely biosolids signature to indicate the presence of Metro District biosolids in soil or streambed sediment from the study area. Antimony, cadmium, cobalt, copper, molybdenum, nickel, nitrogen, phosphorus, selenium, tungsten, vanadium, zinc, detergent degradates, disinfectants, fire retardants, fragrances, pharmaceuticals or their degradates, and plasticizers would be the most likely biosolids signature for groundwater and surface water in the study area. More biosolids-signature components detected and larger concentration differences from untreated materials, baseline, and blank samples indicate more evidence of biosolids presence or effects. Although the inorganic constituent concentrations were relatively large in samples from one monitoring well, the concentrations of organic wastewater compounds in groundwater samples were not correspondingly large. Concentrations of organic wastewater compounds in the groundwater samples from all five monitoring wells were less than the minimum reporting levels with only a few detections. Some of the organic wastewater compounds detected could have anthropogenic sources that are not biosolids. Concentrations of priority parameters in groundwater varied spatially and temporally but generally were less than Colorado regulatory limits. Concentrations of dissolved nitrate, arsenic, and selenium, in addition to chloride, sulfate, total dissolved solids, boron, iron, manganese, and uranium, in samples from some wells exceeded the Colorado standards. Concentrations of dissolved nitrate (three wells), molybdenum (one well), selenium (two wells), and uranium (one well) in shallow groundwater had significant (alpha = 0.05) upward trends in some parts of the study area. The biosolids-signature results indicate that the aquifers intercepted by the five routinely sampled wells likely have received some recharge through treated (biosolids-applied) fields or biosolids-affected ponds. Adverse effects from this biosolids-related recharge range from few (if any) at one well to large and significantly (alpha = 0.05) increasing nitrate concentrations at another well. A statistical evaluation of five paired wheat-grain samples from treated (biosolids-applied) fields and untreated (control) fields did not indicate any evidence that biosolids applications significantly (alpha = 0.05 or 0.10) increased concentration of any of these constituents in wheat grain. The wheat-grain concentrations from this study were similar to those from other studies for fields in North America where no biosolids were applied. The data for the limited crop samples indicate that biosolids applications are not increasing the concentrations of arsenic, cadmium, copper, lead, mercury, molybdenum, nickel, selenium, sulfur, and zinc in mature wheat grain from the study area.

Colorado↗