Search USGSSearch

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 397 records · Page 22Linked to original sources

Simulation of hydrologic conditions and suspended-sediment loads in the San Antonio River Basin downstream from San Antonio, Texas, 2000-12

Suspended sediment in rivers and streams can play an important role in ecological health of rivers and estuaries and consequently is an important issue for water-resource managers. To better understand suspended-sediment loads and transport in a watershed, the U.S. Geological Survey (USGS), in cooperation with the San Antonio River Authority, developed a Hydrological Simulation Program—FORTRAN model to simulate hydrologic conditions and suspended-sediment loads during 2000–12 for four watersheds, which comprise the overall study area in the San Antonio River Basin (hereinafter referred to as the “USGS–2014 model”). The study area consists of approximately 2,150 square miles encompassing parts of Bexar, Guadalupe, Wilson, Karnes, DeWitt, Goliad, Victoria, and Refugio Counties. The USGS–2014 model was calibrated for hydrology and suspended sediment for 2006–12. Overall, model-fit statistics and graphic evaluations from the calibration and testing periods provided multiple lines of evidence indicating that the USGS–2014 model simulations of hydrologic and suspended-sediment conditions were mostly “good” to “very good.” Model simulation results indicated that approximately 1,230 tons per day of suspended sediment exited the study area and were delivered to the Guadalupe River during 2006–12, of which approximately 62 percent originated upstream from the study area. Sample data and simulated model results indicate that most of the suspended-sediment load in the study area consisted of silt- and clay-sized particles (less than 0.0625 millimeters). The Cibolo Creek watershed was the largest contributor of suspended sediment from the study area. For the entire study area, open/developed land and cropland exhibited the highest simulated soil erosion rates; however, the largest contributions of sediment (by land-cover type) were pasture and forest/rangeland/shrubland, which together composed approximately 80 percent of the land cover of the study area and generated about 70 percent of the suspended-sediment load from the study area.

Texas

Lithology, hydrologic characteristics, and water quality of the Arkansas River Valley alluvial aquifer in the vicinity of Van Buren, Arkansas

A study to assess the potential of the Arkansas River Valley alluvial aquifer in the vicinity of Van Buren, Arkansas, as a viable source of public-supply water was conducted by the U.S. Geological Survey in cooperation with the Little Rock, District, U.S. Army Corps of Engineers. An important study component was to identify possible changes in hydrologic conditions following installation of James W. Trimble Lock and Dam 13 (December 1969) on the Arkansas River near the study area. Data were gathered for the study in regard to the lithology, hydrologic characteristics, and water quality of the aquifer. Lithologic information was obtained from drillers’ logs of wells drilled from 1957 through 1959. Water-quality samples were collected from 10 irrigation wells and analyzed for inorganic constituents and pesticides. To evaluate the potential viability of the alluvial aquifer in the Van Buren area, these data were compared to similar stratigraphic, lithologic, and groundwater-quality data from the Arkansas River Valley alluvial aquifer at Dardanelle, Ark., where the aquifer provides a proven, productive, sole-source of public-supply water. Drillers’ logs for 59 wells in the Van Buren study area revealed well depths ranging from 25 to 52 feet (ft), with a mean depth of 42 ft. The thickness of the lower sand/gravel interval serving as the water-producing zone ranged from 5 to 47 ft, with a mean thickness of 29 ft. The presence of gravel was noted in only 4 of 59 well logs available for review from the study area. Percent sand was calculated from well logs in the study area, and these sand percentages were overlain onto an orthophotograph map to examine the areal distribution of sand percentage in relation to geomorphologic features of the flood plain in the study area. The logs denoting the greatest percent sand tend to occur in areas near to the river and on the concave (point bar) side of abandoned channels, while the lower percent sand tends to occur on the convex (channel fill and backswamp deposits) side of the abandoned channels. Comparison of hydrographs from water levels collected between 1957 and 1972 to cumulative departure from mean monthly and mean annual precipitation showed overall good fit and explained the long-term decreasing water levels from the earliest period of record through October 1967, followed by a sharp rise in water levels concurrent with rises in cumulative departure from mean monthly and mean annual precipitation. Hydrographs for four wells ranging from 0.8 to 4.5 miles upstream from the dam and potentially affected by rising river stage were compared to graphs of river stage and cumulative departure from mean monthly precipitation. Water levels for these wells showed minimal discernible effect by rising river stage following dam completion. Periods of increased precipitation compared closely to increases in water level for all hydrographs, regardless of river stage, and periods of no precipitation resulted in declining water levels, although river stage continued to slowly rise during these same periods. The Arkansas River has greater salinity than local groundwater, providing a quantitative tracer for any groundwater recharge originating from the river. Comparison of predam and postdam groundwater-chloride concentrations showed no increase in chloride concentrations after dam installation, which is consistent with hydrologic data. These data suggest that the dominant source of groundwater recharge in the Arkansas River Valley alluvial aquifer is infiltration of precipitation through proximal, coarse channel deposits, with minimal influx of river water. Groundwater-quality data collected from 10 wells in the study area indicated a calcium-bicarbonate water type. No primary drinking-water standards were exceeded for any constituents, and iron and manganese were the only constituents exceeding secondary drinking-water regulations. Six of the 10 well-water samples were analyzed for the presence of pesticides, as row-crop agriculture is the dominant land use in the study area. Six herbicide compounds and one herbicide metabolite were detected at concentrations substantially below those of the Federal primary drinking-water standards and health advisories. The hydrologic and geochemical data gathered for this study provide a qualitative assessment of the potential of the Arkansas River Valley alluvial aquifer as a source of public water supply in the Van Buren area. Results indicate minimal influx of water from the Arkansas River, and recharge to the aquifer appears to be dominantly by infiltration of precipitation through overlying alluvium. If vertical wells are used as a source of public water supply, then several wells will have to be used in combination at relatively low pumping rates and placed in areas with a greater percent sand. Use of a horizontal well configuration near the river to increase production may depend on infiltration of river water to supplement water removed from storage, especially where areas of lower permeability sediments might be encountered within the surrounding alluvium. If a poor hydraulic connection exists between the river and the alluvium, as indicated by this study, then production will depend on ample precipitation and recharge throughout the year and groundwater storage sufficient to prevent declining water levels where pumping rates exceed recharge.

