Search USGS⌕ Search

SEARCH · Search USGS

Results for “Hydrology”

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

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

At least 1,063 records · Page 59Linked to original sources

Hydrologic study at Farm Creek Marsh, Dorchester County, Maryland, from April 2015 to April 2016

In 2015, the U.S. Geological Survey began a 1-year hydrologic study to investigate the extent and cause of inundation at Farm Creek Marsh, in Dorchester County, Maryland. In combination with a tide and precipitation gage, a representative section of the marsh was instrumented with surface-water monitors and shallow groundwater piezometers to capture the spatial and temporal extent of inundation. In addition, water-quality data (major ions and nutrients) were collected to help discern the cause of inundation. Results indicate that during the year-long study, all sites were periodically inundated, ranging from a total of 108 days to the entire study period of 353 days. The depth of inundation was typically between 0 and 0.2 feet (ft) (above land surface), with the exception of large storm events. Less than 0.5 ft of elevation was the difference between a site being inundated during the entire study period of 353 days and a site being inundated for 36 consecutive days out of 108 total days of inundation during the study period. Water-quality data showed a large difference in pH between marsh surface water (6.1 to 6.9 standard pH units) and shallow groundwater (3.0 to 3.6 standard pH units), with differences also observed in concentrations of silica, iron, manganese, and potassium. Collectively, the combination of water-quality, hydrologic, and soils data indicate that inundation is caused by tide and storm events rather than groundwater discharge.

Maryland↗

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

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

Illinois↗

Hydrologic and hydraulic analyses of selected streams in Stark County, Ohio

To update and expand a part of the Federal Emergency Management Agency Flood Insurance Study, the U.S. Geological Survey, the Muskingum Watershed Conservancy District, and the Stark County Commissioners began a cooperative study. The study consisted of hydrologic and hydraulic analyses for selected reaches of 14 streams in Stark County, Ohio: Broad-Monter Creek, Chatham Ditch, East Branch Nimishillen Creek, Fairhope Ditch, Firestone Ditch, Hayden Ditch, Middle Branch Nimishillen Creek, Middle Branch Nimishillen Creek Tributary Number 1, Nimishillen Creek, Reemsnyder Ditch, Sherrick Run, unnamed stream, West Branch Nimishillen Creek, and Zimber Ditch. The study totaled nearly 50 miles of stream reaches. Instantaneous peak streamflows for floods with 10-, 4-, 2-, 1-, and 0.2-percent and 1-percent plus annual exceedance probabilities were estimated using historical streamflow data from the streamgages Nimishillen Creek at North Industry, Ohio (U.S. Geological Survey station number 03118500), and Middle Branch Nimishillen Creek at Canton, Ohio (U.S. Geological Survey station number 03118000), regional flood regression equations, and streamflow urbanization techniques. The annual exceedance probability streamflows were then used in a Hydrologic Engineering Center-River Analysis System step-backwater model to determine water-surface profiles, flood-inundation boundaries for the 10-, 4-, 2-, 1-, and 0.2-percent and 1-percent plus annual exceedance probability floods, and a regulatory floodway along a selected reach of each stream. Model input included DEM-derived cross sections supplemented with field surveys of open channel cross sections and hydraulic structures, field estimates of roughness values, and annual exceedance probability flood estimates from regional regression equations and historical streamflow data. Flood-inundation boundaries were mapped for the 1- and 0.2-percent annual exceedance probability floods and a regulatory floodway for each stream reach.

Ohio↗

Groundwater flow conceptualization of the Pahute Mesa–Oasis Valley Groundwater Basin, Nevada—A synthesis of geologic, hydrologic, hydraulic-property, and tritium data

This report provides a groundwater-flow conceptualization that integrates geologic, hydrologic, hydraulic-property, and radionuclide data in the Pahute Mesa–Oasis Valley (PMOV) groundwater basin, southern Nevada. Groundwater flow in the PMOV basin is of interest because 82 underground nuclear tests were detonated, most near or below the water table. A potentiometric map and nine sets of hydrostratigraphic and hydrologic cross sections supplement the conceptualization. Potentiometric contours indicate that groundwater in the PMOV basin generally flows south-southwest and discharges at Oasis Valley. Groundwater encounters an alternating sequence of low- and high-transmissivity rocks, referred to as dams and pools, respectively, as it moves from east to west across eastern Pahute Mesa. Flow from all Pahute Mesa nuclear tests is to Oasis Valley and is well-constrained by water-level data. Flow converges along a corridor of high transmissivity between Pahute Mesa and Oasis Valley. The location of the lateral PMOV basin boundary is well defined, and this boundary, with a few minor exceptions, represents a no-flow boundary. Some boundary uncertainty exists in the northeastern part of the basin, but potential flow-rate estimates across the northeastern boundary resulting from this uncertainty are small relative to the basin groundwater budget. Recharge in the PMOV basin is derived from episodic pulses of modern water and the diffuse percolation of old water (greater than 1,000 years). Episodic recharge is a minor recharge component observed as a rise in groundwater levels that occurs 3 months to 1 year following a wet winter. Minor amounts of episodic recharge through an unsaturated zone in excess of 1,000 feet (ft) requires preferential flow through faults and fractures. The dominant recharge component is slow, steady, diffuse percolation of old water through the unsaturated zone. A large component of old water recharging the groundwater system is consistent with observations of isotopically light deuterium and oxygen 18 compositions in water from wells on Pahute Mesa and central Oasis Valley. About half the recharge in the PMOV basin is derived from the eastern Pahute Mesa area. The remaining recharge is derived primarily from other highland areas including Timber Mountain, Belted and Kawich Ranges, and Black Mountain. The PMOV groundwater system is nearly steady state, where recharge is balanced by the 5,900 acre-feet per year of natural discharge at Oasis Valley. This assumption is reasonable because the basin is dominated by steady-state conditions, where long-term changes in groundwater storage are minimal. Total groundwater withdrawals from 1963 to 2018 have amounted to less than 10 percent of annual groundwater discharge and less than 0.2 percent of the basin’s groundwater storage. Therefore, present-day (2020) conditions are considered representative of predevelopment (pre-1950) conditions in nearly all areas of the basin. The lower PMOV basin boundary is defined at 4,000 ft below the water table to encompass all underground nuclear tests and tritium plumes. This boundary defines the lower boundary of radionuclide migration. However, nearly all flow and tritium transport occur in the upper 1,600 ft of the saturated zone because a transmissivity-with-depth relation indicates that greater than 90 percent of the transmissivity contributing to groundwater flow occurs within 1,600 ft of the water table. Rocks at deeper depths have low transmissivity because argillic and mineralized alterations plug the fractures. Volcanic rocks form the primary aquifers and confining units in the PMOV basin. Volcanic hydrogeologic units (HGUs) and hydrostratigraphic units (HSUs) have transmissivity distributions that span up to eight orders of magnitude with considerable overlap between distributions. Despite the large overlap between units, mean transmissivities of aquifers are one-to-two orders of magnitude greater than the confining units. However, all volcanic-rock HGUs and HSUs are composite units, meaning that they can function spatially as either an aquifer or confining unit .

