Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Nature, Scientific Data”

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 829 records · Page 46Linked to original sources

Effects of Surface-Water Diversions on Habitat Availability for Native Macrofauna, Northeast Maui, Hawaii

Effects of surface-water diversions on habitat availability for native stream fauna (fish, shrimp, and snails) are described for 21 streams in northeast Maui, Hawaii. Five streams (Waikamoi, Honomanu, Wailuanui, Kopiliula, and Hanawi Streams) were chosen as representative streams for intensive study. On each of the five streams, three representative reaches were selected: (1) immediately upstream of major surface-water diversions, (2) midway to the coast, and (3) near the coast. This study focused on five amphidromous native aquatic species (alamoo, nopili, nakea, opae, and hihiwai) that are abundant in the study area. The Physical Habitat Simulation (PHABSIM) System, which incorporates hydrology, stream morphology and microhabitat preferences to explore relations between streamflow and habitat availability, was used to simulate habitat/discharge relations for various species and life stages, and to provide quantitative habitat comparisons at different streamflows of interest. Hydrologic data, collected over a range of low-flow discharges, were used to calibrate hydraulic models of selected transects across the streams. The models were then used to predict water depth and velocity (expressed as a Froude number) over a range of discharges up to estimates of natural median streamflow. The biological importance of the stream hydraulic attributes was then assessed with the statistically derived suitability criteria for each native species and life stage that were developed as part of this study to produce a relation between discharge and habitat availability. The final output was expressed as a weighted habitat area of streambed for a representative stream reach. PHABSIM model results are presented to show the area of estimated usable bed habitat over a range of streamflows relative to natural conditions. In general, the models show a continuous decrease in habitat for all modeled species as streamflow is decreased from natural conditions. The PHABSIM modeling results from the intensively studied streams were normalized to develop relations between the relative amount of diversion from a stream and the resulting relative change in habitat in the stream. These relations can be used to estimate changes in habitat for diverted streams in the study area that were not intensively studied. The relations indicate that the addition of even a small amount of water to a dry stream has a significant effect on the amount of habitat available. Equations relating stream base-flow changes to habitat changes can be used to provide an estimate of the relative habitat change in the study area streams for which estimates of diverted and natural median base flow have been determined but for which detailed habitat models were not developed. Stream water temperatures, which could have an effect on stream ecology and taro cultivation, were measured in five streams in the study area. In general, the stream temperatures measured at any of the monitoring sites were not elevated enough, based on currently available information, to adversely effect the growth or mortality of native aquatic macrofauna or to cause wetland taro to be susceptible to fungi and associated rotting diseases.

Scientific Investigations Report↗

The potential effects of sodium bicarbonate, a major constituent from coalbed natural gas production, on aquatic life

The production water from coalbed natural gas (CBNG) extraction contains many constituents. The U.S. Environmental Protection Agency has established aquatic life criteria for some of these constituents, and it is therefore possible to evaluate their risk to aquatic life. However, of the major ions associated with produced waters, chloride is the only one with an established aquatic life criterion (U.S. Environmental Protection Agency, 1988). The focus of this research was NaHCO 3 , a compound that is a major constituent of coalbed natural gas produced waters in the Tongue and Powder River Basins. This project included laboratory experiments, field in situ experiments, a field mixing zone study, and a fishery presence/absence assessment. Though this investigation focuses on the Tongue and Powder River Basins, the information is applicable to other watersheds where sodium bicarbonate is a principle component of product water either from CBNG or from traditional or unconventional oil and gas development. These data can also be used to separate effects of saline discharges from those potentially posed by other constituents. Finally, this research effort and the additional collaboration with USGS Water Resources and Mapping, Bureau of Land Management, US Environmental Protection Agency, State of Montana, State of Wyoming, Montana State University, University of Wyoming, and others as part of a Powder River Aquatic Task Group, can be used as a model for successful approaches to studying landscapes with energy development. The laboratory acute toxicity experiments were completed with a suite of organisms, including 7 species of fish, 5 species of invertebrates, and 1 amphibian species. Experiments performed on these multiple species resulted in LC50s that ranged from 1,120 to greater than (>) 8,000 milligrams sodium bicarbonate per liter (mg NaHCO 3 /L) (also defined as 769 to >8,000 milligrams bicarbonate per liter (mg HCO 3 -/L) or total alkalinity expressed as 608 to >4,181 milligrams calcium carbonate per liter (mg CaCO 3 /L)) that varied across species and lifestage within a species. The age at which fish were exposed to NaHCO 3 significantly affected the severity of toxic responses for some organisms. The chronic toxicity of NaHCO 3 was defined in experiments that lasted from 7—60 days post-hatch. For these experiments, sublethal effects such as growth and reproduction, in addition to significant reductions in survival were included in the final determination of effects. Chronic toxicity was observed at concentrations that ranged from 450 to 800mg NaHCO 3 /L (also defined as 430 to 657 mg HCO 3 -/L or total alkalinity expressed as 354 to 539 mg CaCO 3 /L) and the specific concentration depended on the sensitivity of the four species of invertebrates and fish exposed. Sublethal investigations during chronic studies revealed percent decrease in the activity of sodium-potassium adenosine triphosphatase (Na/K ATPase, an enzyme involved in ionoregulation) and the age of the fish at the onset of the decrease may affect the ability of fathead minnow to survive exposures to NaHCO 3 . A database of toxicity evaluations of NaHCO 3 on aquatic life has been constructed. Using these data, sample acute and chronic criteria of 459 and 381 mg NaHCO 3 /L, respectively, can be calculated for the protection of aquatic life. The final derivation and implementation of such criteria is, of course, left to the discretion of the concerned management agencies. A combination of in situ experiments, static-renewal experiments performed simultaneously with in situ experiments, and static renewal experiments performed with site water in the laboratory, demonstrated that untreated coalbed natural gas (CBNG) product water from the Tongue and Powder River Basins reduces survival of fathead minnow and pallid sturgeon. More precisely, the survival of early-lifestage fathead minnow, especially those less than 6-days post hatch (dph), likely is reduced significantly in the field when concentrations of NaHCO 3 rise above 1,500 mg/L. However, age was not a factor for pallid sturgeon and they were sensitive to product water regardless of age. Treatment with the Higgins Loop™ technology and dilution of untreated water increased survival in the laboratory. Both of these situations reduced ammonia in addition to the concentrations of NaHCO 3 . These experiments addressed the acute toxicity of effluent waters being added to the main stem rivers, but did not address issues related to the volumes of water that may be added to the watershed. Mixing zones of the three outfalls studied ranged from approximately 800—1,200 m below the confluence and the areas within these mixing zones with acutely lethal concentrations of NaHCO 3 (as defined by the presence of concentrated dye) are limited. The areas with concentrations of NaHCO 3 more than the concentrations likely to cause significant mortality, and more than the calculated sample water-quality criteria in the Tongue and Powder River Basins appear to be limited to tributaries and parts of mixing zones with considerable additions of untreated effluent.