Arkansas

Hydrogeology and hydrologic conditions of the Ozark Plateaus aquifer system

The hydrogeology and hydrologic characteristics of the Ozark Plateaus aquifer system were characterized as part of ongoing U.S. Geological Survey efforts to assess groundwater availability across the Nation. The need for such a study in the Ozark Plateaus physiographic province (Ozark Plateaus) is highlighted by increasing demand on groundwater resources by the 5.3 million people of the Ozark Plateaus, water-level declines in some areas, and potential impacts of climate change on groundwater availability. The subject study integrates knowledge gained through local investigation within a regional perspective to develop a regional conceptual model of groundwater flow in the Ozark Plateaus aquifer system (Ozark system), a key phase of groundwater availability assessment. The Ozark system extends across much of southern Missouri and northwestern and north-central Arkansas and smaller areas of southeastern Kansas and northeastern Oklahoma. The region is one of the major karst landscapes in the United States, and karst aquifers are predominant in the Ozark system. Groundwater flow is ultimately controlled by aquifer and confining unit lithologies and stratigraphic relations, geologic structure, karst development, and the character of surficial lithologies and regolith mantle. The regolith mantle is a defining element of Ozark Plateaus karst, affecting recharge, karst development, and vulnerability to surface-derived contaminants. Karst development is more advanced—as evidenced by larger springs, hydraulic characteristics, and higher well yields—in the Salem Plateau and in the northern part of the Springfield Plateau (generally north of the Arkansas-Missouri border) as compared with the southern part of the Springfield Plateau in Arkansas, largely due to thinner, less extensive regolith and purer carbonate lithology. Precipitation is the ultimate source of all water to the Ozark system, and the hydrologic budget for the Ozark system includes inputs from recharge, losing-stream sections, and groundwater inflows and losses of water to gaining-stream sections, groundwater withdrawals, and surface-water and groundwater outflows to neighboring systems. Groundwater recharge, estimated by a soil-water-balance model, represents about 24 percent, or 11 inches, of 43.9 inches annual precipitation. Recharge is spatially variable, being greater in the northern Springfield Plateau and Salem Plateau than in the southern Springfield Plateau (generally south of the Arkansas border) because of differences in regolith mantle extent and thickness and carbonate lithology and hydraulic properties. Increased precipitation and decreased agricultural land use during the period 1951 through 2011 increased recharge by approximately 5 percent. Although all Ozark streams have losing, neutral, and gaining sections, they are dominantly gaining and are a net sink for groundwater with nearly 90 percent of groundwater recharge returned to springs and streams. Groundwater pumping is a small but important loss of water in the Ozark system hydrologic budget; water-level declines and local cones of depression have been observed around pumping centers and strong concerns exist over potential effects on stream and spring flow. Data indicate that societal needs for freshwater resources in the Ozark Plateaus will continue to increase and will do so in the context of changing climate and hydrology. Groundwater will continue to be an important part of supporting these societal needs and also local ecosystems. The unique character and hydrogeologic variability across the Ozark system will control how the system responds to future stress. Groundwater of the Ozark system in the northern study area is more dynamic, has greater storage and larger flux, and has greater potential for further development than in the part of the study area south of the Arkansas-Missouri border. Further south in Arkansas, a line exists, roughly defined as 5 miles south of the Springfield Plateau-Boston Mountains boundary, beyond which further extensive municipal or commercial development appears unlikely under current economic and resource-need conditions. A small part of the Ozark system groundwater budget is currently drafted for use, leaving an apparently large component available for further development and use—particularly in the northern Springfield Plateau and Salem Plateau; however, the effects of increased pumping on groundwater’s role in maintaining ecosystems and ecosystem services are not quantitatively well understood, and the close relation between groundwater and surface water highlights the importance of further quantitative assessment.

Arkansas, Kansas, Missouri, Oklahoma

Variability of hydrological droughts in the conterminous United States, 1951 through 2014

Spatial and temporal variability in the frequency, duration, and severity of hydrological droughts across the conterminous United States (CONUS) was examined using monthly mean streamflow measured at 872 sites from 1951 through 2014. Hydrological drought is identified as starting when streamflow falls below the 20th percentile streamflow value for 3 consecutive months and ending when streamflow remains above the 20th percentile streamflow value for 3 consecutive months. Mean drought frequency for all aggregated ecoregions in CONUS is 16 droughts per 100 years. Mean drought duration is 5 months, and mean drought severity is 39 percent on a scale ranging from 0 percent to 100 percent (with 100% being the most severe). Hydrological drought frequency is highest in the Western Mountains aggregated ecoregion and lowest in the Eastern Highlands, Northeast, and Southeast Plains aggregated ecoregions. Hydrological drought frequencies of 17 or more droughts per 100 years were found for the Central Plains, Southeast Coastal Plains, Western Mountains, and Western Xeric aggregated ecoregions. Drought duration and severity indicate spatial variability among the sites, but unlike drought frequency, do not show coherent spatial patterns. A comparison of an older period (1951–82) with a recent period (1983–2014) indicates few sites have statistically significant changes in drought frequency, drought duration, or drought severity at a 95-percent confidence level.

Scientific Investigations Report

Simulations of hydrologic response in the Apalachicola-Chattahoochee-Flint River Basin, Southeastern United States

