USGS Science⌕ Search

SEARCH · USGS Science

Results for “Cooperative Investigations Report”

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

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

At least 1,603 records · Page 89Linked to original sources

Measurement of unsaturated hydraulic properties and evaluation of property-transfer models for deep sedimentary interbeds, Idaho National Laboratory, Idaho

Operations at the Idaho National Laboratory (INL) have the potential to contaminate the underlying Eastern Snake River Plain (ESRP) aquifer. Methods to quantitatively characterize unsaturated flow and recharge to the ESRP aquifer are needed to inform water-resources management decisions at INL. In particular, hydraulic properties are needed to parameterize distributed hydrologic models of unsaturated flow and transport at INL, but these properties are often difficult and costly to obtain for large areas. The unsaturated zone overlying the ESRP aquifer consists of alternating sequences of thick fractured volcanic rocks that can rapidly transmit water flow and thinner sedimentary interbeds that transmit water much more slowly. Consequently, the sedimentary interbeds are of considerable interest because they primarily restrict the vertical movement of water through the unsaturated zone. Previous efforts by the U.S. Geological Survey (USGS) have included extensive laboratory characterization of the sedimentary interbeds and regression analyses to develop property-transfer models, which relate readily available physical properties of the sedimentary interbeds (bulk density, median particle diameter, and uniformity coefficient) to water retention and unsaturated hydraulic conductivity curves. During 2013–14, the USGS, in cooperation with the U.S. Department of Energy, focused on further characterization of the sedimentary interbeds below the future site of the proposed Remote Handled Low-Level Waste (RHLLW) facility, which is intended for the long-term storage of low-level radioactive waste. Twelve core samples from the sedimentary interbeds from a borehole near the proposed facility were collected for laboratory analysis of hydraulic properties, which also allowed further testing of the property-transfer modeling approach. For each core sample, the steady-state centrifuge method was used to measure relations between matric potential, saturation, and conductivity. These laboratory measurements were compared to water-retention and unsaturated hydraulic conductivity parameters estimated using the established property-transfer models. For each core sample obtained, the agreement between measured and estimated hydraulic parameters was evaluated quantitatively using the Pearson correlation coefficient (r). The highest correlation is for saturated hydraulic conductivity ( K sat ) with an r value of 0.922. The saturated water content (q sat ) also exhibits a strong linear correlation with an r value of 0.892. The curve shape parameter (λ) has a value of 0.731, whereas the curve scaling parameter (y o ) has the lowest r value of 0.528. The r values demonstrate that model predictions correspond well to the laboratory measured properties for most parameters, which supports the value of extending this approach for quantifying unsaturated hydraulic properties at various sites throughout INL.

Idaho↗

Main-stem seepage and base-flow recession time constants in the Niobrara National Scenic River Basin, Nebraska, 2016–18

The Niobrara River of northern Nebraska is a valuable water resource that sustains irrigated agriculture and recreation, as well as a diverse ecosystem. Large-quantity withdrawals from the source aquifer system have the potential to reduce the flow into the river and to adversely affect the free-flowing condition of the Niobrara National Scenic River (NSR). Therefore, to understand the magnitude and characteristics of those flows, the U.S. Geological Survey (USGS), in cooperation with the National Park Service, began a study to quantify seepage gains/losses along the eastern half of the Niobrara NSR and to create a map characterizing the base-flow recession time constant (tau) in the Niobrara NSR study area. In 2016, a seepage study was completed to quantify seepage gains/losses along the eastern half of the Niobrara NSR. The seepage study results indicated that the main-stem streamflow on the Niobrara River increases 375 cubic feet per second (ft 3 /s) in the 39.9-mile study reach (river mile 119.3 to river mile 79.4). Although most of the streamflow increases are attributed to tributary inflows (297 ft 3 /s, 79 percent), 78 ft 3 /s are attributed to seepage gains within the reach. Seepage rates in the study reach ranged from 1.41 cubic feet per second per mile ([ft 3 /s]/mi) to 2.56 (ft 3 /s)/mi, with a mean seepage rate of 2 (ft 3 /s)/mi. Tau values were calculated at 10 sites in the Niobrara NSR study area, and kriging geostatistical techniques were used to develop a contour map to estimate tau values at locations where streamflow was not measured. The minimum tau value was 12.1 days at Willow Creek at Atwood Road near Carns, Nebraska (USGS station 06463670), and the maximum value was 45.5 days at Tyler Falls at Fort Niobrara National Wildlife Refuge near Valentine, Nebr. (USGS station 06461150).

Nebraska↗

Water-budgets and recharge-area simulations for the Spring Creek and Nittany Creek Basins and parts of the Spruce Creek Basin, Centre and Huntingdon Counties, Pennsylvania, Water Years 2000–06

