Search USGS⌕ Search

SEARCH · Search USGS

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 847 records · Page 47Linked to original sources

Assessment, water-quality trends, and options for remediation of acidic drainage from abandoned coal mines near Huntsville, Missouri, 2003-2004

Water from abandoned underground coal mines acidifies receiving streams in the Sugar Creek Basin and Mitchell Mine Basin near Huntsville, Missouri. A 4.35-kilometer (2.7-mile) reach of Sugar Creek has been classified as impaired based on Missouri's Water Quality Standards because of small pH values [< (less than) 6.5]. Samples collected from Sugar Creek from July 2003 to June 2004 did not have pH values outside of the specified range of 6.5 to 9.0. However, large concentrations of iron [416 to 2,320 mg/L (milligrams per liter)], manganese (8.36 to 33.5 mg/L), aluminum (0.870 to 428 mg/L), and sulfate (2,990 to 13,700 mg/L) in acidic mine drainage (AMD) from two mine springs as well as small and diffuse seeps were observed to have an effect on water quality in Sugar Creek. Metal and sulfate loads increased and pH decreased immediately downstream from Sugar Creek's confluence with the Calfee Slope and Huntsville Gob drainages that discharge AMD into Sugar Creek. Similar effects were observed in the Mitchell Mine drainage that receives AMD from a large mine spring. Comparisons of water-quality samples from this study and two previous studies by the U.S. Geological Survey in 1987-1988 and the Missouri Department of Natural Resources in 2000-2002 indicate that AMD generation in the Sugar Creek Basin and Mitchell Mine Basin is declining, but the data are insufficient to quantify any trends or time frame. AMD samples from the largest mine spring in the Calfee Slope subbasin indicated a modest but significant increase in median pH from 4.8 to 5.2 using the Wilcoxan rank-sum test (p <0.05) and a decrease in median specific conductance from 5,000 to 3,540 ?S/cm (microsiemens per centimeter at 25 degrees Celsius) during a 17-year period. AMD samples from the largest mine spring in the Mitchell Mine Basin indicated an increase in median pH values from 5.6 to 6.0 and a decrease in median specific conductance from 3,050 to 2,450 ?S/cm during the same period. Remediation of AMD at or near the sites of the three largest mine springs is geochemically feasible based on alkalinity addition rates and increased pH determined by cubitainer experiments and geochemical mixing experiments using the computer model PHREEQCI. Alkalinity values for seven cubitainer experiments conducted to simulate anoxic treatment options exceeded the targeted value for alkalinity [90 mg/L as calcium carbonate (CaCO3)] specified in Missouri's Total Maximum Daily Load program by 18 percent or more, but maximum pH values were between 6.2 and 6.3, which is less than the targeted pH value of 6.5. Treatment of AMD by mixing with stream water or sewage effluent can further increase pH as indicated by geochemical modeling, but will not totally achieve water-quality goals because of limited discharges. A combination of treatments including settling ponds, oxic or anoxic limestone drains, and possibly successive alkalinity producing systems to remediate AMD will likely be required in the Sugar Creek Basin and Mitchell Mine Basin to consistently meet Missouri's Water Quality Standards.

Scientific Investigations Report↗

Temporal changes in nutrient concentrations in the Lower Grand River and selected drainage basins, Missouri and Iowa, during the Mississippi River Basin Healthy Watersheds Initiative (2010–23)

