USGS Science⌕ Search

SEARCH · USGS Science

Results for “Nature Scientific Reports”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 1,657 records · Page 92Linked to original sources

Hexavalent and total chromium at low reporting concentrations in source-water aquifers and surface waters used for public supply in Illinois, 2013

On the basis of their recent review of the human health effects of hexavalent chromium [Cr(VI)] in public drinking water, the U.S. Environmental Protection Agency is considering the need for Federal regulation of Cr(VI). Presently, only total chromium is regulated, at a Maximum Contaminant Level (MCL) of 100 micrograms per liter (µg/L). The occurrence of Cr(VI) in groundwater and surface waters generally is attributed to industrial sources, but can be of natural origin. California’s recently established MCL for Cr(VI) of 10 µg/L illustrates the drinking-water concerns associated with Cr(VI). To improve understanding of the possible impact of a Cr(VI)-specific standard that approximates the California level on the management of Illinois’ public drinking water, the U.S. Geological Survey, in cooperation with the Illinois Environmental Protection Agency, assessed the occurrence and distribution of Cr(VI) in the State’s public-water supplies. During 2013, untreated water samples were collected to be analyzed for Cr(VI) and total chromium [Cr(T)] at 119 water-supply wells and 32 surface-water intakes; also, 32 treated surface-water samples were collected near the point of treatment and 32 near the furthest point of distribution. Public-supply sample sites were selected by a stratified random method. Samples typically were analyzed within 24 hours of collection at reporting limits of 0.02 µg/L for Cr(VI) and 0.1 µg/L for Cr(T). The occurrence of Cr(VI) was compared with selected geophysical, physical, and sampling factors that might more fully explain its distribution and magnitude of concentrations. The maximum concentration of Cr(VI) in groundwater was 2.1 µg/L. Maximum concentrations in untreated and treated surface water were 0.29 µg/L and 2.4 µg/L, respectively. All sample concentrations were below the California MCL; only 35 percent were below that State’s non-enforceable public health goal of 0.02 of µg/L. Cr(VI) was undetected in 43 percent of untreated groundwater samples, with a median of 0.06 µg/L when detected. All but two (94 percent) of untreated surface-water samples had detections. In untreated surface water, the median concentration was 0.09 µg/L, whereas in treated (tap and distributed) water the median was 0.20 µg/L. Surface waters treated with lime for softening typically had the greatest Cr(VI) concentrations (maximum, 2.4 µg/L; median, 1.2 µg/L). The maximum concentration of Cr(T) in groundwater was 1.8 µg/L. Maximum concentrations in untreated and treated surface water were 1.8 µg/L and 2.5 µg/L, respectively. All sample concentrations were below the Federal MCL. Total chromium was detected in 65 percent of untreated groundwater samples, with a median of 0.40 µg/L, when detected. All but one (97 percent) of untreated surface-water samples had detections. In untreated surface water, the median concentration was 0.40 µg/L, whereas in treated (tap and distributed) water the median was 0.30 µg/L. As with Cr(VI), surface waters treated with lime typically had the greatest Cr(T) concentrations. Examination of factors that might account for or be associated with the occurrence of Cr(VI) in public-supply source waters found few clearly evident factors. Associations in frequencies of occurrence and range of concentrations indicate that surface waters and groundwaters of shallow, unconsolidated, unconfined aquifers, particularly alluvial aquifers, are possibly most commonly affected by anthropogenic sources of Cr(VI). Groundwaters of deep (greater than 500 feet) bedrock aquifers, particularly the Cambrian-Ordovician aquifer system, are possibly most commonly affected by geologic sources of Cr(VI). Additional study, with supporting geologic and geochemical data that were not collected in this study, would be necessary to verify these associations. There was a weak positive relation (ρ = 0.23) between concentrations of Cr(VI) and Cr(T) in untreated water samples, with a much stronger positive relation (ρ = 0.86 and ρ = 0.90, respectively) in samples collected soon after treatment and near the endpoint of distribution. The stronger relation and greater similarity between Cr(VI) and Cr(T) concentrations in treated water samples indicate that Cr(VI) represents a greater proportion of the measured concentrations of Cr(T) in treated waters than in untreated waters. The analysis of spikes and other quality-assurance samples indicate uncertainties associated with obtaining or confirming consistently accurate analytical results for Cr(VI) at near the applied reporting limit of 0.02 µg/L.

Illinois↗

Assessment of groundwater trends near Crex Meadows, Wisconsin

Crex Meadows Wildlife Area (Crex) is a 30,000-acre property in Burnett County, Wisconsin. Crex is managed by the Wisconsin Department of Natural Resources (WDNR) with the goal of providing public recreation opportunities while also protecting the quality of native ecological communities and species on the property. The WDNR’s management strategy includes controlling water levels at flowages in Crex using a system of dikes, water control structures, ditches, and a diversion pump. For the past several decades there has been concern among nearby landowners that the water manage-ment strategy at Crex may be contributing to groundwater flooding in adjacent, privately held properties. This issue has been particularly contentious during periods when regional groundwater elevations are already high. This study was conducted in response to those concerns. For the study, a network of 12 monitoring wells was installed in and to the west of Crex. Groundwater elevations were recorded in the wells before, during, and after water-level changes in the western Crex flowages to assess if groundwater elevations to the west of Crex are detectably affected by the flowage water levels. This study successfully collected groundwater elevations in 11 study wells during a 3-month period in 2019 when water elevations in the Dike 6 flowage and Erickson flowage were lowered and then raised. The data logger at a 12th location failed and no data were recorded. The groundwater elevation trends in these study wells were compared with groundwater elevation trends at a regional U.S. Geological Survey well to provide information for determining if changing the flowage elevations had a noticeable response in the study wells west of Crex Meadows. This analysis was done by (1) evaluating study well groundwater elevation trends compared to the regional well, (2) using a scatter plot of study well and regional well data during raising and lowering periods, (3) assessing horizontal hydraulic gradient data during the study period, and (4) assessing the cumulative departure from the mean groundwater elevation for each well. Overall, regional groundwater elevations had a down-ward trend before and during the flowage lowering period and then had an upward trend during the flowage raising period. This pattern was observed in the regional well and in all the study wells adjacent to and several miles from the flowages. The similarity in patterns indicates that precipitation and regional groundwater flow conditions were the dominant drivers of the system during the study period. The scatter plot and cumulative departure from the mean analysis showed that in addition to regional trends, wells 1, 6, and 7 were likely affected by the changes in the flowage water levels. Overall, at least on the timescale of this study, water management at Crex likely did not have detectable effects on wells outside the Crex property. Wells installed on the Crex property including the wells in the lakebeds of the flowages (wells 1 and 7) and possibly well 6 east of the flowages showed what seems to be minor affects due to water management at Crex.

Wisconsin↗

Appendix C: Summary of Major Metallogenic Belts in Northeast Asia (the Russian Far East, Yakutia, Siberia, Transbaikalia, Northern China, Mongolia, South Korea, and Japan)