This report describes the results of a study by the U.S. Geological Survey in cooperation with ClearWater Conservancy and the Pennsylvania Department of Environmental Protection to develop a hydrologic model to simulate a water budget and identify areas of greater than average recharge for the Spring Creek Basin in central Pennsylvania. The model was developed to help policy makers, natural resource managers, and the public better understand and manage the water resources in the region. The Groundwater and Surface-water FLOW model (GSFLOW), which is an integration of the Precipitation-Runoff Modeling System (PRMS) and the Modular Groundwater Flow Model (MODFLOW-NWT), was used to simulate surface water and groundwater in the Spring Creek Basin for water years 2000–06. Because the groundwater and surface-water divides for the Spring Creek Basin do not coincide, the study area includes the Nittany Creek Basin and headwaters of the Spruce Creek Basin. The hydrologic model was developed by the use of a stepwise process: (1) develop and calibrate a PRMS model and steady-state MODFLOW-NWT model; (2) re-calibrate the steady-state MODFLOW-NWT model using potential recharge estimates simulated from the PRMS model, and (3) integrate the PRMS and MODFLOW-NWT models into GSFLOW. The individually calibrated PRMS and MODFLOW-NWT models were used as a starting point for the calibration of the fully coupled GSFLOW model. The GSFLOW model calibration was done by comparing observations and corresponding simulated values of streamflow from 11 streamgages and groundwater levels from 16 wells. The cumulative water budget and individual water budgets for water years 2000–06 were simulated by using GSFLOW. The largest source and sink terms are represented by precipitation and evapotranspiration, respectively. For the period simulated, a net surplus in the water budget was computed where inflows exceeded outflows by about 1.7 billion cubic feet (0.47 inches per year over the basin area); storage increased by about the same amount to balance the budget. The rate and distribution of recharge throughout the Spring Creek, Nittany Creek, and Spruce Creek Basins is variable as a result of the high degree of hydrogeologic heterogeneity and karst features. The greatest amount of recharge was simulated in the carbonate-bedrock valley, near the toe slopes of Nittany and Tussey Mountains, in the Scotia Barrens, and along the area coinciding with the Gatesburg Formation. Runoff extremes were observed for water years 2001 (dry year) and 2004 (wet year). Simulated average recharge rates (water reaching the saturated zone as defined in GSFLOW) for 2001 and 2004 were 5.4 in/yr and 22.0 in/yr, respectively. Areas where simulations show large variations in annual recharge between wet and dry years are the same areas where simulated recharge was large. Those areas where rates of groundwater recharge are much higher than average, and are capable of accepting substantially greater quantities of recharge during wet years, might be considered critical for maintaining the flow of springs, stream base flow, or the source of water to supply wells. The slopes of the Bald Eagle, Tussey, and Nittany Mountains are relatively insensitive to variations in recharge, primarily because of reduced infiltration rates and steep slopes.

Pennsylvania↗

Hydrogeologic framework and simulation of predevelopment groundwater flow, eastern Abu Dhabi Emirate, United Arab Emirates

Groundwater in eastern Abu Dhabi in the United Arab Emirates is an important resource that is widely used for irrigation and domestic supplies in rural areas. The U.S. Geological Survey and the Environment Agency—Abu Dhabi cooperated on an investigation to integrate existing hydrogeologic information and to answer questions about regional groundwater resources in Abu Dhabi by developing a numerical groundwater flow model based on MODFLOW–2005 software. The groundwater flow model developed in this investigation provides an improved understanding of groundwater conditions in the eastern region of the Emirate of Abu Dhabi. The flow model simulates steady-state predevelopment conditions from before the rapid growth of modern pumping in the 1980s and was calibrated with 1,342 groundwater-level observations by use of automated and manual calibration techniques. The calibrated model provides good accuracy, with a mean error of 0.50 meters and a standard error of 5.92 meters for simulated groundwater levels. The results of the regional water budget simulation show that gap recharge, which is groundwater inflow through mountain-front gap alluvium, is the greatest source of water to the aquifer. In the base simulation scenario, gap recharge represents 80 percent of total inflow (119,470 of 149,403 cubic meters per day) and the greatest outflow from the aquifer is from evapotranspiration (93 percent of total outflow). Model scenario and sensitivity results reveal a need for data that more thoroughly and more accurately describe aquifer hydraulic conductivity, inflow to the aquifer from the Oman Mountains, and recharge from precipitation on the piedmont. Additional long-term aquifer pumping test observations would improve understanding of aquifer hydraulic conductivity, which would also improve model accuracy. Future studies can modify the model to understand the effect of land-use change and water use on groundwater supplies and simulate more complex groundwater flow conditions in a predictive mode.

Scientific Investigations Report↗

Characterization of suspended sediment flux and streamflow trends in the Fountain Creek watershed, Colorado, 1998 through 2022