This report describes a cooperative study by the U.S. Geological Survey and Missouri Department of Natural Resources that evaluated temporal changes in total nitrogen (TN) and total phosphorus (TP) concentrations in the Lower Grand River hydrologic unit. The study focused on trends since 2010, when the basin was designated as a priority drainage basin of the Mississippi River Basin Healthy Watersheds Initiative (MRBI). At three local drainage basins within the Lower Grand hydrological unit (MRBI sites), stream nutrient trends were evaluated using flow-adjusted (FA) TN and TP concentrations for water years 2011 through 2023. FATN concentration trends were not statistically significant for any MRBI site. One site (site 2) showed a statistically significant increasing trend in FATP concentration, indicating a possible increase in phosphorus sources in parts of the basin. Overall, streamflow variability appeared to be the dominant factor affecting nutrient concentrations at MRBI sites. At five regional drainage basins, including the Grand River and nearby rivers with data from 1994 through 2023 (long-term sites), annual flow-normalized (FN) TN and TP concentrations were evaluated for trends before (water years 2000–10) and during (water years 2010–23) the MRBI. For water years 2010 through 2023, annual FNTN and FNTP concentrations decreased in the Grand River, as well as in the Nodaway and Chariton Rivers, which were not targeted by the MRBI. The Grand River (site 9) reversed from increasing to decreasing FNTP concentrations after 2010, with a 26-percent reduction. Annual FNTN and FNTP concentrations also decreased at the Missouri River sites. While nutrient reductions in the Grand River may reflect the effects of implemented conservation practices, similar trends in nearby, nontargeted rivers and the absence of strong decreasing trends at MRBI sites suggest that broader regional factors, instead of or in addition to MRBI efforts, may have contributed to nutrient reductions in the Grand River.

Iowa, Missouri↗

Community estimate of global glacier mass changes from 2000 to 2023

Glaciers are indicators of ongoing anthropogenic climate change 1 . Their melting leads to increased local geohazards 2 , and impacts marine 3 and terrestrial 4 , 5 ecosystems, regional freshwater resources 6 , and both global water and energy cycles 7 , 8 . Together with the Greenland and Antarctic ice sheets, glaciers are essential drivers of present 9 , 10 and future 11 , 12 , 13 sea-level rise. Previous assessments of global glacier mass changes have been hampered by spatial and temporal limitations and the heterogeneity of existing data series 14 , 15 , 16 . Here we show in an intercomparison exercise that glaciers worldwide lost 273 ± 16 gigatonnes in mass annually from 2000 to 2023, with an increase of 36 ± 10% from the first (2000–2011) to the second (2012–2023) half of the period. Since 2000, glaciers have lost between 2% and 39% of their ice regionally and about 5% globally. Glacier mass loss is about 18% larger than the loss from the Greenland Ice Sheet and more than twice that from the Antarctic Ice Sheet 17 . Our results arise from a scientific community effort to collect, homogenize, combine and analyse glacier mass changes from in situ and remote-sensing observations. Although our estimates are in agreement with findings from previous assessments 14 , 15 , 16 at a global scale, we found some large regional deviations owing to systematic differences among observation methods. Our results provide a refined baseline for better understanding observational differences and for calibrating model ensembles 12 , 16 , 18 , which will help to narrow projection uncertainty for the twenty-first century 11 , 12 , 18 .

Nature↗

Topographic change detection at select archeological sites in Grand Canyon National Park, Arizona, 2007-2010