A suite of hydrologic models has been developed for the Apalachicola-Chattahoochee-Flint River Basin (ACFB) as part of the National Water Census, a U.S. Geological Survey research program that focuses on developing new water accounting tools and assessing water availability and use at the regional and national scales. Seven hydrologic models were developed using the Precipitation-Runoff Modeling System (PRMS), a deterministic, distributed-parameter, process-based system that simulates the effects of precipitation, temperature, land cover, and water use on basin hydrology. A coarse-resolution PRMS model was developed for the entire ACFB, and six fine-resolution PRMS models were developed for six subbasins of the ACFB. The coarse-resolution model was loosely coupled with a groundwater model to better assess the effects of water use on streamflow in the lower ACFB, a complex geologic setting with karst features. The PRMS coarse-resolution model was used to provide inputs of recharge to the groundwater model, which in turn provide simulations of groundwater flow that were aggregated with PRMS-based simulations of surface runoff and shallow-subsurface flow. Simulations without the effects of water use were developed for each model for at least the calendar years 1982–2012 with longer periods for the Potato Creek subbasin (1942–2012) and the Spring Creek subbasin (1952–2012). Water-use-affected flows were simulated for 2008–12. Water budget simulations showed heterogeneous distributions of precipitation, actual evapotranspiration, recharge, runoff, and storage change across the ACFB. Streamflow volume differences between no-water-use and water-use simulations were largest along the main stem of the Apalachicola and Chattahoochee River Basins, with streamflow percentage differences largest in the upper Chattahoochee and Flint River Basins and Spring Creek in the lower Flint River Basin. Water-use information at a shorter time step and a fully coupled simulation in the lower ACFB may further improve water availability estimates and hydrologic simulations in the basin.

Apalachicola-Chattahoochee-Flint River Basin

Hydrologic conditions and simulation of groundwater and surface water in the Great Dismal Swamp of Virginia and North Carolina

The U.S. Geological Survey (USGS), in cooperation with the U.S Fish and Wildlife Service, has investigated the hydrology of the Great Dismal Swamp (Swamp) National Wildlife Refuge (Refuge) in Virginia and North Carolina and developed a three-dimensional numerical model to simulate groundwater and surface-water hydrology. The model was developed with MODFLOW-NWT, a USGS numerical groundwater flow modeling program, in combination with the Surface-Water Routing Process, a software package that simulates dynamic surface-water flows, water control structure management, and groundwater/surface-water interactions. The steady-state model was calibrated to average spring conditions by using automated parameter estimation software (PEST) to reduce simulation errors and assess model parameter sensitivity. The model was then used to simulate wet and dry climatic conditions and a variety of hypothetical scenarios in which water levels in the Swamp were raised and lowered by simulated management of water control structures. Results of the model simulations indicate that, under average spring conditions, precipitation is the primary water input (92%); surface-water (5%) and groundwater (3%) inflows make up the remainder. The primary outflow (or loss) is evapotranspiration (55%), with surface outflows (about 41%) and groundwater outflow (about 4%) making up the remainder. Simulated adjustment of water control structure weir levels demonstrates that groundwater levels are affected by water levels in adjacent ditches and that surface-water and groundwater levels can be controlled through management of water control structures, allowing the Refuge to better manage fire risks and preserve forested-wetland ecosystems in the Refuge. The 13 water control structures proposed in the simulated scenario representing possible future conditions effectively raised simulated water levels in the northeastern corner of the study area, a goal of the Refuge management. Results of this study demonstrate use of MODFLOW with the Surface-Water Routing Process for simulating water management options in peat wetlands and will help Refuge managers to better understand existing hydrologic conditions, assess the hydrologic effects of planned changes to water control structures, and apply the new simulation tool to guide water management on the Refuge.

North Carolina

Severity and extent of alterations to natural streamflow regimes based on hydrologic metrics in the conterminous United States, 1980–2014

Alteration of the natural streamflow regime by land and water management, such as land-cover change and dams, is associated with aquatic ecosystem degradation. The severity and geographic extent of streamflow alteration at regional and national scales, however, remain largely unquantified. The primary goal of this study is to characterize the severity and extent of alterations to natural streamflow regimes for 1980–2014 based on hydrologic metrics at 3,355 U.S. Geological Survey streamgages in the conterminous United States. Twelve hydrologic metrics with known relevance to aquatic ecosystem health were used to characterize the streamflow regime. Alterations to the 12 hydrologic metrics were quantified by taking ratios of the metrics calculated from observed daily streamflow records divided by the same metrics predicted for natural conditions by random forest statistical models. Some level of streamflow alteration (diminishment or inflation of hydrologic metrics) compared to natural conditions was indicated at about 80 percent of the assessed streamgages across the conterminous United States. The severity of alteration differed among ecoregions because of differences in dominant land and water management practices. Finally, when compared over the period 1980–2014, climate variability generally played a minor role in the alteration of streamflows across the United States when compared to the effects of land and water management.

Conterminous United States

The hydrologic system of the south Florida peninsula—Development and application of the Biscayne and Southern Everglades Coastal Transport (BISECT) model