Montana;Wyoming↗

Dissolved-solids loads discharged from irrigated areas near Manila, Utah, May 2007-October 2012, and relation of loads to selected variables

The Manila/Washam Salinity Project (MWSP) is a cooperative effort by the Natural Resources Conservation Service (NRCS) and local farmers and ranchers to reduce the transport of dissolved solids to Flaming Gorge Reservoir from irrigated agricultural lands near Manila, Utah. To estimate dissolved-solids loads from the MWSP area, discharge and water quality from Birch Spring Draw and other selected outflows and inflows were monitored from May 2007 to October 2012. An average annual May–April streamflow of 5,960 acre-feet discharged from Birch Spring Draw at site BSD-2 to Flaming Gorge Reservoir during 2007–12, containing an average dissolved-solids load of 14,660 tons. An average May–April net dissolved-solids load of 24,300 tons per year discharged from the MWSP area, estimated from the relation between streamflow and dissolved-solids concentration at site BSD-2 and other measured inflows and outflows. The amount of precipitation and water available for irrigation are important factors affecting the dissolved-solids load in outflow from the MWSP area. Net dissolved-solids load discharged from the MWSP area increased with increasing canal streamflow and precipitation measured at Manila during the irrigation season, from May to October, each year. The net tons of dissolved solids discharged from the MWSP area per acre-foot of canal water increased with increasing irrigation season precipitation during May 2007–October 2012. Irrigation improvements began to be implemented in 2007 to reduce dissolved-solids loads discharged from the MWSP area. The theoretical annual net dissolved-solids load where the cumulative NRCS calculated dissolved-solids load reduction is added to the net MWSP dissolved-solids load is what would be expected if there was no irrigation improvement in the area associated with the MWSP. The theoretical data points lie very near the baseline representing the pre-MWSP dissolved-solids load to canal streamflow relation. The proximity of the theoretical data points to the baseline shows that the NRCS calculations of reduction in dissolved-solids load are generally supported by the data collected during this study.

Utah, Wyoming↗

Assessing potential groundwater-level declines from future withdrawals in the Hualapai Valley, northwestern Arizona

A numerical groundwater flow model of the Hualapai Valley Basin in northwestern Arizona was developed to assist water-resource managers in understanding the potential effects of projected groundwater withdrawals on groundwater levels in the basin. The Hualapai Valley Hydrologic Model (HVHM) simulates the hydrologic system for the years 1935 through 2219, including future withdrawal scenarios that simulate large-scale agricultural expansion with and without enhanced groundwater recharge from potential new infiltration basin projects. HVHM is a highly parameterized model (75,586 adjustable parameters) capable of simulating grid-scale variability in aquifer properties (for example, conductivity, specific yield, and specific storage) and system stresses (for instance, natural recharge and groundwater withdrawals). Parameter estimation and uncertainty quantification were performed using an iterative ensemble smoother software (PESTPP-IES) to produce an ensemble of models fit to historical data. Results via the future withdrawal scenario from this ensemble indicate that mean groundwater level will decline at wells in the Kingman subbasin 87 to 128 feet by the year 2050 and 204 to 241 feet by the year 2080. Mean groundwater level is expected to decline at wells in the Hualapai subbasin between 44 and 210 feet by 2050 and between 107 and 350 feet by 2080. The enhanced recharge scenario results show potential for these declines to be partially mitigated in the Kingman subbasin by between 8 and 23 feet in 2050 and between 23 and 43 feet in 2080. The enhanced recharge scenario has no simulated effect on groundwater levels in the Hualapai subbasin. All planned enhanced infiltration projects are located in the Kingman subbasin, which is simulated to become hydraulically disconnected from the Hualapai subbasin owing to groundwater-level declines before 2050. Mean depth to water in the Kingman subbasin as simulated in the future withdrawal scenario will exceed 1,200 feet between the years 2155 and 2214 (median year 2171). In the future withdrawal plus enhanced recharge scenario, mean depth to water in the Kingman subbasin exceeds 1,200 feet between the years 2163 and 2207 (median year 2180), except for one model realization in which the subbasin does not reach an mean depth to water of 1,200 feet by the end of forecast simulation (year 2220). Simulated dewatering of the basin margins reduces scenario pumping rates by as much as 7 percent in 2029 and 12 percent in 2079 below specified rates. Forecasts of groundwater-level declines are based on the reduced simulated pumping rates.

Arizona↗

Completion summary for boreholes USGS 148, 148A, and 149 at the Materials and Fuels Complex, Idaho National Laboratory, Idaho