Human occupation in Grand Canyon, Arizona, dates from at least 11,000 years before present to the modern era. For most of this period, the only evidence of human occupation in this iconic landscape is provided by archeological sites. Because of the dynamic nature of this environment, many archeological sites are subject to relatively rapid topographic change. Quantifying the extent, magnitude, and cause of such change is important for monitoring and managing these archeological sites. Such quantification is necessary to help inform the continuing debate on whether and how controlled releases from Glen Canyon Dam, located immediately upstream of Grand Canyon National Park, are affecting site erosion rates, artifact transport, and archeological resource preservation along the Colorado River in Grand Canyon. Although long-term topographic change resulting from a variety of natural processes is inherent in the Grand Canyon region, continued erosion of archeological sites threatens both the archeological resources and our future ability to study evidence of past cultural habitation. Thus, this subject is of considerable interest to National Park Service managers and other stakeholders in the Glen Canyon Dam Adaptive Management Program. Understanding the causes and effects of archeological site erosion requires a knowledge of several factors, including the location, timing, and magnitude of the changes occurring in relation to archeological resources, the rates of change, and the relative contribution of potential causes. These potential causes include sediment depletion associated with managed flows from Glen Canyon Dam, site-specific weather and overland flow patterns, visitor impacts, and long-term regional climate change. To obtain this information, highly accurate, spatially specific data are needed from sites undergoing change. Using terrestrial lidar techniques, and building upon three previous surveys of archeological sites performed in 2006 and 2007, we collected two new datasets in April and September 2010 and processed and improved upon existing methods to generate high-accuracy (3 to 5 cm vertical change threshold) topographic change-detection maps for 10 survey areas encompassing 9 archeological sites along the Colorado River corridor. We also used terrestrial lidar techniques to investigate several other metrics for studying archeological site stability, including monitoring cultural structures and artifacts and remotely measuring cryptobiotic soil crust areas. Our topographic change results indicate that 9 of 10 survey areas showed signs of either erosion, deposition, or both during the 2007&ndash;2010 time interval and that these changes can be linked to a variety of geomorphic processes, primarily overland flow gullying and aeolian sand transport. In several cases, large (>50 cm) vertical change occurred, and in one case, more than 100 m 3 of sediment was eroded. Further, for all sites monitored throughout the river corridor during this time period, the overall signal was related to erosion rather than deposition. These results highlight the potential for rapid archeological site change in Grand Canyon. Whereas the topographic change results presented herein provide the highest level of change detection yet performed on entire archeological sites in Grand Canyon, additional work in combining these results with site-specific weather, hydrology, and geomorphology data is needed to provide a more thorough understanding of the causes of the documented topographic changes. Linking lidar-derived measurements of topographic changes with these other data sources should provide land managers with a scientific basis for making management decisions regarding archeological resources in Grand Canyon National Park and assist in answering open questions regarding the influence that sediment-depleted flows from Glen Canyon Dam have on archeological site stability.

Arizona↗

Framework for regional synthesis of water-quality data for the glacial aquifer system in the United States

The glacial aquifer system is the largest principal aquifer in aerial extent and ground-water use for public supply in the United States. A principal aquifer is defined as a regionally extensive aquifer or aquifer system that has the potential to be used as a source of potable water (U.S. Geological Survey, 2003). Multiple aquifers often are grouped into large, extensive aquifer systems such as the glacial aquifer system. The glacial aquifer system is considered here to include all unconsolidated aquifers above bedrock north of the line of continental glaciation throughout the country (fig. 1). Total withdrawals from the glacial aquifer system were 3,560 million gallons per day in 2000, which constitutes almost 5 percent of total withdrawals from all aquifers in the United States (Maupin and Barber, 2005). Approximately 41 million people relied on the glacial aquifer for public supply and domestic use in 2000. The U.S. Geological Survey National Water-Quality Assessment (NAWQA) Program began assessing the glacial aquifer system in 1991. The assessment of water-quality data on a regional scale, such as the glacial aquifer system, is coincident with the regional framework established by the Regional Aquifer-System Analysis Program (RASA) (Sun and others, 1997). From 1978 to 1995, the RASA Program systematically evaluated 25 of the Nation's most important groundwater systems including studies in the glacial aquifer system in the northeast, Midwest, and northern Midwest United States. The NAWQA Program is building on the work of the RASA Program to study the water quality of 16 of the most important ground-water systems (Lapham and others, 2005). Over 1,700 water-quality samples have been collected by the NAWQA Program from 1991 to 2004 to assess the glacial aquifer system. This large data set is unique in that the samples have been collected using a consistent sampling protocol, and multiple nested samples. The nested samples address the recently recharged shallow ground water, deeper water from principal aquifers often used for domestic supply, and source water used for public supplies within the glacial aquifer system. Information concerning the NAWQA Program including study unit boundaries is shown in figure 1 (Lapham and others, 2005). A framework for comparison of water quality across the glacial aquifer system has been developed based on two primary characteristics: intrinsic susceptibility and vulnerability. Intrinsic susceptibility, which is a measure of the ease at which water enters and moves through aquifer material, is a characteristic of the aquifer and overlying material and of the hydrologic conditions. Intrinsic susceptibility is independent of the chemical characteristics of the contaminant and its sources. In this way, intrinsic susceptibility assessments do not target specific natural or anthropogenic sources of contamination but instead consider only the physical factors affecting the flow of water to, and through the ground-water resource (Focazio and others, 2002). On a regional scale, intrinsic susceptibility is represented by the spatial distribution of fine- or coarse-grained material at the land surface, and the physical setting of the aquifer system. Vulnerability, which is a function of both intrinsic susceptibility and the proximity and characteristics of contaminant sources, includes consideration of features related to anthropogenic sources of contaminants, such as the character of the upgradient land use (for example, urban, agricultural, undeveloped, and others); as well as features related to natural sources of contaminants, such as the mineralogy of the aquifer material or the geochemical conditions within the aquifer system. The framework helps categorize this large region into areas of similar hydrogeologic characteristics for which water quality can be compared. The purpose of this report is to describe this framework and how it will be used for regional synthesis of water-quality da