The Biscayne and Southern Everglades Coastal Transport (BISECT) model was developed by the U.S. Geological Survey under the Greater Everglades Priority Ecosystem Studies Initiative to evaluate, both separately and in conjunction, the likely effects on surface-water stages and flows, hydroperiod, and groundwater levels and salinity in south Florida of (1) a vertical Biscayne aquifer barrier to maintain higher wetland levels, (2) possible future changes to current water-management practices, and (3) sea-level rise. The BISECT model is a combination of the Tides and Inflows to the Mangrove Everglades (TIME) and Biscayne models of the western and eastern parts of south Florida including Everglades National Park, the southern Miami-Dade urban area, and the Biscayne Bay coast and simulates hydrodynamic surface-water flow and three-dimensional groundwater conditions dynamically for the period 1996–2004 by using the Flow and Transport in a Linked Overland/Aquifer Density-Dependent System (FTLOADDS) simulator. BISECT includes a number of parameter and algorithmic refinements that improve simulation results relative to the TIME and Biscayne models and represents the hydrologic system more explicitly, including (1) improved topographic representations, (2) refined Manning’s friction coefficients, (3) improved evapotranspiration computation through spatially variable albedo, (4) increased vertical aquifer discretization, and (5) extension of the western boundary farther offshore. Sensitivity analyses demonstrate that simulated flows into Long Sound have a different pattern of response to tidal amplitude, wind, and frictional resistance changes than do other coastal streams in the model; flows at Broad River and Lostmans River are most sensitive to tidal amplitude, wind, and frictional resistance changes; and flow to the Everglades coastal streams is substantially affected by surface-water/groundwater interactions in the eastern urban areas. Insight into the hydrologic system came from scenario simulations that represent proposed management actions, such as grouting of the aquifer to prevent seepage from the wetlands and changes to water deliveries proposed by the Comprehensive Everglades Restoration Plan (CERP), and projected sea-level rise. These scenario management changes are considered separately to isolate their specific effects and also in conjunction with sea-level rise. Scenario simulations show that (1) attempts to prevent seepage from the wetlands by grouting the aquifer along the L 31N levee produce minimal effects on surface-water levels; (2) the increased water deliveries proposed in the CERP redistribute flow to the northwestern coastal part of the study area with a minimal reduction to the southeast and a more substantial reduction in flows in the intervening coastal zones, mitigating some sea-level rise effects; (3) sea-level rise has a larger effect on the hydrology (water levels, flow, and salinity) than does CERP restoration; and (4) support for ecological models and hydrologic studies can be provided by applying BISECT to scenarios influenced by climatic and anthropogenic changes or by meteorological variability, such as extreme wet or dry periods.

Florida

An enhanced hydrologic stream network based on the NHDPlus medium resolution dataset

The National Hydrography Dataset Plus, Version 2.1 (NHDPlusV2.1) is an attribute-rich digital stream network for the conterminous United States, serving as a foundational infrastructure for reporting hydrologic information at both regional and national scales. SPAtially Referenced Regressions On Watershed attributes (SPARROW) is a process-based statistical model that relies on a digital hydrologic network like NHDPlusV2.1 to establish spatial relations between quantities of monitored contaminant loads and contaminant sources, accounting for the physical characteristics along flow paths affecting contaminant transport. The U.S. Geological Survey National Water Quality Assessment project adopted and modified the medium-resolution NHDPlusV2.1 network for use as the primary framework supporting SPARROW modeling. This report describes the enhancements made to improve the routing capabilities and the value-added attributes of NHDPlusV2.1 to support modeling and other hydrologic analyses. These enhancements include corrections to inconsistencies in network/routing information, filling in missing attribute values of associated characteristics, accounting of water use affecting flow, new variables useful for interpreting network data, revised flowline attributes such as slope and flow, and incorporation of ancillary spatial data into the network. The resulting dataset containing the enhancements to the network is named E2NHDPlusV2_US. Although the enhancements described in the report were developed for use in SPARROW modeling, the enhancements are expected to be useful for a wide variety of hydrologic studies within the United States.

Conterminous United States

Groundwater hydrology and chemistry of Jamestown Island, Virginia—Potential effects of tides, storm surges, and sea-level rise on archaeological, cultural, and ecological resources

As the site of the first permanent English settlement in North America in 1607, Jamestown Island, Colonial National Historical Park (COLO), Virginia, contains a rich archaeological record that extends from the Paleoindian period (15,000 to 8,000 years ago) through the 20th century. The island is located on the lower James River near the mouth of Chesapeake Bay. Jamestown Island vegetation is dominated by upland forests surrounded by tidal, freshwater-to-oligohaline marshes. Along the Virginia coast, relative sea-level rise was more than 2.5 times the global average during the 20th century. Consequently, the National Park Service (NPS) has identified COLO as one of the 25 national parks most threatened by climate change. Surface waters across the island are hydraulically connected to the laterally continuous Surficial aquifer. The land-surface altitude of the island is low, with two-thirds of the island less than 5 feet (ft) above the North American Vertical Datum of 1988 (NAVD 88). Consequently, sea-level rise, combined with tides and storm surges, threatens the island and its resources as surface-water and groundwater levels rise, saltwater enters the Surficial aquifer, and groundwater chemistry changes. The impact of sea-level rise on the island’s surface-water resources has been well studied, but groundwater effects have been largely ignored. Quantifying the effects of tides, storm surges, and sea-level rise on groundwater levels and chemistry is essential to developing an effective strategy for managing climate-induced changes. The first step in developing a response strategy includes a parkwide general risk assessment for archaeological sites on the island, so that sites can be prioritized for management actions. The U.S. Geological Survey and the NPS began a study in 2015 to develop a long-term groundwater-monitoring program to evaluate this risk and to develop an updated management strategy. The groundwater-monitoring program consists of 45 wells and piezometers in two individual clusters and three transects across the island in different hydrologic and chemical settings. Samples for water quality were collected from the wells and piezometers from October 2015 through September 2018 at variable time intervals. Results of the monitoring identified disparate hydrologic and chemical responses to saltwater intrusion across the island. Specific conductance (an indicator of salinity) of groundwater beneath several marshes responded differently to changes in James River salinity. Groundwater response to changes in James River specific conductance appeared to be controlled by land-surface altitude and slope, differences in lateral and vertical sediment characteristics, distance from surface waters, and the degree of surface water/groundwater connectivity between channels and the aquifer. Groundwater chemistry data from monitoring wells at Black Point, a low-altitude, upland setting, are in contrast with conditions observed in Island House observation wells, a high-altitude, upland setting. Specific conductance (less than 200 microsiemens per centimeter [μS/cm]) and pH (greater than 5.0) of groundwater beneath much of the uplands that characterize the Island House observation wells are typical of groundwater in noncarbonate sedimentary aquifers recharged by precipitation. At Black Point, specific conductance ranged from 2,490 to 15,200 μS/cm, and pH ranged from 3.1 to 6.6 standard units. At the Black Point observation wells, the most saline and dense water was at the water table rather than deeper in the aquifer, causing a density inversion that persisted throughout the study. The density inversion likely resulted from differences in permeability between the shallow clay and fine-grained sands and the deeper coarse-grained sand and gravel. Groundwater with the lowest pH was at the water table. As saline groundwater flows through organic sediment beneath the marshes, bacterial biodegradation of organic matter creates anoxic conditions. Continued biodegradation concomitantly reduces iron-oxide minerals in the sediment and sulfate in saline water. When oxygen is reintroduced into groundwater, iron and sulfur can reoxidize to form sulfuric acid, locally lowering the pH of the water. This report describes the groundwater monitoring network design, rationale for site selection, monitoring approach, and results of monitoring from October 2015 through September 2018. Maps of inundation at selected water-level altitudes are included to identify the risk to archaeological, cultural, and ecological resources. The monitoring results of the hydrology and chemistry data are interpreted, and the different hydrologic and chemical settings are described. The implications of the study results for management decisions are presented, and suggestions for improving the monitoring network are included.