The U.S. Geological Survey evaluated long-term suspended sediment flux and streamflow datasets for temporal trends (monotonic and step trends) at 10 streamgage sites within the Fountain Creek watershed in central Colorado using the Mann-Kendall test (monotonic trend) and the Wilcoxon signed-rank test (step trend). Data were collected in cooperation with Colorado Springs Stormwater Enterprise. In this study, 10 sites with long-term suspended sediment records were evaluated during their operational periods, which ranged from 1998 through 2022. To assess how stream behavior might relate to shifts in suspended sediment transport, the Richards-Baker flashiness index, a measure of flashiness, was evaluated for each site. The Richards-Baker flashiness index was calculated for the same months as streamflow and suspended sediment (April through September) using all available streamflow data for a given site. Additionally, cumulative double-mass curves were developed to define temporal variation in the relation between streamflow and suspended sediment loads. This was completed by assessing differences in slopes before and after an observed break in the double-mass curve plots. Five streamgage sites showed statistically significant ( p <0.05) negative trends for suspended sediment flux, and nonsignificant decreases were indicated at the other five sites. The statistically significant negative trends are distributed across the watershed and include smaller tributaries and the main stem of Fountain Creek closer to its confluence with the Arkansas River. Such a broad distribution of negative trends is most likely an indication of improved water quality in the watershed with regards to suspended sediment. Assessing causes for the decreases in suspended sediment loads is beyond the scope of this report; however, spatially distributed bank erosion control projects, stormflow retention projects in the watershed, and changes in climate and storm patterns could be some of the potential drivers of the suspended sediment load trends in this watershed. The flashiness of the streamflow was assessed and can be dismissed on the basis of this study as a potential driver of trends in suspended sediment flux because the magnitude of changes was found to be negligible at all sites.

Colorado↗

Analytical Results for Municipal Biosolids Samples from a Monitoring Program near Deer Trail, Colorado (U.S.A.), 2007

Since late 1993, the Metro Wastewater Reclamation District of Denver (Metro District), a large wastewater treatment plant in Denver, Colorado, has applied Grade I, Class B biosolids to about 52,000 acres of nonirrigated farmland and rangeland near Deer Trail, Colorado (U.S.A.). In cooperation with the Metro District in 1993, the U.S. Geological Survey (USGS) began monitoring ground water at part of this site. In 1999, the USGS began a more comprehensive monitoring study of the entire site to address stakeholder concerns about the potential chemical effects of biosolids applications to water, soil, and vegetation. This more comprehensive monitoring program recently has been extended through 2010. Monitoring components of the more comprehensive study include biosolids collected at the wastewater treatment plant, soil, crops, dust, alluvial and bedrock ground water, and streambed sediment. Streams at the site are dry most of the year, so samples of streambed sediment deposited after rain were used to indicate surface-water effects. This report will present only analytical results for the biosolids samples collected at the Metro District wastewater treatment plant in Denver and analyzed during 2007. We have presented earlier a compilation of analytical results for the biosolids samples collected and analyzed for 1999 through 2006. More information about the other monitoring components is presented elsewhere in the literature. Priority parameters for biosolids identified by the stakeholders and also regulated by Colorado when used as an agricultural soil amendment include the total concentrations of nine trace elements (arsenic, cadmium, copper, lead, mercury, molybdenum, nickel, selenium, and zinc), plutonium isotopes, and gross alpha and beta activity. Nitrogen and chromium also were priority parameters for ground water and sediment components. In general, the objective of each component of the study was to determine whether concentrations of priority parameters (1) were higher than regulatory limits, (2) were increasing with time, or (3) were significantly higher in biosolids-applied areas than in a similar farmed area where biosolids were not applied. Previous analytical results indicate that the elemental composition of the biosolids from the Denver plant was consistent during 1999-2006 and this consistency continues with the samples for 2007; total concentrations of regulated trace elements remained consistently lower than the regulatory limits for the entire monitoring period. Our previously reported data (1999-2006) and data presented in this report were used to compile an inorganic-chemical biosolids signature that can be contrasted with the geochemical signature for this site. The biosolids signature and an understanding of the geology and hydrology of the site can be used to separate biosolids effects from natural geochemical effects. Elements of particular interest for a biosolids signature include bismuth, copper, silver, mercury, phosphorus, and silver. An alternative method of digestion of biosolids was also recently investigated, and the results are presented in this report. A microwave digestion using only nitric acid at controlled elevated temperature and pressure was tested to replace the much more time-consuming and labor-intensive, traditional four-acid, hotplate method for the preparation of solutions to be analyzed by inductively coupled plasma-mass spectrometry (ICP-MS). Elements of concern determined by ICP-MS following digestion include cadmium, copper, lead, molybdenum, nickel, and zinc. The microwave 'digestion' proved to be a strong acid leach, and it was less efficient at digesting the biosolids samples with consistently lower recoveries (compared to the four-acid digestion value) for most elements, but especially for the elements of concern - copper, nickel, and zinc. Other elements traditionally associated with the silicate or oxide minerals demonstrated low recoveries, especially titaniu

Open-File Report↗

Flood-inundation maps for a 15-mile reach of the Kalamazoo River from Marshall to Battle Creek, Michigan