The major purposes of this chapter are to provide (1) an overview of the regional geology, tectonics, and metallogenesis of Northeast Asia for readers who are unfamiliar with the region, (2) a general scientific introduction to the succeeding chapters of this volume, and (3) an overview of the methodology of metallogenic and tectonic analysis used in this study. We also describe how a high-quality metallogenic and tectonic analysis, including construction of an associated metallogenic-tectonic model will greatly benefit other mineral resource studies, including synthesis of mineral-deposit models; improve prediction of undiscovered mineral deposit as part of a quantitative mineral-resource-assessment studies; assist land-use and mineral-exploration planning; improve interpretations of the origins of host rocks, mineral deposits, and metallogenic belts, and suggest new research. Research on the metallogenesis and tectonics of such major regions as Northeast Asia (eastern Russia, Mongolia, northern China, South Korea, and Japan) and the Circum-North Pacific (the Russian Far East, Alaska, and the Canadian Cordillera) requires a complex methodology including (1) definitions of key terms, (2) compilation of a regional geologic base map that can be interpreted according to modern tectonic concepts and definitions, (3) compilation of a mineral-deposit database that enables a determination of mineral-deposit models and clarification of the relations of deposits to host rocks and tectonic origins, (4) synthesis of a series of mineral-deposit models that characterize the known mineral deposits and inferred undiscovered deposits in the region, (5) compilation of a series of metallogenic-belt belts constructed on the regional geologic base map, and (6) construction of a unified metallogenic and tectonic model. The summary of regional geology and metallogenesis presented here is based on publications of the major international collaborative studies of the metallogenesis and tectonics of Northeast Asia that have been led by the U.S. Geological Survey (USGS). These studies have produced two broad types of publications (1) a series of regional geologic, mineral-deposit, and metallogenic-belt maps, with companion descriptions of the region, and (2) a suite of metallogenic and tectonic analyses of the same region. The study area consists of eastern Russia (most of eastern Siberia and the Russian Far East), Mongolia, northern China, South Korea, Japan, and adjacent offshore areas. The major cooperative agencies are the Russian Academy of Sciences; the Academy of Sciences of the Sakha Republic (Yakutia); VNIIOkeangeologia and Ministry of Natural Resources of the Russian Federation; the Mongolian Academy of Sciences; the Mongolian University of Science and Technology; the Mongolian National University; Jilin University, Changchun, People?s Republic of China, the China Geological Survey; the Korea Institute of Geosciences and Mineral Resources; the Geological Survey of Japan/AIST; the University of Texas, Arlington, and the U.S. Geological Survey (USGS). This study builds on and extends the data and interpretations from a previous project on the Major Mineral Deposits, Metallogenesis, and Tectonics of the Russian Far East, Alaska, and the Canadian Cordillera conducted by the USGS, the Russian Academy of Sciences, the Alaska Division of Geological and Geophysical Surveys, and the Geological Survey of Canada. The major products of this project were summarized by Naumova and others (2006) and are described in appendix A.

Professional Paper↗

Groundwater-flow and land-subsidence model of Antelope Valley, California

Antelope Valley, California, is a topographically closed basin in the western part of the Mojave Desert, about 50 miles northeast of Los Angeles. The Antelope Valley groundwater basin is about 940 square miles and is separated from the northern part of Antelope Valley by faults and low-lying hills. Prior to 1972, groundwater provided more than 90 percent of the total water supply in the valley; since 1972, it has provided between 50 and 90 percent. Most groundwater pumping in the valley occurs in the Antelope Valley groundwater basin, which includes the rapidly growing cities of Lancaster and Palmdale. Groundwater-level declines of more than 270 feet in some parts of the groundwater basin have resulted in an increase in pumping lifts, reduced well efficiency, and land subsidence of more than 6 feet in some areas. Future urban growth and limits on the supply of imported water may increase reliance on groundwater. In 2011, the Los Angeles County Superior Court of California ruled that the Antelope Valley groundwater basin is in overdraft—groundwater extractions are in excess of the Court-defined safe yield of the groundwater basin. The Court determined that the safe yield of the adjudicated area of the basin was 110,000 acre-feet per year (acre-ft/yr). Natural recharge is an important component of total groundwater recharge in Antelope Valley; however, the exact quantity and distribution of natural recharge, primarily in the form of mountain-front recharge, is uncertain, with total estimates ranging from 30,000 to 160,000 acre-ft/yr. Technical experts, retained by parties to the adjudication, used 60,000 acre-ft/yr to estimate the sustainable yield of the basin, and this value was used in this study. In order to better understand the uncertainty associated with natural recharge and to provide a tool to aid in groundwater management, a numerical model of groundwater flow and land subsidence in the Antelope Valley groundwater basin was developed using old and new geohydrologic information. The groundwater-flow system consists of three aquifers: the upper, middle, and lower aquifers. The three aquifers, which were identified on the basis of the hydrologic properties, age, and depth of the unconsolidated deposits, consist of gravel, sand, silt, and clay alluvial deposits and clay and silty clay lacustrine deposits. Prior to groundwater development in the valley, recharge was primarily the infiltration of runoff from the surrounding mountains. Groundwater flowed from the recharge areas to discharge areas around the playas where it discharged from the aquifer system as either evapotranspiration or from springs. Partial barriers to horizontal groundwater flow, such as faults, have been identified in the groundwater basin. Water-level declines owing to groundwater development have eliminated the natural sources of discharge, and pumping for agricultural and urban uses have become the primary source of discharge from the groundwater system. Infiltration of return flow from agricultural irrigation has become an important source of recharge to the aquifer system. The groundwater-flow model of the basin was discretized horizontally into a grid of 130 rows and 118 columns of square cells 1 kilometer (0.621 mile) on a side, and vertically into four layers representing the upper (two layers), middle (one layer), and lower (one layer) aquifers. Faults that were thought to act as horizontal-flow barriers were simulated in the model. The model was calibrated to simulate steady-state conditions, represented by 1915 water levels and transient-state conditions during 1915–95, by using water-level and subsidence data. Initial estimates of the aquifer-system properties and stresses were obtained from a previously published numerical model of the Antelope Valley groundwater basin; estimates also were obtained from recently collected hydrologic data and from results of simulations of groundwater-flow and land-subsidence models of the Edwards Air Force Base area. Some of these initial estimates were modified during model calibration. Groundwater pumpage for agriculture was estimated on the basis of irrigated crop acreage and crop consumptive-use data. Pumpage for public supply, which is metered, was compiled and entered into a database used for this study. Estimated annual agricultural pumpage peaked at 395,000 acre-feet (acre-ft) in 1951 and then declined because of declining agricultural production. Recharge from irrigation return flows was assumed to be 30 percent of agricultural pumpage; delays associated with return flow moving through the unsaturated zone were also simulated. The annual quantity of mountain-front recharge initially was based on estimates from previous studies. The model was calibrated using the PEST software suite; prior information from the area was incorporated through the use of Tikhonov regularization. During model calibration, the estimated mountain-front recharge was reduced from the previous estimate of 30,300 acre-ft/yr to 29,150 acre-ft/yr. Results of the simulations using the calibrated model indicate that simulated groundwater pumpage exceeded recharge in most years, resulting in an estimated cumulative depletion in groundwater storage of 8,700,000 acre-ft during the transient-simulation period (1915–2005). About 15,000,000 acre-ft of cumulative groundwater pumpage was simulated during the transient-simulation period (1915–2005), reaching a maximum rate of about 400,000 acre-ft/yr in 1951. Groundwater pumpage resulted in simulated hydraulic heads declining by more than 150 feet (ft) compared to 1915 conditions in agricultural areas. The decline in hydraulic head in the groundwater basin is the result of this depletion of groundwater storage. In turn, the simulated decline in hydraulic head in the groundwater basin has resulted in the decrease in natural discharge from the basin and has caused compaction of aquitards, resulting in land subsidence. The areal distribution of total simulated land subsidence for 2005, after about 90 years of groundwater development, indicates that land subsidence occurred throughout almost the entire Lancaster subbasin, with a maximum of about 9.4 ft in the central and eastern parts of the subbasin. An important objective of this study was to systematically address the uncertainty in estimates of natural recharge and related aquifer parameters by using the groundwater-flow and land-subsidence model with observational data and expert knowledge. After the model was calibrated to the observations and a reasonable parameter set obtained, the parameter null space—parameter values that do not appreciably affect the model calibration but may have importance for prediction—was identified. The effect of parameter uncertainty on the estimation of mountain-front recharge was addressed using the Null-Space Monte Carlo method. The Pareto trade-off method of visualizing uncertainty was also used to portray the reasonableness of larger natural-recharge rates. Results indicate that the total mountain-front recharge likely ranges between 28,000 and 44,000 acre-ft/yr, which is appreciably less than published estimates of 60,000 acre-ft/yr. Additionally, expected errors associated with agricultural pumpage estimates used in this study were found to have relatively little effect on the estimates of mountain-front recharge, reflecting the difficulty in increasing recharge through manipulation of other components of the water budget. The calibrated model was used to simulate the response of the aquifer to potential future pumping scenarios: (1) no change in the distribution of pumpage, or status quo; (2) redistribution of pumpage; and (3) artificial recharge. All three of these scenarios specify a total pumpage throughout the Antelope Valley of 110,000 acre-ft/yr according to the safe yield value ruled by the Los Angeles County Superior Court of California. This reduction in groundwater pumpage is assumed uniform throughout the basin, based on a 10-percent reduction of the total pumpage in 2005 to achieve the 110,000 acre-ft/yr level. The calibrated Antelope Valley groundwater-flow and land-subsidence model was used to simulate the hydrologic effects of the three groundwater-management scenarios during a 50-year period by using the reduced, temporally constant, pumpage distribution. Results from the first scenario indicated that the total drawdown observed since predevelopment would continue, with values exceeding 325 ft near Palmdale; consequently, land subsidence would also continue, with additional subsidence (since 2005) exceeding 3 ft in the central part of the Lancaster subbasin. The second scenario evaluated redistributing pumpage from areas in the Lancaster subbasin where simulated hydraulic-head declines were the greatest to areas where declines were smallest. Neither a formal optimization algorithm nor water-rights allocations were considered when redistributing the pumpage. Results indicated that hydraulic heads near Palmdale, where the pumpage was reduced, would recover by about 200 ft compared to 2005 conditions, with only 30 ft of additional drawdown in the northwestern part of the Lancaster subbasin, where the pumpage was increased. The magnitude of the simulated additional land subsidence decreased slightly compared to the first, status quo, scenario but land subsidence continued to be simulated throughout most of the northern part of the Lancaster subbasin. The third scenario consisted of two artificial-recharge simulations along the Upper Amargosa Creek channel and at a site located north of Antelope Buttes. Results indicate that applying artificial recharge at these sites would yield continued drawdowns and associated land subsidence. However, the magnitudes of drawdown and subsidence would be smaller than those simulated in the status quo scenario, indicating that artificial-recharge operations in the Antelope Valley could be expected to reduce the magnitude and extent of continued water-level declines and associated land subsidence.