In 2019, the U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy, drilled and constructed boreholes USGS 148A and USGS 149 for stratigraphic framework analyses and long-term groundwater monitoring of the eastern Snake River Plain aquifer at the Idaho National Laboratory (INL) in southeastern Idaho. Initially, boreholes USGS 148A and USGS 149 were continuously cored to allow the USGS and INL subcontractor to collect select geophysical and seismic data and evaluate properties of recovered core material. The USGS geophysical data and descriptions of core material are described in this report; however, data collected by the INL contractor, including seismic data, are not included as part of the report. The unsaturated zone at both borehole locations is relatively thick, depth to water was measured at approximately 663.6 feet (ft) below land surface (BLS) in USGS 148A, and at approximately 654.1 ft BLS at USGS 149. On completion of coring and data collection, both boreholes (USGS 148A and USGS 149) were repurposed as monitoring wells. Well USGS 148A was constructed to a depth of 759 ft BLS and instrumented with a dedicated submersible pump and measurement line; well USGS 149 was constructed to a depth of 974 ft BLS and instrumented with a multilevel monitoring system (Westbay TM ). Geophysical data, collected by the USGS, were used to characterize the subsurface geology and aquifer conditions. Natural gamma log measurements were used to assess sediment-layer thickness and location. Neutron and gamma-gamma source logs were used to confirm fractured and vesicular basalt identified for aquifer testing and multilevel monitoring well zone testing. Acoustic televiewer logs, collected for well USGS 149, were used to identify fractures and assess groundwater movement when compared with neutron measurements. Furthermore, gyroscopic deviation measurements were used to measure horizontal and vertical displacement for the constructed boreholes USGS 148A and USGS 149. A single-well aquifer test was done in well USGS 148A during November 6–7, 2019, to provide estimates of transmissivity and hydraulic conductivity. Estimates for transmissivity and hydraulic conductivity were 6.34×10 3 feet squared per day and 3.17 feet per day, respectively. The aquifer test was run overnight (21.3 hours) and measured drawdown was relatively small (0.09 ft) at sustained pumping rates ranging from 15.7 to 16.1 gallons per minute. The transmissivity estimates for well USGS 148A were slightly lower than those determined from previous aquifer tests for wells near the Materials and Fuels Complex, but well within range of other aquifer tests done at the INL. Water-quality samples, collected from well USGS 148A and from four zones in well USGS 149, were analyzed for cations, anions, metals, nutrients, volatile organic compounds, stable isotopes, and radionuclides. Water samples for most of the inorganic constituents showed similar chemistry in USGS 148A and all four zones in USGS 149. Water samples for stable isotopes of oxygen and hydrogen indicated some possible influence of irrigation on the water quality. Nitrate plus nitrite concentrations indicated influence from anthropogenic sources. The volatile organic compound and radiochemical data indicated that wastewater disposal practices at the Materials and Fuels Complex or from drilling had no detectable influence on these wells.

Idaho↗

Estimation of low-flow statistics at ungaged sites on streams in the Lower Hudson River Basin, New York, from data in geographic information systems

A variety of watershed properties available in 2015 from geographic information systems were tested in regression equations to estimate two commonly used statistical indices of the low flow of streams, namely the lowest flows averaged over 7 consecutive days that have a 1 in 10 and a 1 in 2 chance of not being exceeded in any given year (7-day, 10-year and 7-day, 2-year low flows). The equations were based on streamflow measurements in 51 watersheds in the Lower Hudson River Basin of New York during the years 1958–1978, when the number of streamflow measurement sites on unregulated streams was substantially greater than in subsequent years. These low-flow indices are chiefly a function of the area of surficial sand and gravel in the watershed; more precisely, 7-day, 10-year and 7-day, 2-year low flows both increase in proportion to the area of sand and gravel deposited by glacial meltwater, whereas 7-day, 2-year low flows also increase in proportion to the area of postglacial alluvium. Both low-flow statistics are also functions of mean annual runoff (a measure of net water input to the watershed from precipitation) and area of swamps and poorly drained soils in or adjacent to surficial sand and gravel (where groundwater recharge is unlikely and riparian water loss to evapotranspiration is substantial). Small but significant refinements in estimation accuracy resulted from the inclusion of two indices of stream geometry, channel slope and length, in the regression equations. Most of the regression analysis was undertaken with the ordinary least squares method, but four equations were replicated by using weighted least squares to provide a more realistic appraisal of the precision of low-flow estimates. The most accurate estimation equations tested in this study explain nearly 84 and 87 percent of the variation in 7-day, 10-year and 7-day, 2-year low flows, respectively, with standard errors of 0.032 and 0.050 cubic feet per second per square mile. The equations use natural values of streamflow and watershed properties; logarithmic transformations yielded less accurate equations inconsistent with some conceptualized relationships.

New York↗

Simulation of groundwater and surface-water flow in the upper Deschutes Basin, Oregon

This report describes a hydrologic model for the upper Deschutes Basin in central Oregon developed using the U.S. Geological Survey (USGS) integrated Groundwater and Surface-Water Flow model (GSFLOW). The upper Deschutes Basin, which drains much of the eastern side of the Cascade Range in Oregon, is underlain by large areas of permeable volcanic rock. That permeability, in combination with the large annual precipitation at high elevations, results in a substantial regional aquifer system and a stream system that is heavily groundwater dominated. The upper Deschutes Basin is also an area of expanding population and increasing water demand for public supply and agriculture. Surface water was largely developed for agricultural use by the mid-20th century, and is closed to additional appropriations. Consequently, water users look to groundwater to satisfy the growing demand. The well‑documented connection between groundwater and the stream system, and the institutional and legal restrictions on streamflow depletion by wells, resulted in the Oregon Water Resources Department (OWRD) instituting a process whereby additional groundwater pumping can be permitted only if the effects to streams are mitigated, for example, by reducing permitted surface-water diversions. Implementing such a program requires understanding of the spatial and temporal distribution of effects to streams from groundwater pumping. A groundwater model developed in the early 2000s by the USGS and OWRD has been used to provide insights into the distribution of streamflow depletion by wells, but lacks spatial resolution in sensitive headwaters and spring areas. The integrated model developed for this project, based largely on the earlier model, has a much finer grid spacing allowing resolution of sensitive headwater streams and important spring areas, and simulates a more complete set of surface processes as well as runoff and groundwater flow. In addition, the integrated model includes improved representation of subsurface geology and explicitly simulates the effects of hydrologically important fault zones not included in the previous model. The upper Deschutes Basin GSFLOW model was calibrated using an iterative trial and error approach using measured water-level elevations (water levels) from 800 wells, 144 of which have time series of 10 or more measurements. Streamflow was calibrated using data from 21 gage locations. At 14 locations where measured flows are heavily influenced by reservoir operations and irrigation diversions, so called “ naturalized ” flows, with the effects of reservoirs and diversion removed, developed by the Bureau of Reclamation, were used for calibration. Surface energy and moisture processes such as solar radiation, snow accumulation and melting, and evapotranspiration were calibrated using national datasets as well as data from long-term measurement sites in the basin. The calibrated Deschutes GSFLOW model requires daily precipitation, minimum and maximum air temperature data, and monthly data describing groundwater pumping and artificial recharge from leaking irrigation canals (which are a significant source of groundwater recharge). The calibrated model simulates the geographic distribution of hydraulic head over the 5,000 ft range measured in the basin, with a median absolute residual of about 53 ft. Temporal variations in head resulting from climate cycles, pumping, and canal leakage are well simulated over the model area. Simulated daily streamflow matches gaged flows or calculated naturalized flows for streams including the Crooked and Metolius Rivers, and lower parts of the mainstem Deschutes River. Seasonal patterns of runoff are less well fit in some upper basin streams. Annual water balances of streamflow are good over most of the model domain. Model fit and overall capabilities are appropriate for the objectives of the project. The integrated model results confirm findings from other studies and models indicating that most streamflow in the upper Deschutes Basin comes directly from groundwater discharge. The integrated model provides additional insights about the components of streamflow including direct groundwater discharge to streams, interflow, groundwater discharge to the land surface (Dunnian flow), and direct runoff (Hortonian flow). The new model provides improved capability for exploring the timing and distribution of streamflow capture by wells, and the hydrologic response to changes in other external stresses such as canal operation, irrigation, and drought. Because the model uses basic meteorological data as the primary input; and simulates surface energy and moisture balances, groundwater recharge and flow, and all components of streamflow; it is well suited for exploring the hydrologic response to climate change, although no such simulations are included in this report. The model was developed as a tool for future application; however, example simulations are provided in this report. In the example simulations, the model is used to explore the influence of well location and geologic structure on stream capture by pumping wells. Wells were simulated at three locations within a 12-mi area close to known groundwater discharge areas and crossed by a regional fault zone. Simulations indicate that the magnitude and timing of stream capture from pumping is largely controlled by the geographic location of the wells, but that faults can have a large influence on the propagation of pumping stresses.