Digital flood-inundation maps for a 15-mile reach of the Kalamazoo River from Marshall to Battle Creek, Michigan, were created by the U.S. Geological Survey (USGS) in cooperation with the U.S. Environmental Protection Agency to help guide remediation efforts following a crude-oil spill on July 25, 2010. The spill happened on Talmadge Creek, a tributary of the Kalamazoo River near Marshall, during a flood. The floodwaters transported the spilled oil down the Kalamazoo River and deposited oil in impoundments and on the surfaces of islands and flood plains. Six flood-inundation maps were constructed corresponding to the flood stage (884.09 feet) coincident with the oil spill on July 25, 2010, as well as for floods with annual exceedance probabilities of 0.2, 1, 2, 4, and 10 percent. Streamflow at the USGS streamgage at Marshall, Michigan (USGS site ID 04103500), was used to calculate the flood probabilities. From August 13 to 18, 2010, 35 channel cross sections, 17 bridges and 1 dam were surveyed. These data were used to construct a water-surface profile for the July 25, 2010, flood by use of a one-dimensional step-backwater model. The calibrated model was used to estimate water-surface profiles for other flood probabilities. The resulting six flood-inundation maps were created with a geographic information system by combining flood profiles with a 1.2-foot vertical and 10-foot horizontal resolution digital elevation model derived from Light Detection and Ranging data.

Michigan↗

Hydrologic balance, water quality, chemical-mass balance, and geochemical modeling of hyperalkaline ponds at Big Marsh, Chicago, Illinois, 2016–17

Hyperalkaline (pH greater than 12) ponds and groundwater exist at Big Marsh near Lake Calumet, Chicago, Illinois, a site used by the steel industry during the mid-1900s to deposit steel- and iron-making waste, in particular, slag. The hyperalkaline ponds may pose a hazard to human health and the environment. The U.S. Geological Survey (USGS), in cooperation with the Environmental Protection Agency (EPA) and in collaboration with the City of Chicago’s Park District, completed a study to evaluate the hydrologic balance, water quality, and chemical-mass balance of hyperalkaline ponds at Big Marsh and geochemical modeling used to evaluate remediation options for water quality at the site based on data collected in 2016–17. Synoptic measurements of surface-water and groundwater elevations were used to determine flow directions and to enable a preliminary estimate of the hydrologic balance for the ponds. Water-quality samples also were collected and analyzed for selected constituents including major anions and cations, nutrients, metals, and trace elements. The results of the water-quality analyses were used to develop a geochemical model to evaluate concentrations, factors affecting pH, and the state of equilibrium between surface waters and atmospheric carbon dioxide. The geochemical model was used to evaluate remediation scenarios using riprap, spillways, or active aeration. The results indicate that active aeration will decrease the pH to near 7.5 in about 8 hours, the fastest rate of the scenarios. Passive aeration, such as riprap or spillways, also can be effective at decreasing the pH in about 45 hours, but spatial obstacles limit their implementation. Seasonal variations in temperature also affect the rate of equilibration, where colder temperatures may have a lower pH than warmer temperatures and may affect the timing and frequency of remediation.

Illinois↗

Approaches for assessing long-term annual yields of highway and urban runoff in selected areas of California with the Stochastic Empirical Loading and Dilution Model (SELDM)

The California Department of Transportation, commonly known as CalTrans, and other municipal separate storm sewer system permittees in California as well as other State departments of transportation nationwide need information about potential loads and yields (loads per unit area) of constituents of concern in stormwater runoff and discharges from stormwater best management practices (BMPs). Although its National Pollution Discharge Elimination System stormwater permit is focused on areas subject to total maximum daily load (TMDL) regulations, CalTrans builds and maintains BMPs to minimize the adverse effects of roadway runoff on receiving waters throughout the State. This report describes approaches used by the U.S. Geological Survey in cooperation with CalTrans for using the Stochastic Empirical Loading and Dilution Model (SELDM) to assess long-term annual yields of highway and urban runoff in selected areas of California. In this study, a series of regional and local yields were simulated to provide statewide planning-level estimates and more refined TMDL-specific yield values. SELDM was used to analyze 368 State roadway and urban runoff yields for 53 runoff quality constituents. The analyses included 222 random-seed analyses, 60 regional State roadway-runoff analyses, 24 regional urban roadway-runoff analyses, and 62 focused TMDL-area analyses. This report describes approaches and statistics used to analyze available hydrologic and runoff quality data in all analyses. Results for all analyses are provided in the model archive, but only a selected subset of results are presented as examples in this report. State roadway runoff, urban runoff, and BMP discharge yields for total suspended solids, total nitrogen, total phosphorus, and total zinc were selected as examples because they are widespread constituents of concern with substantial amounts of State roadway and urban runoff monitoring data. In this report, a hypothetical basin was specified by using available geographic information to demonstrate use of the State roadway and urban runoff yields to estimate long-term annual stormwater loads from developed areas. Application of these yields to the hypothetical basin indicates that although State-roadway yields may be higher than urban-runoff yields for some constituents, State-roadway loads may be a small proportion of total stormwater loads because State roadways themselves are a small fraction of the total impervious area in such basins. Although application of results from this study may have considerable uncertainty for any particular stormwater outfall, the study does provide robust estimates to support basin-scale runoff-load analyses in California. These analyses also provide estimates for the 12 U.S. Environmental Protection Agency level III ecoregions that are completely or partially within the boundaries of the State of California.