California↗

The Connecticut Streamflow and Sustainable Water Use Estimator—A decision-support tool to estimate water availability at ungaged stream locations in Connecticut

Freshwater streams in Connecticut are subject to many competing demands, including public water supply; agricultural, commercial, and industrial water use; and ecosystem and habitat needs. In recent years, drought has further stressed Connecticut’s water resources. To sustainably allocate and manage water resources among these competing uses, Federal, State, and local water-resource managers require data and modeling tools to estimate the water availability at a variety of temporal and spatial scales for planning purposes. The Connecticut Streamflow and Sustainable Water Use Estimator (CT SSWUE), developed by the U.S. Geological Survey in cooperation with the Connecticut Department of Energy and Environmental Protection, is a decision-support tool for estimating daily unaltered streamflow and sustainable water use at ungaged sites in Connecticut. The CT SSWUE estimates unaltered daily mean streamflow and water-use-adjusted streamflow for the period from October 1, 1960, to September 30, 2015, and the monthly sustainable net withdrawal at ungaged sites in Connecticut. Unaltered streamflow is the estimated daily mean streamflow in a drainage basin in the absence of any water withdrawals or wastewater discharges and with minimal human development. Sustainable net withdrawal is the maximum net withdrawal (withdrawal minus wastewater discharges) that can be drawn from a basin without critically depleting the water available through natural streamflow patterns. Sustainable net withdrawal is defined for this study as the difference between the unaltered daily mean streamflow and a user-defined target minimum streamflow. Weighted least squares and Tobit regression techniques were used to develop equations for estimating streamflow at ungaged sites at 19 streamflow quantiles with exceedance probabilities ranging from 0.005 to 99.995 percent. Regressions were based on streamflow quantiles and basin characteristics from 36 reference streamgages in and around Connecticut. Four basin characteristics—drainage area, mean of the soil permeability, mean of the average annual precipitation, and ratio of the length of streams that overlay sand and gravel deposits to the total length of streams in the basin—are used as explanatory variables in the equations. At an ungaged site, interpolation between the streamflow quantiles estimated from the regression equations produces a continuous flow-duration curve. A time series of daily mean streamflow at an ungaged site is then estimated by assuming that for each day, the streamflow quantile occurs on the same date at both a reference streamgage and the ungaged site. In a remove-one cross validation, estimated unaltered daily mean streamflow agreed well with observed values at reference streamgages, with a few exceptions. Nash Sutcliffe efficiency ranged from −0.43 to 0.97 with a median value of 0.88. The normalized root-mean-square error ranged from 16.6 to 120.4 percent with a median value of 34.5 percent. An empirical method for estimating 95-percent prediction intervals for unaltered daily and monthly mean streamflow was developed and tested by using the cross-validation data. Prediction intervals for unaltered daily mean streamflow at the cross-validation reference streamgages performed well in most cases. Gaged streamflow values from the cross-validation data fell within the prediction intervals a median 96.6 percent of the time for daily mean time series and 93.9 percent of the time for monthly mean time series. The CT SSWUE computes water-use-adjusted streamflow using spatially referenced water-use information provided by the Connecticut Department of Energy and Environmental Protection. Available water-use information included permitted and registered water withdrawals and permitted wastewater discharges during 1998 to 2015 for the Thames River Basin and central coastal drainage basins. Water-use information was incorporated into the U.S. Geological Survey StreamStats web application for Connecticut and can be used for computing water-use-adjusted streamflow and sustainable net withdrawal at selected points of interest. Altered daily streamflow is computed by applying average daily withdrawals and wastewater discharges to the water balance equation. Average daily surface water withdrawals and wastewater discharges are applied directly to the daily water balance equation. Time-lagged alterations on streamflow from groundwater withdrawals or wastewater discharges are estimated by using a response-coefficient method developed from results of previously published, calibrated groundwater models.

Connecticut↗

Simulation of the ground-water-flow system in the Kalamazoo County area, Michigan

A ground-water-flow model was developed to investigate the ground-water resources of Kalamazoo County. Ground water is widely used as a source of water for drinking and industry in Kalamazoo County and the surrounding area. Additionally, lakes and streams are valued for their recreational and aesthetic uses. Stresses on the ground-water system, both natural and human-induced, have raised concerns about the long-term availability of ground water for people to use and for replenishment of lakes and streams. Potential changes in these stresses, including withdrawals and recharge, were simulated using a ground-water-flow model. Simulations included steady-state conditions (in which stresses remained constant and changes in storage were not included) and transient conditions (in which stresses changed in seasonal and monthly time scales and storage within the system was included). Steady-state simulations were used to investigate the long-term effects on water levels and streamflow of a reduction in recharge or an increase in pumping to projected 2010 withdrawal rates, withdrawal and application of water for irrigation, and a reduction in recharge in urban areas caused by impervious surfaces. Transient simulations were used to investigate changes in withdrawals to match seasonal and monthly patterns under various recharge conditions, and the potential effects of the use of water for irrigation over the summer months. With a reduction in recharge, simulated water levels declined over most of the model area in Kalamazoo County; with an increase in pumping, water levels declined primarily near pumping centers. Because withdrawals by wells intercept water that would have discharged possibly to a stream or lake, model simulations indicated that streamflow was reduced with increased withdrawals. With withdrawal and consumption of water for irrigation, simulated water levels declined. Assuming a reduction in recharge due to urbanization, water levels declined and flow to streams was reduced based on steady-state simulation results. Transient results indicated a reduction of water levels with the simulated use of water for irrigation over the summer months. Generally the transient simulation with recharge only in the winter provided the best fit to observed water levels collected during synoptic water-level measurements in some wells and to the trends observed in water levels for other wells. Analysis of the regional hydrologic budgets provides an increased understanding of water movement within the ground-water-flow system in Kalamazoo County. Budgets for the steady-state simulations indicated that with reduced recharge, less water was available for streamflow and less water left the model area through the model boundaries. Similarly, with an increase in pumping rates, less water was available to enter streams and become streamflow. When recharge was assumed to remain constant and when it was allowed to vary throughout the year, the amount of water that entered storage was greater than that which left storage. However, when recharge was distributed through October?May only or when recharge rates were reduced from October to May, the amount of water that entered storage was less than that which left storage. Thus, on the basis of model simulations, with reduced recharge or increased withdrawals, water must come from storage, rivers, or from ground-flow-system boundaries to meet withdrawal demands.