Oregon↗

Stream-channel and watershed delineations and basin-characteristic measurements using lidar elevation data for small drainage basins within the Des Moines Lobe landform region in Iowa

Basin-characteristic measurements related to stream length, stream slope, stream density, and stream order have been identified as significant variables for estimation of flood, flow-duration, and low-flow discharges in Iowa. The placement of channel initiation points, however, has always been a matter of individual interpretation, leading to differences in stream definitions between analysts. This study investigated five different methods to define stream initiation using 3-meter light detection and ranging (lidar) digital elevation models (DEMs) data for 17 streamgages with drainage areas less than 50 square miles within the Des Moines Lobe landform region in north-central Iowa. Each DEM was hydrologically enforced and the five stream initiation methods were used to define channel initiation points and the downstream flow paths. The five different methods to define stream initiation were tested side-by-side for three watershed delineations: (1) the total drainage-area delineation, (2) an effective drainage-area delineation of basins based on a 2-percent annual exceedance probability (AEP) 12-hour rainfall, and (3) an effective drainage-area delineation based on a 20-percent AEP 12-hour rainfall. Generalized least squares regression analysis was used to develop a set of equations for sites in the Des Moines Lobe landform region for estimating discharges for ungaged stream sites with 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent AEPs. A total of 17 streamgages were included in the development of the regression equations. In addition, geographic information system software was used to measure 58 selected basin-characteristics for each streamgage. Results of the regression analyses of the 15 lidar datasets indicate that the datasets that produce regional regression equations (RREs) with the best overall predictive accuracy are the National Hydrographic Dataset, Iowa Department of Natural Resources, and profile curvature of 0.5 stream initiation methods combined with the 20-percent AEP 12-hour rainfall watershed delineation method. These RREs have a mean average standard error of prediction (SEP) for 4-, 2-, and 1-percent AEP discharges of 53.9 percent and a mean SEP for all eight AEPs of 55.5 percent. Compared to the RREs developed in this study using the basin characteristics from the U.S. Geological Survey StreamStats application, the lidar basin characteristics provide better overall predictive accuracy.

Iowa↗

Response in the water quality of Delavan Lake, Wisconsin, to changes in phosphorus loading—Setting new goals for loading from its drainage basin

During 1989–92, an extensive rehabilitation project was completed in and around Delavan Lake, Wisconsin, to improve the lake’s water quality. However, in 2016, the lake was listed by the Wisconsin Department of Natural Resources as impaired for excessive algal growth (high chlorophyll a concentrations), and high phosphorus input was listed as its likely cause. In addition, the recent (2017–21) mean summer water clarity (as measured with a Secchi disk) was shallower than the goal set by the community (3.0 meters). Based primarily on flow and water-quality data collected in Jackson Creek, which is the main tributary of the lake, the mean annual phosphorus loading to the lake during water years (WYs) 2017–21 was 6,570 kilograms per year (kg/yr), and 306 kg/yr came from uncontrollable sources (atmospheric deposition and groundwater). Phosphorus loading during these years was about 48 percent higher than the long-term mean loading from WY 1984 to WY 2021. Based on results from Canfield-Bachmann phosphorus models, Carlson trophic state index relations, and the Jones and Bachmann chlorophyll a relation, external phosphorus loading would need to be decreased from 6,570 to 5,270 kg/yr (a 21-percent reduction in the potentially controllable external phosphorus load from the base period of WYs 2017–21) for chlorophyll a concentrations greater than 20 micrograms per liter to be detected no more than 5.0 percent of the time (the Wisconsin Department of Natural Resources criterion for chlorophyll a impairment for the lake). Based on Carlson trophic state index relations, external loading would need to be decreased from 6,570 to 4,380 kg/yr (a 35-percent reduction in the potentially controllable external phosphorus load) for summer mean Secchi depths to increase to 3.0 meters. Therefore, for Delavan Lake to reach the water-quality criteria for impairment and the goals for all three water-quality constituents, a 35-percent reduction in the potentially controllable phosphorus load is needed, which equates to a reduction in total phosphorus loading from 6,570 to 4,380 kg/yr. A 35-percent reduction in phosphorus loading to improve the water quality of Delavan Lake is less than the 49-percent reduction in phosphorus loading required for the area near Delavan Lake to improve the water quality of the Rock River and its tributaries indicated in the Rock River total maximum daily load.

Wisconsin↗