Nevada↗

Hydrologic and hydraulic analyses of the Grand River, Red Cedar River, and Sycamore Creek near Lansing, Michigan

The U.S. Geological Survey (USGS) completed hydrologic and hydraulic analyses for selected reaches of the Grand River, Red Cedar River, and Sycamore Creek near Lansing, Michigan, in cooperation with the city of Lansing. The study comprised a 3.1-mile reach of the Grand River, a 30.3-mile reach of the Red Cedar River, and a 12.0-mile reach of Sycamore Creek. The information produced from the study can be used to update and expand an existing Federal Emergency Management Agency Flood Insurance Study for Ingham County, Mich. Historical streamflow data from USGS streamgages on Grand River at Lansing, Mich. (station number 04113000); Red Cedar River at East Lansing, Mich. (station number 04112500); Red Cedar River near Williamston, Mich. (station number 04111379); and Sycamore Creek at Holt Road near Holt, Mich. (station number 04112850) were used to estimate instantaneous peak streamflows for floods with 10-, 4-, 2-, 1-, and 0.2-percent annual exceedance probabilities (AEPs) and a “1-percent plus” AEP. The Hydrologic Engineering Center’s River Analysis System step-backwater model was used to determine water-surface elevation profiles for the 10-, 4-, 2-, 1-, and 0.2-percent AEP floods, the 1-percent plus AEP flood, and a regulatory floodway for each stream reach. The hydraulic models were calibrated based on stage-streamflow ratings at USGS streamgages. Flood-inundation boundaries for the 1- and 0.2-percent annual exceedance probability floods and regulatory floodway were created for each stream.

Michigan↗

Hydrologic and hydraulic analyses of selected streams near the city of Rittman in Wayne and Medina Counties, Ohio

The U.S. Geological Survey, in cooperation with the Muskingum Watershed Conservancy District and the city of Rittman, Ohio, did a study to provide data to update and expand parts of two Federal Emergency Management Agency Flood Insurance Studies. The study consisted of hydrologic and hydraulic analyses for selected reaches of four streams (Chippewa Creek, Little Chippewa Creek, Styx River, and the unnamed tributary to Styx River) near the city of Rittman in Wayne and Medina Counties, Ohio. The study covered 36.2 miles of stream reaches. Instantaneous peak streamflows for floods with 10-, 4-, 2-, 1-, and 0.2-percent and 1-percent plus annual exceedance probabilities were estimated using historical streamflow data from three U.S. Geological Survey streamgages and regional flood-frequency regression equations. The flood-frequency estimates were then used in a Hydrologic Engineering Center River Analysis System step-backwater model to determine water-surface profiles; flood-inundation boundaries for the 10-, 4-, 2-, 1-, and 0.2-percent and 1-percent plus annual exceedance probabilities; and a regulatory floodway for the study reaches. Model inputs included cross sections derived from a digital elevation model supplemented with field surveys of open-channel cross sections and hydraulic structures, field estimates of Manning’s roughness values, and flood estimates determined from regional regression equations and historical streamflow data. Flood-inundation boundaries were mapped for each stream reach for the 1- and 0.2-percent annual exceedance probability floods and a regulatory floodway. All data used in the creation of the flood-inundation boundaries are available through a U.S. Geological Survey data release (Ostheimer, 2021) and will be submitted to the Federal Emergency Management Agency for inclusion in updated Flood Insurance Studies for Wayne and Medina Counties.

Ohio↗

Hydrologic and water-quality conditions in the Cedar River alluvial aquifer, Linn County, Iowa, 1990–2019

Alluvial aquifers in Iowa have more wells with nitrate exceeding drinking-water standards than other aquifers; are susceptible to contamination by organic contaminants; and have high concentrations of naturally occurring iron and manganese in depositional areas that contain abundant organic matter. The U.S. Geological Survey, in cooperation with the City of Cedar Rapids, Iowa, studied the Cedar River alluvial aquifer in Linn County, Iowa, from 1990 to 2019 to understand the effect of municipal pumping on spatial and temporal hydrologic and water-quality variability. The Cedar River alluvial aquifer is the source of water for the city of Cedar Rapids, Iowa. Withdrawal of large quantities of water for municipal and industrial supply has altered the normal flow of water in the alluvial aquifer. Pumping induces flow from the Cedar River and the underlying bedrock aquifer into the alluvial aquifer. Water quality in the alluvial aquifer varies along the Cedar River. Changes in nitrate, ammonia, manganese, and iron in the alluvial aquifer are seen as the upstream free-flowing reach of the Cedar River transitions to a partially regulated downstream reach, likely because of differences in reduction-oxidation conditions in the aquifer, which are controlled by infiltration from the Cedar River under normal conditions and when wells are being pumped. Nitrate, normally found in oxygenated environments, had the highest concentrations in the most upstream wells in the Seminole well field and the lowest concentrations in the most downstream wells in the East well field. In contrast, ammonia, manganese, and iron, normally found in greatest abundance in anoxic (reducing) conditions, had the greatest concentrations in the most downstream wells. Additionally, dissolved nitrate plus nitrite nitrogen concentrations in wells were substantially less and manganese concentrations were greater in production wells near backwater wetlands in contrast to wells near the Cedar River. Temporal variability in water quality in the alluvial aquifer was driven by pumping that increased flow from the Cedar River into the alluvial aquifer and ultimately led to changes in reduction-oxidation conditions of the aquifer. Increasing dissolved nitrate plus nitrite nitrogen concentrations in the Cedar River from 1990 to 2019 were mirrored in the alluvial aquifer. Anoxic conditions are prevalent in the alluvial aquifer next to the Cedar River when the aquifer is not under pumping stress. However, production well pumping caused induced infiltration of oxygenated river water into the aquifer resulting in increased dissolved nitrate plus nitrite nitrogen concentrations and pesticides and decreased naturally occurring dissolved iron and manganese. Hydrologic and water-quality conditions in the Cedar River alluvial aquifer from 1990 to 2019 provide baseline conditions needed to evaluate the effects of current and future nutrient reduction efforts and land-use changes in the Cedar River Basin on water quality of the Cedar River alluvial aquifer and its source water, the Cedar River. This summary and analysis provide information that can assist the City of Cedar Rapids Utilities Water Department in managing groundwater resources, and provides information that could be used develop a groundwater-quality model to characterize variability over larger areas of the alluvial aquifer, allowing water providers to plan for future water needs of their users.

Iowa↗

Hydrologic budget of the Harney Basin groundwater system, southeastern Oregon