Scientific Investigations Report↗

Method for estimating potential wetland extent by utilizing streamflow statistics and flood-inundation mapping techniques: Pilot study for land along the Wabash River near Terre Haute, Indiana

Potential wetland extents were estimated for a 14-mile reach of the Wabash River near Terre Haute, Indiana. This pilot study was completed by the U.S. Geological Survey in cooperation with the U.S. Department of Agriculture, Natural Resources Conservation Service (NRCS). The study showed that potential wetland extents can be estimated by analyzing streamflow statistics with the available streamgage data, calculating the approximate water-surface elevation along the river, and generating maps by use of flood-inundation mapping techniques. Planning successful restorations for Wetland Reserve Program (WRP) easements requires a determination of areas that show evidence of being in a zone prone to sustained or frequent flooding. Zone determinations of this type are used by WRP planners to define the actively inundated area and make decisions on restoration-practice installation. According to WRP planning guidelines, a site needs to show evidence of being in an "inundation zone" that is prone to sustained or frequent flooding for a period of 7 consecutive days at least once every 2 years on average in order to meet the planning criteria for determining a wetland for a restoration in agricultural land. By calculating the annual highest 7-consecutive-day mean discharge with a 2-year recurrence interval (7MQ2) at a streamgage on the basis of available streamflow data, one can determine the water-surface elevation corresponding to the calculated flow that defines the estimated inundation zone along the river. By using the estimated water-surface elevation ("inundation elevation") along the river, an approximate extent of potential wetland for a restoration in agricultural land can be mapped. As part of the pilot study, a set of maps representing the estimated potential wetland extents was generated in a geographic information system (GIS) application by combining (1) a digital water-surface plane representing the surface of inundation elevation that sloped in the downstream direction of flow and (2) land-surface elevation data. These map products from the pilot study will aid the NRCS and its partners with the onsite inundation-zone verification in agricultural land for a potential restoration and will assist in determining at what elevation to plant hardwood trees for increased survivability on ground above frequently flooded terraces.

Indiana↗

Comparison of Hydrologic and Water-Quality Characteristics of Two Native Tallgrass Prairie Streams with Agricultural Streams in Missouri and Kansas

