SEARCH · Geology Reports

Results for “Water Information Series Report”

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

SourceResearch connectionFrom yearThrough year

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

At least 721 records · Page 40Linked to original sources

USGS2025

Second Integrated Hydro-Terrestrial Modeling (IHTM 2.0) workshop: USGCRP federal agency perspectives

U.S. federal leadership, mobilizing close interagency coordination across sectors, is key to tackling compound climate and human impacts on our nation’s water resources. Recognizing this need, federal agencies and academic collaborators conducted a series of workshops to advance Integrated Hydro-Terrestrial Modeling (IHTM). IHTM is focused on supporting a multiscale framework to accelerate research insights, better integrate operational and planning perspectives, and bridge national-to-regional capabilities to address major interdependent societal water challenges. This framework leverages the capabilities of different agencies and institutional partners across sectors to advance a shared vision for use-inspired modeling and open science. The IHTM conceptual framework has developed over several community workshops in recent years: ● In September 2019, a community coordinating group on IHTM hosted an interagency workshop (“IHTM 1.0”) with support from NSF, DOE, and USGS. The resulting report identified the need for agencies to invest resources in designing pilot use cases, setting the stage for an integrated national IHTM capability. ● In November 2020, the U.S. Global Change Research Program (USGCRP) held a virtual Coastal-IHTM (C-IHTM) workshop in collaboration with the MultiSector Dynamics Community of Practice, which engages researchers across universities and national labs. Building on IHTM 1.0, this workshop aimed to identify coastal modeling capabilities and appropriate coastal applications of the IHTM framework. The workshop report detailed needs for open science, geographical or topical use cases, integrated modeling frameworks, and linking communities of practice. ● From October 31 to November 2, 2023, USGCRP hosted a hybrid IHTM 2.0 community workshop that brought together 160 scientists and managers from nine federal departments/agencies (DHS-FEMA, DOC-NOAA, DoD-USACE, DOE, DOI [U.S. Bureau of Reclamation (USBR), USGS], EPA, NASA, NSF, USDA [ARS, FS, NRCS]) and 32 non-federal academic, non-profit, and private institutions and consortiums. The workshop was organized by a federal Interagency Steering Committee with representatives from DOE, USGS, NSF, NASA, and NOAA, as well as a Scientific Organizing Committee comprising experts from academia and additional federal agency technical experts. The IHTM 2.0 workshop focused on moving the IHTM 1.0 concepts to action by making advances in designing national and regional integrated modeling experiments in five geographic testbeds. These modeling testbeds seek to facilitate collaborative water resources research and assessment that responds to societal needs, extending from basic research to resource management operations, including the assessment of water quantity, quality, and use at national and regional scales. Ultimately, we envision that the IHTM testbeds will improve our ability to provide water resources assessment scenarios to decision-makers and the public. This document summarizes USGCRP and member agencies’ perspectives of the IHTM 2.0 workshop. Sections 2 and 3 build on the work of the IHTM 2.0 workshop Scientific Organizing Committee to outline the importance of IHTM, summarize the workshop, and propose near- and medium-term IHTM priorities. Section 4 presents select ongoing and potential interagency activities that could help realize the IHTM 2.0 vision and align with three of USGCRP’s strategic goals for this decade: Advancing Science, Informing Decisions, and Engaging the Nation. Finally, Section 5 summarizes this report and offers information about how to connect with USGCRP’s IHTM activities. USGCRP and its member agencies are well-positioned to coordinate federal efforts to enhance IHTM capabilities in the coming years.

Conference Paper↗

USGS2015

Flood-Inundation Maps for the North River in Colrain, Charlemont, and Shelburne, Massachusetts, From the Confluence of the East and West Branch North Rivers to the Deerfield River