Groundwater-level declines and limited quantitative knowledge of the groundwater-flow system in the Harney Basin prompted a cooperative study between the U.S. Geological Survey and the Oregon Water Resources Department to evaluate the groundwater-flow system and budget. This report provides a hydrologic budget of the Harney Basin groundwater system that includes separate groundwater budgets for upland and lowland areas to avoid double counting water that recharges in the uplands, discharges to streams and springs in the uplands, flows downstream to the lowlands, and recharges the lowland groundwater system. Lowlands generally represent the conterminous valleys within the center of the basin, including floodplains of the major streams and uplands represent all other areas in the basin. The upland groundwater budget is minimally affected by groundwater development and generally represents the budget of the natural system. In upland areas during 1982–2016, mean-annual recharge totaled 288,000 acre-feet (acre-ft) and mean-annual discharge totaled 239,000 acre-ft, resulting in a net recharge of 49,000 acre-ft. Upland groundwater recharge occurs as infiltration of precipitation and snowmelt and was estimated using the USGS Soil-Water-Balance model calibrated to estimates of runoff, evapotranspiration (ET), base flow, and snow-water equivalent. Groundwater discharge to streams is the predominant discharge mechanism in upland areas and was estimated as 225,000 acre-feet per year (acre-ft/yr) during 1982–2016 using hydrograph separation and summer low-flow estimates in streamgaged watersheds and a linear relation between estimated streamflow and base flow in ungaged watersheds. The remaining upland discharge occurs through springs (14,000 acre-ft/yr) that either emerge downgradient of locations where groundwater discharge to streams was estimated or are routed to irrigated areas. Spring discharge was estimated as a compilation of current and historical measurements. The net upland recharge, which is 17 percent of total upland recharge, ultimately recharges lowland areas as groundwater flow from uplands to lowlands. The lowland groundwater budget for the Harney Basin represents a combination of natural conditions and human activity as more than 99 percent of groundwater development has occurred either inside or within 2 miles of the lowland boundary. In lowland areas during 1982–2016, mean annual groundwater recharge totaled 173,000 acre-ft and groundwater discharge totaled 283,000 acre-ft, indicating discharge exceeded recharge by more than 60 percent. Excluding groundwater pumping, the lowland groundwater budget is more in balance with a mean annual recharge of 165,000 acre-ft and a mean annual discharge of 131,000 acre-ft during 1982–2016. The 23-percent difference between non-pumping recharge and discharge mostly represents the cumulative uncertainty in the estimates of the various groundwater budget components but also likely includes a small reduction in natural groundwater discharge captured by pumping. Lowland groundwater is predominantly recharged by infiltration of surface water (116,000 acre-ft/yr) through streams, floodwater, and irrigation, with a lesser amount as groundwater inflow from uplands and minimal recharge beneath Malheur and Harney Lakes. Recharge from streams and floodwater (natural and irrigation) was estimated using a balance of measured and estimated surface-water inflow to and outflow from lowland areas including streamflow, springflow, and ET where a portion of surface-water inflow to lowland areas is comprised of upland discharge to streams and springs. Groundwater ET (119,000 acre-ft/yr) is the predominant natural discharge mechanism in lowland areas and was estimated as the mean from two remote-sensing based approaches incorporating groundwater ET measurements from other similar basins and 23 years (1987–2015) of Landsat imagery. Discharge of lowland groundwater into Malheur and Harney Lakes is about 700 acre-ft/yr and is represented in groundwater ET estimates. The remaining natural groundwater discharge from lowland areas issues from Sodhouse Spring (8,900 acre-ft/yr) and as groundwater flow to the Malheur River Basin through Virginia Valley (3,100 acre-ft/yr). The relatively large amount of groundwater discharged to springs in Warm Springs Valley (25,000 acre-ft/yr) is accounted for in groundwater ET estimates. Natural groundwater discharge in lowland areas of the Harney Basin has remained relatively constant during the last 80 years based on comparisons with estimates north of Malheur Lake and west of Harney Lake published in the 1930s. Annual net amount of groundwater pumped (pumpage) from the Harney Basin during 2017–18 averaged 144,000 acre-ft. The net value is the difference between pumpage (about 152,000 acre-ft/yr) and reinfiltration of groundwater pumped for irrigation and non-irrigation purposes (about 8,000 acre-ft/yr). Net pumpage was estimated in concurrent studies that compiled groundwater-use data and coupled reported groundwater pumpage data from wells with remote-sensing-based ET estimates from groundwater-irrigated fields. Total pumpage for irrigation has increased from about 54,000 acre-ft/yr during 1991–92 to 145,000 acre-ft/yr during 2017–18. Presently, pumpage is greatest in the lowland region north of Malheur Lake (81,000 acre-ft/yr), with lesser amounts to the north and northwest of Harney Lake (41,000 acre-ft/yr) and to the south and east of Malheur Lake (22,000 acre-ft/yr). During this study, mean annual lowland groundwater discharge (including pumpage) exceeded mean annual recharge, indicating that the lowland hydrologic budget is out of balance. Net groundwater pumpage during 2017–18 is similar to groundwater discharge from all other sources in the lowlands and is four times the imbalance between non-pumping lowland recharge and discharge (34,000 acre-ft/yr). Declining groundwater levels at depth across many parts of the Harney Basin lowlands indicate that pumpage is depleting aquifer storage and is likely capturing a small amount of natural groundwater discharge to springs and ET in some lowland areas. If pumping continues, aquifer storage depletion will continue until the capture rate of natural discharge to springs and ET is equal to the pumping rate. If groundwater development occurs in upland areas and reduces either the streamflow or groundwater inflow to lowland areas, the deficit in the lowland water budget will increase.

Oregon↗

Groundwater hydrology in the area of Savannah and Gunstocker Creeks in northeastern Hamilton, southern Meigs, and northwestern Bradley Counties, Tennessee, 2007–09