This report presents the results of a study by the U.S. Geological Survey, in cooperation with the Missouri Department of Natural Resources, to analyze and compare hydrologic and water-quality characteristics of tallgrass prairie and agricultural basins located within the historical distribution of tallgrass prairie in Missouri and Kansas. Streamflow and water-quality data from two remnant, tallgrass prairie basins (East Drywood Creek at Prairie State Park, Missouri, and Kings Creek near Manhattan, Kansas) were compared to similar data from agricultural basins in Missouri and Kansas. Prairie streams, especially Kings Creek in eastern Kansas, received a higher percentage of base flow and a lower percentage of direct runoff than similar-sized agricultural streams in the region. A larger contribution of direct runoff from the agricultural streams made them much flashier than prairie streams. During 22 years of record, the Kings Creek base-flow component averaged 66 percent of total flow, but base flow was only 16 to 26 percent of flows at agricultural sites of various record periods. The large base-flow component likely is the result of greater infiltration of precipitation in prairie soils and the resulting greater contribution of groundwater to streamflow. The 1- and 3-day annual maximum flows were significantly greater at three agricultural sites than at Kings Creek. The effects of flashier agricultural streams on native aquatic biota are unknown, but may be an important factor in the sustainability of some native aquatic species. There were no significant differences in the distribution of dissolved-oxygen concentrations at prairie and agricultural sites, and some samples from most sites fell below the 5 milligrams per liter Missouri and Kansas standard for the protection of aquatic life. More than 10 percent of samples from the East Drywood Creek prairie stream were less than this standard. These data indicate low dissolved-oxygen concentrations during summer low-flow periods may be a natural phenomenon for small prairie streams in the Osage Plains. Nutrient concentrations including total nitrogen, ammonia, nitrate, and total phosphorus were significantly less in base-flow and runoff samples from prairie streams than from agricultural streams. The total nitrogen concentration at all sites other than one of two prairie sampling sites were, on occasion, above the U.S. Environmental Protection Agency recommended criterion for total nitrogen for the prevention of nutrient enrichment, and typically were above this recommended criterion in runoff samples at all sites. Nitrate and total phosphorus concentrations in samples from the prairie streams generally were below the U.S. Environmental Protection Agency recommended nutrient criteria in base-flow and runoff samples, whereas samples from agricultural sites generally were below the criteria in base-flow samples and generally above in runoff samples. The lower concentrations of nutrient species in prairie streams is likely because prairies are not fertilized like agricultural basins and prairie basins are able to retain nutrients better than agricultural basins. This retention is enhanced by increased infiltration of precipitation into the prairie soils, decreased surface runoff, and likely less erosion than in agricultural basins. Streamflow in the small native prairie streams had more days of zero flow and lower streamflow yields than similar-sized agricultural streams. The prairie streams were at zero flow about 50 percent of the time, and the agricultural streams were at zero flow 25 to 35 percent of the time. Characteristics of the prairie basins that could account for the greater periods of zero flow and lower yields when compared to agricultural streams include greater infiltration, greater interception and evapotranspiration, shallower soils, and possible greater seepage losses in the prairie basins. Another difference between the prairie and agricultural strea

Scientific Investigations Report↗

Water quality of streams draining abandoned and reclaimed mined lands in the Kantishna Hills area, Denali National Park and Preserve, Alaska, 2008–11

The Kantishna Hills are an area of low elevation mountains in the northwest part of Denali National Park and Preserve, Alaska. Streams draining the Kantishna Hills are clearwater streams that support several species of fish and are derived from rain, snowmelt, and subsurface aquifers. However, the water quality of many of these streams has been degraded by mining. Past mining practices generated acid mine drainage and excessive sediment loads that affected water quality and aquatic habitat. Because recovery through natural processes is limited owing to a short growing season, several reclamation projects have been implemented on several streams in the Kantishna Hills region. To assess the current water quality of streams in the Kantishna Hills area and to determine if reclamation efforts have improved water quality, a cooperative study between the U.S. Geological Survey and the National Park Service was undertaken during 2008-11. High levels of turbidity, an indicator of high concentrations of suspended sediment, were documented in water-quality data collected in the mid-1980s when mining was active. Mining ceased in 1985 and water-quality data collected during this study indicate that levels of turbidity have declined significantly. Turbidity levels generally were less than 2 Formazin Nephelometric Units and suspended sediment concentrations generally were less than 1 milligram per liter during the current study. Daily turbidity data at Rock Creek, an unmined stream, and at Caribou Creek, a mined stream, documented nearly identical patterns of turbidity in 2009, indicating that reclamation as well as natural revegetation in mined streams has improved water quality. Specific conductance and concentrations of dissolved solids and major ions were highest from streams that had been mined. Most of these streams flow into Moose Creek, which functions as an integrator stream, and dilutes the specific conductance and ion concentrations. Calcium and magnesium are the dominant cations, and bicarbonate and sulfate are the dominant anions. Water samples indicate that the water from Rock Creek, Moose Creek, Slate Creek, and Eldorado Creek is a calcium bicarbonate-type water. The remaining sites are a calcium sulfate type water. U.S. Environmental Protection Agency guidelines for arsenic and antimony in drinking water were exceeded in water at Slate Creek and Eureka Creek. Concentrations of arsenic, cadmium, chromium, copper, lead, nickel, and zinc in streambed sediments at many sites exceed sediment quality guideline thresholds that could be toxic to aquatic life. However, assessment of these concentrations, along with the level of organic carbon detected in the sediment, indicate that only concentrations of arsenic and chromium may be toxic to aquatic life at many sites. In 2008 and 2009, 104 macroinvertebrate taxa and 164 algae taxa were identified from samples collected from seven sites. Of the macroinvertebrates, 86 percent were insects and most of the algae consisted of diatoms. Based on the National Community Index, Rock Creek, a reference site, and Caribou Creek, and a mined stream that had undergone some reclamation, exhibited the best overall stream conditions; whereas Slate Creek and Friday Creek, two small streams that were mined extensively, exhibited the worst stream conditions. A non-metric multi-dimensional scaling analysis of the macroinvertebrate and algae data showed a distinct grouping between the 2008 and 2009 samples, likely because of differences between a wet, cool summer in 2008 and a dry, warm summer in 2009.