Virginia

Hydrologic change in the St. Louis River Basin from iron mining on the Mesabi Iron Range, northeastern Minnesota

This study compares the results of two regional steady-state U.S. Geological Survey Modular Three-Dimensional Finite-Difference Ground-Water Flow (MODFLOW) models constructed to quantify the hydrologic changes in the St. Louis River Basin from iron mining on the Mesabi Iron Range in northeastern Minnesota. The U.S. Geological Survey collaborated in this study with bands of the Minnesota Chippewa Tribe, and the Minnesota Pollution Control Agency to inform management decisions about aquatic resources in the St. Louis River Basin. A model constructed and calibrated to represent average 1995–2015 mining conditions produced regional groundwater heads and flows. A pre-mining scenario model was constructed from this mining model but had the land and bedrock surfaces restored to pre-mining topographies and had modeled mining features (mine pits, tailings basins, waste-rock piles, and mining-disturbed areas) eliminated to represent general pre-mining stratigraphy and hydrogeology. Many of the features important to the hydrology of this mining area (like individual mine pits) are difficult to represent in groundwater models and required the use of modeling tools to indirectly account for their effects. The difference between the results of these two models represents mining’s effects on the hydrology in the Mesabi Iron Range area of the St Louis River Basin. The mining and pre-mining regional models also can provide boundary conditions and initial properties for future local or site-specific groundwater-flow models in the area. Total groundwater flow through the mining model is 171 million cubic feet per day. Areal recharge is the largest source of groundwater (78 and 81 percent of total groundwater flow in the mining and pre-mining scenario models, respectively). Seepage from streams and lakes provides another 17 percent of the total groundwater flow through both models. Water leaves aquifers through seepage to streams (discharge as base flow, 43 percent in both models) and areal seepage to the land surface (surface seepage), for example to wetlands (45 and 49 percent, mining and pre-mining scenario models respectively). Comparison of the results from the mining and pre-mining scenario models shows that iron mining has produced measurable hydrologic changes in the St. Louis River Basin, but that most of those changes and the highest magnitude changes occur near the mining features. Flow changes to and from surface-water bodies like streams and wetlands were analyzed in detail because of their importance in sustaining surface waters and aquatic life. Overall, groundwater flow in the mining model was 3.62 million cubic feet per day (2.2 percent) greater than total pre-mining model groundwater flow. This was caused by an increase in recharge from tailings basins and a decrease in discharge from surface seepage. Groundwater discharge to mine pits was the largest change in groundwater flows between the models (a change representing 2.8 percent of total pre-mining model groundwater flow). Net recharge to groundwater from tailings basins (2.4 percent), net decrease in surface seepage from groundwater (2.7 percent), and net increase in seepage to streams (1.0 percent) were all in this same range of total pre-mining model groundwater flow. Groundwater lost through mine-pit withdrawals was nearly offset by groundwater gained through recharge from tailings basins. However, because losses and gains occurred in different areas, the effect of mining can have more substantial effects on local areas than the model-wide averages represent.

Minnesota

Evaluating drivers of hydrology, water quality, and benthic macroinvertebrates in streams of Fairfax County, Virginia, 2007–18