The U.S. Geological Survey, in cooperation with the Savannah Valley Utility District, evaluated the groundwater hydrology of the Valley and Ridge carbonate rock aquifer in northeastern Hamilton, southern Meigs, and northwestern Bradley Counties, Tennessee, from 2007 through 2009. The evaluation included, and built on, the results of test drilling conducted in the area in 1974 to determine the potential for groundwater as a source of public supply for the utility and the results of an investigation conducted to define recharge areas for wells used by groundwater-source public-supply water systems throughout Hamilton County in the early 1990s. Groundwater-level data collected from wells open to the aquifer in the study area were used to prepare potentiometric-surface maps for fall 1992, spring and fall 1993, summer 2008, and spring 2009 conditions. Two primary groundwater basins were delineated from the maps—the larger of which coincides with the watershed of Savannah Creek in the southern part of the study area and the smaller of which coincides with the watershed of Gunstocker Creek in the northern part of the study area. Both basins are characterized by potentiometric surfaces that contain a central area of low-altitude groundwater levels and low gradients relative to the basin margins that reflect the orientation of enhanced permeability along dissolution-enlarged features that have developed parallel to strike in the aquifer. The recharge area of the Savannah Creek groundwater basin is estimated to be about 31 square miles, and the recharge area of the Gunstocker Creek groundwater basin is estimated to be about 17 square miles. Recharge to the aquifer in the Savannah Creek and Gunstocker Creek groundwater basins primarily occurs in the uplands area along White Oak Mountain in the eastern part of the study area and along the western boundaries of the basins. Groundwater flows toward the potentiometric lows in each basin, discharging as base flow to the streams and to springs locally. Groundwater withdrawals for public supply by the utility influence the potentiometric low in the north-central part of the Savannah Creek groundwater basin and disrupt flow in the creek and nearby Anderson Spring, particularly during the summer and fall seasons. No large groundwater withdrawals currently occur in the Gunstocker Creek basin, but there is potential for groundwater supply development in the basin. A conceptual model of the groundwater hydrology of the area developed from the evaluation indicates that Chickamauga Lake is the base-level control on groundwater discharge from the Savannah Creek and Gunstocker Creek basins and that lake stage affects the potentiometric surfaces and groundwater discharge in the most downgradient parts of the basins as a result of inferred hydraulic connection between the aquifer and the lake. The model also infers that captured surface water from sections of Savannah Creek and the Hiwassee River that are embayed by the lake could recharge the aquifer and serve as a source of water withdrawn by wells in each basin if the potentiometric surfaces were lowered to altitudes less than the stage of the lake, particularly under potential future groundwater-development scenarios in the Gunstocker Creek basin. Geochemical analysis of samples collected from six wells for the study indicate that groundwater in the Valley and Ridge aquifer in the area generally is a calcium-magnesium-bicarbonate type, and although the water generally is hard, it is suitable for most uses. Trace-element concentrations were less than primary drinking-water criteria in all the samples. Results of the investigation indicate that options are available for additional groundwater withdrawal in the study area. Water-level data collected since 1975 at the Savannah Valley Utility District Smith Road well site indicate that some additional amount of groundwater is available for withdrawal from the aquifer in the Savannah Creek groundwater basin. The potentiometric low within the Gunstocker Creek groundwater basin indicates that an area with enhanced permeability is present as a northeastern counterpart to the potentiometric low within the Savannah Creek basin. Because the Gunstocker Creek basin is about one-half the total area of the Savannah Creek basin, a commensurate decrease in available groundwater storage is likely. Furthermore, groundwater withdrawal locations in the Gunstocker Creek basin would be closer to—and possibly connected hydraulically to—the Hiwassee River, thus increasing the potential for induced surface-water recharge in the basin if sustained drawdown from pumping lowered groundwater levels to altitudes less than the stage of the river.

Tennessee↗

Hydrologic and geochemical characterization of the Petaluma River watershed, Sonoma County, California

Executive Summary The objectives of the study are to (1) develop an updated assessment of the hydrogeology and geochemistry of the Petaluma valley watershed (PVW) and (2) develop an integrated hydrologic model for the PVW. The purpose of this report is to describe the conceptual model of the hydrologic, hydrogeologic, and water-quality characteristics of the PVW and a numerical groundwater-flow model of PVW.

California↗

Update and recalibration of the Rio Grande Transboundary Integrated Hydrologic Model, New Mexico and Texas, United States, and northern Chihuahua, Mexico

The Rio Grande Transboundary Integrated Hydrologic Model (RGTIHM) was developed through an interagency effort between the U.S. Geological Survey and the Bureau of Reclamation to provide a tool for analyzing the hydrologic system response to the historical evolution of water use and potential changes in water supplies and demands in the Hatch Valley (also known as Rincon Valley in the study area) and Mesilla Basin, New Mexico and Texas, United States, and northern Chihuahua, Mexico. Reclamation operates the Rio Grande Project (RGP) to store and deliver surface water for irrigation and municipal use within the study area and in the El Paso Valley south of the El Paso Narrows. Biases in the RGTIHM’s simulation of streamflow and aquifer storage depletion and the availability of new estimates of historical agricultural consumptive use in the study area initiated an update and recalibration of the RGTIHM. In addition to the new estimates of historical agricultural consumptive use, updates were made to more accurately represent the natural system and included adjustments to the initial groundwater levels; streamflow rating tables; Rio Grande, canal, and drain streambed elevations; tributary streambed elevations; surface-water inflows and diversions; RGP surface-water deliveries and canal waste; on-farm efficiency; the routing of surface-water runoff within the MODFLOW Farm Process; and general head boundaries used to simulate interbasin groundwater flow. Model settings, including the assignment of hydraulic conductivity and storage properties to model layers and the MODFLOW solver package, were adjusted to improve numerical stability, and the model was recalibrated to better simulate the natural system. The updated and recalibrated RGTIHM demonstrates a robust ability to simulate the spatially and temporally variable measurements, estimates, or reports of hydraulic head, surface-water flows, agricultural pumping, RGP surface-water deliveries and canal waste, and decadal aquifer storage changes, with improvements over the previous version of the model.

Chihuahua, New Mexico, Texas↗

Evaluation of potential stresses and hydrologic conditions driving water-level fluctuations in well ER-5-3-2, Frenchman Flat, southern Nevada