Michigan↗

Assessment of an in-channel redistribution technique for large woody debris management in Locust Creek, Linn County, Missouri

The U.S. Geological Survey, in cooperation with the Missouri Department of Conservation and Missouri Department of Natural Resources, completed a study to assess a mechanical redistribution technique used for the management of large woody debris (LWD) jams in Locust Creek within Pershing State Park and Fountain Grove Conservation Area, Linn County, Missouri. Extensive LWD jams were treated from 1996 to 2009 using a low-impact technique in which LWD from the jams was redistributed to reopen the channel and to mimic the natural geomorphic process of channel migration and adjustment to an obstruction. The scope of the study included the comparison of selected channel geometry characteristics and bed material particle-size distribution in seven LWD treatment reaches with that of adjacent untreated reaches (unaffected by LWD accumulations) of Locust Creek. A comparison of 1996 and 2015 survey cross sections in treated and untreated reaches and photograph documentation were used to assess channel geomorphic change and the stability of redistributed LWD. The physical characteristics of LWD within jams present in the study reach during 2015–16 also were documented. Based on the general lack of differences in channel metrics between treated and untreated reaches, it can be concluded that the mechanical redistribution technique has been an effective treatment of extensive LWD jams in Locust Creek. Channel alterations, including aggradation, streamflow piracy, and diversions, have resulted in temporal and spatial changes in the Locust Creek channel that may affect future applications of the redistribution technique in Pershing State Park. The redistribution technique was used to effectively manage LWD in Locust Creek at a potentially lower financial cost and reduced environmental disturbance than the complete removal of LWD. A comparison of four channel metrics (bankfull cross-sectional area, channel width, streamflow capacity, and width-depth ratio) for individual treatment reaches with adjacent untreated reaches indicated no statistically significant difference in most comparisons. Where statistically significant differences in channel metrics were determined between individual reaches, the channel metrics in treatment reaches were significantly less than adjacent untreated reaches in some comparisons, and significantly greater than adjacent untreated reaches in others. Without immediate posttreatment cross sections in treated and untreated reaches for comparison, it is impossible to say with certainty that a lack of significant differences in channel metrics is a result of posttreatment channel adjustment or, conversely, that any significant differences that remain are a result of the treatment of LWD. Characteristics of LWD in accumulations sampled within the study area in 2015 indicate that most sampled pieces were in the 1–2 foot diameter size class, the 5–16 foot length class, and the advanced decay class. Most of documented LWD pieces were loose and not buried, about 20 percent on average had a root wad attached, and about 6.5 percent on average were sawn logs. Most of sampled material was less than one-half of the bankfull channel width, indicating it was easily transportable, and the advanced decay class of material entering the study area indicated that it was likely sourced from outside of Pershing State Park. Redistributed LWD associated with treatment seems to be intact in the 1996 treated reaches from direct observation and from inference because there was net channel aggradation between 1996 and 2015 in comparison surveys. The change in channel area resulting from aggradation in time (1996 to 2015) in treated and untreated reaches exceeded the differences in channel characteristics between the treated and untreated channels in 2015 surveys.

Missouri↗

Hydrogeochemical investigation of the Standard Mine Vicinity, Upper Elk Creek Basin, Colorado

Ground- and surface-water samples were collected in the vicinity of the Standard Mine in west-central Colorado in order to characterize the local ground-water flow system, determine metal concentrations in local ground water, and better understand factors controlling the discharge of metal-rich waters from the mine. The sampling program included a one-time sampling of springs, mine adits, and exploration pits in Elk Basin and Redwell Basin; repeated sampling throughout one year of Standard Mine Level 1 discharge and Elk Creek near its confluence with Coal Creek; and a one-time sampling of underground sites in Levels 3 and 5 of the Standard Mine. Samples were analyzed for major ions and trace elements, stable isotopes of hydrogen (2H/1H) and oxygen (18O/16O), strontium isotopes, and tritium and dissolved noble gases (including helium isotopes) for tritium/helium-3 age dating. No clear correlations were observed between natural ground-water discharge locations and map-scale faults and lithology. Surface observations and the location of ground-water discharge suggest that simple topography, rather than large-scale geologic features, primarily controls the occurrence and flow of shallow ground water in Elk Basin. Discrete inflows from cross faults or other features were not observed in Levels 3 and 5 of the Standard Mine. Instead, water entered the mine as relatively persistent dripping from gouge and breccia within the Standard fault, which both tunnels follow. Therefore, the Standard fault itself is probably the main pathway of ground-water flow from the shallow subsurface to the mine workings. Low pH (as low as 3.2) and elevated concentrations of zinc, lead, cadmium, copper, and manganese (commonly exceeding water-quality standards for Elk Creek) were measured in samples located within or immediately downgradient of areas where sulfides are abundant, including the Standard fault, the Elk Lode portal, and the breccia pipe in Redwell Basin. Concentrations of these metals were typically low and pH values were circumneutral at surrounding locations. Metal concentrations in samples collected from underground workings in the Standard Mine were also generally higher than in samples collected at aboveground sites located outside of sulfide-rich areas. Metal concentrations in discharge from the Level 1 tunnel were among the highest measured in Elk Basin. All of these observations suggest that sulfide-rich mineralized rock is the primary control on dissolved metal concentrations and pH in ground water in the Standard Mine vicinity. Waste-rock piles apparently exert another major control on metal concentrations and pH; the lowest pH and highest metal concentrations typically are found in discharge from waste-rock piles. Concentrations of several chemical constituents along with strontium isotope data indicate that none of the sampled waters could have been the primary source of metals in discharge from Level 1. Therefore, this study did not identify the primary source location for metals in Level 1 discharge. Possible sources must be located below Levels 3 and 5 or farther back into the mountainside than the ends of Levels 3 and 5. Apparent tritium/helium-3 ground-water ages ranged from 0 to 9 yr, and a considerable majority were <1 yr. Tritium data and computed initial tritium values (measured tritium plus measured tritiogenic helium-3) suggest that much of the ground water in the Standard Mine vicinity was weeks to months old rather than years old. Tritium, d2H, and d18O data from water entering into and discharging from the Standard Mine displayed spatial and temporal patterns indicating that these tracers were influenced by seasonal variations in their concentration in precipitation. The tracer data therefore suggest that ground water entering into and discharging from the Standard Mine was largely composed of water <1 yr old. Pronounced seasonal variations in geochemistry in Level 1 discharge also are consistent with short r

Colorado↗

Development and analysis of regional curves for streams in the non-urban valley and ridge physiographic province, Maryland, Virginia, and West Virginia

Regression relations for bankfull stream characteristics based on drainage area (often called 'regional curves') are used in natural stream channel design to verify field determinations of bankfull discharge and stream channel characteristics. Bankfull stream characteristics were assessed for stream reaches at 41 streamflow-gaging stations in the Valley and Ridge Physiographic Province in Maryland, Virginia, and West Virginia. Data collected included bankfull cross-sectional geometry, flood plain geometry, and longitudinal profile data. In addition, particle-size distributions of streambed material were determined and data on basin characteristics were compiled for each reach. Regional curves were developed for bankfull cross-sectional area, width, and discharge with R2 values of 0.95, 0.89, 0.87, and 0.91, respectively. Examination of the regional curves residuals indicates that there is more variability in bankfull cross-sectional area, width, and discharge for smaller streams than for larger streams. In contrast, there is more variability for bankfull mean depth for larger streams than for smaller streams. Geographic analysis of regional curve residuals indicated that there were no further subdivisions within the Valley and Ridge Physiographic Province in the three-state study area for which individual sets of regional curves should be developed. In addition, two separate sets of regional curves were developed with data from the 41 sites to examine potential differences in the relations between the southern (n = 9) and central (n = 32) sections of the province. There were differences in slope and intercept between the two bankfull discharge test relations and a difference in intercept for the width test relations at the 95-percent confidence level. However, the results of this analysis were inconclusive and therefore one set of regional curves for the study area is presented in this report. The regional curves were compared to regression models developed from similar data collected in the Pennsylvania and Maryland portions of the province. No statistical difference in the slope or intercept of regression lines of the three data sets was detected for any of the four bankfull parameters at the 95-percent confidence level. Basin characteristics such as percentage of basin forested (percent forested) and percentage of basin underlain by carbonate bedrock (percent carbonate) were analyzed to evaluate variability among regression points. Multivariate regression relations including explanatory terms for percent carbonate and drainage area produced higher R2 values than the regional curves for bankfull cross-sectional area (R2 = 0.95), bankfull width (R2 = 0.92), and bankfull discharge (R2 = 0.93). There was no improvement for the bankfull mean depth relation from adding the additional term. Inclusion of the other basin characteristics in multivariate relations did not improve the regression models. Regression models developed for the 1.5-year discharge for all streamflow-gaging stations with peak discharge data throughout Virginia (n = 486) and throughout the Valley and Ridge Physiographic Province in Virginia (n = 147) were compared to the regional curve relating bankfull discharge to drainage area. A similar trend in decreasing variability with increasing drainage area was observed for the 1.5-year discharge for all stations in Virginia . This indicates that the change in variability observed in the discharge regional curve likely would exist with a larger data set. There was no statistical difference at the 95-percent confidence level between regression relations for the southern section of the province (n = 40) and the central section (n = 107). This finding supports maintaining only one set of regional curves for the study area. Not all of the variability in the regional curves is explained by drainage area alone. Causes of the remaining variability likely vary among study sites. Users of the regional curves developed in this study are cautioned that, because of inherent variability in hydrologic data, a thorough analysis of the fluvial system and the drainage basin is necessary when evaluating the bankfull characteristics of a stream.