In 2007, the U.S. Geological Survey partnered with Fairfax County, Virginia, to establish a long-term water-resources monitoring program to evaluate the hydrology, water quality, and ecology of Fairfax County streams and the watershed-scale effects of management practices. Fairfax County uses a variety of management practices, policies, and programs to protect and restore its water resources, but the effects of such strategies are not well understood. This report used streamflow, water-quality, and ecological monitoring data collected from 20 Fairfax County watersheds from 2007 through 2018 to assess the effects of management practices, landscape factors, and climatic conditions on observed nutrient, sediment, salinity, and benthic-macroinvertebrate community responses. Urbanization, climatic variability, and an increase in management practices occurred within Fairfax County during the study period. Impervious cover, housing units, wastewater infrastructure, and (or) stormwater infrastructure increased in most study watersheds. Climatic conditions varied among study years; countywide estimates of average-annual air temperature differed by about 3 degrees Celsius, and total precipitation ranged from about 34 to 63 inches per year. The effects of the management practices, implemented to reduce nitrogen, phosphorus, and (or) sediment loads, are considered in this study. These management practices primarily consist of stormwater retrofits and stream restorations; however, stream restorations account for most of the financial investment and expected load reductions. Management practices were implemented in half of the study watersheds, and most practices were installed and reductions credited late in the study period. Changes in hydrologic response during storm events were evaluated over the study period because many management practices that were implemented were designed to achieve nutrient and sediment reductions by slowing or intercepting runoff. The average number and length of storm events was mostly unchanged throughout the monitoring network. Four watersheds with 10 years of streamflow data showed a mixture of trends in stormflow peak, volume, and rate-of-change. Event-mean nutrient and sediment concentrations from these watersheds were evaluated during storm events and generally showed increases in total phosphorus (TP) and suspended sediment and reductions or no changes in total nitrogen (TN). Landscape inputs of nitrogen and phosphorus and the percentage of inputs delivered to streams were estimated for the study watersheds. Estimated phosphorus from fertilizer and nitrogen from atmospheric deposition represented large nutrient inputs in most watersheds; amounts of other nonpoint sources varied based on land use. Estimated nitrogen inputs declined throughout Fairfax County and in most study watersheds from 2008 through 2018; in comparison, phosphorus input changes were relatively small. Most nonpoint-nutrient inputs were retained on the landscape and did not reach streams, with slightly more nitrogen retention than phosphorus, on average. Retention rates were lower for years with more precipitation and streamflow. After adjusting for streamflow, TN and TP loads were generally higher for years with more nutrient inputs. Calculated as a function of flow-adjusted loads, TP retention declined at most stations from 2009 through 2018, in comparison, TN retention was relatively unchanged. Landscape and climatic conditions affected spatial differences and changes in Fairfax County stream conditions from 2009 through 2018. TN concentrations were higher and increases over time were larger in watersheds with elevated septic-system density. TP concentrations were higher in watersheds with more turfgrass; concentrations were lower, but had larger increases over time, in watersheds with deeper soils. Suspended-sediment concentrations were higher in watersheds with greater stream densities. Specific conductance was higher in watersheds with more developed land use and shallower soils. Benthic-macroinvertebrate index of biotic integrity (IBI) scores were lower in watersheds with high road density and had larger increases over time in bigger, more developed watersheds. Annual variability in TN and TP concentrations and benthic-macroinvertebrate IBI scores was affected by precipitation; annual variability in suspended sediment concentrations and specific conductance was affected by air temperature. After accounting for influences from landscape and climatic conditions, expected management-practice effects were not consistently observed in monitored stream responses. These effects were assessed by comparing expected management-practice load reductions with the timing, direction, and magnitude of changes in storm-event hydrology, nutrient and sediment loads, median-annual water-quality conditions, and benthic-macroinvertebrate IBI scores. An important consideration for future investigations of management-practice effects is how to control for water-quality and ecological variability caused by geologic properties, the urban environment, precipitation, and (or) air temperature. The interpretation of management-practice effects in this report was likely influenced by a combination of factors, including (1) the amount, timing, and location of management-practice implementation; (2) unmeasured landscape and climatic factors; (3) uncertain management-practice expectations; (4) hydrologic variability; and (5) analytical assumptions. Through continued data-collection efforts, particularly after management practices have been completed, many of these factors may become less influential in the future.

Virginia

Hydrology, water-quality, and watershed characteristics in 15 watersheds in Gwinnett County, Georgia, water years 2002–20

The U.S. Geological Survey, in cooperation with Gwinnett County Department of Water Resources, established the Long-Term Trend Monitoring program in 1996 to monitor and analyze the hydrologic and water-quality conditions in Gwinnett County, Georgia. Gwinnett County is a suburban to urban area northeast of the city of Atlanta in north-central Georgia. The monitoring program currently consists of 15 watersheds ranging in size from 1.3 to about 161 square miles. This report synthesizes watershed characteristics and hydrologic and water-quality monitoring data collected for water years (WYs) 2002–20. The 15 study watersheds were characterized for land-surface elevations, average land-surface slopes, septic densities, sanitary sewer densities, and detention pond areas. Temporal patterns in watershed characteristics were determined for land cover (2001–19), percent imperviousness (2000–20), population density (2000–20), and building density (1950–2022). In 2001, most of the watersheds had at least 45 percent of their land cover composed of developed land cover groups, and by 2019, at least 59 percent of each watershed was developed. Land cover changes occurred most rapidly between 2004 and 2008 at most watersheds. Percent imperviousness in the study watersheds varied substantially and ranged from 14.75 to 55.13 percent in 2019. Precipitation and runoff were quantified at all study watersheds for WYs 2002–20, and the hydrologic cycle was evaluated both annually and seasonally. Several 1-year or longer droughts occurred during this period. Study area precipitation averaged 51.5 inches per year and runoff averaged 22.5 inches per year. Variations in annual runoff were largely determined by annual precipitation but were also dependent upon watershed storage. Runoff varied seasonally because of high evapotranspiration rates in the summer and changes in base flow associated with seasonal changes in watershed storage. Fifty-one percent of runoff in the study area occurred as base flow. Watersheds with higher imperviousness had higher stormflows because of increased surface runoff and lower base flows because of reduced infiltration that recharges watershed storage. Turbidity, water temperature, and specific conductance were continuously measured at each study site. These constituents varied seasonally, diurnally, and with streamflow. A minimum of two base-flow and six stormflow samples were collected per year at each watershed and were analyzed for 21 water-quality constituents (water temperature, laboratory specific conductance, pH, and turbidity, biochemical and chemical oxygen demand, suspended sediments, nutrients, base cations, trace metals, and total dissolved solids). Concentrations of most particulate constituents were approximately one-half or more orders of magnitude higher in stormflow samples than in base-flow samples. Total copper and zinc stormflow concentrations exceeded the national recommended aquatic life criteria for acute conditions to varying degrees. Annual loads and yields were estimated for 12 constituents (which include suspended sediments, nutrients, base cations, trace metals, and total dissolved solids) using a surrogate regression model approach and the Beale load estimator. Loads were typically higher for years with higher runoff. The proportional range of annual loads for total suspended solids, suspended-sediment concentrations, total phosphorus, and total lead, however, were 3.2 to 4.8 times larger than for annual runoff. Higher-than-expected annual sediment loads occurred in the years that also had some of the highest peak flows during the period, indicating that large storms are responsible for much of the sediment transport. Large development projects in proximity to streams also were related to years with high sediment loads. Yields from the Crooked Creek and North Fork Peachtree Creek watersheds were typically among the highest for 8 of the 12 constituents. These watersheds had the two highest amounts of developed medium plus high intensity land cover and the two highest percentages of imperviousness. Moderate to strong correlations were identified between seven of the constituent yields and the percentage of developed medium and high intensity land cover groups. Temporal trends in concentrations and loads were identified for 140 of the 300 possible watershed-time period-constituent combinations. There were substantially more negative than positive temporal trends identified during WYs 2003–10, whereas the number of negative and positive temporal trends were similar during WYs 2010–20. Measures of sediment transport had the most negative temporal trends. A few watersheds had consistent trends across several constituents; however, these trends did not appear to be associated with temporal changes in development or imperviousness. This study provides a thorough assessment of watershed characteristics, hydrology, and water-quality conditions and trends for the 15 study watersheds and can be used to identify possible factors that affect runoff and water quality and determine changes in water-quality conditions. Watershed managers can use these data and analyses to inform management decisions regarding the designated uses of streams, minimization of flooding, protection of aquatic habitats, and optimization of the effectiveness of best management practices.