Well ER-5-3-2 is part of a well network designed to monitor long-term water levels and radionuclide concentrations downgradient from underground nuclear tests that occurred in Frenchman Flat, an area of the U.S. Department of Energy Nevada National Security Site in southern Nevada. Interpretation of monitoring records for well ER-5-3-2 was confounded by previously unexplained water-level fluctuations in the well hydrograph. This study integrated geologic, hydrologic, and water-chemistry data to evaluate potential stresses and hydrologic conditions that likely affected the well ER-5-3-2 hydrograph. Numerical groundwater models were applied to evaluate four model scenarios: (1) wellbore leakage without recharge, (2) wellbore leakage with recharge, (3) equilibration to vertical heterogeneities between shallow (low transmissivity) and deep (higher transmissivity) carbonate zones, and (4) equilibration to lateral heterogeneities in carbonate rocks. Meteoric recharge was not the cause of the 21-foot (ft) water-level rise in well ER-5-3-2 from 2001 to 2011 or the 4-ft decline from 2012 to 2016. Based on observed water-level fluctuations in nearby wells, the water-level rise and decline from recharge for these periods was less than 3 and 1 ft, respectively. The lateral-heterogeneity scenario is based on the assumption that the 21-ft water-level rise from 2001 to 2011 was a natural water-level reequilibration following the pumping-induced depressurization of a large volume of high transmissivity and low-storage carbonate rock that is surrounded by low transmissivity and high-storage carbonate rock. The lateral-heterogeneity scenario was discounted because simulated water levels cannot match the well ER-5-3-2 hydrograph. Underground nuclear testing and temperature effects were discounted based on hydraulic connections and water-temperature data. Wellbore-leakage scenarios are based on the assumption that the water-level rise was sustained from leakage rates required to cause a localized mounding in the carbonate system near well ER-5-3-2, where the carbonate transmissivity is 530 square feet per day. Even though simulated and measured water levels compare favorably for scenarios of wellbore leakage with and without recharge, large volumes (178–184 million gallons) of groundwater from volcanic rocks would be required to leak into the carbonate system, which is not supported by water-chemistry data. An alternative conceptualization of wellbore leakage is based on the assumption that the 21-ft water-level rise from 2001 to 2011 was sustained by the hydraulic disconnection of well ER-5-3-2 from the carbonate system. The disconnection occurred several months after a constant-rate test in well ER-5-3-2 when carbonate rocks were hydraulically disconnected from the well by either (1) the shifting of sloughed fill in the open hole or (2) the encrusting of carbonate precipitate in the well screen. The hydraulic disconnection effectively sealed the well and caused a 21-ft water-level rise from wellbore leakage during 2001–11. In this case, total wellbore leakage from 2001 to 2011 was about 50 gallons. The 4-ft water-level decline from 2012 to 2016 was conceptualized to have occurred from the slow breaking of the seal and reconnection of the well to the carbonate system. This alternative conceptualization of wellbore leakage was consistent with water-chemistry analyses because the computed wellbore leakage (50 gallons) was small relative to purged volumes (30,000–40,000 gallons) for sampling, and the water chemistry would not be expected to change. The shallow-deep carbonate scenario provided another explanation for the well ER-5-3-2 hydrograph. This scenario is based on the assumption that well-construction effects and vertical heterogeneity of the carbonate system explain the ER-5-3-2 water-level trend. Well-construction effects are attributed to a temporary clogging of the open interval below the well screen that was opened during pumping events, which affected the hydraulic connection of deep transmissive carbonate rocks to the wellbore. The 21-ft water-level rise from 2001 to 2011 was a natural equilibration to shallow, low-transmissivity carbonate rocks during a period when the lower open interval was clogged. The 4-ft decline from 2012 to 2016 represents equilibration between the shallow and deep intervals, because of a partial unclogging of the connection between the two intervals. The low water levels from 2016 to 2021 resulted from pumping for sampling and an unclogging of the open interval so that the low head in the deep carbonate dominated the water level. Despite potential well-construction effects, from either a wellbore leakage or shallow-deep carbonate scenario, samples collected from well ER-5-3-2 are representative of the carbonate system.

Nevada↗

Groundwater hydrology, groundwater and surface-water interactions, aquifer testing, and groundwater-flow simulations for the Fountain Creek alluvial aquifer, near Colorado Springs, Colorado, 2018–20

From 2018 through 2020, the U.S. Geological Survey, in cooperation with the Air Force Civil Engineering Center, conducted an integrated study of the Fountain Creek alluvial aquifer located near Colorado Springs, Colorado. The objective of the study was to characterize hydrologic conditions for the alluvial aquifer pertinent to the potential for transport of solutes. Specific goals of this report were to characterize the groundwater hydrology of the area, to quantify groundwater and surface-water interactions, to estimate hydraulic properties of the aquifer using aquifer testing, and to complete numerical simulations of groundwater flow. Synoptic groundwater-level elevation measurements completed throughout this study, and as part of other U.S. Geological Survey programs between 1994 and 2020, indicate groundwater-level elevations fluctuate on annual and interannual timeframes. Groundwater-level fluctuations likely were caused by temporally variable groundwater recharge and discharge components in the area, with many wells showing maximum groundwater-level elevations during the winter months (November through March). From an interannual perspective, groundwater-level fluctuations appear to have reached maximum values during 2000 to 2003, decreased during 2003 to 2006, and remained relatively constant since that time, with the exception of several wells which have displayed rising groundwater-level elevations since 2018. Spatial evaluation of groundwater-level elevations indicates groundwater flow is generally from northeast to southwest within the vicinity of several alluvial paleochannels occurring along the northeastern margin of the aquifer. Within the center of the aquifer along Fountain Creek, groundwater flow is generally from north to south, approximately paralleling surface-water flow. To quantitatively understand the potential effect of groundwater recharge and groundwater pumping on fluctuations in groundwater-level elevation, a statistical transfer-function-noise model was applied. Results of the statistical model indicate throughout most of the aquifer, fluctuations were primarily the result of recharge seasonality. In the main stem of the aquifer where groundwater pumping wells were more concentrated, however, groundwater-level elevation fluctuations were also attributable to groundwater pumping through time. Three-dimensional evaluation of the aquifer geometry near Fountain Creek was combined with synoptic streamflow measurement and accounting of stream gains and losses to evaluate groundwater and surface-water interactions in the study area. Streamflow gain or loss calculations indicate Fountain Creek both gains from and loses flow to the alluvial aquifer, and gaining or losing reaches of the stream may be partially controlled by the depth to bedrock near the stream. Reaches with streamflow gains tend to coincide with areas where the estimated depth to bedrock is decreasing, meaning the alluvial aquifer is likely thinning in these areas and groundwater-flow paths may be converging and discharging groundwater to the stream. Losing reaches tended to coincide with locally greater depth to bedrock where the alluvial aquifer is likely thicker and has greater storage potential for surface water lost from Fountain Creek. Results of aquifer testing indicate hydraulic conductivity, estimated from slug tests and single-well pumping tests, ranged from 0.32 to 1,410 feet per day (ft/d) and 4.13 to 664 ft/d, respectively. These results are similar to the range of values from previous aquifer tests in the study area. Hydraulic conductivities from aquifer testing for this study were generally greater than the estimates of previous slug tests and had a mean value less than the estimates from previous pumping tests. Spatial evaluation of aquifer testing results indicates hydraulic conductivity tends to be greater in the main stem of the alluvial aquifer and lower in paleochannels upgradient from the main stem of the aquifer. The spatial variation in hydraulic conductivity may be attributed to the geomorphologic processes that formed the alluvial aquifer. Compacted sediment in the paleochannels has not been potentially transported sufficient distance to cause grain-size sorting, resulting in a poorly sorted deposit and lower hydraulic conductivities. In the central portion of the alluvial aquifer, near Fountain Creek, the sediments have been transported farther from their source areas and are likely better sorted, removing finer grained sediments that would cause lower hydraulic conductivity. A numerical groundwater-flow model was calibrated for the Fountain Creek alluvial aquifer for 2000–19 to simulate water-budget components, groundwater-flow directions, and groundwater-flow paths. The model simulated precipitation recharge, groundwater and surface-water interactions, evapotranspiration, high-volume groundwater pumping by pumping wells, and external inflows and outflows occurring along the boundaries of the alluvial aquifer. Model calibration was completed using manual and automated approaches, the latter of which assisted in quantifying model results sensitivity to input parameters. The calibrated model corresponds well with groundwater-level elevation observations, with a mean residual (observed minus simulated groundwater-level elevation) equal to −0.60 feet. Simulated groundwater base flow to streams was typically within 10 percent of base flow estimated by independent methods. Groundwater and surface-water interactions represented the largest water-budget components of the aquifer, with the second largest groundwater discharge component coming from pumping wells. Groundwater and surface-water interactions represent both the largest gain and loss terms in the water budget, because these interactions differ spatially, meaning in some areas of the model domain groundwater is being recharged by streams, whereas in other areas, groundwater is discharged to streams. Estimates of advective groundwater-flow paths indicate pumping wells may capture groundwater recharged from losing streams and groundwater that flows into the main stem of the alluvial aquifer from paleochannels.