Maryland, Virginia, West Virginia↗

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

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

Ohio↗

Questa baseline and pre-mining ground-water quality investigation. 19. Leaching characteristics of composited materials from mine waste-rock piles and naturally altered areas near Questa, New Mexico

The goal of this study is to compare and contrast the leachability of metals and the acidity from individual mine waste-rock piles and natural erosional scars in the study area near Questa, New Mexico. Surficial multi-increment (composite) samples less than 2 millimeters in diameter from five waste-rock piles, nine erosional-scar areas, a less-altered site, and a tailings slurry-pipe sample were analyzed for bulk chemistry and mineralogy and subjected to two back-to-back leaching procedures. The first leaching procedure, the U.S. Geological Survey Field Leach Test (FLT), is a short-duration leach (5-minute shaking and 10-minute settling) and is intended to leach readily soluble materials. The FLT was immediately followed by an 18-hour, end-over-end rotation leaching procedure. Comparison of results from the back-to-back leaching procedures can provide information about reactions that may take place upon migration of leachates through changing geochemical conditions (for example, pH changes), both within the waste-rock and scar materials and away from the source materials. For the scar leachates, the concentrations of leachable metals varied substantially between the scar areas sampled. The scar leachates have low pH (pH 3.2-4.1). Under these low-pH conditions, cationic metals are solubilized and mobile, but anionic species, such as molybdenum, are less soluble and less mobile. Generally, metal concentrations in the waste-rock leachates did not exceed the upper range of those metal concentrations in the erosional-scar leachates. One exception is molybdenum, which is notably higher in the waste-rock leachates compared with the scar leachates. Most of the waste-rock leachates were at least mildly acidic (pH 3.0-6.2). The pH values in the waste-rock leachates span a large pH range that includes some pH-dependent solubility and metal-attenuation reactions. An increase in pH with leaching time and agitation indicates that there is pH-buffering capacity in some of the waste-rock piles. As pH increased in the waste-pile leachates, concentrations of several metals decreased with increasing time and agitation. Similar pH-dependent reactions may take place upon migration of the leachates in the waste-rock piles. Bulk chemistry, mineralogy, and leachate sulfur-isotope data indicate that the Capulin and Sugar Shack West waste-rock piles are compositionally different from the younger Sugar Shack South, Sugar Shack Middle, and Old Sulphur Gulch piles. The Capulin and Sugar Shack West piles have the lowest-pH leachates (pH 3.0-4.1) of the waste-pile samples, and the source material for the Capulin and Sugar Shack West piles appears to be similar to the source material for the erosional-scar areas. Calcite dissolution, in addition to gypsum dissolution, appears to produce the calcium and sulfate concentrations in leachates from the Sugar Shack South, Sugar Shack Middle, and Old Sulphur Gulch piles.

Scientific Investigations Report↗

Design and analysis of a natural-gradient ground-water tracer test in a freshwater tidal wetland, West Branch Canal Creek, Aberdeen Proving Ground, Maryland

A natural-gradient ground-water tracer test was designed and conducted in a tidal freshwater wetland at West Branch Canal Creek, Aberdeen Proving Ground, Maryland. The objectives of the test were to characterize solute transport at the site, obtain data to more accurately determine the ground-water velocity in the upper wetland sediments, and to compare a conservative, ionic tracer (bromide) to a volatile tracer (sulfur hexafluoride) to ascertain whether volatilization could be an important process in attenuating volatile organic compounds in the ground water. The tracer test was conducted within the upper peat unit of a layer of wetland sediments that also includes a lower clayey unit; the combined layer overlies an aquifer. The area selected for the test was thought to have an above-average rate of ground-water discharge based on ground-water head distributions and near-surface detections of volatile organic compounds measured in previous studies. Because ground-water velocities in the wetland sediments were expected to be slow compared to the underlying aquifer, the test was designed to be conducted on a small scale. Ninety-seven ?-inch-diameter inverted-screen stainless-steel piezometers were installed in a cylindrical array within approximately 25 cubic feet (2.3 cubic meters) of wetland sediments, in an area with a vertically upward hydraulic gradient. Fluorescein dye was used to qualitatively evaluate the hydrologic integrity of the tracer array before the start of the tracer test, including verifying the absence of hydraulic short-circuiting due to nonnatural vertical conduits potentially created during piezometer installation. Bromide and sulfur hexafluoride tracers (0.139 liter of solution containing 100,000 milligrams per liter of bromide ion and 23.3 milligrams per liter of sulfur hexafluoride) were co-injected and monitored to generate a dataset that could be used to evaluate solute transport in three dimensions. Piezometers were sampled 2 to 15 times each, from July 1998 through September 1999, to assess background conditions and monitor tracer movement. During the test, 644 samples were analyzed for fluorescein, 617 samples were analyzed for bromide with an ion-selective electrode, 213 samples were analyzed for bromide with colorimetric methods, and 603 samples were analyzed for sulfur hexafluoride, including samples collected prior to tracer injection to determine background concentrations. Additional samples were analyzed for volatile organic compounds (96 samples) and methane (37 samples) to determine the distribution of these contaminants and the extent of methanogenic conditions within the tracer array; however, these data were not used for the analysis of the test. During the tracer test, the fluorescein dye, bromide, and sulfur hexafluoride were transported predominantly in the upward direction, although all three tracers also moved outward in all directions from the injection point, and it is likely that some tracer mass moved beyond the lateral edges of the array. An analysis of the tracer-test data was performed through the use of breakthrough curves and isoconcentration contour plots. Results show that movement of the fluorescein dye, a non-conservative tracer, was retarded compared to the other two tracers, likely as a result of sorption onto the wetland sediments. Suspected loss of tracer mass along the lateral edges of the array prevented a straightforward quantitative analysis of tracer transport and ground-water velocity from the bromide and sulfur-hexafluoride data. In addition, the initial density of the bromide/sulfur hexafluoride solution (calculated to be 1.097 grams per milli2 Ground-Water Tracer Test, West Branch Canal Creek, Aberdeen Proving Ground, MD liter) could have caused the solution to sink below the injection point before undergoing dilution and moving back up into the array. For these reasons, the data analysis in this report was performed largely through qualitative method

Scientific Investigations Report↗

Overview of the ARkStorm scenario