California↗

Numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units in the Rapid City area, South Dakota

The city of Rapid City and other water users in the Rapid City area obtain water supplies from the Minnelusa and Madison aquifers, which are contained in the Minnelusa and Madison hydrogeologic units. A numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units in the Rapid City area was developed to synthesize estimates of water-budget components and hydraulic properties, and to provide a tool to analyze the effect of additional stress on water-level altitudes within the aquifers and on discharge to springs. This report, prepared in cooperation with the city of Rapid City, documents a numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units for the 1,000-square-mile study area that includes Rapid City and the surrounding area. Water-table conditions generally exist in outcrop areas of the Minnelusa and Madison hydrogeologic units, which form generally concentric rings that surround the Precambrian core of the uplifted Black Hills. Confined conditions exist east of the water-table areas in the study area. The Minnelusa hydrogeologic unit is 375 to 800 feet (ft) thick in the study area with the more permeable upper part containing predominantly sandstone and the less permeable lower part containing more shale and limestone than the upper part. Shale units in the lower part generally impede flow between the Minnelusa hydrogeologic unit and the underlying Madison hydrogeologic unit; however, fracturing and weathering may result in hydraulic connections in some areas. The Madison hydrogeologic unit is composed of limestone and dolomite that is about 250 to 610 ft thick in the study area, and the upper part contains substantial secondary permeability from solution openings and fractures. Recharge to the Minnelusa and Madison hydrogeologic units is from streamflow loss where streams cross the outcrop and from infiltration of precipitation on the outcrops (areal recharge). MODFLOW-2000, a finite-difference groundwater-flow model, was used to simulate flow in the Minnelusa and Madison hydrogeologic units with five layers. Layer 1 represented the fractured sandstone layers in the upper 250 ft of the Minnelusa hydrogeologic unit, and layer 2 represented the lower part of the Minnelusa hydrogeologic unit. Layer 3 represented the upper 150 ft of the Madison hydrogeologic unit, and layer 4 represented the less permeable lower part. Layer 5 represented an approximation of the underlying Deadwood aquifer to simulate upward flow to the Madison hydrogeologic unit. The finite-difference grid, oriented 23 degrees counterclockwise, included 221 rows and 169 columns with a square cell size of 492.1 ft in the detailed study area that surrounded Rapid City. The northern and southern boundaries for layers 1-4 were represented as no-flow boundaries, and the boundary on the east was represented with head-dependent flow cells. Streamflow recharge was represented with specified-flow cells, and areal recharge to layers 1-4 was represented with a specified-flux boundary. Calibration of the model was accomplished by two simulations: (1) steady-state simulation of average conditions for water years 1988-97 and (2) transient simulations of water years 1988-97 divided into twenty 6-month stress periods. Flow-system components represented in the model include recharge, discharge, and hydraulic properties. The steady-state streamflow recharge rate was 42.2 cubic feet per second (ft3/s), and transient streamflow recharge rates ranged from 14.1 to 102.2 ft3/s. The steady-state areal recharge rate was 20.9 ft3/s, and transient areal recharge rates ranged from 1.1 to 98.4 ft3/s. The upward flow rate from the Deadwood aquifer to the Madison hydrogeologic unit was 6.3 ft3/s. Discharge included springflow, water use, flow to overlying units, and regional outflow. The estimated steady-state springflow of 32.8 ft3/s from seven springs was similar to the simulated springflow of 31.6 ft3/s, which included 20.5 ft3

South Dakota↗

Mercury Loads in the South River and Simulation of Mercury Total Maximum Daily Loads (TMDLs) for the South River, South Fork Shenandoah River, and Shenandoah River: Shenandoah Valley, Virginia

Due to elevated levels of methylmercury in fish, three streams in the Shenandoah Valley of Virginia have been placed on the State's 303d list of contaminated waters. These streams, the South River, the South Fork Shenandoah River, and parts of the Shenandoah River, are downstream from the city of Waynesboro, where mercury waste was discharged from 1929-1950 at an industrial site. To evaluate mercury contamination in fish, this total maximum daily load (TMDL) study was performed in a cooperative effort between the U.S. Geological Survey, the Virginia Department of Environmental Quality, and the U.S. Environmental Protection Agency. The investigation focused on the South River watershed, a headwater of the South Fork Shenandoah River, and extrapolated findings to the other affected downstream rivers. A numerical model of the watershed, based on Hydrological Simulation Program-FORTRAN (HSPF) software, was developed to simulate flows of water, sediment, and total mercury. Results from the investigation and numerical model indicate that contaminated flood-plain soils along the riverbank are the largest source of mercury to the river. Mercury associated with sediment accounts for 96 percent of the annual downstream mercury load (181 of 189 kilograms per year) at the mouth of the South River. Atmospherically deposited mercury contributes a smaller load (less than 1 percent) as do point sources, including current discharge from the historic industrial source area. In order to determine how reductions of mercury loading to the stream could reduce methylmercury concentrations in fish tissue below the U.S. Environmental Protection Agency criterion of 0.3 milligrams per kilogram, multiple scenarios were simulated. Bioaccumulation of mercury was expressed with a site-specific exponential relation between aqueous total mercury and methylmercury in smallmouth bass, the indicator fish species. Simulations indicate that if mercury loading were to decrease by 98.9 percent from 189 to 2 kilograms per year, fish tissue methylmercury concentrations would drop below 0.3 milligrams per kilogram. Based on the simulations, the estimated maximum load of total mercury that can enter the South River without causing fish tissue methylmercury concentrations to rise above 0.3 milligrams per kilogram is 2.03 kilograms per year for the South River, and 4.12 and 6.06 kilograms per year for the South Fork Shenandoah River and Shenandoah River, respectively.