Colorado↗

Hydrologic analysis of an earthen embankment dam in southern Westchester County, New York

In 2001, the New York City Department of Environmental Protection installed 25 wells on the southern embankment of the Hillview Reservoir in Westchester County in an unsuccessful attempt to locate the source of a large seep (seep A) that began flowing continuously in 1999. In 2005, the U.S. Geological Survey began a cooperative study with the NYCDEP to characterize the hydrology of the local groundwater system and identify potential sources of seep A and other seeps on the embankment. At least two groundwater-flow zones—one shallow and the other deep—overlie the bedrock at the Hillview Reservoir in southern Westchester County, New York. Analyses of slug tests of wells drilled into the southern embankment of the reservoir were used to determine the three-dimensional distribution of hydraulic conductivity of the embankment materials. The wells with the minimum and maximum hydraulic conductivity values are in the deep saturated zone on the southern embankment, where hydraulic conductivity ranges from 0.0012 to 2 feet per day. Hydraulic conductivity ranges from 0.0026 to 1 foot per day in the shallow saturated zone and from 0.021 to 0.27 foot per day in the toe of the embankment. A hydraulic conductivity of 0.016 foot per day was determined for one toe well partially screened in the crystalline-bedrock aquifer. In 2005, the U.S. Geological Survey began a cooperative study with New York City Department of Environmental Protection to characterize the local groundwater-flow system and identify potential sources of seeps on the southern embankment of the Hillview Reservoir in southern Westchester County, New York. Long-term hydrologic data indicated that water levels trended downward in 29 of 41 sites, including the reservoir basin that was monitored during the 12-year study period; data from a National Weather Service precipitation gage at Central Park indicated annual precipitation also trended downward during the same 12-year period. Of the seven wells in which water levels trended upward during the study, two of the wells are on the west side of the southern embankment, proximal to a major water supply conduit, whereas the five remaining wells are screened in the toe. These data indicate an increasing hydrostatic pressure within the deep system and the toe of the dam, which could result in future seeps on the southern embankment near these wells. Results of 11 suspended-sediment samples collected from seeps along the southern embankment at 234.1- and 221.6-feet elevation, and another drainage outflow point between 2007 and 2015 indicate a poor correlation between suspended-sediment concentration and discharge. From the flowing seep at 234.1 feet, suspended-sediment concentrations ranged from 1 milligram per liter at a flow of 2.6 gallons per minute (that is, 1 milligram per 0.26 gallons) during March 2008 to 16 milligrams per liter at 12 gallons per minute during July 2014. At about 12 gallons per minute discharge, suspended-sediment concentration from samples collected at that seep during different sampling events, ranged from 3 to 16 milligrams per liter. From the seep at 221.6 feet elevation, the suspended sediment concentration was 2 milligrams per liter at a discharge of 3.4 gallons per minute and 2 milligrams per liter at a discharge of 1.1 gallons per minute. Only one sample was collected at the drainage outflow point, for which the suspended sediment concentration was 2 milligrams per liter at a discharge of 2.4 gallons per minute. Anomalously high-water levels were recorded in deep-system wells between June 5, 2013, and January 14, 2014. The period for the increase and the decrease back to more typical water-level elevations occurred rapidly during a 13-hour period in each instance. The sudden and rapid changes, in addition to the spatial distribution of magnitude of water-level response indicate that leaky water infrastructure was the source of recharge to the affected wells. A major water supply conduit was drained for repairs between July 7 and 10, 2010. The seeps indicated an immediate response and a substantial hydraulic connection to the water supply conduit. Approximately 10.5 hours after the water supply conduit was drained, flow from a seep on the southern embankment decreased from about 20 gallons per minute to less than 1 gallon per minute. This seep is located at about the same elevation and within the vicinity of the water supply conduit. A travel-time of about 10.5 hours from the source to the seep at 234.1 feet elevation was estimated from the dewatering timeline. During the 3-month shutdown of the water supply conduit, the previously flowing seeps remained dry until precipitation resulted in discharge of about 0.7 gallon per minute at the higher elevation seep, indicating a minor contribution from precipitation to the total seepage discharge. Discharge from the seeps resumed almost immediately coincident with the refilling of the water supply conduit, supporting the hydraulic connection observations during the drainage stage. In addition, during the refilling of the water supply conduit on September 21, 2010, a new seep (I) was observed on the southern embankment. Discharge from this new seep remained relatively constant until it became inaccessible under construction stone from subsequent embankment repairs by the New York City Department of Environmental Protection. Precipitation after the refilling stage of the shutdown seemed to have induced a rise in water levels in the toe wells and an increase in discharge from the seep at 234.1 feet elevation. The post shutdown discharge was less than 12 gallons per minute, compared to a discharge of about 20 gallons per minute before the repairs. The lower discharge rate measured during the period of historically higher discharge rates for the fall season indicates that the repair of the major water supply conduit may have contributed to a reduced discharge from the seeps. There were no definitive responses to the shutdown in any of the wells near the major water supply conduit. The more transmissive deep system of the southern embankment near the major water supply conduit and its associated infrastructure seems to be the preferential flow path for leaking infrastructure. The wells screened in this system showed a response during the deep system anomaly and have some of the highest hydraulic conductivities of the tested wells. All the seeps are in the elevation range of the deep system from approximately the crystalline bedrock surface around 200 feet elevation to the contact between the deep and shallow saturated zones of the reservoir at about 250 feet elevation.

New York↗

A comparison of contemporary and historical hydrology and water quality in the foothills and coastal plain of the Arctic National Wildlife Refuge, Arctic Slope, northern Alaska

The Arctic National Wildlife Refuge is a unique landscape in northern Alaska with limited water resources, substantial biodiversity of rare and threatened species, as well as oil and gas resources. The region has unique hydrology related to perennial springs, and the formation of large aufeis fields—sheets of ice that grow in the river channels where water reaches the surface in the winter and freezes. This work aims to update our understanding of water resources and water quality in the springs, streams, rivers, and lakes of this region, returning to sites sampled by the U.S. Geological Survey in the 1970s. We resampled eight streams, four springs, and six lakes for hydrological metrics, water quality, and macroinvertebrates, and recalculated flood-frequency metrics for rivers using updated data and modern techniques. Aufeis field melt rates were also assessed for the past several decades. Although the available data preclude trend determinations in most cases, our analysis and comparison to the historical sampling indicates an increase in dissolved ions for streams and springs, faster and earlier aufeis melt, and similar macroinvertebrate populations.