The U.S. Geological Survey, Multi Hazards Demonstration Project (MHDP) uses hazards science to improve resiliency of communities to natural disasters including earthquakes, tsunamis, wildfires, landslides, floods and coastal erosion. The project engages emergency planners, businesses, universities, government agencies, and others in preparing for major natural disasters. The project also helps to set research goals and provides decision-making information for loss reduction and improved resiliency. The first public product of the MHDP was the ShakeOut Earthquake Scenario published in May 2008. This detailed depiction of a hypothetical magnitude 7.8 earthquake on the San Andreas Fault in southern California served as the centerpiece of the largest earthquake drill in United States history, involving over 5,000 emergency responders and the participation of over 5.5 million citizens. This document summarizes the next major public project for MHDP, a winter storm scenario called ARkStorm (for Atmospheric River 1,000). Experts have designed a large, scientifically realistic meteorological event followed by an examination of the secondary hazards (for example, landslides and flooding), physical damages to the built environment, and social and economic consequences. The hypothetical storm depicted here would strike the U.S. West Coast and be similar to the intense California winter storms of 1861 and 1862 that left the central valley of California impassible. The storm is estimated to produce precipitation that in many places exceeds levels only experienced on average once every 500 to 1,000 years. Extensive flooding results. In many cases flooding overwhelms the state’s flood-protection system, which is typically designed to resist 100- to 200-year runoffs. The Central Valley experiences hypothetical flooding 300 miles long and 20 or more miles wide. Serious flooding also occurs in Orange County, Los Angeles County, San Diego, the San Francisco Bay area, and other coastal communities. Windspeeds in some places reach 125 miles per hour, hurricane-force winds. Across wider areas of the state, winds reach 60 miles per hour. Hundreds of landslides damage roads, highways, and homes. Property damage exceeds $ 300 billion, most from flooding. Demand surge (an increase in labor rates and other repair costs after major natural disasters) could increase property losses by 20 percent. Agricultural losses and other costs to repair lifelines, dewater (drain) flooded islands, and repair damage from landslides, brings the total direct property loss to nearly $ 400 billion, of which $ 20 to $ 30 billion would be recoverable through public and commercial insurance. Power, water, sewer, and other lifelines experience damage that takes weeks or months to restore. Flooding evacuation could involve 1.5 million residents in the inland region and delta counties. Business interruption costs reach $ 325 billion in addition to the $ 400 billion property repair costs, meaning that an ARkStorm could cost on the order of $ 725 billion, which is nearly 3 times the loss deemed to be realistic by the ShakeOut authors for a severe southern California earthquake, an event with roughly the same annual occurrence probability. The ARkStorm has several public policy implications: (1) An ARkStorm raises serious questions about the ability of existing federal, state, and local disaster planning to handle a disaster of this magnitude. (2) A core policy issue raised is whether to pay now to mitigate, or pay a lot more later for recovery. (3) Innovative financing solutions are likely to be needed to avoid fiscal crisis and adequately fund response and recovery costs from a similar, real, disaster. (4) Responders and government managers at all levels could be encouraged to conduct risk assessments, and devise the full spectrum of exercises, to exercise ability of their plans to address a similar event. (5) ARkStorm can be a reference point for application of Federal Emergency Management Agency (FEMA) and California Emergency Management Agency guidance connecting federal, state and local natural hazards mapping and mitigation planning under the National Flood Insurance Plan and Disaster Mitigation Act of 2000. (6) Common messages to educate the public about the risk of such an extreme disaster as the ARkStorm scenario could be developed and consistently communicated to facilitate policy formulation and transformation. These impacts were estimated by a team of 117 scientists, engineers, public-policy experts, insurance experts, and employees of the affected lifelines. In many aspects the ARkStorm produced new science, such as the model of coastal inundation. The products of the ARkStorm are intended for use by emergency planners, utility operators, policymakers, and others to inform preparedness plans and to enhance resiliency.

California↗

Status of Water Levels and Selected Water-Quality Conditions in the Sparta-Memphis Aquifer in Arkansas and the Status of Water Levels in the Sparta Aquifer in Louisiana, Spring 2005

The U.S. Geological Survey in cooperation with the Arkansas Natural Resources Commission, the Arkansas Geological Commission, and the Louisiana Department of Transportation and Development has monitored water levels in the Sparta Sand of Claiborne Group and Memphis Sand of Claiborne Group since the 1920's. Ground-water withdrawals have increased while water levels have declined since monitoring was initiated. This report has been produced to describe ground-water levels in the aquifers in the Sparta Sand and Memphis Sand and provide information for the management of this valuable resource. The 2005 potentiometric-surface map of the aquifers in the Sparta Sand and Memphis Sand was constructed using water-level data collected in 333 wells in Arkansas and 120 wells in Louisiana during the spring of 2005. The highest water-level altitude measured in Arkansas was 327 feet above National Geodetic Vertical Datum of 1929 located in Grant County in the outcrop at the western boundary of the study area; the lowest water-level altitude was 189 feet below National Geodetic Vertical Datum of 1929 in Union County. The highest water-level altitude measured in Louisiana was 246 feet above National Geodetic Vertical Datum of 1929 located in Bossier Parish in the outcrop area near the western boundary of the study area; the lowest water-level altitude was 226 feet below National Geodetic Vertical Datum of 1929 in central Ouachita Parish. Three large depressions centered in Columbia, Jefferson, and Union Counties in Arkansas are the result of large withdrawals for industrial and public supplies. In Louisiana, three major pumping centers are in Ouachita, Jackson, and Lincoln Parishes. Water withdrawals from these major pumping centers primarily is used for industrial and public-supply purposes. Withdrawals from Ouachita and Lincoln Parishes and Union County, Arkansas, primarily for industrial purposes, have caused the resulting cones of depression to coalesce so that the -40 foot potentiometric contour encircles the three pumping centers. Seven smaller depressions are evident on the 2005 Sparta-Memphis potentiometric-surface map located in Webster and Winn Parishes, Louisiana, and Calhoun, Cleveland, western Columbia, Desha, and Lafayette Counties, Arkansas. The depression in Calhoun County initially was shown in the 1996-1997 potentiometric surface. The depression in Desha County initially was shown in the 1999 potentiometric surface. The depressions in Webster and Winn Parishes were shown as early as 1975. The depressions in Cleveland, western Columbia, and Lafayette Counties initially were shown in the 2003 potentiometric surface. A map of differences in water-level measurements between 2001 and 2005 was constructed using the difference between water-level measurements from 294 wells in Arkansas and 29 wells in Louisiana. The difference in water levels between 2001 and 2005 ranged from -30.1 to 44.6 feet. The largest rise of 44.6 feet in water level measured was in Union County in Arkansas. The largest decline of 30.1 feet in water level measured was in Columbia County in Arkansas. Areas with a general rise in water levels in Arkansas are shown in Arkansas, Columbia, Craighead, Jefferson, Prairie, and the western half of Union Counties. The area around west-central Union County had rises as much as 44.6 feet, with seven wells showing a rise of 20 feet or greater, which is an annual rise of 5 feet or greater. Areas in Arkansas with a general decline in water level are shown in western Bradley, eastern Calhoun, Cleveland, Cross, Desha, Drew, Lafayette, Lee, Lincoln, Lonoke, Poinsett, and the eastern half of Union Counties. In Louisiana, the water-level difference map showed a general rise in water levels in northern Claiborne, northern Webster, and northwestern Union Parishes mainly because of a decrease in industrial withdrawals in southern Arkansas, particularly Union County. Another rise in water level was indicated in western

Scientific Investigations Report↗

Geophysical bed sediment characterization of the Androscoggin River from the former Chlor-Alkali Facility Superfund Site, Berlin, New Hampshire, to the state border with Maine, August 2009