Georgia

Updating and recalibrating the integrated Santa Rosa Plain Hydrologic Model to assess stream depletion and to simulate future climate and management scenarios in Santa Rosa, Sonoma County, California

The Santa Rosa Plain Hydrologic Model (SRPHM) was developed and published in 2014 through a collaboration between the U.S. Geological Survey (USGS) and Sonoma Water to analyze the hydrologic system in the Santa Rosa Plain watershed, help meet the increasing demand for fresh water, and prepare for future uncertainties in water resources. The original model simulated hydrological conditions and water use from water years 1975 to 2010. Recently (2023), the USGS, in cooperation with Sonoma Water and the California State Water Resources Control Board, updated the SRPHM model to extend its simulation period to the end of the 2018 calendar year, incorporate new estimates of rural and agricultural water use, and use efficient input format for climate variables. The updated model was recalibrated, and evaluation of the new model calibration is included in this report. This report presents the results of comparing the hydraulic heads, streamflow, and groundwater budget simulated by the updated model with those generated by the original model and observed data. The main difference in the simulated budget between the original and updated SRPHM is the estimates of agricultural pumping, rural domestic pumping, and return flow generated from rural water use that was not simulated in the original model. The revised agricultural pumping is simulated using the agricultural package, which constrains pumping to available groundwater. The use of the agricultural package leads to a more realistic estimation of agricultural water use, with revised agricultural pumping being one-third less than that in the original model. The revised rural pumping is about half of the pumping in the original model because of using detailed parcel data to estimate population density in rural areas instead of coarse census tracts. Overall, average total inflows for water years 2006–10 simulated by the updated model were about 2 percent less than the original model, and the average total updated outflows were nearly 5 percent less than the original model. The updated model was then used to generate stream depletion maps, simulate climate change scenarios during 2019–99, and simulate water rights allocation using the Model for Decision Support in Integrated River Basin Management (MODSIM). The results from simulating eight future climate scenarios indicated either an increase in groundwater storage or no significant change in the next 80 years, along with an increase in recharge, an increase in actual evapotranspiration in six out of eight climate projections, and an increase in surface runoff. The increases in the simulated future groundwater storage, recharge, evapotranspiration, and runoff in most climate projections are mainly driven by the projected increase in precipitation in most of the future climate scenarios. The updated model also was used to test a pilot case study demonstrating water-resource allocation among different users with different water rights using the integrated MODSIM-Groundwater and Surface-Water Flow Model (GSFLOW) platform. The updated SRPHM serves as a valuable tool for analyzing historical and future hydrologic conditions in the Santa Rosa Plain watershed and preparing for future uncertainties.

California

Development of a hydrologic flow model of the Zuñi Salt Lake and surrounding areas, west-central New Mexico

The terminal Zuñi Salt Lake is in a maar in west-​central New Mexico and contains hypersaline water that has long been used by Native Americans for religious purposes and the collection of salt. The U.S. Geological Survey (USGS), in cooperation with the Zuni Tribe of the Zuni Reservation, New Mexico, and Bureau of Reclamation, developed the Zuñi Salt Lake hydrologic flow model to simulate the steady-state conditions of the lake and surrounding groundwater-flow system. The model is a tool that can be used to analyze the potential hydrologic responses to resource development in the region under different water-use scenarios, which can support negotiations that could help protect future use of the lake. The USGS modular hydrologic model, MODFLOW 6, was used to simulate the hydrologic system of the lake and surrounding areas. To define the layering of the model, a three-dimensional hydrogeologic framework was constructed, which included seven informal hydrogeologic units and the locations of faults that truncate and offset the units. Calibration of the model was performed using the iterative ensemble smoother included in a parameter estimation software suite (PEST++). The iterative ensemble smoother approach resulted in a posterior parameter ensemble of model realizations with parameter and simulated values that show similar success in matching observations. Groundwater flux to the Zuñi Salt Lake simulated by using the posterior parameter distribution ranged from 26 to 730 acre-feet per year with a mean, median, and standard deviation of 530, 570, and 170 acre-feet per year, respectively. The relative contribution of groundwater discharge to the lake was simulated using a particle-tracking model for MODFLOW (MODPATH) to be from, in order of decreasing contribution, the Mesaverde, Cenozoic, Mancos, and Dakota hydrogeologic units. Simulations of future conditions using the ensemble of models can provide an empirical estimate of forecast uncertainty without substantial additional analysis.

New Mexico

User guide to Exploration and Graphics for RivEr Trends (EGRET) and dataRetrieval: R packages for hydrologic data

Evaluating long-term changes in river conditions (water quality and discharge) is an important use of hydrologic data. To carry out such evaluations, the hydrologist needs tools to facilitate several key steps in the process: acquiring the data records from a variety of sources, structuring it in ways that facilitate the analysis, processing the data with routines that extract information about changes that may be happening, and displaying findings with graphical techniques. A pair of tightly linked R packages, called dataRetrieval and EGRET (Exploration and Graphics for RivEr Trends), have been developed for carrying out each of these steps in an integrated manner. They are designed to easily accept data from three sources: U.S. Geological Survey hydrologic data, U.S. Environmental Protection Agency (EPA) STORET data, and user-supplied flat files. The dataRetrieval package not only serves as a “front end” to the EGRET package, it can also be used to easily download many types of hydrologic data and organize it in ways that facilitate many other hydrologic applications. The EGRET package has components oriented towards the description of long-term changes in streamflow statistics (high flow, average flow, and low flow) as well as changes in water quality. For the water-quality analysis, it uses Weighted Regressions on Time, Discharge and Season (WRTDS) to describe long-term trends in both concentration and flux. EGRET also creates a wide range of graphical presentations of the water-quality data and of the WRTDS results. This report serves as a user guide to these two R packages, providing detailed guidance on installation and use of the software, documentation of the analysis methods used, as well as guidance on some of the kinds of questions and approaches that the software can facilitate.