A series of 10 digital flood-inundation maps were developed for a 3.3-mile reach of the North River in Colrain, Charlemont, and Shelburne, Massachusetts, by the U.S. Geological Survey in cooperation with the Federal Emergency Management Agency. The coverage of the maps extends from the confluence of the East and West Branch North Rivers to the Deerfield River. Peak-flow estimates at the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities were computed for the reach from updated flood-frequency analyses. These peak flows were routed through a one-dimensional step-backwater hydraulic model to obtain the corresponding peak water-surface elevations and to place the tropical storm Irene flood of August 28, 2011, into historical context. The hydraulic model was calibrated by using the current [2015] stage-discharge relation at the U.S. Geological Survey streamgage North River at Shattuckville, MA (station number 01169000), and from documented high-water marks from the tropical storm Irene flood, which had a peak flow with approximately a 0.2-percent annual exceedance probability. A hydraulic model was used to compute water-surface profiles for 10 flood stages referenced to the streamgage and ranging from 6.6 feet (ft; 464.5 ft North American Vertical Datum of 1988 [which is approximately bankfull]) to 18.3 ft (476.2 ft North American Vertical Datum of 1988 [which is the stage of the 0.2-percent annual exceedance probability peak flow and exceeds the maximum recorded water level at the streamgage and the National Weather Service major flood stage of 13.0 ft]. The mapped stages of 6.6 to 18.3 ft were selected to match the stages of flows for bankfull; the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities; and an incremental stage of 17.0 ft. The simulated water-surface profiles were combined with a geographic information system digital elevation model derived from light detection and ranging (lidar) data with a 0.5-ft vertical accuracy to create a set of flood-inundation maps. The availability of the flood-inundation maps, combined with information regarding near-real-time stage from the U.S. Geological Survey North River at Shattuckville, MA streamgage can provide emergency management personnel and residents with information that is critical for flood response activities, such as evacuations and road closures, and postflood recovery efforts. The flood-inundation maps are nonregulatory, but provide Federal, State, and local agencies and the public with estimates of the potential extent of flooding during selected peak-flow events. Introduction

Massachusetts↗

USGS2015

Median nitrate concentrations in groundwater in the New Jersey Highlands Region estimated using regression models and land-surface characteristics

Nitrate-concentration data are used in conjunction with land-use and land-cover data to estimate median nitrate concentrations in groundwater underlying the New Jersey (NJ) Highlands Region. Sources of data on nitrate in 19,670 groundwater samples are from the U.S. Geological Survey (USGS) National Water Information System (NWIS) and the NJ Private Well Testing Act (PWTA). In a study conducted by the USGS, in cooperation with the New Jersey Department of Environmental Protection, logistic regression was used to relate measured nitrate concentrations to five explanatory variables (percent urban and agricultural land use, septic-system density, total length of streams, and number of known contaminated sites) quantified in 610-meter-square grid cells. A method for calculating the median concentrations of nitrate from a series of logistic regression models was developed. Two calibration and two validation procedures showed that the logistic-regression-based method can estimate groundwater-nitrate concentrations in the Highlands Region accurately to within 0.1 milligram per liter as nitrogen (mg/L as N). Limitations of the logistic-regression-based method include the inability to select a logistic model with exactly 0.5 probability of exceeding the threshold value and lack of an algorithm to directly calculate the median value. Quantile regression was evaluated as a suitable alternative and was slightly less accurate than the logistic-regression method in estimating median groundwater nitrate concentrations in the Highlands Region. Multiple-linear regression with log-transformed nitrate-concentration data and the same five explanatory values was less accurate than either logistic or quantile regression in estimating median nitrate concentrations. On the basis of 4,516 2000 x 2000 foot grid cells that contain wells with data stored in NWIS and the PWTA database, the estimated median nitrate concentration for the entire Highlands Region is about 1.25 mg/L as N, and estimated median concentrations range from about 1.05 to 1.78 mg/L as N among 11 smaller administratively defined areas within the Highlands Region that vary in percentages of urban land use, agricultural land use, and septic-system density. The Kaplan-Meier method of estimating summary statistics from left-censored data was applied in order to include nondetects (left-censored data) in median nitrate-concentration calculations. Median concentrations also were determined using three alternative methods of handling nondetects. Treatment of the 23 percent of samples that were nondetects had little effect on estimated median nitrate concentrations because method detection limits were mostly less than median values.

New Jersey↗

USGS2023

Lake Ontario April prey fish survey results and Alewife assessment, 2023

The April bottom trawl survey and Alewife Alosa pseudoharengus population assessment provides science to inform Lake Ontario fisheries management. The 2023 survey included 215 trawls in the main lake and embayments, and sampled depths from 6.5 to 252 m (21-833 ft). The survey captured 1,012,178 fish from 32 species with a total weight of 12,136 kg (26,700 lbs.). Alewife were 92% of the catch by number while Rainbow Smelt, Osmerus mordax , Deepwater Sculpin, Myoxocephalus thompsonii , and Round Goby, Neogobius melanostomus , comprised 3%, 3%, and 1% of the catch, respectively. To improve the accuracy of prey fish biomass and density estimates we reanalyzed trawl sensor data from each of three participating survey vessels and created vessel-specific relationships predicting how bottom trawl bottom contact time, wing width, and area-swept varies with depth. Total Alewife biomass increased in 2023 due to growth and survival of the abundant 2020 year class (now age-3) and an abundant 2022 year class (age-1). The 2023 mean Alewife biomass (81.1 kg·ha -1 ) was the largest since whole lake sampling began in 2016 and was the ninth largest value observed in the modern time series (1997-2023, maximum value in 2000 = 91.8 kg·ha -1 ). The 2023 Alewife density (6795 n·ha -1 ) was the greatest density observed in the modern time series. These high biomass and density values are due to above average Alewife reproductive success in 2020 and 2022. Simulation modeling suggests the 2024 and 2025 Alewife biomass index may be substantially higher than the 2023 observations. In 2023, the Rainbow Smelt biomass index increased relative to the 2022 index, as did the biomass index for Cisco, Coregonus artedi . In contrast, Emerald Shiner Notropis atherinoides and Threespine Stickleback Gasterosteus aculeatus , biomass values continue to be low (< 0.01 kg·ha -1 ). Three Bloater Coregonus hoyi , were captured during the 2023 survey. Hydroacoustic sampling conducted during the bottom trawl survey estimated prey fish densities in pelagic habitats not sampled by the bottom trawl (3 m below the surface to 3 m above the lake bottom) and these densities were hundreds to thousands of times lower than bottom trawl-based densities. These results support the idea that, in April, when the warmest water is on the lake bottom, Alewife and most other pelagic prey fish are near the lake bottom and can be effectively sampled with bottom trawling.

Lake Ontario↗

USGS2017

Response of currents and water quality to changes in dam operations in Hoover Reservoir, Columbus, Ohio, August 24–28, 2015

Hoover Reservoir, an important drinking water supply for the City of Columbus, Ohio, has been the source of a series of taste and odor problems in treated drinking water during the past few years. These taste and odor problems were caused by the compounds geosmin and 2-methylisoborneol, which are thought to have been related to cyanobacteria blooms. In an effort to reduce the phosphorus available for cyanobacteria blooms at fall turnover, the City of Columbus began experimenting with the dam’s selective withdrawal system to remove excess phosphorus in the hypolimnion, which is released from bottom sediments during summer anoxic conditions. The U.S. Geological Survey completed two synoptic survey campaigns to assess distributions of water quality and water velocity in the lower part of Hoover Reservoir to provide information on the changes to reservoir dynamics caused by changing dam operations. One campaign (campaign 1) was done while water was being withdrawn from the reservoir through the dam’s middle gate and the other (campaign 2) while water was being withdrawn through the dam’s lower gate. Velocities were measured using an acoustic Doppler current profiler, and water-quality parameters were measured using an autonomous underwater vehicle equipped with water-quality sensors. Along with the water-quality and water-velocity data, meteorological, inflow and outflow discharges, and independent water-quality data were compiled to monitor changes in other parameters that affect reservoir behavior. Monthly nutrient data, collected by the City of Columbus, were also analyzed for trends in concentration during periods of expected stratification. Based on the results of the two campaigns, when compared to withdrawing water through the middle gate, withdrawing water through the lower gate seemed to increase shear-driven mixing across the thermocline, which resulted in an increase in the depth of the epilimnion throughout the lower part of Hoover Reservoir. The observations from this study, if repeatable and driven primarily by changes in gate operations, can inform nutrient management strategies for Hoover Reservoir. Increased mixing across the thermocline may potentially supply nutrients from the hypolimnion to algae in the epilimnion. Although operation of the lower gate has the potential to export nutrients from the hypolimnion (where the concentrations of nutrients have typically been higher during summer months) through two mechanisms (direct withdrawal and mixing into the epilimnion), supply of nutrients to the epilimnion through enhanced mixing could lead to a short-term increase in algal populations. Therefore, further study is recommended to (1) test the repeatability of the results of gate changes on water-quality distributions and circulation patterns in lower Hoover Reservoir, (2) identify the immediate effect of gate changes on nutrient concentrations in the water column, and (3) identify the best management practices to reduce the nutrient storage in the hypolimnion of Hoover Reservoir without increasing the potential for nutrient transport to the highly productive epilimnion.

Ohio↗

USGS2006

Proceedings of the U.S. Geological Survey Sixth Biennial Geographic Information Science Workshop, Denver, Colorado, April 24-28, 2006

The U.S. Geological Survey's (USGS) Sixth Biennial Geographic Information Science Workshop April 24 - 28, 2006, at the Denver Federal Center in Denver, Colorado, provides a unique opportunity for multi-disciplinary Geographic Information Systems (GIS) and associated scientific professionals to share, learn, present, and discuss a wide variety of geospatial-related topics. Information is exchanged through a series of plenary sessions, hands-on technical workshops, user and commercial vendor demonstrations, lecture and poster sessions, and specialty meetings. Workshop attendance is limited to Department of Interior (DOI) and USGS employees and contractors; however, plenary, hands-on, and lecture sessions have presenters from other Federal agencies, numerous commercial vendors, universities, and several consortia. Over 175 participants are expected to attend the Workshop. Several prominent speakers are featured at this Workshop. Monday evening Star Guest Speaker and National Aeronautics and Space Administration (NASA) Astronaut Captain Dominic Gorie will talk about his experiences as a veteran of three space flights and over 32 days in space, including the NASA Space Shuttle Radar Topography Mission that mapped more than 47 million miles of the Earth's land surface. Selected as an astronaut candidate by NASA in December 1994, Captain Gorie is currently Chief of the Astronaut Shuttle Branch. Monday evening also features a town hall meeting with Geographic Information Office (GIO) leaders Karen Siderelis, Kevin Gallagher, Bob Pierce, Steve Guptill, Mark DeMulder, John Mahoney, and Mark Negri, who will discuss changes and activities within the GIO in an open discussion format. Tuesday plenary sessions feature keynote speaker Dr. P. Patrick Leahy, Acting USGS Director. Dr. Leahy holds undergraduate and graduate degrees in geology (1968) and geophysics (1970) from Boston College. In 1979, he received his doctorate in geology from Rensselaer Polytechnic Institute, where he specialized in regional ground-water studies and hydraulics. Dr. Leahy has been with the USGS since 1974, and has held various technical and managerial positions, including Associate Director for Geology and Chief of the National Water-Quality Assessment (NAWQA) Program. Dr. Leahy will be discussing a broad range of topics including Bureau accomplishments, initiatives, and budgets. The purpose of this proceedings volume is to serve as an activity reference for Workshop attendees as well as an archive of technical abstracts submitted, presented, and discussed at the Workshop. Author, co-author, and presenter names, affiliations, and contact information are listed with presentation titles along with submitted abstracts. Some hands-on sessions are offered twice. In these instances, abstracts submitted for publication are presented in the proceedings on both days they are offered. All acronyms used in these proceedings are explained in the text of each abstract. The term 'ArcGIS' refers to an integrated collection of GIS software products produced by Environmental Systems Research Institute, Inc. (ESRI).

Scientific Investigations Report↗

USGS2015

A model for evaluating stream temperature response to climate change in Wisconsin

Expected climatic changes in air temperature and precipitation patterns across the State of Wisconsin may alter future stream temperature and flow regimes. As a consequence of flow and temperature changes, the composition and distribution of fish species assemblages are expected to change. In an effort to gain a better understanding of how climatic changes may affect stream temperature, an approach was developed to predict and project daily summertime stream temperature under current and future climate conditions for 94,341 stream kilometers across Wisconsin. The approach uses a combination of static landscape characteristics and dynamic time-series climatic variables as input for an Artificial Neural Network (ANN) Model integrated with a Soil-Water-Balance (SWB) Model. Future climate scenarios are based on output from downscaled General Circulation Models (GCMs). The SWB model provided a means to estimate the temporal variability in groundwater recharge and provided a mechanism to evaluate the effect of changing air temperature and precipitation on groundwater recharge and soil moisture. The Integrated Soil-Water-Balance and Artificial Neural Network version 1 (SWB-ANNv1) Model was used to simulate daily summertime stream temperature under current (1990–2008) climate and explained 76 percent of the variation in the daily mean based on validation at 67 independent sites. Results were summarized as July mean water temperature, and individual stream segments were classified by thermal class (cold, cold transition, warm transition, and warm) for comparison of current (1990–2008) with future climate conditions. Integrating the SWB Model with the ANN Model provided a mechanism by which downscaled global or regional climate model results could be used to estimate the potential effects of climate change on future stream temperature on a daily time step. To address future climate scenarios, statistically downscaled air temperature and precipitation projections from 10 GCMs and 2 time periods were used with the SWB-ANNv1 Model to project future stream temperature. Projections of future stream temperatures at mid- (2046–65) and late- (2081–2100) 21st century showed the July mean water temperature increasing for all stream segments with about 80 percent of stream kilometers increasing by 1 to 2 degrees Celsius (°C) by mid-century and about 99 percent increasing by 1 to 3 °C by late-century. Projected changes in stream temperatures also affected changes in thermal classes with a loss in the total amount of cold-water, cold-transition, and warm-transition thermal habitat and a gain in warm-water and very warm thermal habitat for both mid- and late-21st century time periods. The greatest losses occurred for cold-water streams and the greatest gains for warm-water streams, with a contraction of cold-water streams in the Driftless Area of western and southern Wisconsin and an expansion of warm-water streams across northern Wisconsin. Results of this study suggest that such changes will affect the composition of fish assemblages, with a loss of suitable habitat for cold-water fishes and gain in suitable habitat for warm-water fishes. In the end, these projected changes in thermal habitat attributable to climate may result in a net loss of fisheries, because many warm-water species may be unable to colonize habitats formerly occupied by cold-water species because of other habitat limitations (e.g., stream size, gradient). Although projected stream temperatures may vary greatly, depending on the emissions scenario and models used, the results presented in this report represent one possibility. The relative change in stream temperature can provide useful information for planning for potential climate impacts to aquatic ecosystems. Model results can be used to help identify vulnerabilities of streams to climate change, guide stream surveys and thermal classifications, prioritize the allocation of scarce financial resources, identify approaches to climate adaptation to best protect and enhance resiliency in stream thermal habitat, and provide information to make quantitative assessments of statewide stream resources.

Wisconsin↗

USGS2010

Simulation of streamflow, evapotranspiration, and groundwater recharge in the lower San Antonio River Watershed, South-Central Texas, 2000-2007

The U.S. Geological Survey (USGS), in cooperation with the San Antonio River Authority, the Evergreen Underground Water Conservation District, and the Goliad County Groundwater Conservation District, configured, calibrated, and tested a watershed model for a study area consisting of about 2,150 square miles of the lower San Antonio River watershed in Bexar, Guadalupe, Wilson, Karnes, DeWitt, Goliad, Victoria, and Refugio Counties in south-central Texas. The model simulates streamflow, evapotranspiration (ET), and groundwater recharge using rainfall, potential ET, and upstream discharge data obtained from National Weather Service meteorological stations and USGS streamflow-gaging stations. Additional time-series inputs to the model include wastewater treatment-plant discharges, withdrawals for cropland irrigation, and estimated inflows from springs. Model simulations of streamflow, ET, and groundwater recharge were done for 2000-2007. Because of the complexity of the study area, the lower San Antonio River watershed was divided into four subwatersheds; separate HSPF models were developed for each subwatershed. Simulation of the overall study area involved running simulations of the three upstream models, then running the downstream model. The surficial geology was simplified as nine contiguous water-budget zones to meet model computational limitations and also to define zones for which ET, recharge, and other water-budget information would be output by the model. The model was calibrated and tested using streamflow data from 10 streamflow-gaging stations; additionally, simulated ET was compared with measured ET from a meteorological station west of the study area. The model calibration is considered very good; streamflow volumes were calibrated to within 10 percent of measured streamflow volumes. During 2000-2007, the estimated annual mean rainfall for the water-budget zones ranged from 33.7 to 38.5 inches per year; the estimated annual mean rainfall for the entire watershed was 34.3 inches. Using the HSPF model it was estimated that for 2000-2007, less than 10 percent of the annual mean rainfall on the study watershed exited the watershed as streamflow, whereas about 82 percent, or an average of 28.2 inches per year, exited the watershed as ET. Estimated annual mean groundwater recharge for the entire study area was 3.0 inches, or about 9 percent of annual mean rainfall. Estimated annual mean recharge was largest in water-budget zone 3, the zone where the Carrizo Sand outcrops. In water-budget zone 3, the estimated annual mean recharge was 5.1 inches or about 15 percent of annual mean rainfall. Estimated annual mean recharge was smallest in water-budget zone 6, about 1.1 inches or about 3 percent of annual mean rainfall. The Cibolo Creek subwatershed and the subwatershed of the San Antonio River upstream from Cibolo Creek had the largest and smallest basin yields, about 4.8 inches and 1.2 inches, respectively. Estimated annual ET and annual recharge generally increased with increasing annual rainfall. Also, ET was larger in zones 8 and 9, the most downstream zones in the watershed. Model limitations include possible errors related to model conceptualization and parameter variability, lack of data to quantify certain model inputs, and measurement errors. Uncertainty regarding the degree to which available rainfall data represent actual rainfall is potentially the most serious source of measurement error.

Texas↗

USGS2010

Simulation of Streamflow, Evapotranspiration, and Groundwater Recharge in the Lower San Antonio River Watershed, South-Central Texas, 2000-2007

Open-File Report↗

USGS2005

Studies of the Chesapeake Bay impact structure - Introduction and discussion

The late Eocene Chesapeake Bay impact structure on the Atlantic margin of Virginia is the largest known impact crater in the United States, and it may be the Earth's best preserved example of a large impact crater that formed on a predominantly siliciclastic continental shelf. The 85-kilometer-wide (53-milewide) crater also coincides with a region of saline ground water. It has a profound influence on ground-water quality and flow in an area of urban growth. The USGS-NASA Langley corehole at Hampton , Va., is the first in a series of new coreholes being drilled in the crater, and it is the first corehole to penetrate the entire crater-fill section and uppermost crystalline basement rock. The Langley corehole is located in the southwestern part of the crater's annular trough. A comprehensive effort to understand the crater's materials, architecture, geologic history , and formative processes, as well as its influence on ground water, includes the drilling of coreholes accompanied by high-resolution seismic-reflection and seismic-refraction surveys, audio-magnetotelluric surveys, and related multidisciplinary research. The studies of the core presented in this volume provide detailed information on the outer part of the crater, including the crystalline basement , the overlying impact -modified and impact -generated sediments (physical geology, paleontology, shocked minerals, and crystalline ejecta), and the upper Eocene to Quaternary postimpact sedimentary section (stratigraphy, paleontology, and paleoenvironments). The USGS-NASA Langley corehole has a total depth below land surface of 635.1 meters (m; 2,083.8 feet (ft)). The deepest unit in the corehole is the Neoproterozoic Langley Granite. The top of this granite at 626.3 m (2,054.7 ft) depth is overlain by 390.6 m (1,281.6 ft) of impact -modified and impact -generated siliciclastic sediments. These crater-fill materials are preserved beneath a 235.6-m-thick (773.12-ft-thick) blanket of postimpact sediments. A high-resolution seismic-reflection and seismic-refraction profile that crosses the Langley drill site is tied to the core by borehole geophysical logs, and it reveals the details of extensional collapse structures in the western annular trough. Electrical cross sections based on audio-magnetotelluric (AMT) soundings image a nearly vertical zone of high resistivity at the outer margin of the annular trough, possibly indicating fresh ground water at that location, and they show impedance trends that match the curvature of the structure . They also image the subsurface contact between conductive sediments and resistive crystalline basement , showing that the depth to crystalline basement is relatively constant in the western part of the annular trough. Chemical and isotopic data indicate that saline ground water of the Virginia inland saltwater wedge or bulge is a mixture of freshwater and seawater, and evidence for a mixing zone at the crater's outer margin supports the concept of differential flushing of residual seawater to create the bulge. Ground-water brine in the central part of the crater was produced by evaporation, and brine production from the heat of the impact is at least theoretically possible.

Chesapeake Bay↗

USGS2013

National wildlife refuge visitor survey 2012--Individual refuge results

The National Wildlife Refuge System (Refuge System), established in 1903 and managed by the U.S. Fish and Wildlife Service (Service), is the leading network of protected lands and waters in the world dedicated to the conservation of fish, wildlife and their habitats. There are 560 national wildlife refuges and 38 wetland management districts nationwide, encompassing more than 150 million acres. The Refuge System attracts nearly 45 million visitors annually, including 34.8 million people who observe and photograph wildlife, 9.6 million who hunt and fish, and nearly 675,000 teachers and students who use refuges as outdoor classrooms. Understanding visitor perceptions of refuges and characterizing their experiences on refuges are critical elements of managing these lands and meeting the goals of the Refuge System. The Service collaborated with the U.S. Geological Survey to conduct a national survey of visitors regarding their experiences on national wildlife refuges. The purpose of the survey was to better understand visitor experiences and trip characteristics, to gauge visitors’ levels of satisfaction with existing recreational opportunities, and to garner feedback to inform the design of programs and facilities. The survey results will inform performance, planning, budget, and communications goals. Results will also inform Comprehensive Conservation Plans (CCPs), visitor services, and transportation planning processes. This Data Series consists of 25 separate data files. Each file describes the results of the survey for an individual refuge and contains the following information: • Introduction: An overview of the Refuge System and the goals of the national surveying effort. • Methods: The procedures for the national surveying effort, including selecting refuges, developing the survey instrument, contacting visitors, and guidance for interpreting the results.• Refuge Description: A brief description of the refuge location, acreage, purpose, recreational activities, and visitation statistics, including a map (where available) and refuge website link. • Sampling at This Refuge: The sampling periods, locations, and response rate for the refuge. • Selected Survey Results: Key findings for the refuge, including: o Visitor and trip characteristics o Visitor spending in the local communities o Visitors opinions about the refuge o Visitor opinions about National Wildlife Refuge System topics • Conclusion • References Cited • Survey Frequencies (Appendix A): The survey instrument with frequency results for the refuge. • Visitor Comments (Appendix B): The verbatim responses to the open-ended survey questions for the refuge.Individual-refuge results for the 53 participating refuges in the 2010-2011 national effort are available at http://pubs.usgs.gov/ds/643/ as part of USGS Data Series 643. Combined results for the 53 participating refuges in the 2010-2011 national effort are available at http://pubs.usgs.gov/ds/685/ as part of USGS Data Series 685.

Data Series↗

USGS2004

Conceptualization and simulation of the Edwards aquifer, San Antonio region, Texas

A new numerical ground-water-flow model (Edwards aquifer model) that incorporates important components of the latest information and plausible conceptualization of the Edwards aquifer was developed. The model includes both the San Antonio and Barton Springs segments of the Edwards aquifer in the San Antonio region, Texas, and was calibrated for steady-state (1939–46) and transient (1947–2000) conditions, excluding Travis County. Transient simulations were conducted using monthly recharge and pumpage (withdrawal) data. The model incorporates conduits simulated as continuously connected (other than being separated in eastern Uvalde and southwestern Medina Counties), one-cell-wide (1,320 feet) zones with very large hydraulic-conductivity values (as much as 300,000 feet per day). The locations of the conduits were based on a number of factors, including major potentiometric-surface troughs in the aquifer, the presence of sinking streams, geochemical information, and geologic structures (for example, faults and grabens). The simulated directions of flow in the Edwards aquifer model are most strongly influenced by the presence of simulated conduits and barrier faults. The simulated flow in the Edwards aquifer is influenced by the locations of the simulated conduits, which tend to facilitate flow. The simulated subregional flow directions generally are toward the nearest conduit and subsequently along the conduits from the recharge zone into the confined zone and toward the major springs. Structures simulated in the Edwards aquifer model influencing groundwater flow that tend to restrict flow are barrier faults. The influence of simulated barrier faults on flow directions is most evident in northern Medina County. A water budget is an accounting of inflow to, outflow from, and storage change in the aquifer. For the Edwards aquifer model steady-state simulation, recharge (from seepage losses from streams and infiltration of rainfall) accounts for 93.5 percent of the sources of water to the Edwards aquifer, and inflow through the northern and northwestern model boundaries contributes 6.5 percent. The largest discharges are springflow (73.7 percent) and ground-water withdrawals by wells (25.7 percent). The principal source of water to the Edwards aquifer for the Edwards aquifer model transient simulation was recharge, constituting about 60 percent of the sources of water (excluding change in storage) to the Edwards aquifer during 1956, a drought period, and about 97 percent of the sources (excluding change in storage) during 1975, a period of above-normal rainfall and recharge. The principal discharges from the Edwards aquifer for the transient simulation were springflow and withdrawals by wells. During 1956, representing drought conditions, the change in storage (net water released from storage) was much greater than recharge, accounting for 75.9 percent of the total flow compared to 14.5 percent for recharge. Conversely, during 1975, representing above-normal rainfall and recharge conditions, recharge constituted 79.9 percent of the total flow, compared to 7.1 percent for the change in storage (net water added to storage). A series of sensitivity tests was made to ascertain how the model results were affected by variations greater than and less than the calibrated values of input data. Simulated hydraulic heads in the Edwards aquifer model were most sensitive to recharge, withdrawals, hydraulic conductivity of the conduit segments, and specific yield and were comparatively insensitive to spring-orifice conductance, northern boundary inflow, and specific storage. Simulated springflow in the Edwards aquifer model was most sensitive to recharge, withdrawals, hydraulic conductivity of the conduit segments, specific yield, and increases in northern boundary inflow and was comparatively insensitive to spring-orifice conductance and specific storage.

Texas↗

USGS2013

Hydrology and water quality of Shell Lake, Washburn County, Wisconsin, with special emphasis on the effects of diversion and changes in water level on the water quality of a shallow terminal lake

Shell Lake is a relatively shallow terminal lake (tributaries but no outlets) in northwestern Wisconsin that has experienced approximately 10 feet (ft) of water-level fluctuation over more than 70 years of record and extensive flooding of nearshore areas starting in the early 2000s. The City of Shell Lake (City) received a permit from the Wisconsin Department of Natural Resources in 2002 to divert water from the lake to a nearby river in order to lower water levels and reduce flooding. Previous studies suggested that water-level fluctuations were driven by long-term cycles in precipitation, evaporation, and runoff, although questions about the lake’s connection with the groundwater system remained. The permit required that the City evaluate assumptions about lake/groundwater interactions made in previous studies and evaluate the effects of the water diversion on water levels in Shell Lake and other nearby lakes. Therefore, a cooperative study between the City and U.S. Geological Survey (USGS) was initiated to improve the understanding of the hydrogeology of the area and evaluate potential effects of the diversion on water levels in Shell Lake, the surrounding groundwater system, and nearby lakes. Concerns over deteriorating water quality in the lake, possibly associated with changes in water level, prompted an additional cooperative project between the City and the USGS to evaluate efeffects of changes in nutrient loading associated with changes in water levels on the water quality of Shell Lake. Numerical models were used to evaluate how the hydrology and water quality responded to diversion of water from the lake and historical changes in the watershed. The groundwater-flow model MODFLOW was used to simulate groundwater movement in the area around Shell Lake, including groundwater/surface-water interactions. Simulated results from the MODFLOW model indicate that groundwater flows generally northward in the area around Shell Lake, with flow locally converging toward the lake. Total groundwater inflow to Shell Lake is small (approximately 5 percent of the water budget) compared with water entering the lake from precipitation (83 percent) and surface-water runoff (13 percent). The MODFLOW model also was used to simulate average annual hydrologic conditions from 1949 to 2009, including effects of the removal of 3 billion gallons of water during 2003–5. The maximum decline in simulated average annual water levels for Shell Lake due to the diversion alone was 3.3 ft at the end of the diversion process in 2005. Model simulations also indicate that although water level continued to decline through 2009 in response to local weather patterns (local drought), the effects of the diversion decreased after the diversion ceased; that is, after 4 years of recovery (2006–9), drawdown attributable to the diversion alone decreased by about 0.6 ft because of increased groundwater inflow and decreased lake-water outflow to groundwater caused by the artificially lower lake level. A delayed response in drawdown of less than 0.5 ft was transmitted through the groundwater-flow system to upgradient lakes. This relatively small effect on upgradient lakes is attributed in part to extensive layers of shallow clay that limit lake/groundwater interaction in the area. Data collected in the lake indicated that Shell Lake is polymictic (characterized by frequent deep mixing) and that its productivity is limited by the amount of phosphorus in the lake. The lake was typically classified as oligotrophic-mesotrophic in June, mesotrophic in July, and mesotrophic-eutrophic in August. In polymictic lakes like Shell Lake, phosphorus released from the sediments is not trapped near the bottom of the lake but is intermittently released to the shallow water, resulting in deteriorating water quality as summer progresses. Because the productivity of Shell Lake is limited by phosphorus, the sources of phosphorus to the lake were quantified, and the response in water quality to changes in phosphorus inputs were evaluated by means of eutrophication models. During 2009, the total input of phosphorus to Shell Lake was 1,730 pounds (lb), of which 1,320 lb came from external sources (76 percent) and 414 lb came from internal loading from sediments in the lake (24 percent). The largest external source was from surface-water runoff, which delivered about 52 percent of the total phosphorus load compared with about 13 percent of the water input. The second largest source was from precipitation (wetfall and dryfall), which delivered 19 percent of the load compared to about 83 percent of the water input. Contributions from septic systems and groundwater accounted for about 3 and 2 percent, respectively. Increased runoff raises water levels in the lake but does not necessarily increase phosphorus loading because phosphorus concentrations in the tributaries decline during increased flow, possibly because of shorter retention times in upstream wetlands. Phosphorus loading to the lake in 2009 represented what occurred after a series of dry years; therefore, this information was combined with data from 2011, a wet year, to estimate phosphorus loading during a range of hydrologic conditions by estimating loading from each component of the phosphorus budget for each year from 1949 to 2011. Comparisons of historical water-quality records with historical water levels and applications of a hydrodynamic model (Dynamic Lake Model, DLM) and empirical eutrophication models were used to understand how changes in water level and the coinciding changes in phosphorus loading affect the water quality of Shell Lake. DLM simulations indicate that large changes in water level (approximately 10 ft) affect the persistence of stratification in the lake. During periods with low water levels, the lake is a well-mixed, polymictic system, with water quality degrading slightly as summer progresses. During periods with high water levels, the lake is more stratified, and phosphorus from internal loading is trapped in the hypolimnion and released later in summer, which results in more extreme seasonality in water quality and better clarity in early summer. Results of eutrophication model simulations using a range in external phosphorus inputs illustrate how water quality in Shell Lake (phosphorus and chlorophyll a concentrations and Secchi depths) responds to changes in external phosphorus loading. Results indicate that a 50-percent reduction in external loading from that measured in 2009 would be required to change phosphorus concentrations from 0.018 milligram per liter (mg/L) (measured in 2009) to 0.012 mg/L (estimated for the mid-1800s from analysis of diatoms in sediment cores). Such reductions in phosphorus loading cannot be accomplished by targeting septic systems or internal loading alone because septic systems contribute only about 3 percent of the phosphorus input to the lake, and internal loading from the sediments of Shell Lake contributes only about 25 percent of phosphorus input. Complete elimination of phosphorus from septic systems and internal loading would decrease the phosphorus concentrations in the lake by 0.003–0.004 mg/L. Therefore, reducing phosphorus concentration in the lake more than by 0.004 mg/L requires decreasing phosphorus loading from surface-water contributions, primarily runoff to the lake. Reconstructed changes in water quality from 1860 to 2010, based on changes in the diatom communities archived in the sediments and eutrophication model simulations, suggest that anthropogenic changes in the watershed (sawmill construction in 1881; the establishment of the village of Shell Lake; and land-use changes in the 1920s, including increased agriculture) had a much larger effect on water quality than the natural changes associated with fluctuations in water level. Although the effects of natural changes in water level on water quality appear to be small, changes in water level do have a modest effect on water quality, primarily manifested as small improvements during higher water levels. Fluctuations in water level, however, have a larger effect on the seasonality of water-quality patterns, with better water quality, especially increased Secchi depths, in early summer during years with high water levels.

Wisconsin↗

USGS2010

Seasonal and Spatial Distribution of Freshwater Flow and Salinity in the Ten Thousand Islands Estuary, Florida, 2007-2009

The watershed of the Ten Thousand Islands (TTI) estuary has been substantially altered through the construction of canals and roads for the Southern Golden Gate Estates (SGGE), Barron River Canal, and U.S. 41 (Tamiami Trail). Two restoration projects designed to improve freshwater delivery to the estuary are the Picayune Strand Restoration Project, which includes the Southern Golden Gate Estates, and the Tamiami Trail Culverts Project; both are part of the Comprehensive Everglades Restoration Plan. To address hydrologic information needs critical for monitoring the effects of these restoration projects, the U.S. Geological Survey initiated a study in October 2006 to characterize freshwater outflows from the rivers, internal circulation and mixing within the estuary, and surface-water exchange between the estuary and Gulf of Mexico. The effort is conducted in cooperation with the South Florida Water Management District and complemented by monitoring performed by the Rookery Bay National Estuarine Research Reserve. Surface salinity was measured during moving boat surveys using a flow-through system that operated at planing speeds averaging 20 miles per hour. The data were logged every 10 seconds by a data recorder that simultaneously logged location information from a Global Positioning System. The major rivers, bays, and nearshore Gulf of Mexico region of the TTI area were surveyed in approximately 5 hours by two boats traversing about 200 total miles. Salinity and coordinate data were processed using inverse distance weighted interpolation to create salinity contour maps of the entire TTI region. Ten maps were created from salinity surveys performed between May 2007 and May 2009 and illustrate the dry season, transitional, and wet season salinity patterns of the estuarine rivers, inner bays, mangrove islands, and Gulf of Mexico boundary. The effects of anthropogenic activities are indicated by exceptionally low salinities associated with point discharge into the estuary from the Faka Union Canal and Barron River during the wet season. Low salinities in Faka Union Bay may cause reduced diversity and density of submerged aquatic vegetation, fish, and benthic organisms compared with neighboring Fakahatchee Bay. The Faka Union Canal System reduced the size of the watershed for the western TTI estuary, resulting in increased wet season salinities compared to those for the eastern TTI estuary, the watershed of which is composed of the relatively pristine Fakahatchee Strand Preserve State Park. Minimal river discharge and high evaporation caused hypersaline conditions to develop throughout the entire TTI region during the dry season. The 2007-2008 drought and passage of Tropical Storm Fay on August 18-19, 2008, demonstrated the effects of seasonal rainfall on salinity patterns, with substantially higher salinities observed during the 2007 wet season compared to those for the 2008 wet season. The salinity maps, coupled with data from the monitoring stations, provide baseline information of seasonal and spatial distribution of freshwater flow and salinity in the TTI estuary, and a means of monitoring the effects of restoration in improving freshwater delivery to the estuary.

Data Series↗

USGS2009

The Farm Process Version 2 (FMP2) for MODFLOW-2005 - Modifications and Upgrades to FMP1

The ability to dynamically simulate the integrated supply-and-demand components of irrigated agricultural is needed to thoroughly understand the interrelation between surface water and groundwater flow in areas where the water-use by vegetation is an important component of the water budget. To meet this need, the computer program Farm Process (FMP1) was updated and refined for use with the U.S. Geological Survey's MODFLOW-2005 groundwater-flow model, and is referred to as MF2005-FMP2. The updated program allows the simulation, analysis, and management of nearly all components of human and natural water use. MF2005-FMP2 represents a complete hydrologic model that fully links the movement and use of groundwater, surface water, and imported water for water consumption of irrigated agriculture, but also of urban use, and of natural vegetation. Supply and demand components of water use are analyzed under demand-driven and supply-constrained conditions. From large- to small-scale settings, the MF2005-FMP2 has the unique set of capabilities to simulate and analyze historical, present, and future conditions. MF2005-FMP2 facilitates the analysis of agricultural water use where little data is available for pumpage, land use, or agricultural information. The features presented in this new version of FMP2 along with the linkages to the Streamflow Routing (SFR), Multi-Node Well (MNW), and Unsaturated Zone Flow (UZF) Packages prevents mass loss to an open system and helps to account for 'all of the water everywhere and all of the time'. The first version, FMP1 for MODFLOW-2000, is limited to (a) transpiration uptake from unsaturated root zones, (b) on-farm efficiency defined solely by farm and not by crop type, (c) a simulation of water use and returnflows related only to irrigated agriculture and not also to non-irrigated vegetation, (d) a definition of consumptive use as potential crop evapotranspiration, (e) percolation being instantly recharged to the uppermost active aquifer, (f) automatic routing of returnflow from runoff either to reaches of tributary stream segments adjacent to a farm or to one reach nearest to the farm's lowest elevation, (g) farm-well pumping from cell locations regardless of whether an irrigation requirement from these cells exists or not, and (h) specified non-routed water transfers from an undefined source outside the model domain. All of these limitations are overcome in MF2005-FMP2. The new features include (a) simulation of transpiration uptake from variably saturated, fully saturated, or ponded root zones (for example, for crops like rice or riparian vegetation), (b) definition of on-farm efficiency not only by farm but also by crop, (c) simulation of water use and returnflow from non-irrigated vegetation (for example, rain-fed agriculture or native vegetation), (d) use of crop coefficients and reference evapotranspiration, (e) simulation of the delay between percolation from farms through the unsaturated zone and recharge into the uppermost active aquifer by linking FMP2 to the UZF Package, (f) an option to manually control the routing of returnflow from farm runoff to streams, (g) an option to limit pumping to wells located only in cells where an irrigation requirement exists, and (h) simulation of water transfers to farms from a series of well fields (for example, recovery well field of an aquifer-storage-and-recovery system, ASR). In addition to the output of an economic budget for each farm between irrigation demand and supply ('Farm Demand and Supply Budget' in FMP1), a new output option called 'Farm Budget' was created for FMP2, which allows the user to track all physical flows into and out of a water accounting unit at all times. Such a unit can represent individual farms, farming districts, natural areas, or urban areas. The example model demonstrates the application of MF2005-FMP2 with delayed recharge through an unsaturated zone, rejected infiltration in a riparian area, changes in de

Techniques and Methods↗

USGS2026

Wind River subbasin restoration: Annual report of U.S. Geological Survey activities January 2024 through December 2024

We sampled juvenile wild Oncorhynchus mykiss (Steelhead Trout) in headwater streams of the Wind River, Washington, to characterize population attributes and investigate life-history metrics, particularly migratory patterns, and early life-stage survival. Additionally, this sampling provides an index of abundance of non-native Brook Trout Salvelinus confluentus . We used passive integrated transponder (PIT) tagging and a series of instream PIT-tag interrogation systems (PTISs) to track juveniles and adults. The Wind River subbasin is considered a wild Steelhead refuge by Washington Department of Fish and Wildlife (WDFW). No hatchery Steelhead Trout have been released in the Wind River subbasin since 1997, and hatchery adults are estimated at less than one percent of spawners in most years. Over 25 years of Steelhead Trout status and trend monitoring and research in the subbasin is contributing to understanding of population response to numerous restoration actions in the subbasin, including removal of Hemlock Dam from Trout Creek in 2009, which had an inadequate adult ladder and contributed to increased water temperatures, to near lethal levels for salmonids. Data from our study, and companion work by WDFW, are contributing to the Columbia Basin Fish and Wildlife Program (2009) Research, Monitoring, and Evaluation (RM&E) Strategy of Fish Population Status Monitoring. Specifically, this work addresses the sub-strategies of 1) Assessing the Status and Trends of Diversity of Natural Origin Fish Populations and Uncertainties Research regarding differing life histories of a wild Steelhead Trout population, 2) Assessing the Status and Trend of Adult Natural Origin Fish Populations, and 3) Monitoring and Evaluating the Effectiveness of Tributary Habitat Actions Relative to Environmental, Physical, or Biological Performance Objectives. During summer and fall 2024, we PIT-tagged 2,075 Steelhead Trout parr (age-0 and age-1), in the Trout Creek and upper Wind River watersheds. Age-0 and age-1 parr were at high abundance due to strong spawner numbers in 2023 and 2024 (Estimate of 932 in 2023; and a moderate number of spawners in 2024 though an estimate is not available at this time; 24-year median = 619; Charlie Cochran, WDFW Fish Biologist, personal commun., 2025). An additional 100 parr (age-0 and age-2 or older) were tagged to provide fish for estimating detection efficiencies at PTISs. Steelhead Trout parr were recaptured and detected through repeat headwater sampling, smolt trapping, instream PTISs and Columbia River mainstem PIT-tag detection. We maintained and upgraded six instream PTISs to detect PIT-tagged Steelhead Trout parr, smolts, and adults, providing data for population assessments, and life-cycle research. Detection data from PIT-tagged adult Steelhead Trout at our four primary PTISs allow assessment of adult escapement to tributary watersheds within the Wind River subbasin. Most adult Steelhead detected were tagged by WDFW at an adult trap in the lower Wind River and some as smolts at screw traps. Adult Steelhead Trout detection efficiency estimates at our primary PTIS in Trout Creek have been greater than 99 percent during eight of the past ten years and have exceeded 97% at our primary PTIS in the Wind River during nine of the past ten years. Adult escapement estimates to tributary watersheds are helping evaluate the efficacy of the 2009 removal of Hemlock Dam from rkm 2.0 of Trout Creek. Detections at the instream PTISs have shown trends of age-0 and age-1 Steelhead Trout parr emigration from natal areas during summer and fall, in addition to the expected movement of parr and smolts in spring. Our data suggest that most fish from a cohort that migrate downstream will do so at age-1 for additional rearing in the subbasin. We have estimated that from 15% to 56% of parr tagged as age-0 fish in headwater areas make downstream migrations at age-1 for additional rearing (median = 29%). We have estimated that from 1% to 27% of Steelhead Trout parr, tagged as age-1 fish, make downstream migrations during fall (median = 3.5%): this is especially pronounced in the upper Wind River portion of the watershed (median 10%). These findings raise questions about parr rearing habitat and whether migrations are density- or habitat-quality driven, and answers to such questions are long-term goals of this study as active and passive habitat restoration occurs. Repeat sampling at sites in the subbasin within and between years has enabled assessment of juvenile Steelhead Trout growth patterns. Growth rates (relative change in weight) of age-0 PIT-tagged parr during summer have been similar across the subbasin, though slightly lower in the Trout Creek watershed (Trout Creek sites, 0.57 - 1.25%/day, Wind River sites = 0.94 – 2.48%/day). Summer growth rates have been lower for age-1 parr in the Trout Creek watershed than the upper Wind River watershed (Trout Creek sites, = 0.03 - 0.14%/day, Wind River sites = 0.05 – 0.19%/day). Yearly relative growth was similar across the subbasin for both age-0 and age-1 tagged parr. During 2024 range of mean daily relative change in weight from the previous year for age-0 tagged fish was: Trout Creek sites = 0.96 – 1.41%/day, Wind River sites = 0.80 – 1.30 %/day; and for age1 tagged fish was: Trout Creek sites = 0.37 - 0.50 %/day, Wind River sites = no age-1 tagged recaptures. Non-native Brook Trout are present in the subbasin, chiefly the Trout Creek watershed. Mean percent-of-catch that is Brook Trout, at four sample sites in Trout Creek, has declined from the period 1998 – 2003 to the period 2011 – 2024 (range of decrease = 19 - 82%). Percent-of-catch (of 64 samples, 52 < 10%; 3 > 20%) and number of Brook Trout at the Trout Creek sites from 2011 through 2024, though variable, have generally declined. Habitat restoration evaluation efforts are critical to ensure efficient use of money and resources but require long-term monitoring of population metrics. Detailed viable salmonid population and life-history data, such as that provided by PIT-tagging and instream PTIS networks inform fisheries policy and management and enable assessment of long-term effects of habitat restoration actions such as the removal of Hemlock Dam on Trout Creek and proposed major instream habitat restoration in the upper Wind River.