The former Chlor-Alkali Facility in Berlin, New Hampshire, was listed on the U.S. Environmental Protection Agency National Priorities List in 2005 as a Superfund site. The Chlor-Alkali Facility lies on the east bank of the Androscoggin River. Elemental mercury currently discharges from that bank into the Androscoggin River. The nature, extent, and the speciation of mercury and the production of methyl mercury contamination in the adjacent Androscoggin River is the subject of continuing investigations. The U.S. Geological Survey, in cooperation with Region I of the U.S. Environmental Protection Agency, used geophysical methods to determine the distribution, thickness, and physical properties of sediments in the Androscoggin River channel at a small area of an upstream reference reach and downstream from the site to the New Hampshire&ndash;Maine State border. Separate reaches of the Androscoggin River in the study area were surveyed with surface geophysical methods including ground-penetrating radar and step-frequency electromagnetics. Results were processed to assess sediment characteristics including grain size, electrical conductivity, and pore-water specific conductance. Specific conductance measured during surface- and pore-water sampling was used to help interpret the results of the geophysical surveys. The electrical resistivity of sediment samples was measured in the laboratory with intact pore water for comparison with survey results. In some instances, anthropogenic features and land uses, such as roads and power lines affected the detection of riverbed properties using geophysical methods; when this occurred, the data were removed. Through combining results, detailed riverbed sediment characterizations were made. Results from ground-penetrating radar surveys were used to image and measure the depth to the riverbed, depth to buried riverbeds, riverbed thickness and to interpret material-type variations in terms of relative grain size. Fifty two percent of the riverbed in the study area was covered with gravel and finer sediments. The electrically resistive river water and sediment in this study area were conducive to the penetration of the ground-penetrating radar and step-frequency electromagnetic signals and allowed for effective sediment characterization by geophysical methods. The reach between the former Chlor-Alkali Facility and the Riverside Dam, had small areas of fine sediment (estimated 11 percent of riverbed area), found on the upstream left bank and the downstream right bank, with an electromagnetic conductivity (31.4 millisiemens per meter (mS/m) maximum) that was higher than the upstream reference reach. The greatest electromagnetic conductivity (195 mS/m), pore-water specific conductance (324 mS/m) and lab measured sediment conductivity of (76.8 mS/m, measured with a direct-current resistivity test box) in the study were measured approximately 1 mile (mi) downstream of the site from a sandbar on the left bank. Reaches adjacent to and within 2 mi downstream from the site had elevated electromagnetic conductivity despite having lower estimated percentages of riverbed area covered in sediment (11, 25, and 61 percent, respectively) than the reference reach (97). Typically finer grained sediment with similar mineralogy will be more conductive. The Shelburne Reservoir is approximately 8 mi downstream from the site had the second greatest pore-water specific conductance measured, 45.8 mS/m. Many of the locations with the largest step-frequency electromagnetic values have not been sampled for pore water and sediment.

New Hampshire↗

Real-time piscicide tracking using Rhodamine WT dye for support of application, transport, and deactivation strategies in riverine environments

Piscicide applications in riverine environments are complicated by the advection and dispersion of the piscicide by the flowing water. Proper deactivation of the fish toxin is required outside of the treatment reach to ensure that there is minimal collateral damage to fisheries downstream or in connecting and adjacent water bodies. In urban settings and highly managed waterways, further complications arise from the influence of industrial intakes and outfalls, stormwater outfalls, lock and dam operations, and general unsteady flow conditions. These complications affect the local hydrodynamics and ultimately the transport and fate of the piscicide. This report presents two techniques using Rhodamine WT dye for real-time tracking of a piscicide plume—or any passive contaminant—in rivers and waterways in natural and urban settings. Passive contaminants are those that are present in such low concentration that there is no effect (such as buoyancy) on the fluid dynamics of the receiving water body. These methods, when combined with data logging and archiving, allow for visualization and documentation of the application and deactivation process. Real-time tracking and documentation of rotenone applications in rivers and urban waterways was accomplished by encasing the rotenone plume in a plume of Rhodamine WT dye and using vessel-mounted submersible fluorometers together with acoustic Doppler current profilers (ADCP) and global positioning system (GPS) receivers to track the dye and map the water currents responsible for advection and dispersion. In this study, two methods were used to track rotenone plumes: (1) simultaneous injection of dye with rotenone and (2) delineation of the upstream and downstream boundaries of the treatment zone with dye. All data were logged and displayed on a shipboard laptop computer, so that survey personnel provided real-time feedback about the extent of the rotenone plume to rotenone application and deactivation personnel. Further, these strategies facilitate adjustment of rotenone application and deactivation strategies in real time if necessary based on the observed advection and dispersion of the rotenone plume. Two large-scale and complex applications of rotenone in the Chicago Area Waterway System (CAWS) in 2009 and 2010 to combat invasive Asian carp are documented in this report. The application in Chicago Sanitary and Ship Canal (CSSC) in December 2009 involved more than 1,800 gallons of rotenone injected at multiple stations through a 6.2-mile reach of the canal near Lockport, Illinois. The rotenone plume was encased in Rhodamine WT dye so that two survey boats provided real-time feedback to shore personnel regarding the plume extent as it advected downstream. Real-time tracking of the rotenone was essential in this large-scale application because of the multistage injection strategy and the numerous deactivation points required to minimize collateral damage to fisheries in surrounding and receiving water bodies. All timing of application and deactivation operations relied on dye tracking. A second application of rotenone in May 2010 to the Little Calumet River near O’Brien Lock and Dam (Illinois) provided another opportunity for dye-tracking support operations; however, application and deactivation strategies were designed considering zero-flow conditions within the reach of interest. Therefore, dye was injected at the upstream and downstream boundaries of the rotenone application reach and was used to track movement of water in and out of a treatment reach, allowing proper deactivation to occur and avoiding unnecessary damage to fisheries downstream. The data collected during the real-time tracking operations for both applications allowed full documentation of the rotenone treatment for archival purposes and provided information for future applications. The methods presented in this report for real-time tracking and documentation of piscicide applications in riverine environments worked exceptionally well and allowed the multiagency Asian Carp Rapid Response Workgroup to carry out large-scale rotenone applications in urban waterways in an environmentally responsible manner with minimal collateral damage to fisheries outside the treatment reach. Traveltime information extracted from the boat-mounted and fixed-position fluorometers agrees well with empirical predictions from a preliminary dye study (mock rotenone injection) on this system completed in November 2009 on the CSSC and with previously published methods for estimating traveltimes of the peak, leading edge, and trailing edge of the plume. Although the rotenone application strategy called for zero-flow conditions on the Little Calumet River in 2010, downstream advection of treated water did occur, and dye tracing combined with velocity mapping allowed this advection to be documented and exposed the unique hydrodynamics and mixing within this reach. The large volumes of data collected during the operations allow documentation and visualization of the rotenone applications, thus providing feedback to planners and archival of the treatments for future reference. The methods developed in this report are directly transferrable to piscicide applications in water bodies in other locations, including rivers, ponds, or lakes, and can be used for real-time tracking of any passive contaminant that may enter a water body.

Illinois↗

Characterization of water quality and biological communities, Fish Creek, Teton County, Wyoming, 2007-2011