Alaska↗

Predicted nitrate and arsenic concentrations in basin-fill aquifers of the Southwestern United States

The National Water-Quality Assessment (NAWQA) Program of the U.S. Geological Survey (USGS) is conducting a regional analysis of water quality in the principal aquifer systems across the United States. The Southwest Principal Aquifers (SWPA) study is building a better understanding of the susceptibility and vulnerability of basin-fill aquifers in the region to groundwater contamination by synthesizing baseline knowledge of groundwater-quality conditions in 16 basins previously studied by the NAWQA Program. The improved understanding of aquifer susceptibility and vulnerability to contamination is assisting in the development of tools that water managers can use to assess and protect the quality of groundwater resources. Human-health concerns and economic considerations associated with meeting drinking-water standards motivated a study of the vulnerability of basin-fill aquifers to nitrate con­tamination and arsenic enrichment in the southwestern United States. Statistical models were developed by using the random forest classifier algorithm to predict concentrations of nitrate and arsenic across a model grid that represents about 190,600 square miles of basin-fill aquifers in parts of Arizona, California, Colorado, Nevada, New Mexico, and Utah. The statistical models, referred to as classifiers, reflect natural and human-related factors that affect aquifer vulnerability to contamina­tion and relate nitrate and arsenic concentrations to explana­tory variables representing local- and basin-scale measures of source, aquifer susceptibility, and geochemical conditions. The classifiers were unbiased and fit the observed data well, and misclassifications were primarily due to statistical sampling error in the training datasets. The classifiers were designed to predict concentrations to be in one of six classes for nitrate, and one of seven classes for arsenic. Each classification scheme allowed for identification of areas with concentrations that were equal to or exceeding the U.S. Environmental Protection Agency drinking-water standard. Whereas 2.4 percent of the area underlain by basin-fill aquifers in the study area was predicted to equal or exceed this standard for nitrate (10 milligrams per liter as N; mg/L), 42.7 percent was predicted to equal or exceed the standard for arsenic (10 micrograms per liter; μg/L). Areas predicted to equal or exceed the drinking-water standard for nitrate include basins in central Arizona near Phoenix; the San Joaquin, Inland, and San Jacinto basins of California; and the San Luis Valley of Colorado. Much of the area predicted to equal or exceed the drinking-water standard for arsenic is within a belt of basins along the western portion of the Basin and Range Physiographic Province in Nevada, California, and Arizona. Predicted nitrate and arsenic concentrations are substantially lower than the drinking-water standards in much of the study area—about 93.0 percent of the area underlain by basin-fill aquifers was less than one-half the standard for nitrate (5.0 mg/L), and 50.2 percent was less than one-half the standard for arsenic (5.0 μg/L).

Arizona, California, Colorado, Nevada, New Mexico,↗

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&mdash;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&trade; 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&mdash;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&ndash;April streamflow of 5,960 acre-feet discharged from Birch Spring Draw at site BSD-2 to Flaming Gorge Reservoir during 2007&ndash;12, containing an average dissolved-solids load of 14,660 tons. An average May&ndash;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&ndash;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↗