Scientific Investigations Report↗

Analysis of factors affecting plume remediation in a sole-source aquifer system, southeastern Nassau County, New York

Several plumes of dissolved, chlorinated solvents, including trichloroethylene, have been identified in a sole-source aquifer near the former Northrop Grumman Bethpage Facility and Naval Weapons Industrial Reserve Plant sites in southeastern Nassau County, New York. Past investigations have documented that the groundwater contamination originated from this industrial area and now extends to the south, in the direction of groundwater flow. The intermixed plumes are commonly referred to as the “Navy Grumman groundwater plume.” Detailed groundwater-flow modeling was needed for the New York State Department of Environmental Conservation (NYSDEC) to evaluate design options necessary for the construction, operation, optimization, maintenance, and monitoring of a groundwater extraction and treatment cleanup plan selected in a December 2019 Amended Record of Decision by the NYSDEC to comprehensively address these plumes. Consequently, the NYSDEC began a cooperative study with the U.S. Geological Survey in 2020 to better understand the local hydrogeologic framework using two independent approaches to characterize aquifer heterogeneity and update an existing regional groundwater-flow model to provide transient boundary conditions for new inset groundwater-flow models of the plume area. We developed these detailed inset models for the two independent aquifer characterizations using history-matching techniques coupled with a novel approach to risk-based management optimization of the remedial design. We also used the updated regional model to assess this optimized groundwater extraction and treatment design for potential saltwater intrusion. The ensembles of parameters resulting from history matching provided a platform with which to evaluate capture by water-supply and remedial wells using particle-tracking techniques. Using the ensemble to select a risk stance, we performed multiobjective optimization to identify various configurations of remedial pumping that are consistent with external constraints and that favor potentially competing objectives. Multiple solutions provide tradeoffs that NYSDEC can consider. In general, pumping redistribution may help to prevent further contamination migration downgradient. These and other study results are intended to support decisions for the remedial design focused on the local area encompassing the full extent of the Navy Grumman groundwater plume.

New York↗

Sediment transport, terrace incision, and sediment chemistry in subbasins of differing land use, Stinking Water Creek Basin near Rangely, Colorado

In many areas of western Colorado, Late Cretaceous Mancos Shale is a source of sediment, salinity, and selenium to surface water. Natural weathering processes and human activity can expose and move sediments from the Mancos Shale, leading to increased salinity and selenium concentrations in surface water. The U.S. Geological Survey, in cooperation with the Bureau of Land Management, studied four subbasins of Stinking Water Creek near Rangely, Colorado, to provide information regarding landscapes on Mancos Shale that have undergone land-use change. Rangeland grazing is present in all four subbasins. Subbasins A1 and A2 have land use related to energy development and subbasins B1 and B2 do not. Aerial photographs, sediment chemistry, optically stimulated luminescence (OSL) ages, cross-section surveys, and Watershed Erosion Prediction Project (WEPP) modeling were used to characterize the history of land-use change by estimating stream channel incision rates and sediment, salinity, and selenium loads. All subbasins had upper and lower terraces and a modern streambed channel, and these surfaces were dated using OSL. In subbasins A1, A2, and B1, the streambed age was older than the terraces, indicating the channels are incising and widening into the older alluvium. In subbasin B2, however, the streambed age was much younger than the terraces, indicating aggradation is the main process. Differences in channel morphology were not always related to land use but to the subbasin position (upstream compared to downstream). Channels in subbasins A1 and A2 were generally wider and shallower than channels in subbasins B1 and B2. Landscape disturbances related to energy development in subbasins A1 and A2 could have caused channel evolution to progress from incision to widening processes.

Colorado↗

Effectiveness of rapid 'ōhi'a death management strategies at a focal disease outbreak on Hawai'i Island