Evaluation of the groundwater-flow model for the Ohio River alluvial aquifer near Carrollton, Kentucky, updated to conditions in September 2010

The Ohio River alluvial aquifer near Carrollton, Ky., is an important water resource for the cities of Carrollton and Ghent, as well as for several industries in the area. The groundwater of the aquifer is the primary source of drinking water in the region and a highly valued natural resource that attracts various water-dependent industries because of its quantity and quality. This report evaluates the performance of a numerical model of the groundwater-flow system in the Ohio River alluvial aquifer near Carrollton, Ky., published by the U.S. Geological Survey in 1999. The original model simulated conditions in November 1995 and was updated to simulate groundwater conditions estimated for September 2010. The files from the calibrated steady-state model of November 1995 conditions were imported into MODFLOW-2005 to update the model to conditions in September 2010. The model input files modified as part of this update were the well and recharge files. The design of the updated model and other input files are the same as the original model. The ability of the updated model to match hydrologic conditions for September 2010 was evaluated by comparing water levels measured in wells to those computed by the model. Water-level measurements were available for 48 wells in September 2010. Overall, the updated model underestimated the water levels at 36 of the 48 measured wells. The average difference between measured water levels and model-computed water levels was 3.4 feet and the maximum difference was 10.9 feet. The root-mean-square error of the simulation was 4.45 for all 48 measured water levels. The updated steady-state model could be improved by introducing more accurate and site-specific estimates of selected field parameters, refined model geometry, and additional numerical methods. Collection of field data to better estimate hydraulic parameters, together with continued review of available data and information from area well operators, could provide the model with revised estimates of conductance values for the riverbed and valley wall, hydraulic conductivities for the model layer, and target water levels for future simulations. Additional model layers, a redesigned model grid, and revised boundary conditions could provide a better framework for more accurate simulations. Additional numerical methods would identify possible parameter estimates and determine parameter sensitivities.

Kentucky↗

Ground-Water System in the Chimacum Creek Basin and Surface Water/Ground Water Interaction in Chimacum and Tarboo Creeks and the Big and Little Quilcene Rivers, Eastern Jefferson County, Washington

A detailed study of the ground-water system in the unconsolidated glacial deposits in the Chimacum Creek Basin and the interactions between surface water and ground water in four main drainage basins was conducted in eastern Jefferson County, Washington. The study will assist local watershed planners in assessing the status of the water resources and the potential effects of ground-water development on surface-water systems. A new surficial geologic map of the Chimacum Creek Basin and a series of hydrogeologic sections were developed by incorporating LIDAR imagery, existing map sources, and drillers' logs from 110 inventoried wells. The hydrogeologic framework outlined in the study will help characterize the occurrence of ground water in the unconsolidated glacial deposits and how it interacts with the surface-water system. Water levels measured throughout the study show that the altitude of the water table parallels the surface topography and ranges from 0 to 400 feet above the North American Vertical Datum of 1988 across the basin, and seasonal variations in precipitation due to natural cycles generally are on the order of 2 to 3 feet. Synoptic stream-discharge measurements and instream mini-piezometers and piezometers with nested temperature sensors provided additional data to refine the positions of gaining and losing reaches and delineate seasonal variations. Chimacum Creek generally gains water from the shallow ground-water system, except near the community of Chimacum where localized losses occur. In the lower portions of Chimacum Creek, gaining conditions dominate in the summer when creek stages are low and ground-water levels are high, and losing conditions dominate in the winter when creek stages are high relative to ground-water levels. In the Quilcene Bay area, three drainage basins were studied specifically to assess surface water/ground water interactions. The upper reaches of Tarboo Creek generally gain water from the shallow ground-water system throughout most of the year and the lower reaches have little or no gains. The Big Quilcene River generally gains water from the shallow ground-water system after it emerges from a bedrock canyon and loses water from the town of Quilcene to the mouth of the river in Quilcene Bay. The Little Quilcene River generally loses water to the shallow ground-water system, although two localized areas were found to have gaining conditions. The Big Quilcene and Little Quilcene Rivers incur significant losses on the alluvial plain at the head of Quilcene Bay. Each of the creeks examined had a unique pattern of gaining and losing reaches, owing to the hydraulic conductivity of the streambed material and the relative altitude of the surrounding water table. Although the magnitudes of gains and losses varied seasonally, the spatial distribution did not vary greatly, suggesting that patterns of gains and losses in surface-water systems depend greatly on the geology underlying the streambed.

Scientific Investigations Report↗

Annotated bibliography of scientific research relevant to oil and gas reclamation best management practices in the western United States, published from 1969 through 2020