Fish Creek, an approximately 25-kilometer-long tributary to Snake River, is located in Teton County in western Wyoming near the town of Wilson. Fish Creek is an important water body because it is used for irrigation, fishing, and recreation and adds scenic value to the Jackson Hole properties it runs through. Public concern about nuisance growths of aquatic plants in Fish Creek has been increasing since the early 2000s. To address these concerns, the U.S. Geological Survey conducted a study in cooperation with the Teton Conservation District to characterize the hydrology, water quality, and biologic communities of Fish Creek during 2007–11. The hydrology of Fish Creek is strongly affected by groundwater contributions from the area known as the Snake River west bank, which lies east of Fish Creek and west of Snake River. Because of this continuous groundwater discharge to the creek, land-use activities in the west bank area can affect the groundwater quality. Evaluation of nitrate isotopes and dissolved-nitrate concentrations in groundwater during the study indicated that nitrate was entering Fish Creek from groundwater, and that the source of nitrate was commonly a septic/sewage effluent or manure source, or multiple sources, potentially including artificial nitrogen fertilizers, natural soil organic matter, and mixtures of sources. Concentrations of dissolved nitrate and orthophosphate, which are key nutrients for growth of aquatic plants, generally were low in Fish Creek and occasionally were less than reporting levels (not detected). One potential reason for the low nutrient concentrations is that nutrients were being consumed by aquatic plant life that increases during the summer growing season, as a result of the seasonal increase in temperature and larger number of daylight hours. Several aspects of Fish Creek’s hydrology contribute to higher productivity and biovolume of aquatic plants in Fish Creek than typically observed in streams of its size in Wyoming. Especially in the winter, the proportionately large, continuous gain of groundwater into Fish Creek in the perennial section keeps most of the creek free of ice. Because sunlight can still reach the streambed in Fish Creek and the water is still flowing, aquatic plants continue to photosynthesize in the winter, albeit at a lower level of productivity. Additionally, the cobble and large gravel substrate in Fish Creek provides excellent attachment points for aquatic plants, and when combined with Fish Creek’s channel stability allows rapid growth of aquatic plants once conditions allow during the spring. The aquatic plant community of Fish Creek was different than most streams in Wyoming in that it contains many different macrophytes—including macroalgae such as long streamers of Cladophora , aquatic vascular plants, and moss; most other streams in the state contain predominantly algae. From the banks of Fish Creek, the bottom of the stream sometimes appeared to be a solid green carpet. A shift was observed from higher amounts of microalgae in April/May to higher amounts macrophytes in August and October, and differences in the relative abundance of microalgae and macrophytes were statistically significant between seasons. Differences in dissolved-nitrate concentrations and in the nitrogen-to-phosphorus ratio were significantly different between seasons, as concentrations of dissolved nitrate decreased from April/May to August and October. It is likely that dissolved-nitrate concentrations in Fish Creek were lower in August and October because macrophytes were quickly utilizing the nutrient, and a negative correlation between macro-phytes and nitrate was found. Macroinvertebrates also were sampled because of their role as indicators of water quality and their documented responses to perturbation such as degradation of water quality and habitat. Statistically significant seasonal differences were noted in the macroinvertebrate community. Taxa richness and relative abundance of Ephemeroptera, Plecoptera, and Trichoptera, which tend to be intolerant of water-quality degradation, decreased from April/May to August; the same time period saw a corresponding increase in Diptera and noninsects, particularly Oligochaeta (worms) that are more tolerant. Seasonal changes in macroinvertebrate functional feeding groups were significantly different. The relative abundance of gatherer-collector and scraper feeding groups decreased from April/May to August, accompanied by an increase in filterer-collector and shredders feeding groups. Seasonal changes in feeding groups might be due to the seasonal shift in aquatic plant communities, as indicated by comparison with other streams in the area that had fewer aquatic macrophytes than Fish Creek. Statistical tests of macroinvertebrate metrics indicated few differences between years or biological sampling sites on Fish Creek, although the site farthest upstream sometimes was different not only in terms of macroinvertebrates but also in streamflow, water quality, and aquatic plants. Potential effects of contributions of additional nutrients to the Fish Creek ecosystem beyond the conditions sampled during the study period are not known. However, because virtually all of the detectable dissolved nitrate commonly was consumed by aquatic plants in August (leaving dissolved nitrate less than the reporting level in water samples), it is possible that increased nutrient contributions could cause increased growth of aquatic plants. Additional long-term monitoring of the stream, with concurrent data analysis and interpretation would be needed to determine the effects of additional nutrients on the aquatic plant community and on higher levels of the food chain.

Wyoming↗

Comparison of two regression-based approaches for determining nutrient and sediment fluxes and trends in the Chesapeake Bay watershed

Nutrient and sediment fluxes and changes in fluxes over time are key indicators that water resource managers can use to assess the progress being made in improving the structure and function of the Chesapeake Bay ecosystem. The U.S. Geological Survey collects annual nutrient (nitrogen and phosphorus) and sediment flux data and computes trends that describe the extent to which water-quality conditions are changing within the major Chesapeake Bay tributaries. Two regression-based approaches were compared for estimating annual nutrient and sediment fluxes and for characterizing how these annual fluxes are changing over time. The two regression models compared are the traditionally used ESTIMATOR and the newly developed Weighted Regression on Time, Discharge, and Season (WRTDS). The model comparison focused on answering three questions: (1) What are the differences between the functional form and construction of each model? (2) Which model produces estimates of flux with the greatest accuracy and least amount of bias? (3) How different would the historical estimates of annual flux be if WRTDS had been used instead of ESTIMATOR? One additional point of comparison between the two models is how each model determines trends in annual flux once the year-to-year variations in discharge have been determined. All comparisons were made using total nitrogen, nitrate, total phosphorus, orthophosphorus, and suspended-sediment concentration data collected at the nine U.S. Geological Survey River Input Monitoring stations located on the Susquehanna, Potomac, James, Rappahannock, Appomattox, Pamunkey, Mattaponi, Patuxent, and Choptank Rivers in the Chesapeake Bay watershed. Two model characteristics that uniquely distinguish ESTIMATOR and WRTDS are the fundamental model form and the determination of model coefficients. ESTIMATOR and WRTDS both predict water-quality constituent concentration by developing a linear relation between the natural logarithm of observed constituent concentration and three explanatory variables—the natural log of discharge, time, and season. ESTIMATOR uses two additional explanatory variables—the square of the log of discharge and time-squared. Both models determine coefficients for variables for a series of estimation windows. ESTIMATOR establishes variable coefficients for a series of 9-year moving windows; all observed constituent concentration data within the 9-year window are used to establish each coefficient. Conversely, WRTDS establishes variable coefficients for each combination of discharge and time using only observed concentration data that are similar in time, season, and discharge to the day being estimated. As a result of these distinguishing characteristics, ESTIMATOR reproduces concentration-discharge relations that are closely approximated by a quadratic or linear function with respect to both the log of discharge and time. Conversely, the linear model form of WRTDS coupled with extensive model windowing for each combination of discharge and time allows WRTDS to reproduce observed concentration-discharge relations that are more sinuous in form. Another distinction between ESTIMATOR and WRTDS is the reporting of uncertainty associated with the model estimates of flux and trend. ESTIMATOR quantifies the standard error of prediction associated with the determination of flux and trends. The standard error of prediction enables the determination of the 95-percent confidence intervals for flux and trend as well as the ability to test whether the reported trend is significantly different from zero (where zero equals no trend). Conversely, WRTDS is unable to propagate error through the many (over 5,000) models for unique combinations of flow and time to determine a total standard error. As a result, WRTDS flux estimates are not reported with confidence intervals and a level of significance is not determined for flow-normalized fluxes. The differences between ESTIMATOR and WRTDS, with regard to model form and determination of model coefficients, have an influence on the determination of nutrient and sediment fluxes and associated changes in flux over time as a result of management activities. The comparison between the model estimates of flux and trend was made for combinations of five water-quality constituents at nine River Input Monitoring stations. The major findings with regard to nutrient and sediment fluxes are as follows: (1)WRTDS produced estimates of flux for all combinations that were more accurate, based on reduction in root mean squared error, than flux estimates from ESTIMATOR; (2) for 67 percent of the combinations, WRTDS and ESTIMATOR both produced estimates of flux that were minimally biased compared to observed fluxes(flux bias = tendency to over or underpredict flux observations); however, for 33 percent of the combinations, WRTDS produced estimates of flux that were considerably less biased (by at least 10 percent) than flux estimates from ESTIMATOR; (3) the average percent difference in annual fluxes generated by ESTIMATOR and WRTDS was less than 10 percent at 80 percent of the combinations; and (4) the greatest differences related to flux bias and annual fluxes all occurred for combinations where the pattern in observed concentration-discharge relation was sinuous (two points of inflection) rather than linear or quadratic (zero or one point of inflection). The major findings with regard to trends are as follows: (1) both models produce water-quality trends that have factored in the year-to-year variations in flow; (2) trends in water-quality condition are represented by ESTIMATOR as a trend in flow-adjusted concentration and by WRTDS as a flow normalized flux; (3) for 67 percent of the combinations with trend estimates, the WRTDS trends in flow-normalized flux are in the same direction and magnitude to the ESTIMATOR trends in flow-adjusted concentration, and at the remaining 33 percent the differences in trend magnitude and direction are related to fundamental differences between concentration and flux; and (4) the majority (85 percent) of the total nitrogen, nitrate, and orthophosphorus combinations exhibited long-term (1985 to 2010) trends in WRTDS flow-normalized flux that indicate improvement or reduction in associated flux and the majority (83 percent) of the total phosphorus (from 1985 to 2010) and suspended sediment (from 2001 to 2010) combinations exhibited trends in WRTDS flow-normalized flux that indicate degradation or increases in the flux delivered.

Delaware, Maryland, New York, Pennsylvania, Virgin↗