Washington↗

USGS2005

Classification of hydrogeologic areas and hydrogeologic flow systems in the basin and range physiographic province, southwestern United States

The hydrogeology of the Basin and Range Physiographic Province in parts of Arizona, California, New Mexico, Utah, and most of Nevada was classified at basin and larger scales to facilitate information transfer and to provide a synthesis of results from many previous hydrologic investigations. A conceptual model for the spatial hierarchy of the hydrogeology was developed for the Basin and Range Physiographic Province and consists, in order of increasing spatial scale, of hydrogeologic components, hydrogeologic areas, hydrogeologic flow systems, and hydrogeologic regions. This hierarchy formed a framework for hydrogeologic classification. Hydrogeologic areas consist of coincident ground-water and surface-water basins and were delineated on the basis of existing sets of basin boundaries that were used in past investigations by State and Federal government agencies. Within the study area, 344 hydrogeologic areas were identified and delineated. This set of basins not only provides a framework for the classification developed in this report, but also has value for regional and subregional purposes of inventory, study, analysis, and planning throughout the Basin and Range Physiographic Province. The fact that nearly all of the province is delineated by the hydrogeologic areas makes this set well suited to support regional-scale investigations. Hydrogeologic areas are conceptualized as a control volume consisting of three hydrogeologic components: the soils and streams, basin fill, and consolidated rocks. The soils and streams hydrogeologic component consists of all surface-water bodies and soils extending to the bottom of the plant root zone. The basin-fill hydrogeologic component consists of unconsolidated and semiconsolidated sediment deposited in the structural basin. The consolidated-rocks hydrogeologic component consists of the crystalline and sedimentary rocks that form the mountain blocks and basement rock of the structural basin. Hydrogeologic areas were classified into 19 groups through a cluster analysis of 8 characteristics of each area's hydrologic system. Six characteristics represented the inflows and outflows of water through the soils and streams, basin fill, and consolidated rocks, and can be used to determine the hydrogeologic area's position in a hydrogeologic flow system. Source-, link-, and sink-type hydrogeologic areas have outflow but not inflow, inflow and outflow, and inflow but not outflow, respectively, through one or more of the three hydrogeologic components. Isolated hydrogeologic areas have no inflow or outflow through any of the three hydrogeologic components. The remaining two characteristics are indexes that represent natural recharge and discharge processes and anthropogenic recharge and discharge processes occurring in the hydrogeologic area. Of the 19 groups of hydrogeologic areas, 1 consisted of predominantly isolated-type hydrogeologic areas, 7 consisted of source-type hydrogeologic areas, 9 consisted of link-type hydrogeologic areas, and 2 consisted of sink-type hydrogeologic areas. Groups comprising the source-, link-, and sink-type hydrogeologic areas can be distinguished between each other on the basis of the hydrogeologic component(s) through which interbasin flow occurs, as well as typical values for the two indexes. Conceptual models of the hydrologic systems of a representative hydrogeologic area for each group were developed to help distinguish groups and to synthesize the variation in hydrogeologic systems in the Basin and Range Physiographic Province. Hydrogeologic flow systems consist of either a single isolated hydrogeologic area or a series of multiple hydrogeologic areas that are hydraulically connected through interbasin flows. A total of 54 hydrogeologic flow systems were identified and classified into 9 groups. One group consisted of single isolated hydrogeologic areas. The remaining eight groups consisted of multiple hydrogeologic areas and were distinguished o