The ongoing spread of rapid ‘ōhi‘a death (ROD) in the Hawaiian Islands threatens the long-term sustainability of ‘ōhi‘a lehua (Metrosideros polymorpha) forests throughout the state. First identified in the Puna district of Hawai‘i Island in 2014, the disease caused by the novel fungi Ceratocystis lukuohia and Ceratocystis huliohia has now spread island-wide and was recently detected on Kaua‘i, O‘ahu, and Maui. The leading hypothesis for the spread of ROD is through airborne ambrosia beetle frass particles that contain viable Ceratocystis propagules, thus management efforts focus on containing this frass. At the time of this study (2017–2018), the Waipunalei site was the northernmost outbreak of ROD on Hawai‘i Island. The focal nature of the outbreak and accessibility of the location provided the opportunity to monitor the effectiveness of two types of proposed management methods to reduce the airborne spread of potentially infective ambrosia beetle frass: tree felling and insecticide treatments. We placed 23 passive environmental samplers (PES), which monitored for airborne frass and wood particles containing C. lukuohia and C. huliohia in a grid that spanned the outbreak area over 22 weeks. Cross-vane panel traps with 1:1 methanol:ethanol lures were attached to nine of the PES to document wood-boring ambrosia and cerambycid beetle populations during the latter three months of the study. Monitoring with PES began three weeks before management and continued for one month after the last infected trees were felled. Glass microscope slides from the 23 PES were examined for airborne ambrosia beetle frass and wood particles by microscopy. DNA was extracted from the slides and tested by qPCR (quantitative polymerase chain reaction) for C. lukuohia and C. huliohia. We also investigated the correlation of beetle gallery counts with tree height and tested the efficacy of Bifen I/T insecticide (active ingredient: bifenthrin 7.9%) for preventing beetle attacks on the cut surface of ‘ōhi‘a bolts (tree stem sections). Beetle trapping data revealed that the area supports a diverse community of wood-boring beetles, some of which likely attack ‘ōhi‘a and may facilitate the spread of ROD. The number of beetle galleries on felled ‘ōhi‘a trees decreased linearly as tree height increased. We also observed significantly fewer beetle attacks on Bifen I/T treated ‘ōhi‘a bolts than non-treated bolts, but gallery formation nearly ceased in both treated and control bolts by week three. Ceratocystis lukuohia DNA was detected twenty-six times and C. huliohia was detected five times in the PES throughout this study. DNA detections were correlated to frass and wood counts, and the number of felled trees were correlated to wood particle counts but not frass counts. Both the timing and distribution of detections across the sampling grid indicate that tree felling may have reduced airborne detections of Ceratocystis DNA soon after tree felling was completed. A subsequent increase in detections after tree felling ceased may indicate that incomplete removal of infected trees and the appearance of new infections in previously asymptomatic trees could have allowed airborne detections of potentially infectious fungal propagules to once again increase.

Hawaii↗

A guide to safe field operations

Most functions of the U.S. Geological Survey (USGS), Water Resources Division (WRD) require employees to participate in numerous field activities ranging from routine meetings with cooperators, other federal and public officials, and private citizens to potentially hazardous assignments, such as making flood measurements and scuba diving to service underwater instruments. It is paramount that each employee be aware of safety procedures and operational policies of the WRD to ensure that (1) their activities avoid or minimize personal injury to the employee, coworkers, or anyone in the vicinity of the field activity, and (2) their conduct does not infringe on the personal or property rights of any individual or organization. The purpose of the guide is to familiarize employees with the operational and safety procedures expected to be followed by each employee as a representative of the WRD. It is also intended as a training tool for all new employees and a document to be reviewed by each employee before undertaking a field assignment. It includes general procedures that are standard and applicable to all field operations, such as communication, vehicle operation, and adequate preparation for anticipated weather conditions. It also includes a discussion of specific procedures and safety considerations for most of the routine field assignments undertaken by hydrologists and hydrologic technicians of the WRD. The guide is not intended to be a technical handbook outlining step-by-step procedures for performing specific tasks or a comprehensive discussion of every possible activity that may be undertaken by a USGS employee. Employees are referred to the Techniques for Water-Resources Investigations (TWRI) series for specific technical procedures and to the U.S. Geological Survey Safety and Environmental Health Handbook 445-1-H (USGS, August 1989), USGS Occupational Hazards and Safety Procedures Handbook 445-2-H (December 1993), the WRD notebook on Safety Policy and Guidance Memoranda, and other references for procedures and safety issues related to nonroutine activities, such as operations on large vessels and aircraft.

Open-File Report↗

Trends in lake chemistry in response to atmospheric deposition and climate in selected Class I wilderness areas in Colorado, Idaho, Utah, and Wyoming, 1993-2009

