USGS Science⌕ Search

SEARCH · USGS Science

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

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Trends in streamflow and concentrations and flux of nutrients and total suspended solids in the Upper White River at Muncie, near Nora, and near Centerton, Indiana

The U.S. Geological Survey (USGS), in cooperation with The Nature Conservancy, completed a study to estimate and assess trends in streamflow and annual mean concentrations and flux of nutrients (nitrate plus nitrite, total Kjeldahl nitrogen, and total phosphorus) and total suspended solids at three USGS streamgages (hereafter referred to as “study gages”) on the Upper White River at Muncie (USGS station 03347000), near Nora (USGS station 03351000), and near Centerton (USGS station 03354000), Indiana. Water-quality data used in the analyses were collected by several agencies between calendar years 1991 and 2017, and streamflow (discharge) data were collected by the USGS. For most of the water-quality constituents, there were suitable data to facilitate an analysis of the 26-year period extending from calendar years 1991 to 2017 (water years 1992 to 2017); however, shorter analytical periods were necessary for total Kjeldahl nitrogen for the study gages at Muncie and near Centerton and for total suspended solids for the study gage near Centerton. Temporal trends in streamflows at the study gages for the period extending from water years 1978 to 2017 were assessed using Exploration and Graphics for RivEr Trends (EGRET) and Mann-Kendall and Pettitt tests. With just one exception, the annual maximum and mean daily streamflows and the annual minimum 7-day mean streamflows at the study gages demonstrated upward trends (increasing streamflows) in the EGRET analyses. The exception was the annual 7-day minimum streamflow at the study gage near Nora, which indicated no trend. Mann-Kendall tests also indicated that the average trend for the annual maximum daily, annual mean daily, and annual 7-day minimum streamflow statistics between water years 1978 and 2017 was upward at each of the study gages; however, only the trends in the annual mean daily streamflows at the study gage at Muncie and the annual maximum daily streamflows at the study gages near Nora and near Centerton were statistically significant at a 0.05 probability level. The Pettitt tests indicated that a statistically significant step trend (abrupt change) in annual mean daily streamflows occurred at each of the study gages around water year 2001. The seasonal distributions of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen concentrations at the study gages were evaluated to identify patterns and other distinguishing characteristics by examining boxplots of concentrations as a function of month of the year. Seasonal distributions of nitrate plus nitrite concentrations and total suspended solids concentrations differed from each other but were generally similar among the three study gages for a given constituent. Median concentrations of nitrate plus nitrite were highest during the January–June months, whereas median concentrations of total suspended solids were highest during June and July. Seasonal distributions of total phosphorus concentrations were similar at the study gages near Nora and near Centerton, but the seasonal distribution was noticeably different at the study gage at Muncie, which had monthly median concentrations that were substantially lower than at the two downstream study gages (near Nora and near Centerton). The seasonal distribution of total Kjeldahl nitrogen concentrations differed in pattern among the three study gages; however, in general, some of the higher monthly median total Kjeldahl nitrogen concentrations at each study gage were associated with the late spring and summer periods. The Weighted Regressions on Time, Discharge, and Season (WRTDS) method implemented in EGRET was used to estimate water-year annual mean daily concentrations and flux of nutrients and total suspended solids, as well as estimates of concentrations and flux that were “normalized” to remove the effect of year-to-year variation in streamflow. The approximate coefficients of determination for the WRTDS regression models ranged from a high of 0.82 for total phosphorus for the study gage near Centerton to a low of 0.19 for nitrate plus nitrite for the study gage near Nora. Loads and yields of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen were estimated for analytical periods consisting of the longest periods of concurrent record at the three study gages. Loads of each of the constituents increased sequentially from the most upstream study gage to the most downstream study gage; however, the same was not true for yields. The highest yields of total suspended solids, total phosphorus, and total Kjeldahl nitrogen occurred at the most upstream study gage (at Muncie); however, the highest yield of nitrate plus nitrite occurred at the most downstream study gage (near Centerton). WRTDS bootstrap tests were used to assess the magnitude, direction, and likelihood of changes in annual flow-normalized mean daily concentrations and flux of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen at the study gages between water years 1997 and 2017. Changes in flow-normalized concentrations and flux of the constituents between water years 1997 and 2017 were mostly downward (decreasing). The exceptions were likely to highly likely upward (increasing) changes in (1) flow-normalized annual mean daily concentration and annual flux for total suspended solids and total phosphorus at the study gage at Muncie, (2) flow-normalized annual mean daily total phosphorus concentration at the study gage near Centerton, (3) flow-normalized annual flux of total phosphorus at the study gage near Centerton, and (4) flow-normalized annual mean daily nitrate plus nitrite concentration at the study gage near Centerton. Although an upward change in flow-normalized nitrate plus nitrite concentrations was likely at the study gage near Centerton, flow-normalized annual flux of nitrate plus nitrite at that study gage was determined to have a highly likely downward change. EGRET and Exploration and Graphics for RivEr Trends Confidence Intervals (EGRETci) analyses can be used to improve our understanding of how concentrations and flux change as functions of time and streamflow, as well as provide information on how the relations between streamflow and constituent concentrations have changed within the calendar year between any 2 years included in the analyses. Examples of those uses, illustrating changes between calendar years 1992 and 2017, were given for total suspended solids concentrations at the study gage near Nora and for nitrate plus nitrite concentrations at the study gage near Centerton.