Techniques and Methods

Simulation of the shallow hydrologic system in the vicinity of Middle Genesee Lake, Wisconsin, using analytic elements and parameter estimation

Middle Genesee Lake is a ground-water flow-through lake located in a developing area in southeastern Wisconsin. Because the lake is in good connection with the shallow ground-water system, hydrologic stresses to the shallow ground-water system could adversely affect the lake system. In order to assess the effects of potential stresses on the lake, a study was completed by the U.S. Geological Survey, in cooperation with the Middle Genesee Lake Management District. The objective of the study was to identify areas that contribute ground water to the lake and estimate the hydrologic budget of the lake and hydraulic parameters affecting ground-water flow. A two-dimensional, steady-state analytic element model of the lake and surrounding area was developed using the computer code GFLOW. A parameter estimation model, UCODE, was used to optimize the calibration to measured water levels and streamflow. The calibrated model was used to evaluate the effect of three hypothetical stress scenarios on the stage of Middle Genesee Lake; the simulations were linked to UCODE, which formally incorporated parameter uncertainty into 95-percent confidence intervals around the simulated value. The scenarios included: (1) pumping from upgradient irrigation wells, (2) pumping from Lower Genesee Lake to lower lake levels, and (3) reduction in recharge resulting from development. The results of the simulations demonstrated that lake levels could be affected by hydrologic stresses in the shallow hydrologic system, with effects ranging from a 2.7 feet decline in lake stage resulting from pumping in Lower Genesee Lake to a 0.1 feet decline in lake stage from development in part of the upgradient recharge area. The range of lake stage decline increased when parameter uncertainty was included, from a decline of 3.1 feet for pumping from Lower Genesee Lake to no reduction in lake stage for the development in the recharge area. Whereas these simulated effects are within the natural variation in lake stage, they represent a systematic reduction of ground-water flow to the lake. Therefore, these hypothetical stresses are expected to establish a new, lower, baseline lake stage over which the natural variation due to climatic effects are added and subtracted.

Wisconsin

Hydrologic budgets for the Madison and Minnelusa aquifers, Black Hills of South Dakota and Wyoming, water years 1987-96

The Madison and Minnelusa aquifers are two of the most important aquifers in the Black Hills area of South Dakota and Wyoming. Quantification and evaluation of various hydrologic budget components are important for managing and understanding these aquifers. Hydrologic budgets are developed for two scenarios, including an overall budget for the entire study area and more detailed budgets for subareas. Budgets generally are combined for the Madison and Minnelusa aquifers because most budget components cannot be quantified individually for the aquifers. An average hydrologic budget for the entire study area is computed for water years 1987-96, for which change in storage is approximately equal to zero. Annual estimates of budget components are included in detailed budgets for nine subareas, which consider periods of decreasing storage (1987-92) and increasing storage (1993-96). Inflow components include recharge, leakage from adjacent aquifers, and ground-water inflows across the study area boundary. Outflows include springflow (headwater and artesian), well withdrawals, leakage to adjacent aquifers, and ground-water outflow across the study area boundary. Leakage, ground-water inflows, and ground-water outflows are difficult to quantify and cannot be distinguished from one another. Thus, net ground-water flow, which includes these components, is calculated as a residual, using estimates for the other budget components. For the overall budget for water years 1987-96, net ground-water outflow from the study area is computed as 100 ft3/s (cubic feet per second). Estimates of average combined budget components for the Madison and Minnelusa aquifers are: 395 ft3/s for recharge, 78 ft3/s for headwater springflow, 189 ft3/s for artesian springflow, and 28 ft3/s for well withdrawals. Hydrologic budgets also are quantified for nine subareas for periods of decreasing storage (1987-92) and increasing storage (1993-96), with changes in storage assumed equal but opposite. Common subareas are identified for the Madison and Minnelusa aquifers, and previous components from the overall budget generally are distributed over the subareas. Estimates of net ground-water flow for the two aquifers are computed, with net ground-water outflow exceeding inflow for most subareas. Outflows range from 5.9 ft3/s in the area east of Rapid City to 48.6 ft3/s along the southwestern flanks of the Black Hills. Net groundwater inflow exceeds outflow for two subareas where the discharge of large artesian springs exceeds estimated recharge within the subareas. More detailed subarea budgets also are developed, which include estimates of flow components for the individual aquifers at specific flow zones. The net outflows and inflows from the preliminary subarea budgets are used to estimate transmissivity of flow across specific flow zones based on Darcy?s Law. For estimation purposes, it is assumed that transmissivities of the Madison and Minnelusa aquifers are equal in any particular flow zone. The resulting transmissivity estimates range from 90 ft2/d to about 7,400 ft2/d, which is similar to values reported by previous investigators. The highest transmissivity estimates are for areas in the northern and southwestern parts of the study area, and the lowest transmissivity estimates are along the eastern study area boundary. Evaluation of subarea budgets provides confidence in budget components developed for the overall budget, especially regarding precipitation recharge, which is particularly difficult to estimate. Recharge estimates are consistently compatible with other budget components, including artesian springflow, which is a dominant component in many subareas. Calculated storage changes for subareas also are consistent with other budget components, specifically artesian springflow and net ground-water flow, and also are consistent with water-level fluctuations for observation wells. Ground-water budgets and flowpaths are especially complex i

Water-Resources Investigations Report