Alaska↗

Geomorphic change, hydrology, and hydraulics of Caulks Creek, Wildwood, Missouri

Caulks Creek is a small stream that flows through the city of Wildwood in western St. Louis County, Missouri. The U.S. Geological Survey, in cooperation with the city of Wildwood, has documented historical and recent geomorphic change along Caulks Creek, simulated the hydrologic and hydraulic response of Caulks Creek to a variety of design storm scenarios, and simulated bank retreat resulting from fluvial erosion and mass failure processes. Six study reaches were selected for monitoring geomorphic change based on known locations of erosion issues documented by the city of Wildwood. Recent short-term rates and patterns of geomorphic change in the study reaches, with a focus on bank retreat, were determined from repeat terrestrial light detection and ranging surveys and field observations of the six study reaches. Historical aerial photographs of the study reaches were analyzed to determine long-term rates of bank retreat and channel widening. In general, rapid bank retreat and widening was observed at the outer banks of meander bends and where banks are unforested. Short-term bank retreat varied substantially within individual study reaches, across the study area, and during the study period from no change to as much as 16 feet of retreat between consecutive surveys (5 to 8 months). The field surveys and visual observations indicated that bank retreat occurs episodically owing to a combination of fluvial erosion and mass failure processes, as well as freezing and thawing cycles. Long-term rates of bank retreat ranged from 0.6 to 4.4 feet per year. Hydrologic and hydraulic models of Caulks Creek were used to quantify the peak, volume, and timing of the flow response and the spatial distribution of hydraulic drivers of erosion (velocity and shear stress) along Caulks Creek for design storm scenarios that represent current (as of this publication) and projected future climate. The projected climate conditions resulted in higher peak flows compared to current conditions, including 6 to 21 percent for the year 2050 and 10 to 42 percent for the year 2099 at the downstream end of the study area. Additionally, for a given design storm, projected climate change is predicted to result in faster flows with greater shear stress, as well as more within-stream variability in velocity and shear stress. Many factors affect the velocity and shear stress at a given location, but in general, somewhat fast velocities and high shear stresses tended to occur where the channel is relatively narrow and straight. The velocity and shear stress in the study reaches (known areas of widening and bank retreat) were not particularly high, at least in part owing to the relatively large widths and high sinuosity of the present-day channel in these reaches. The potential mitigating effect of adding runoff storage to the basin also was examined for a selection of design storm scenarios. Additional runoff storage was more effective at mitigating peak flows and total runoff volumes for higher-frequency, lower-intensity storms than for lower-frequency, higher-intensity storms. The additional storage also resulted in an overall reduction in velocity (by as much as 28 percent) and shear stress (by as much as 40 percent) in the study area. However, the effect of the additional storage on peak flows, total runoff volumes, velocity, and shear stress decreased with distance downstream through the study area. For the simulated scenarios, added runoff storage was effective at mitigating the increases in peak flows, total runoff volumes, velocity, and shear stress caused by projected climate change. Lastly, the bank stability and toe erosion model (BSTEM) was used to predict bank erosion and potential bank failure surfaces at five locations along Caulks Creek for a selection of design storm scenarios. The lower-frequency, higher-magnitude design storms resulted in more bank retreat than the higher-frequency, lower-magnitude design storms, though the magnitude of the difference was site dependent. Although scenarios with additional storage were not directly simulated in BSTEM, it is likely that the additional storage would result in reduced bank retreat compared to the same design storm with existing storage, based on the hydraulic modeling results for scenarios with added runoff storage.

Missouri↗

Groundwater hydrology, groundwater and surface-water interactions, water quality, and groundwater-flow simulations for the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado, 2017–19

In 2017, the U.S. Geological Survey, in cooperation with the Upper Arkansas Water Conservancy District, began a study to provide a comprehensive analysis of the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado. The study included collection of data pertaining to groundwater hydrology, groundwater and surface-water interactions, and water quality in the alluvial aquifer. In addition to providing foundational information on the hydrology of the alluvial aquifer, a numerical groundwater-flow model was developed to estimate the potential effects of additional storage of groundwater in the alluvial aquifer. Groundwater-level elevation data from 30 wells were used to estimate groundwater-flow directions in the alluvial aquifer, which were generally from the southwest to northeast, away from the Sangre de Cristo Mountains and towards perennial streams in the center of the valley. Although some seasonal variation was apparent in groundwater-level elevation records, no statistically significant seasonal trends were indicated. Statistically significant long-term trends were indicated in groundwater-level elevation records for 8 of the 30 wells, and of these wells with statistically significant trends, all but 1 indicated a negative trend of groundwater-level elevations. Spatial evaluation of wells with statistically significant negative groundwater-level elevation trends showed many are in areas of denser well drilling for domestic or other uses, indicating increasing groundwater use could potentially be causing groundwater-level elevation declines. There were instances of wells with no statistically significant groundwater-level elevation trends also located in areas of greater density of well completions. Additional investigations may be necessary to more fully characterize the processes responsible for negative groundwater-level elevation trends. Streamflow gain or loss calculations were completed for low flow in 2017–19 and for high flow in 2018 in nine reaches of streams within the study area. Stream reaches of the upper Texas Creek, upper Grape Creek, upper-middle Grape Creek, and Taylor Creek displayed consistent streamflow loss in each period from 2017 to 2019. These stream reaches represent long-term sources of recharge to the alluvial aquifer. Streamflow gain or loss varies through time in other stream reaches (lower Texas Creek, lower-middle Grape Creek, lower Grape Creek below Westcliffe, and lower Grape Creek above DeWeese Reservoir). The temporally variable behavior indicates these stream reaches may be sources of groundwater recharge or areas of groundwater discharge, likely depending on temporal dynamics between the elevation of the water table and the stream. Water-quality samples were collected from 10 groundwater wells and 10 stream sites during September through November 2019. All groundwater and stream samples were analyzed for major and trace elements and stable isotopes of water. A subset of groundwater samples was also analyzed for the environmental tracers sulfur hexafluoride, tritium, and noble gases. Comparison of water-quality results to U.S. Environmental Protection Agency drinking water-quality standards indicated no constituents exceeded primary standards for human health. Spatial evaluation of water quality indicated the concentrations of various constituents are likely controlled by groundwater and surface-water interactions and by spatial variability in bedrock geology underlying the alluvial aquifer. Specifically, streams shown to gain from groundwater had water chemistry constituent compositions similar to groundwater, whereas streams exiting the Sangre de Cristo Mountains tended to have compositions consistent with snowmelt. Groundwater geochemistry appeared to be partially controlled by oxidation-reduction processes and by proximity to igneous rocks in the Wet Mountains. Environmental tracers used to estimate groundwater age indicated all sampled groundwater contained tracers representing modern recharge (approximately less than 65 years old) but mixing of premodern recharge (approximately more than 65 years old) also occurs. Spatial evaluation of environmental tracers indicated large faults may be conduits for upwelling of older groundwater. No trends were observed in groundwater age with well depth, indicating all sampled wells are located within the zone of active groundwater flow. The presence of modern groundwater in wells with statistically significant negative groundwater-level elevation trends indicates groundwater storage depletions may be partially offset by capture of modern recharge. Repeated sampling of groundwater age would be necessary, however, to determine if any trends in groundwater age exist, which may indicate changing groundwater recharge, storage, or discharge. Additional investigations could also consider quantifying groundwater age in deeper wells to more fully define the depth of active groundwater flow. A numerical groundwater-flow model was developed to estimate components of the water budget, simulate groundwater and surface-water interactions, and evaluate the potential effects of aquifer storage and recovery. Simulated groundwater-level elevations from the calibrated groundwater-flow model are similar to the observed pattern of groundwater-level elevations with higher elevations in the western part of the study area along the Sangre de Cristo Mountains. Simulated water-budget components indicate most of the recharge to the alluvial aquifer is derived from streamflow losses, which is consistent with observations of losing streams along the mountain front. The largest groundwater discharge component of the alluvial aquifer was to streams in the center of the valley, where observations of stream gain or loss indicated the predominance of gaining conditions. Comparison of groundwater and surface-water interactions between the calibrated groundwater-flow model for 2000-19 (the base-case model) and a simulation including additional recharge, representing potential aquifer storage and recovery operations, indicated the additional recharge distributed throughout the area had minimal effects on streamflow in the nearby Grape Creek. An analysis of subregional groundwater budgets showed approximately 54 percent of the additional recharge flowed back to nearby Grape Creek, and the other 46 percent was distributed laterally into adjacent cells in the alluvial aquifer. The comparison of simulations and subregional water budget show the additional recharge did not substantially alter groundwater-level elevations or basin wide groundwater storage. Although the analysis of additional recharge provided in the numerical groundwater-flow model considers only one of many possible recharge scenarios, the model provides a useful tool that could be modified for various scenarios to understand potential effects of managed aquifer recharge.