In 2010, the U.S. Geological Survey, in cooperation with the U.S. Department of Agriculture Forest Service, Air Resource Management, began a study to evaluate long-term trends in lake-water chemistry for 64 high-elevation lakes in selected Class I wilderness areas in Colorado, Idaho, Utah, and Wyoming during 1993 to 2009. The purpose of this report is to describe trends in the chemical composition of these high-elevation lakes. Trends in emissions, atmospheric deposition, and climate variables (air temperature and precipitation amount) are evaluated over a similar period of record to determine likely drivers of changing lake chemistry. Sulfate concentrations in precipitation decreased over the past two decades at high-elevation monitoring stations in the Rocky Mountain region. The trend in deposition chemistry is consistent with regional declines in sulfur dioxide emissions resulting from installation of emission controls at large stationary sources. Trends in nitrogen deposition were not as widespread as those for sulfate. About one-half of monitoring stations showed increases in ammonium concentrations, but few showed significant changes in nitrate concentrations. Trends in nitrogen deposition appear to be inconsistent with available emission inventories, which indicate modest declines in nitrogen emissions in the Rocky Mountain region since the mid-1990s. This discrepancy may reflect uncertainties in emission inventories or changes in atmospheric transformations of nitrogen species that may be affecting deposition processes. Analysis of long-term climate records indicates that average annual mean air temperature minimums have increased from 0.57 to 0.75 &deg;C per decade in mountain areas of the region with warming trends being more pronounced in Colorado. Trends in annual precipitation were not evident over the period 1990 to 2006, although wetter than average years during 1995 to 1997 and drier years during 2001 to 2004 caused a notable decline in precipitation in the middle part of the record.

Colorado;Idaho;Utah;Wyoming↗

Regression Equations for Monthly and Annual Mean and Selected Percentile Streamflows for Ungaged Rivers in Maine

In an effort to delineate hydrologic conditions in Maine, the U.S. Geological Survey, in cooperation with the Maine Department of Transportation, used streamflow data to develop dependent variables for 130 regression equations for estimating monthly and annual mean and 1, 5, 10, 25, 50, 75, 90, 95, and 99 percentile streamflows for ungaged, unregulated rivers in Maine. Daily streamflow data from 24 rural unregulated basins with drainage areas between 14.9 and 1,419 square miles in Maine and northern New Hampshire were used in the derivation of the equations. Streamflow data collected from October 1, 1982, through September 30, 2012, were used to derive the dependent variables for this study to represent current [2015] hydrologic conditions in Maine and northern New Hampshire. Weighted least squares regression techniques were used to derive the final coefficients and measures of uncertainty for the regression equations. Eight basin characteristics serve as the explanatory variables: drainage area, distance from the coast, mean and maximum basin elevation, mean basin slope, mean basin percentage of hydrologic soil group A, fraction of sand and gravel aquifers, and percentage of open water. The largest average errors of prediction are associated with regression equations for the lowest streamflows derived for months during which the lowest streamflows of the year occur (such as the 5 and 1 monthly percentiles for August and September). The regression equations have been derived on the basis of streamflow and basin characteristics data for unregulated, rural drainage basins without substantial streamflow or drainage modifications (for example, diversions and (or) regulation by dams or reservoirs, tile drainage, irrigation, channelization, and impervious paved surfaces), therefore using the equations for regulated or urbanized basins with substantial streamflow or drainage modifications will yield results of unknown error. Input basin characteristics derived using techniques or datasets other than those documented in this report or using values outside the ranges used to develop these regression equations also will yield results of unknown error.

Maine↗

Analysis of groundwater and surface water in areas of isoxaflutole application, Tuscola and Kalamazoo Counties, Michigan

The herbicide 5-cyclopropyl-4-(2-methylsulfonyl-4-trifluoromethylbenzoyl) isoxazole, also known as isoxaflutole (IXF), was conditionally approved for use on corn in Michigan in 2015. The fate of IXF and its degradates in different environmental settings and the processes by which these compounds move to groundwater or to surface-water bodies have been previously studied, but little information about possible persistence and buildup of this herbicide and its degradates in Michigan’s groundwater and surface water is available. Therefore, from 2015 to 2020, the U.S. Geological Survey, in cooperation with the Michigan Department of Agriculture and Rural Development, studied IXF and two of its degradates in two locations where IXF was applied. IXF and its degradates were rarely detected in shallow groundwater downgradient from IXF applications. In contrast, one or more of the three target IXF compounds were detected in 40 percent of surface-water samples. The degradates 1-(2-methylsulfonyl-4-trifluoromethylphenyl)-2-cyano-3-cyclopropyl propan-1-dione), also known as diketonitrile isoxaflutole (DKN), and 2-methylsulfonyl-4-(trifluoromethyl) benzoic acid (BAA), the benzoic acid analogue of IXF, were detected more frequently than the parent compound. At surface-water sites, DKN and BAA reached maximum concentrations within about the first 5 weeks after IXF application after rainfall-runoff events. Concentrations subsequent to post-application maxima decreased through time approximately following first-order, exponential loss kinetics. Carryover of DKN and BAA, from an application year to the following spring, occurred at several surface-water sites, and springtime concentrations were typically 1–5 percent of maximum, post-application concentrations. Virtually no detections were recorded later in the growing season during non-application years. Results indicate rapid loss of IXF and very few detections of the parent compound from the study areas. The degradates DKN and BAA were more frequently detected in surface runoff up to about 1 year after IXF application, but no evidence of longer-term accumulation was found.

Michigan↗