Integrating recent scientific knowledge into management decisions supports effective natural resource management and can lead to better resource outcomes. However, finding and accessing scientific knowledge can be time consuming and costly. To assist in this process, the U.S. Geological Survey has created a series of annotated bibliographies on topics of management concern for lands in the western United States (U.S.). Oil and gas development on public lands is a long-standing and substantial component of local and regional economies and has expanded in recent decades, particularly on public lands in the western U.S. This development is associated with extensive networks of pipelines, roads, and processing facilities, across which reclamation is Federally mandated following initial well pad development (“interim” reclamation) and once resource extraction is complete (“final” reclamation). Reclamation is critical for recovering ecological services to energy-affected lands, including vegetation productivity, wildlife habitat, water and air quality, and soil stability (for example, resistance to wind and water erosion). However, reclamation of oil and gas affected lands in the western U.S. has proved challenging due to an array of regulatory and environmental factors, such as minimally developed soils, short growing seasons, herbivory, high winds, invasive species, rugged terrain, and in particular, arid climates associated with low total precipitation, high evapotranspiration rates, and highly variable precipitation patterns. We compiled and summarized journal articles, government reports, technical reports, proceedings, and theses and dissertations relevant to oil and gas reclamation. We constrained our search to products published on or before December 31, 2020 but did not limit our search by a starting date; the earliest product resulting from this effort was published in March 1969. Second, we manually scanned the last 15 years (2005-2020) of tables of contents in journals, bibliographies, and proceedings of which we were aware would contain articles highly relevant to this bibliography. We carried out the search for these products through multiple means: (1) performing a structured search of two reference databases, (2) examining articles published since 2005 in highly relevant scientific journals and conference proceedings, and (3) reviewing additional material suggested by authors of products identified in steps 1 and 2. Our search was intentionally broad in order to identify as much relevant work as possible, much of which is professionally applied and tested within the industry of oil and gas reclamation, but which remains unpublished in scientific journals. We refined the initial list of products by removing: (1) duplicates, (2) products not written in English, (3) products that were not relevant to the arid ecosystems of western North America, (4) products that were not released as research, data products, or review articles in journals or as formal scientific reports, and (5) products with data which were not relevant to reclamation of oil and gas-affected lands, or for which the study did not present new data, findings, or syntheses relevant to reclamation of oil and gas-affected lands. We summarized each product using a consistent structure (background, objectives, methods, location, findings, and implications) and assigned standardized management topics to each. Management topics are intended to aid online searching within the bibliography and are described in more detail in the Methods Section of this report; they include what type of disturbance the product addresses (well pads, mining, pipelines), what aspect of oil and gas reclamation they pertain to (practices, standards, monitoring), what type of data are present in the product (for instance soil or vegetation recovery data), and an indication if the product were from a source other than a published, peer-reviewed outlet (such as dissertations or unpublished professional reports – these are identified as grey literature). The review process for this annotated bibliography included an initial internal colleague review of each summary, requesting input on each summary from an author of the original product, and a formal peer-review. Our initial searches resulted in 3,197 total products, of which 290 met our criteria for inclusion. “Reclamation Practices” is by far the management topic most addressed, followed by “Reclamation Monitoring,” for example, products assessing what and how monitoring methods are used to track and measure reclamation outcome. This document may be accessed at https://doi.org/10.3133/ofr20231068 or from the U.S. Geological Survey Publication Warehouse ( https://pubs.usgs.gov/ ). The 1-page product summaries herein will also be used to create a bibliography at https://apps.usgs.gov/science-for-resource-managers that includes links to each original product, where available, and in which subject matter will be searchable by topic, location, and year. The studies compiled and summarized here may inform planning and management actions that seek to reclaim landscapes across the western U.S. which have been affected by oil and gas development.

Open-File Report↗

Geohydrology, geochemistry, and numerical simulation of groundwater flow and land subsidence in the Bicycle Basin, Fort Irwin National Training Center, California

Groundwater pumping from Bicycle Groundwater Basin (referred to as Bicycle Basin) in the Fort Irwin National Training Center, California, began in 1967. From 1967 to December 2010, about 46,000 acre-feet of water had been pumped from the basin and transported to the Irwin Basin. During this time, not only did water levels in the basin decline by as much as 100 feet, the quality of the groundwater pumped from the basin also deteriorated in some wells. The U.S. Geological Survey collected geohydrologic data from existing wells, test holes, and 16 additional monitoring wells installed at 6 sites in Bicycle Basin during 1992–2011 to determine the quantity and quality of groundwater available in the basin. Geophysical surveys, including electrical, gravity, and seismic refraction surveys, were completed to help determine the geometry of the structural basin, delineate depths to the interface between Quaternary and Tertiary rocks, map the depth to the water table, and used to develop a geohydrologic framework and groundwater-flow model for Bicycle Basin. Water samples were used to determine the groundwater quality in the basin and to delineate potential sources of poor-quality groundwater. Analysis of stable isotopes of oxygen and hydrogen in groundwater indicated that presentday precipitation is not a major source of recharge to the basin. Tritium and carbon-14 data indicated that most of the groundwater in the basin was recharged prior to 1952 and had an apparent age of 15,625–39,350 years. Natural recharge to the basin was not sufficient to replenish the groundwater pumped from the basin. Interferograms from synthetic aperture radar data (InSAR), analyzed to evaluate land-surface subsidence between 1993 and 2010, showed 0.23 to 1.1 feet of subsidence during this period near one production well north of Bicycle Lake (dry) playa. A groundwater-flow model of Bicycle Basin was developed and calibrated using groundwater levels for 1964– 2010, and a subsidence model using land-surface deformation data for 1993–2010. Between January 1967 and December 2010, the simulated total recharge from precipitation runoff and underflow from adjacent basins was about 5,100 acre-feet and pumpage from the Bicycle Basin was about 47,000 acrefeet of water. Total outflows exceeded natural recharge during this period, resulting in a net loss of about 42,100 acre-feet of groundwater storage in the basin. The Fort Irwin National Training Center is considering various groundwater-management options in the Bicycle Basin. The groundwater-flow model was used to (1) evaluate changes in groundwater levels and subsidence with the addition of capture and recharge of simulated runoff in retention basins (scenario 1) for predevelopment through 2010; (2) simulate a base case (scenario 2) for reference; and (3) compare projections of alternative future pumping strategies for 2011–60 (scenarios 3–5). Model results from the runoff-capture simulation (scenario 1) indicated that total recharge, including runoff captured using retention basins, locally increased water levels, which partially offset, but did not mitigate, groundwater depletion associated with pumping. Groundwater-storage depletion in scenario 1 was about 14 percent less than without runoff capture. Simulated-drawdown results in model layer 1 in the eastern part of the basin indicated that, because of the captured runoff, simulated heads were as much as 100 feet higher in December 2010 than prior to the onset of development in 1967. In contrast, simulated drawdown for model without runoff capture indicated that, without captured runoff, simulated heads for December 2010 in this area were 80–90 feet lower than during the predevelopment period. Subsidence was mitigated slightly in scenario 1 compared to without runoff capture; the largest decrease in subsidence at observation sites was about 0.07 feet.

California↗

Groundwater-quality monitoring program in Chester County, Pennsylvania, 1980-2008