Indiana↗

Hydrogeology and Aquifer Storage and Recovery Performance in the Upper Floridan Aquifer, Southern Florida

Well construction, hydraulic well test, ambient water-quality, and cycle test data were inventoried and compiled for 30 aquifer storage and recovery facilities constructed in the Floridan aquifer system in southern Florida. Most of the facilities are operated by local municipalities or counties in coastal areas, but five sites are currently being evaluated as part of the Comprehensive Everglades Restoration Plan. The relative performance of all sites with adequate cycle test data was determined, and compared with four hydrogeologic and design factors that may affect recovery efficiency. Testing or operational cycles include recharge, storage, and recovery periods that each last days or months. Cycle test data calculations were made including the potable water (chloride concentration of less than 250 milligrams per liter) recovery efficiency per cycle, total recovery efficiency per cycle, and cumulative potable water recovery efficiencies for all of the cycles at each site. The potable water recovery efficiency is the percentage of the total amount of potable water recharged for each cycle that is recovered; potable water recovery efficiency calculations (per cycle and cumulative) were the primary measures used to evaluate site performance in this study. Total recovery efficiency, which is the percent recovery at the end of each cycle, however, can be substantially higher and is the performance measure normally used in the operation of water-treatment plants. The Upper Floridan aquifer of the Floridan aquifer system currently is being used, or planned for use, at 29 of the aquifer storage and recovery sites. The Upper Floridan aquifer is continuous throughout southern Florida, and its overlying confinement is generally good; however, the aquifer contains brackish to saline ground water that can greatly affect freshwater storage and recovery due to dispersive mixing within the aquifer. The hydrogeology of the Upper Floridan varies in southern Florida; confinement between flow zones is better in southwestern Florida than in southeastern Florida. Vertical hydraulic conductivity in the upper part of the aquifer also may be higher in southeastern Florida because of unconformities present at formation contacts within the aquifer that may be better developed in this area. Recovery efficiencies per cycle varied widely. Eight sites had recovery efficiencies of less than about 10 percent for the first cycle, and three of these sites had not yet achieved recoveries exceeding 10 percent, even after three to five cycles. The highest recovery efficiency achieved per cycle was 94 percent. Three southeastern coastal sites and two southwestern coastal sites have achieved potable water recoveries per cycle exceeding 60 percent. One of the southeastern coastal sites and both of the southwestern coastal sites achieved good recoveries, even with long storage periods (from 174 to 191 days). The high recovery efficiencies for some cycles apparently resulted from water banking?an operational approach whereby an initial cycle with a large recharge volume of water is followed by cycles with much smaller recharge volume. This practice flushes out the aquifer around the well and builds up a buffer zone that can maintain high recovery efficiency in the subsequent cycles. The relative performance of all sites with adequate cycle test data was determined. Performance was arbitrarily grouped into ?high? (greater than 40 percent), ?medium? (between 20 and 40 percent), and ?low? (less than 20 percent) categories based primarily on their cumulative recovery efficiency for the first seven cycles, or projected to seven cycles if fewer cycles were conducted. The ratings of three sites, considered to be borderline, were modified using the overall recharge rate derived from the cumulative recharge volumes. A higher overall recharge rate (greater than 300 million gallons per year) can improve recovery efficiency because of the water-bankin

Scientific Investigations Report↗