Professional Paper↗

USGS2020

Hydrogeology and groundwater geochemistry of till confining units and confined aquifers in glacial deposits near Litchfield, Cromwell, Akeley, and Olivia, Minnesota, 2014–18

Confined (or buried) aquifers of glacial origin overlain by till confining units provide drinking water to hundreds of thousands of Minnesota residents. The sustainability of these groundwater resources is not well understood because hydraulic properties of till that control vertical groundwater fluxes (leakage) to underlying aquifers are largely unknown. The U.S. Geological Survey, Iowa State University, Minnesota Geological Survey, and Minnesota Department of Health investigated hydraulic properties and groundwater flow through till confining units using field studies and heuristic MODFLOW simulations. Till confining units in the following late-Wisconsinan stratigraphic units (with locations in parentheses) were characterized: Des Moines lobe till of the New Ulm Formation (Litchfield, Minnesota), Superior lobe till of the Cromwell and Aitkin Formations (Cromwell, Minn.), and Wadena lobe till of the Hewitt Formation (hydrogeology field camp [HFC] near Akeley, Minn.). Pre-Illinoian till of the Good Thunder formation (Olivia, Minn.) was also characterized. Hydraulic and geochemical field data were collected from sediment cores and a series of five piezometer nests. Each nest consisted of five to eight piezometers screened at short vertical intervals in hydrostratigraphic units including (if present) surficial aquifers, till confining units, confined/buried aquifers, and underlying bedrock. Till hydraulic conductivity was estimated from slug tests (horizontal [ K h ]) and constant-rate aquifer tests in the confined aquifer (vertical [ K v ]). Travel times through the till were evaluated with Darcy’s law and stable isotope concentrations. A series of heuristic MODFLOW simulations were used to evaluate groundwater fluxes through till across the range of till hydraulic properties and pumping rates observed at the field sites. The field data demonstrated variability in hydraulic properties between and within till stratigraphic units horizontally and vertically. The variability in hydraulic properties within and between sites resulted in substantial differences in groundwater flux through till. A conceptual understanding that emerges from the vertical till profiles is that they are not homogeneous hydrostratigraphic units with uniform properties; rather, each vertical sequence is a heterogeneous mixture of glacial sediment with differing abilities to transmit water. Till thicknesses varied from 60 to 166 feet, and till textures ranged from a sandy loam (Hewitt Formation, HFC site) to a silt loam/clay loam (Good Thunder formation, Olivia site). Till Kh varied by one to three orders of magnitude within each piezometer nest. Four piezometer nests had downward hydraulic gradients ranging from 0.04 to 0.56, and one nest had a slight upward hydraulic gradient of 0.02. The Cromwell, HFC, and Litchfield 1 sites were examples of “leaky” tills with high Kv (0.001 to 1.1 feet per day [ft/d]) and geometric mean Kh (0.03 to 0.07 ft/d) and extensive vertical hydraulic connectivity between the confined aquifer and the overlying till. Estimated groundwater travel times through these sites ranged from 1 to 81 years, and two of these sites had tritium throughout their till profiles. The tills at the other two sites, Olivia and Litchfield 2, were effective confining units that had low Kv (0.001 to 0.0005 ft/d) and geometric mean Kh (0.0002 to 0.004 ft/d). The till piezometers at these sites had no drawdown response to short-term (up to 10 hours for Olivia and up to 5 days for Litchfield) high-capacity pumping from the confined aquifer. Estimated groundwater travel times through the tills at these sites ranged from 165 to nearly 1,800 years, and tritium was only detected in the upper one-third of these till profiles. Across all sites, the till vertical anisotropy (ratio of K h to K v ) ranged by four orders of magnitude from 0.05 at the Cromwell nest to 70 at the Litchfield 1 nest. Stable isotopes of oxygen and hydrogen indicate that groundwater throughout all five till profiles is younger than the last glacial advance into Minnesota at about 11,000 years ago. The heuristic modeling demonstrated that, for understanding sustainability of groundwater pumping from confined aquifers, knowledge of till hydraulic properties is just as important as knowledge of aquifer hydraulic properties. Substantial differences in groundwater fluxes into and through till were observed across hydrogeologic settings representative of the field sites. Over long periods of time (hundreds of years), pumping-induced hydraulic gradients are established in confined aquifer systems and, even in low hydraulic conductivity tills, these pumping-induced hydraulic gradients increase leakage into and through till compared to ambient conditions. In conclusion, groundwater flowing vertically downward through till confining units (leakage) replenishes water pumped from confined aquifers. Till hydraulic properties, such as those presented in this report, provide important information that can be used to quantify leakage rates through till. Till hydraulic properties are variable over short distances and profoundly affect leakage rates, demonstrating the importance of site-specific till hydraulic data for evaluating the sustainability of groundwater withdrawals from confined aquifers.

Minnesota↗