The U.S. Geological Survey in cooperation with the Chester County Water Resources Authority and the Chester County Health Department began a groundwater-quality monitoring program in 1980 in Chester County, Pa., where a large percentage of the population relies on wells for drinking-water supply. This report documents the program and serves as a reference for data collected through the program from 1980 through 2008. The initial focus of the program was to collect data on groundwater quality near suspected localized sources of contamination, such as uncontrolled landfills and suspected industrial wastes, to determine if contaminants were present that might pose a health risk to those using the groundwater. Subsequently, the program was expanded to address the effects of widely distributed contaminant sources associated with agricultural and residential land uses on groundwater quality and to document naturally occurring constituents, such as radium, radon, and arsenic, that are potential hazards in drinking water. Since 2000, base-flow stream samples have been collected in addition to well-water and spring samples in a few small drainage areas to investigate the relation between groundwater quality measured in well samples and streams. The program has primarily consisted of spatial assessment with limited temporal data collected on groundwater quality. Most data were collected through the monitoring program for reconnaissance purposes to identify and locate groundwater-quality problems and generally were not intended for rigorous statistical analyses that might determine land-use or geochemical factors affecting groundwater quality in space or through time. Results of the program found several contaminants associated with various land uses and human activities in groundwater in Chester County. Volatile organic compounds (such as trichloroethylene) were measured in groundwater near suspected localized contaminant sources in concentrations that exceeded drinking-water standards. Groundwater in some agricultural areas had concentrations of nitrate and some pesticides that exceeded drinking-water standards. Elevated concentrations of chloride were measured near salt storage areas and highways. Formaldehyde was detected in groundwater near cemeteries. In residential areas with on-site wastewater disposal, effects on groundwater quality included elevated nitrate concentrations and low concentrations of volatile organic compounds and wastewater compounds, such as antibiotics and detergents. Base-flow samples indicated that groundwater discharge to streams carried contaminants such as nitrate, pesticides, wastewater compounds, and other contaminants. Radionuclides, including radium-226, radium-228, radium-224, and radon-222, and gross alpha-particle activity were measured in groundwater at levels above established and proposed drinking-water standards in some geologic units, particularly in quartzite and quartzite schists. Arsenic concentrations above drinking-water standards were measured in a few samples and were most likely to occur in groundwater in the shales and sandstones in the northern part of the county. Other potential natural hazards, such as lead from aquifer materials or leached from plumbing because of pH, were present in concentrations above drinking-water standards infrequently (less than 10 percent of samples). Limited temporal sampling suggested that chloride concentrations in groundwater increased in the county since the program began in 1980 through 2008, reflecting increasing population and urbanization in that period.

Scientific Investigations Report↗

Source contributions to suspended sediment and particulate selenium export from the Loutsenhizer Arroyo and Sunflower Drain watersheds in Colorado

Selenium in aquatic ecosystems of the lower Gunnison River Basin in Colorado is affecting the recovery of populations of endangered, native fish species. Dietary exposure is the primary pathway for bioaccumulation of selenium in fish, and particulate selenium can be consumed directly by fish or by the invertebrates on which fish feed. Although selenium can be incorporated into particulate matter via biogeochemical processes, particulate selenium can also enter aquatic ecosystems of the lower Gunnison River Basin from sediments derived from the selenium-rich Mancos Shale. The U.S. Geological Survey, in cooperation with the Colorado Water Conservation Board, conducted this study during 2018–19 to identify sources of selenium-rich suspended sediments from two watersheds underlain by Mancos Shale: Loutsenhizer Arroyo and Sunflower Drain, which is a locally known agricultural drainage near the municipality of Delta, Colorado. A multipronged approach (fieldwork, laboratory work, and computer modeling) referred to as “sediment fingerprinting” was used to evaluate sources of suspended sediments in the streams flowing out of the two studied watersheds. Four potential source types for suspended sediments were identified and sampled (using soil plugs) within the watersheds: rangelands, agricultural fields, arroyo walls, and streambanks. The sediment fingerprinting approach used elemental concentrations and naturally occurring fallout radionuclides as tracers to apportion percent contributions from the four source types of suspended sediments found in streamflow from both watersheds. To determine the dominant sources of suspended sediment in streamflow from both watersheds, a mathematical “unmixing” model was used. Unmixing models apportion source percentages to samples of material in which those sources are mixed. These models used elemental and isotopic data in the suspended sediments to unmix them into proportional contributions from source types. The results indicated that arroyo walls and streambanks generally dominated as sources of the suspended sediment. Arroyo walls and streambanks were channel-adjacent sources, with sediments mobilized by water flowing within the stream channel. These sources accounted for greater than 50 percent of suspended sediment in all but one sample and accounted for 100 percent of suspended sediment in 5 of the 11 samples collected. Rangeland and agricultural field sources were located in uplands outside of stream channels and were detected more often during the non-irrigation season. Rangeland and agricultural field sources each were found in 5 of the 11 samples collected. Concentrations of selenium in sediment-source samples were comparatively greater in streambanks and lower in rangelands, with agricultural fields and arroyo walls being intermediate. As a result, source apportionments for particulate selenium skewed towards sources adjacent to stream channels more than for suspended sediments. Water imports for irrigation have changed the hydrology of the watersheds, and a notable fraction of imported water passes through the watersheds rapidly. The rapid flowthrough water during the irrigation season likely contributes heavily to sediment erosion and transport in Loutsenhizer Arroyo and Sunflower Drain, particularly from channel-adjacent sources of sediment. Decreases in irrigation season streamflow, at least in Loutsenhizer Arroyo, may have decreased sediment erosion and transport during the 2018–20 irrigation seasons compared to the 2015–17 seasons.

Colorado↗

Water-balance techniques for determining available soil-water storage for selected sandy and clay soil study sites in Cass County, North Dakota, 2016–17

The U.S. Geological Survey, in cooperation with the U.S. Department of Agriculture Natural Resources Conservation Service, collected field and remotely sensed data on precipitation, evapotranspiration (ET), and soil-water content to determine available soil-water storage (AWS) at six study sites on sandy and clay soils in Cass County, North Dakota. Data were collected at all the study sites from May 1–October 31, 2016, and from May 1–October 24, 2017. Estimated daily AWS was determined using daily meteorological and potential evapotranspiration (PET) data obtained from various climate stations, and estimated monthly AWS was determined using monthly meteorological and PET data and monthly ET data determined using the Operational Simplified Surface Energy Balance model. AWS during 2016 and 2017 was determined at daily and monthly time steps because of data availability and to assess results using varying time steps. Comparisons of measured and estimated daily values of AWS at the Brewer Lake site indicated poor agreement during May–October 2016 and May–October 2017. Comparisons of measured and estimated daily values of AWS at the Embden East and Embden West sites indicated poor and fair agreement respectively. At the Lynchburg Crop and Lynchburg Grass sites, comparisons of measured and estimated daily values of AWS indicated fair and good relations, respectively, even with the possible effects of soil cracks. Mean estimated values of daily runoff plus soil percolation for the four sandy soil sites indicated that a maximum of about 19 percent of the estimated runoff plus soil percolation could be considered runoff and that the remaining 81 percent could be considered soil percolation, and for the two clay soil sites about 13 percent of the runoff plus soil percolation could have been considered runoff and about 87 percent could have been considered soil percolation. Results indicated little difference between using monthly PET or monthly ET in water-balance equations to estimate monthly AWS for the grouped sandy soil sites, and only slightly better results were obtained using monthly PET than monthly ET to estimate monthly AWS for the grouped clay soil study sites. Overall, the monthly water-balance models did not perform as well as the daily water-balance models for determining AWS at the six study sites. Additional data collection from a longer-period study and adjustments to the models may improve results from the monthly water-balance techniques.