Colorado↗

Potential water-quality and hydrology stressors on freshwater mussels with development of environmental DNA assays for selected mussels and macroinvertebrates in Big Darby Creek Basin, Ohio, 2020–22

The richness and abundance of freshwater mussels in the Big Darby Creek Basin has declined in recent decades, according to survey results published by the Ohio Biological Survey. In October 2016, a major mussel die-off of undetermined cause reportedly affected over 50 miles of Big Darby Creek; however, fishes and other wildlife were not noticeably impacted. Pollution, habitat destruction, climate change, and hydrologic modification have all been theorized as potential reasons for the widespread declines in freshwater mussel populations in North America. To better understand potential stressors to mussels and other aquatic organisms in the Big Darby Creek Basin, the U.S. Geological Survey, in cooperation with the Ohio Water Development Authority, evaluated water quality and temporal changes in hydrology at selected locations. In addition, environmental deoxyribonucleic acid (eDNA) quantitative polymerase chain reaction (qPCR) assays were developed to detect the presence of selected mussels and macroinvertebrates using stream water. Time-weighted average concentrations of pesticides, organic wastewater compounds (OWCs), and polycyclic aromatic hydrocarbons (PAHs) were determined for selected locations within the Big Darby Creek Basin. Passive samplers designed to mimic the respiratory exposure of aquatic organisms and the bioconcentration of organic contaminants into their fatty tissues were deployed three times annually at three sites within the Big Darby Creek Basin in 2020 and 2021. Analyses were done for 204 pesticide compounds, 38 OWCs, and 33 PAHs. Of the 204 pesticide compounds, 70 were detected in at least one sample; 30 were detected in all samples. Herbicides and herbicide degradates were the pesticides most frequently detected and also had some of the highest concentrations of the pesticides detected in this study. Three herbicides (atrazine, ametryn, and metribuzin) were detected in at least 88 percent of samples and two fungicides (azoxystrobin and propiconazole) were detected in all samples. Of the 38 OWCs, 24 were detected in at least one sample; however, only one ( N , N -diethyltoluamide [DEET]) was detected in all samples. Of the 33 PAHs, 29 were detected in at least one sample; 12 were detected in all samples. A continuous water-quality monitor was operated seasonally on Big Darby Creek above Georgesville, Ohio, from 2020 to 2022. Dissolved oxygen concentrations generally followed a daily cycle, peaking in early evening and troughing around sunrise. There were occasional 24-hour swings in dissolved oxygen concentration that had a range exceeding 10 milligrams per liter. However, dissolved oxygen concentrations never fell below Ohio’s aquatic life criteria for warmwater habitats (outside of mixing zones) of 4.0 milligrams per liter as an instantaneous minimum and 5.0 milligrams per liter as a minimum 24-hour average. The Ohio water-quality criteria for temperatures are 29.4 degrees Celsius as an instantaneous maximum and 27.8 degrees Celsius as a 24-hour average maximum. In 2020, there were 10 days when the maximum instantaneous value for temperature was exceeded and 3 consecutive days when the maximum 24-hour average temperature was exceeded. Streamflow time-series data from three gaging stations within the Big Darby Creek Basin were evaluated for trends in annual flow statistics and daily nonexceedance probabilities over time. In general, the evaluation of streamflow conditions at the Big Darby Creek gage (with 97 years of record) indicated that streamflow changed between water years 1922 and 2021. During that time span, flows in general increased, the number of high-flow pulses became more frequent, and low-flow pulses and extreme low-flow periods became less frequent. The only strong indication of trends over time in annual flow statistics for the relatively short records for the other two gages (on Little Darby Creek, with 25 years of record, and Hellbranch Run, with 29 years of record) was that as time went on, reversals between rising and falling periods became more frequent. The U.S. Geological Survey Ohio Water Microbiology Laboratory developed eDNA qPCR assays to detect Epioblasma rangiana (northern riffleshell mussels), Chimarra obscura (a species of caddisfly), Maccaffertium pulchellum (a species of mayfly), and optimized a preexisting eDNA qPCR assay to detect for Ptychobranchus fasciolaris (kidneyshell mussels). The assays were validated by using environmental sampling methods. Assay sensitivity was established by determining the limits of detection and quantification. Water samples were collected at 12 sites in the Big Darby Creek Basin between 2020 and 2022 and analyzed for eDNA with the qPCR assays developed for this study.

Ohio↗