North Dakota↗

Trends in environmental, anthropogenic, and water-quality characteristics in the upper White River Basin, Indiana

The U.S. Geological Survey (USGS), in cooperation with The Nature Conservancy, undertook a study to update and extend results from a previous study (Koltun, 2019, https://doi.org/10.3133/sir20195119 ), using data from 3 additional years and newer estimation methods. Koltun (2019) assessed trends in streamflow, precipitation, and estimated annual mean concentrations and flux of nitrate plus nitrite, total Kjeldahl nitrogen, total phosphorus, and total suspended solids (TSS) for USGS streamflow gages on the upper White River at Muncie, near Nora, and near Centerton, Indiana. Annual mean and maximum daily streamflows had statistically significant upward trends at all study gages between water years 1978 and 2020. An abrupt increase in streamflow occurred around water year 2001. Annual total precipitation at the Indianapolis International Airport increased between calendar years 1932 and 2020 at an average rate of 0.089 inches per year. The current study assessed the magnitude, direction, and likelihood of change in flow-normalized concentrations and flux of TSS, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen between water years 1997 and 2019. With two exceptions, concentration and flux changes that were statistically significant in Koltun (2019, https://doi.org/10.3133/sir20195119 ), which reported changes between water years 1997 and 2017, still have the same statistically significant change directions. The reliability of the current trend result for TSS is uncertain because of a large gap in the TSS record for the Centerton gage. For each constituent, spatial patterns were examined in the sampled distribution of nutrient and TSS concentration data from 20 mainstem, tributary, and distributary locations in the upper White River Basin. The largest median concentrations of TSS, total phosphorus, and total Kjeldahl nitrogen were associated with mainstem upper White River sites downstream from Indianapolis. The median total phosphorus and total Kjeldahl nitrogen concentrations were elevated relative to bracketing upstream/downstream mainstem sites at the upper White River site immediately downstream from Muncie. Data on several anthropogenic factors that could influence the concentrations and fluxes of nutrients and TSS were gathered and analyzed to better understand the factors’ spatial and temporal variations. Those anthropogenic factors included population, land cover, cropping and operational tillage practices, fertilizer application, and upgrades to wastewater treatment systems and delivery processes.

Indiana↗

Simulation of groundwater flow in the regional aquifer system on Long Island, New York, for pumping and recharge conditions in 2005–15

A three-dimensional groundwater-flow model was developed for the aquifer system of Long Island, New York, to evaluate (1) responses of the hydrologic system to changes in natural and anthropogenic hydraulic stresses, (2) the subsurface distribution of groundwater age, and (3) the regional-scale distribution of groundwater travel times and the source of water to fresh surface waters and coastal receiving waters. The model also provides the groundwater flow components used to define model boundaries for possible inset models used for local-scale analyses. The three-dimensional, groundwater flow model developed for this investigation uses the numerical code MODFLOW–NWT to represent steady-state conditions for average groundwater pumping and aquifer recharge for 2005–15. The particle-tracking algorithm MODPATH, which simulates advective transport in the aquifer, was used to estimate groundwater age, delineate the areas at the water table that contribute recharge to coastal and freshwater bodies, and estimate total travel times of water from the water table to discharge locations. A three-dimensional, 1-meter (3.3-foot) topobathymetric model was used to determine land-surface altitudes for the island and seabed altitudes for the surrounding coastal waters. The mapped extents and surface altitudes of major geologic units were compiled and used to develop a three-dimensional hydrogeologic framework of the aquifer system, including aquifers and confining units. Lithologic data from deep boreholes and previous aquifer-test results were used to estimate the three-dimensional distribution of hydraulic conductivity in principal aquifers. Natural recharge from precipitation was estimated for 2005–15 using a modified Thornthwaite-Mather methodology as implemented in a soil-water balance model. Components of anthropogenic recharge—wastewater return flow, storm water inflow, and inflow from leaky infrastructure—also were estimated for 2005–15. Groundwater withdrawals for various sources, including public water supply, industrial, remediation, and agricultural, were compiled or estimated for the same period. These data were incorporated into the model development to represent the aquifer system geometry, boundaries, and initial hydraulic properties of the regional aquifers and confining units within the Long Island aquifer system. Average hydraulic conditions—water levels and streamflows—for 2005–15 were estimated using existing data from the U.S. Geological Survey National Water Information System database. Model inputs were adjusted to best match average hydrologic conditions using inverse methods as implemented in the parameter-estimating software PEST. The calibrated model was used to simulate average hydrologic conditions in the aquifer system for 2005–15. About 656 cubic feet per second of water was withdrawn on average annually for 2005–15 for water supply and an average of about 349 cubic feet per second of water recharged the aquifer annually from return flow and leaky infrastructure. Parts of New York City have drawdowns exceeding 25 feet, mostly because of urbanization and associated large decreases in recharge rates. Large areas in the western part of the island have drawdowns exceeding 10 feet, mostly from large groundwater withdrawals and sewering, which largely eliminates wastewater return flow. Water-table altitudes in eastern parts of the island increased by more than 2 feet in some areas as a result of wastewater return flow in unsewered areas and changes in land use. Changes in streamflows show a similar pattern as water-table altitudes. Streamflows generally decrease in western parts of the island where there are large drawdowns and increase in eastern parts of the island where water-table altitudes increase.

New York↗