Search USGS⌕ Search

SEARCH · Search USGS

Results for “Scientific Investigations Map”

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

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

At least 1,351 records · Page 75Linked to original sources

Bend-scale geomorphic classification and assessment of the Lower Missouri River from Sioux City, Iowa, to the Mississippi River for application to pallid sturgeon management

Management actions intended to increase growth and survival of pallid sturgeon ( Scaphirhynchus albus ) age-0 larvae on the Lower Missouri River require a comprehensive understanding of the geomorphic habitat template of the river. The study described here had two objectives relating to where channel-reconfiguration projects should be located to optimize effectiveness. The first objective was to develop a bend-scale (that is, at the scale of individual bends, defined as “cross-over to cross-over”) geomorphic classification of the Lower Missouri River to help in the design of monitoring and evaluation of such projects. The second objective was to explore whether geomorphic variables could provide insight into varying capacities of bends to intercept drifting larvae. The bend-scale classification was based on geomorphic and engineering variables for 257 bends from Sioux City, Iowa, to the confluence with the Mississippi River near St. Louis, Missouri. We used k-means clustering to identify groupings of bends that shared the same characteristics. Separate 3-, 4-, and 6-cluster classifications were developed and mapped. The three classifications are nested in a hierarchical structure. We also explored capacities of bends to intercept larvae through evaluation of linear models that predicted persistent sand area or catch per unit effort (CPUE) of age-0 sturgeon as a function of the same geomorphic variables used in the classification. All highly ranked models that predict persistent sand area contained mean channel width and standard deviation of channel width as significant variables. Some top-ranked models also included contributions of channel sinuosity and density of navigation structures. The sand-area prediction models have r-squared values of 0.648–0.674. In contrast, the highest-ranking CPUE models have r-squared values of 0.011–0.170, indicating much more uncertainty for the biological response variable. Whereas the persistent sand model documents that physical processes of transport and accumulation are systematic and predictable, the poor performance of the CPUE models indicate that additional processes will need to be considered to predict biological transport and accumulation.

Iowa, Missouri↗

Statewide assessment of karst aquifers in New York with an inventory of closed-depression and focused-recharge features

Karst is a landscape formed from the dissolution of soluble rock or rock containing minerals that are easily dissolved from within the rock. The landscape is characterized by sinkholes, caves, losing streams, springs, and underground drainage systems, which rapidly move water through the karst. The two forms of karst in New York State include carbonate karst, which forms in carbonate rock (limestone, marble, and dolostone), and evaporite karst, which forms in rock that contains the evaporite minerals gypsum and halite. Past and recent studies of karst across the State have shown that areas of focused recharge in karstic carbonate rock allow contaminants to enter aquifer systems with little attenuation. Focused areas of recharge need to be identified to help prevent such contamination from sources on or adjacent to the karst. The New York State Departments of Environmental Conservation and Health are collaborating with the agricultural community to make farmers and farm-planning advisors more aware of karst and how to manage daily farming activities to reduce their impact on surface water and groundwater resources, especially in karst areas. There is also a need to make regulators, planners, and the general public aware of New York’s karst resources and to properly protect and manage these resources to protect the quality of groundwater and surface water that can flow into, through, and from karst bedrock. Using publicly available geospatial data, karst bedrock and closed depressions over or near karst rock were identified across New York. Carbonate, evaporite, and marble geologic units were selected from a statewide 1:250,000-scale bedrock geology dataset. The selected geologic units were intersected with 7.5-minute quadrangle maps to define the study area. The U.S. Geological Survey has compiled an inventory of closed depressions from statewide digital contour data, scanned 7.5-minute topographic maps known as a digital raster graphics, and light detection and ranging (lidar) digital elevation models. Analysis of the data resulted in the identification of 5,023 closed depressions statewide. The inventory was conducted to eliminate duplication of results from analysis of the three data sources. A series of overlay analyses was conducted using the closed depressions and thematic data known to be key factors in determining the probability of a closed depression contributing to focused groundwater recharge; the thematic data include bedrock geology, soil type, soil infiltration rate, and land cover. Though the extent of karst development is important in understanding the interaction between surface water and groundwater in karst terrains, some of the worst cases of groundwater contamination in karst can occur where only minor karst features might be present. The presence of karst—be it a short section of a solutioned fracture or an extensive cave system—requires careful consideration, forward-looking environmental planning, and consistent water-quality protection to preserve New York State’s water resources.

New York↗

Mapping karst groundwater flow paths and delineating recharge areas for springs in the Little Sequatchie and Pryor Cove watersheds, Tennessee

The Little Sequatchie River and Pryor Cove Branch, in southern Tennessee, drain the eastern escarpment of the Cumberland Plateau to the Sequatchie River near the southern end of the Sequatchie Valley. The Little Sequatchie River is the largest tributary to the Sequatchie River by drainage area, covering over 120 square miles. The hydrology of the two drainage areas has been largely altered by karst processes, which has caused the majority of the streams to sink at the contact between the Mississippian Pennington Formation and the underlying Mississippian Bangor Limestone. The U.S. Geological Survey, in cooperation with the U.S. Fish and Wildlife Service and Tennessee Department of Environment and Conservation, initiated a study in 2021 to map the karst groundwater pathways in both watersheds in order to delineate recharge areas for several springs. One of these springs, Sequatchie Cave, represents a significant habitat for two Species of Greatest Conservation Need, the Glyphopsyche sequatchie (Sequatchie caddisfly) and the federally endangered Marstonia ogmorhaphe (royal marstonia). Springs and springflow-dominated streams in the Little Sequatchie River valley and Pryor Cove also provide water for agricultural practices and serve as a drinking water source for nearby communities. During the study, a total of 25 dye injections were conducted over eight rounds from January 2022 through March 2023. Dye traces from these injections helped to delineate recharge areas for six major springs, ranging from 7.3 to 65.2 square miles in area. The majority of the dye traces remained subsurface (from sinkpoint to recovery site) for long distances, with karst groundwater travelling nearly 8 miles before resurfacing. The dye traces also had rapid traveltimes, often travelling hundreds to thousands of feet per hour. The goal of this project was to provide scientific data related to karst groundwater pathways and spring recharge areas to aid State and Federal agencies in making informed decisions to protect and preserve this unique and vulnerable karst system.

Tennessee↗

An update of the Death Valley regional groundwater flow system transient model, Nevada and California

Since the original publication of the Death Valley regional groundwater flow system (DVRFS) numerical model in 2004, more information on the regional groundwater flow system in the form of new data and interpretations has been compiled. Cooperators such as the Bureau of Land Management, National Park Service, U.S. Fish and Wildlife Service, the Department of Energy, and Nye County, Nevada, recognized a need to update the existing regional numerical model to maintain its viability as a groundwater management tool for regional stakeholders. The existing DVRFS numerical flow model was converted to MODFLOW-2005, updated with the latest available data, and recalibrated. Five main data sets were revised: (1) recharge from precipitation varying in time and space, (2) pumping data, (3) water-level observations, (4) an updated regional potentiometric map, and (5) a revision to the digital hydrogeologic framework model. The resulting DVRFS version 2.0 (v. 2.0) numerical flow model simulates groundwater flow conditions for the Death Valley region from 1913 to 2003 to correspond to the time frame for the most recently published (2008) water-use data. The DVRFS v 2.0 model was calibrated by using the Tikhonov regularization functionality in the parameter estimation and predictive uncertainty software PEST. In order to assess the accuracy of the numerical flow model in simulating regional flow, the fit of simulated to target values (consisting of hydraulic heads and flows, including evapotranspiration and spring discharge, flow across the model boundary, and interbasin flow; the regional water budget; values of parameter estimates; and sensitivities) was evaluated. This evaluation showed that DVRFS v. 2.0 simulates conditions similar to DVRFS v. 1.0. Comparisons of the target values with simulated values also indicate that they match reasonably well and in some cases (boundary flows and discharge) significantly better than in DVRFS v. 1.0.

California, Nevada↗

Grammar to graph—An approach for semantic transformation of annotations to triples

Data annotation is the process of labeling data to show the outcome that a related data model should predict. In this study, annotation data were transformed into semantic graph triples, mainly for use with the Resource Description Framework (RDF), a type of entity-relationship-attribute data model for graph databases. The transformation of annotation data to semantic graph triples provides complex linguistic meaning with data handling advantages such as reduced data storage needs, improved logical specification of relations between objects, and reusable classes and properties that support logic and inference. A grammar-based framework in graph form supports user questions and queries. The words defining approximately 334 topographic feature types compiled by the U.S. Geological Survey were tokenized as units of analysis and grouped by part of speech. Their dependency relations were identified for this study using natural language processing libraries. Dependency concepts are used as structured semantic relations among part-of-speech classes. Tokens, units equivalent to words, form instances of classes and were quantified within a tabular output format using PostgreSQL data storage software. Table data were logically aligned as triples following a mapping file and stored with an ontology file using Ontop virtual triplestore software. A grammar ontology schema for the data was synchronized to match queries whose results validated the graph’s structure. The text analysis produced 8 part-of-speech classes of content words for object representations and 4 classes of function words for operational applications. Dependency relations formed 27 ontology properties for topographic subgraph structures. Token occurrences shaped overall ontology salience and formed a lexicon of syntactic terms for subgraph objects and properties. The schema ontology of class and property population shapes formed the lexicon of English terms. SPARQL Protocol and RDF Query Language (SPARQL) was used with the lexicon to conform data to RDF guidelines. This study confirms the hypothesis that although linguistic logic varies from description logic, its approximation applies to ontology design. Property and query use case patterns extracted from the analysis support queries concerning complex topographic relations and patterns normally embedded within text definitions. The method used in this study could be applied to text forms in other domains, such as survey notes.

Scientific Investigations Report↗

Hydrogeologic conceptualization of the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory, Edwards Air Force Base, Antelope Valley, California, 2017–21

Research and testing within the Air Force Research Laboratory at Edwards Air Force Base in the southwestern Mojave Desert has been conducted since the 1950s. This testing requires large quantities of groundwater for cooling and operations involving hazardous constituents that have entered the groundwater. To better understand the potential movement of contaminants through the groundwater system, the U.S. Geological Survey entered into a cooperative agreement with the U.S. Air Force Civil Engineer Center with the purpose of (1) gaining a more detailed understanding of the hydraulic connection between the unconsolidated basin-fill sediment and the underlying weathered and crystalline bedrock and (2) investigating faults and their potential effect on groundwater flow. This study was specific to the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory. The refined hydrogeologic characterization was derived by evaluating surface geologic maps, a regional gravity model, new subsurface geophysical surveys, lithology data from boreholes, and groundwater-level data. Refinement of the hydrogeologic characterization of the study area involved developing a map of basin-fill sediment thickness, estimating the geometry and depth to bedrock, determining the presence of weathered bedrock material, and identifying the location of faults and subsurface structures. The thickness of basin-fill sediment was derived from a regional gravity model. Basin-fill sediment ranges from less than 20 feet (ft) thick to about 150 ft thick and increases substantially toward a deep structural basin in the northwestern part of the study area. The top of bedrock estimated from the new surface geophysical surveys is generally deeper than the top of bedrock from the regional gravity model. These differences in interpreted bedrock tops indicate that there may be a zone of weathered bedrock in the subsurface, overlying crystalline bedrock. Additional subsurface data evaluating the geometry and thickness of weathered bedrock are necessary for a comprehensive evaluation of weathered bedrock in the study area. Previously unmapped extensions of the Leuhman Fault, Spring Fault Zone, and Boron Gate Fault were identified in the new surface geophysical surveys. Evaluation of groundwater flow between bedrock and basin-fill sediment, and the effects of faults on groundwater flow, was done using groundwater-level data from 2020. The groundwater table was in crystalline bedrock throughout most of the Arroyos and Northeast Air Force Research Laboratory subareas, with a few exceptions where it was in basin-fill sediment and weathered bedrock. Basin-fill sediment, and likely weathered bedrock, was partially saturated in parts of the study area where basin-fill sediment was as thick as 150 ft. Weathered bedrock material may be present near the Spring Fault Zone, and if so, was likely saturated in 2020. The groundwater-level data near and within Spring Fault Zone indicate that the Spring Fault Zone acts as a vertical-along-strike conduit for, and as an across-strike impediment to, groundwater flow. In addition, steep groundwater-level gradients in the southeastern part of the study area indicate that a subsurface structure affecting groundwater flow could be related to the Leuhman Fault. Sufficient groundwater-level data were not available to conclusively determine if the Arroyos, Rich, and Boron Gate Faults, and their inferred extensions from the geophysical surveys, affect the movement or direction of groundwater flow.

California↗

Geologic resource evaluation of Pu‘uhonua O Hōnaunau National Historical Park, Hawai‘i. Part II: Benthic habitat mapping

In cooperation with the U.S. National Park Service (NPS), the U.S. Geological Survey (USGS) has mapped the underwater environment in and adjacent to three parks along the Kona coast on the island of Hawai‘i. This report is the second of two produced for the NPS on the geologic resource evaluation of of Pu‘uhonua O Hōnaunau National Historical Park (PUHO) and presents the results of benthic habitat mapping of the offshore waters for PUHO. See Part I (Richmond and others, 2006) for an overview of the regional geology, local volcanics, and a detailed description of coastal landforms in the park. Pu‘uhonua O Hōnaunau National Historical Park boundaries extend only to the mean high tide line and do not officially include the marine environment. However, impacts downslope of any development in the park are of concern to management. The area mapped for this report extends from Hōnaunau Bay, around Pu‘uhonua Point, to Ki‘ilae Bay and the south park boundary and from the shoreline to depths of approximately 40 m (130 ft), where the shelf drops off to a sand-covered bottom. Pu‘uhonua O Hōnaunau National Historical Park lies within the nearly 850-km2 Ki‘ilae watershed, which begins at the crest of Mauna Loa. The Ki‘ilae Watercourse runs through the southern area of the park and empties into Ki‘ilae Bay, but only during periods of extreme rainfall. The waters of Keone‘ele Cove, the ancient royal canoe landing at PUHO, while not formally under NPS jurisdiction, are managed by the park under an agreement with the State of Hawaii. This small embayment is a known haven for sea turtles, which are often found sunning themselves on the near- shore volcanic platform. Impacts to this area include frequent visits by scuba divers and snorkelers to Hōnaunau Bay and a small boat ramp located just to the north of Keone‘ele Cove.

Hawai'i↗

Evaluation of stream capture related to groundwater pumping, Lower Humboldt River Basin, Nevada

The Humboldt River Basin is the only river basin that is contained entirely within the State of Nevada. The effect of groundwater pumping on the Humboldt River is not well understood. Tools are needed to determine stream capture and manage groundwater pumping in the Humboldt River Basin. The objective of this study is to estimate capture and storage change caused by groundwater withdrawals in the lower Humboldt River Basin that can provide the Nevada State Engineer with data and information needed to manage groundwater and surface-water resources. A numerical groundwater flow model was developed for the purpose of estimating stream capture from pre-2016 and future pumping as well as for any location of potential future pumping within the lower Humboldt River Basin. This model was developed using MODFLOW-NWT to represent the lower Humboldt River Basin hydrologic system, including Humboldt River; Rye Patch Reservoir; groundwater evapotranspiration; pumping from municipal, agricultural, mining, and domestic wells; and agricultural drains. Aquifer properties were calibrated using results from numerous single- and multi-well aquifer tests (Nadler, 2020) and through the process of model calibration. Historical capture was estimated for 1960–2016 and predictive capture for the system was projected 100 years into the future (2017–2116) based on historical pumping patterns. Stream capture and drain capture are relatively low for the historical and predictive periods. During the historical period, increased pumping during dry years caused increased connections with capture sources and less water sourced to wells from aquifer storage. Storage and groundwater levels generally recovered during subsequent wet years. Overall, storage change has been the main source of water to wells in the lower Humboldt River Basin, followed by groundwater evapotranspiration capture. During the predictive period, pumping is projected to remain constant and capture 9 percent of stream water after 100 years. Capture and storage change maps were created to visualize spatial variability in potential capture and storage change through time and to provide a database of results that can be used to manage groundwater and surface-water resources. These maps show that potential stream capture would be a minor source of water to wells located across most of the simulated area, except for locations close to the Humboldt River and Rye Patch Reservoir. Drains also would be a minor potential source of water to wells except for those directly adjacent to the drains. In general, the potential supply of water to wells is storage-dominated and over time groundwater evapotranspiration-dominated in the agricultural area. Capture difference maps were generated to visualize where potential capture results might have greater limitations associated with nonlinear flow processes, such as head-dependent boundary conditions. Higher capture differences indicate larger capture map bias and therefore greater capture map uncertainty due to the inability of capture maps to account for nonlinear flow processes. Stream capture differences are highest directly adjacent to the river but are otherwise minimal. Drain capture differences are highest in the region of the agricultural drain network but are otherwise minimal. The Humboldt River, Rye Patch Reservoir, and drains introduce very little nonlinearity to the model, and their associated capture map bias is minimal. Potential groundwater evapotranspiration capture introduces a fair amount of nonlinearity to the model and has the potential to result in significant, localized groundwater evapotranspiration capture map bias over time. Groundwater evapotranspiration capture differences are the result of higher pumping rates lowering the water table below the root zone faster than lower pumping rates and essentially removing groundwater evapotranspiration as a potential source of capture faster than lower pumping rates. Wells that can no longer source their supply through groundwater evapotranspiration capture then generally source more of their water from storage. Thus, storage change bias increases over time as well. Capture prediction uncertainty due to parameter estimation was evaluated using a covariance matrix adaptation-evolution strategy. One hundred Monte Carlo realizations of model parameters were applied to the model to assess capture uncertainty at 13 grid cell locations within the model domain. In general, results indicated that greater capture uncertainty for a given source (river, drains, or evapotranspiration) is associated with proximity of a pumping well to that source. The magnitude of maximum capture fraction uncertainties after 100 years of pumping for stream capture, drain capture, groundwater evapotranspiration capture, and storage change were plus or minus (±) 0.17, ±0.10, ±0.20, and ±0.22, respectively.

Nevada↗

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↗

Hydrogeologic investigations of the Miocene Nogales Formation in the Nogales Area, Upper Santa Cruz Basin, Arizona

Hydrogeologic investigations were conducted to evaluate the groundwater resource potential for the Miocene Nogales Formation in the Nogales area, southern Arizona. Results indicate that parts of the formation may provide new, deeper sources of groundwater for the area. Geologic mapping determined the hydrogeologic framework of the formation by defining lithologic, mineralogic, and stratigraphic characteristics; identifying potential aquifers and confining units; and mapping faults and fractures which likely influence groundwater flow. Geophysical modeling was used to determine the basin geometry and thickness of the Nogales Formation and younger alluvial aquifers and to identify target areas (deep subbasins) which may prove to be productive aquifers. Volcaniclastic sandstone samples from the formation were analyzed for porosity, bulk density, saturated hydraulic conductivity, and fabric. Effective porosity ranges from 16 to 42 percent, bulk density from 1.6 to 2.47 grams per cubic centimeter, and saturated hydraulic conductivity (SHC) from 4 to 57 centimeters per day (4.9×10 -5 to 6.7×10 -4 centimeters per second). Thin sections show that sandstone framework grains consist of quartz, feldspar, biotite, hornblende, pumice, volcanic glass, and opaque minerals. The matrix in most samples consists of pumice fragments, and some contain predominantly silt and clay. Samples with a mostly silt and clay matrix have lower porosity and SHC compared to samples with mostly pumice, which have higher and wider ranges of porosity and SHC. Pore space in the Nogales Formation sediments includes moldic, intercrystalline, and fracture porosity. Some intercrystalline pore space is partially filled with calcite cement. About one third of the samples contain fractures, which correspond to fractures noted in outcrops in all members of the formation. Scanning electron microscope (SEM) and x-ray diffraction (XRD) analyses indicate that most of the samples contained the zeolite clinoptilolite and mixed-layer clay. X-ray diffraction analyses verified clinoptilolite as the only zeolite in Nogales Formation samples; they also verified the presence of smectite and illite clay and some kaolinite. Samples which contain greater amounts of clinoptilolite and lesser amounts of smectite have high porosity and SHC in narrow ranges. However, samples with abundant smectite and lesser amounts of clinoptilolite span the entire ranges of porosity and SHC for the formation. All members of the Nogales Formation are fractured and faulted as a result of Tertiary Basin and Range extensional deformation, which was broadly contemporaneous with deposition of the formation. These structures may have significant influence on groundwater flow in the upper Santa Cruz basin because, although many of the sediments in the formation have characteristics indicating they may be productive aquifers based only on porous-media flow, fracturing in these sediments may further enhance permeability and groundwater flow in these basin-fill aquifers by orders of magnitude.

Arizona↗

Sediment distribution and hydrologic conditions of the Potomac aquifer in Virginia and parts of Maryland and North Carolina

Sediments of the heavily used Potomac aquifer broadly contrast across major structural features of the Atlantic Coastal Plain Physiographic Province in eastern Virginia and adjacent parts of Maryland and North Carolina. Thicknesses and relative dominance of the highly interbedded fluvial sediments vary regionally. Vertical intervals in boreholes of coarse-grained sediment commonly targeted for completion of water-supply wells are thickest and most widespread across the central and southern parts of the Virginia Coastal Plain. Designated as the Norfolk arch depositional subarea, the entire sediment thickness here functions hydraulically as a single interconnected aquifer. By contrast, coarse-grained sediment intervals are thinner and less widespread across the northern part of the Virginia Coastal Plain and into southern Maryland, designated as the Salisbury embayment depositional subarea. Fine-grained intervals that are generally avoided for completion of water-supply wells are increasingly thick and widespread northward. Fine-grained intervals collectively as thick as several hundred feet comprise two continuous confining units that hydraulically separate three vertically spaced subaquifers. The subaquifers are continuous northward but merge southward into the single undivided Potomac aquifer. Lastly, far southeastern Virginia and northeastern North Carolina are designated as the Albemarle embayment depositional subarea, where both coarse- and fine-grained intervals are of only moderate thickness. The entire sediment thickness functions hydraulically as a single interconnected aquifer. A substantial hydrologic separation from overlying aquifers is imposed by the upper Cenomanian confining unit. Potomac aquifer sediments were deposited by a fluvial depositional complex spanning the Virginia Coastal Plain approximately 100 to 145 million years ago. Westward, persistently uplifted granite and gneiss source rocks sustained a supply of coarse-grained sand and gravel. Immature, high-gradient braided streams deposited longitudinal bars and channel fills across the Norfolk arch subarea. By contrast, across the Salisbury and Albemarle embayment subareas, mature, medium- to low-gradient meandering streams deposited medium- to coarse-grained channel fills and point bars segregated from fine-grained overbank deposits. The Virginia depositional complex merged northward across the Salisbury embayment subarea with another complex in Maryland. Here, additional sediments were received from schist source rocks that underwent three cycles of initial uplift and rapid erosion followed by crustal stability and erosional leveling. Because of the predominance of coarse-grained sediments, transmissivity, hydraulic conductivity, and regional velocities of lateral flow through the Potomac aquifer are greatest across the Norfolk arch depositional subarea, but decrease progressively northward with increasingly fine-grained sediments. Confining units hydraulically separate the Potomac aquifer from overlying aquifers, as indicated by large vertical hydraulic gradients. By contrast, most of the Potomac aquifer internally functions hydraulically as a single interconnected aquifer, as indicated by uniformly small vertical gradients. Most fine-grained sediments within the aquifer do not hydraulically separate overlying and underlying coarse-grained sediments. Across the Salisbury embayment depositional subarea, however, hydraulic separation among the vertically spaced subaquifers is imposed by the intervening confining units. The Potomac aquifer is the largest and most heavily used source of groundwater in the Virginia Coastal Plain. Water-level declines as great as 200 feet create the potential for saltwater intrusion. Conventional stratigraphic correlation has been generally ineffective at accurately characterizing complexly distributed fluvial sediments that compose the Potomac aquifer. Consequently, the aquifer’s internal hydraulic connectivity and overall hydrologic function have not been well understood. Water-supply planning and development efforts have been hampered, and interpretations of regulatory criteria for allowable water-level declines have been ambiguous. An investigation undertaken during 2010–11 by the U.S. Geological Survey, in cooperation with the Virginia Department of Environmental Quality, provides a comprehensive regional description of the spatial distribution of Potomac aquifer sediments and their relation to hydrologic conditions. Altitudes and thicknesses of 2,725 vertical sediment intervals represent the spatial distribution of Potomac aquifer sediments in the Virginia Coastal Plain and adjacent parts of Maryland and North Carolina. Sediment intervals are designated as either dominantly coarse or fine grained and were determined by interpretation of geophysical logs and ancillary information from 456 boreholes. Sediment-interval and borehole summary statistical data indicate regional trends in sediment lithology and stratigraphic continuity, upon which three structurally based and hydrologically distinct sediment depositional subareas are designated. Broad patterns of sediment deposition over time are inferred from published sediment pollen-age data. Discrepancies in previously drawn hydrostratigraphic relations between southeastern Virginia and northeastern North Carolina are partly resolved based on borehole geophysical logs and a recently documented geologic map and corehole. A conceptual model theorizes the depositional history of the sediments and geologically accounts for their distribution. Documented pumping tests of the Potomac aquifer at 197 locations produced 336 values of transmissivity and 127 values of storativity. Based on effective aquifer thicknesses, 296 values of sediment hydraulic conductivity and 113 values of sediment specific storage are calculated. Vertical hydraulic gradients are calculated from 9,479 pairs of water levels measured between November 17, 1953, and October 4, 2011, in 129 closely spaced pairs of wells. Borehole sediment-interval and related data provide a means to achieve high yielding production wells in the Potomac aquifer by site-specific targeting of drilling operations toward water-bearing coarse-grained sand and gravel. Advance knowledge of the potential of different parts of the aquifer also aids in planning optimal groundwater-development areas. Depositional subareas further provide a possible context for resource management. Current (2013) regulatory limits on water-level declines are relative to top surfaces of subdivided upper, middle, and lower Potomac aquifers across the entire Virginia Coastal Plain, but have the potential to exceed the same limit relative to a single undivided Potomac aquifer. By contrast, designation of the sediments as a single aquifer in the Norfolk arch and Albemarle embayment subareas—and as a series of vertically spaced subaquifers and intervening confining units in the Salisbury embayment subarea—best reflects understanding of the Potomac aquifer and can avoid the potential for excessive water-level declines. Simulation modeling to evaluate effects of groundwater withdrawals could be designed similarly, including vertical discretization and (or) zonation of the Potomac aquifer based on depositional subareas and a geostatistical distribution of aquifer properties derived from borehole sediment-interval data. Further resource-management information needs extend beyond the developed part of the Potomac aquifer, particularly across the Northern Neck and Middle Peninsula where only the shallowest part of the aquifer is known, and include structural aspects such as faults, basement bedrock, and the Chesapeake Bay impact crater.

Maryland, North Carolina, Virginia↗

Three-Dimensional Geologic Model of Complex Fault Structures in the Upper Seco Creek Area, Medina and Uvalde Counties, South-Central Texas

This multimedia report shows and describes digital three-dimensional faulted geologic surfaces and volumes of the lithologic units of the Edwards aquifer in the upper Seco Creek area of Medina and Uvalde Counties in south-central Texas. This geologic framework model was produced using (1) geologic maps and interpretations of depositional environments and paleogeography; (2) lithologic descriptions, interpretations, and geophysical logs from 31 drill holes; (3) rock core and detailed lithologic descriptions from one drill hole; (4) helicopter electromagnetic geophysical data; and (5) known major and minor faults in the study area. These faults were used because of their individual and collective effects on the continuity of the aquifer-forming units in the Edwards Group. Data and information were compared and validated with each other and reflect the complex relationships of structures in the Seco Creek area of the Balcones fault zone. This geologic framework model can be used as a tool to visually explore and study geologic structures within the Seco Creek area of the Balcones fault zone and to show the connectivity of hydrologic units of high and low permeability between and across faults. The software can be used to display other data and information, such as drill-hole data, on this geologic framework model in three-dimensional space.

Scientific Investigations Report↗

Apparent Resistivity and Estimated Interaction Potential of Surface Water and Groundwater along Selected Canals and Streams in the Elkhorn-Loup Model Study Area, North-Central Nebraska, 2006-07

In 2005, the State of Nebraska adopted new legislation that in part requires local Natural Resources Districts to include the effect of groundwater use on surface-water systems in their groundwater management plan. In response the U.S. Geological Survey, in cooperation with the Upper Elkhorn, Lower Elkhorn, Upper Loup, Lower Loup, Middle Niobrara, Lower Niobrara, Lewis and Clark, and Lower Platte North Natural Resources Districts, did a study during 2006-07 to investigate the surface-water and groundwater interaction within a 79,800-square-kilometer area in north-central Nebraska. To determine how streambed materials affect surface-water and groundwater interaction, surface geophysical and lithologic data were integrated at four sites to characterize the hydrogeologic conditions within the study area. Frequency-domain electromagnetic and waterborne direct- current resistivity profiles were collected to map the near-surface hydrogeologic conditions along sections of Ainsworth Canal near Ainsworth, Nebraska; Mirdan and Geranium Canals near Ord, Nebraska; North Loup River near Ord, Nebraska; and Middle Loup River near Thedford, Nebraska. Lithologic data were collected from test holes at each site to aid interpretation of the geophysical data. Geostatistical analysis incorporating the spatial variability of resistivity was used to account for the effect of lithologic heterogeneity on effective hydraulic permeability. The geostatistical analysis and lithologic data descriptions were used to make an interpretation of the hydrogeologic system and derive estimates of surface-water/groundwater interaction potential within the canals and streambeds. The estimated interaction potential at the Ainsworth Canal site and the Mirdan and Geranium Canal site is generally low to moderately low. The sediment textures at nearby test holes typically were silt and clay and fine-to-medium sand. The apparent resistivity values for these sites ranged from 2 to 120 ohm-meters. The vertical and horizontal variability of the apparent resistivity data were consistently low. Low resistive variability indicates little lithologic heterogeneity for either canal site. The surface-water/groundwater interaction-potential estimates are in agreement with the narrow frequency distribution of resistivity, low apparent resistivities, low spatial heterogeneity, and test-hole grain-size ranges. The estimated surface-water/groundwater interaction potential at the North Loup and Middle Loup River sites is moderate to moderately high. The sediment textures at nearby test holes were predominantly fine, medium, and coarse sand with some silt and silty to sandy clay. The apparent resistivity values for these sites ranged from 34 to 1,338 ohm-meters. The vertical variability of the resistivity data was moderately high. The horizontal variability at these sites is low to moderately low. The higher resistive variability at these sites indicates generally greater lithologic heterogeneity than at either the Ainsworth Canal site or the Mirdan and Geranium Canal site. The surface-water/groundwater interaction-potential estimates are in agreement with the generally moderate to high apparent resistivity, the greater spatial heterogeneity, and the variable lithologic texture. A higher interaction potential as compared to the canal sites is expected because of the higher subsurface resistivity and greater lithologic heterogeneity.

Scientific Investigations Report↗

Evaluation of a Single-Beam Sonar System to Map Seagrass at Two Sites in Northern Puget Sound, Washington

Seagrass at two sites in northern Puget Sound, Possession Point and nearby Browns Bay, was mapped using both a single-beam sonar and underwater video camera. The acoustic and underwater video data were compared to evaluate the accuracy of acoustic estimates of seagrass cover. The accuracy of the acoustic method was calculated for three classifications of seagrass observed in underwater video: bare (no seagrass), patchy seagrass, and continuous seagrass. Acoustic and underwater video methods agreed in 92 percent and 74 percent of observations made in bare and continuous areas, respectively. However, in patchy seagrass, the agreement between acoustic and underwater video was poor (43 percent). The poor agreement between the two methods in areas with patchy seagrass is likely because the two instruments were not precisely colocated. The distribution of seagrass at the two sites differed both in overall percent vegetated and in the distribution of percent cover versus depth. On the basis of acoustic data, seagrass inhabited 0.29 km2 (19 percent of total area) at Possession Point and 0.043 km2 (5 percent of total area) at the Browns Bay study site. The depth distribution at the two sites was markedly different. Whereas the majority of seagrass at Possession Point occurred between -0.5 and -1.5 m MLLW, most seagrass at Browns Bay occurred at a greater depth, between -2.25 and -3.5 m MLLW. Further investigation of the anthropogenic and natural factors causing these differences in distribution is needed.

Scientific Investigations Report↗

Mercury in precipitation in Indiana, January 2004–December 2005

Mercury in precipitation was monitored during 2004–2005 at five locations in Indiana as part of the National Atmospheric Deposition Program–Mercury Deposition Network (NADP–MDN). Monitoring stations were operated at Roush Lake near Huntington, Clifty Falls State Park near Madison, Fort Harrison State Park near Indianapolis, Monroe County Regional Airport near Bloomington, and Indiana Dunes National Lakeshore near Porter. At these monitoring stations, precipitation amounts were measured continuously and weekly samples were collected for analysis of mercury by methods achieving detection limits as low as 0.05 ng/L (nanograms per liter). Wet deposition was computed as the product of mercury concentration and precipitation. The data were analyzed for seasonal patterns, temporal trends, and geographic differences. In the 2 years, 520 weekly samples were collected at the 5 monitoring stations and 448 of these samples had sufficient precipitation to compute mercury wet deposition. The 2-year mean mercury concentration at the five monitoring stations (normalized to the sample volume) was 10.6 ng/L. As a reference for comparison, the total mercury concentration in 41 percent of the samples analyzed was greater than the statewide Indiana water-quality standard for mercury (12 ng/L, protecting aquatic life) and 99 percent of the concentrations exceeded the most conservative Indiana water-quality criterion (1.3 ng/L, protecting wild mammals and birds). The normalized annual mercury concentration at Clifty Falls in 2004 was the fourth highest in the NADP–MDN in eastern North America that year. In 2005, the mercury concentrations at Clifty Falls and Indiana Dunes were the ninth highest in the NADP–MDN in eastern North America. At the five monitoring stations during the study period, the mean weekly total mercury deposition was 0.208 µg/m 2 (micrograms per square meter) and mean annual total mercury deposition was 10.8 µg/m 2 . The annual mercury deposition at Clifty Falls in 2004 and 2005 was in the top 25 percent of the NADP–MDN stations in eastern North America. Mercury concentrations and deposition varied at the five monitoring stations during 2004–2005. Mercury concentrations in wet-deposition samples ranged from 1.2 to 116.6 ng/L and weekly mercury deposition ranged from 0.002 to 1.74 µg/m 2 . Data from weekly samples exhibited seasonal patterns. During April through September, total mercury concentrations and deposition were higher than the median for all samples. Annual precipitation at four of the five monitoring stations was within 10 percent of normal both years, with the exception of Indiana Dunes, where precipitation was 23 percent below normal in 2005. Episodes of high mercury deposition, which were the top 10 percent of weekly mercury deposition at the five monitoring stations, contributed 39 percent of all mercury deposition during 2004–2005. Mercury deposition more than 1.04 µg/m 2 (5 times the mean weekly deposition) was recorded for 12 samples. These episodes of highest mercury deposition were recorded at all five monitoring stations, but the most (7 of 12) were at Clifty Falls and contributed 34.4 percent of the total deposition at that station during 2004–2005. Weekly samples with high mercury deposition may help to explain the differences in annual mercury deposition among the five monitoring stations in Indiana. A statistical evaluation of the monitoring data for 2001–2005 indicated several statistically significant temporal trends. A statewide (5-station) decrease (p = 0.007) in mercury deposition and a statewide decrease (p = 0.059) in mercury concentration were shown. Decreases in mercury deposition (p = 0.061 and p = 0.083) were observed at Roush Lake and Bloomington. A statistically significant trend was not observed for precipitation at the five monitoring stations during this 5-year period. A potential explanation for part of the statewide decrease in mercury concentration and mercury deposition was a 28 percent decrease in the total estimated annual mercury emissions in Indiana between 2002 and 2005. Mercury deposition statistically was correlated most closely to precipitation in the 448 samples, 2004–2005, and this relation was demonstrated by statewide maps of annual precipitation and annual mercury deposition based on precipitation data from 127 National Weather Service Cooperative Observer Program stations. However, one area in southeastern Indiana in the vicinity of Clifty Falls exhibited high mercury deposition that might be related more to mercury concentration than to precipitation. This is because areas with the same range of precipitation as southeastern Indiana were mapped with less mercury deposition. Other data demonstrate a geographic difference for mercury in precipitation in the vicinity of the Clifty Falls monitoring station. The weekly mercury concentrations at Clifty Falls were statistically higher than concentrations at Roush Lake, Fort Harrison, and Bloomington. Clifty Falls data ranked highest among the five monitoring stations for mercury concentration and mercury deposition, 2004–2005, and in the previous 3 years. Episodes of high mercury deposition were recorded most often at Clifty Falls in 2004–2005 and in the previous 3 years. Statistical trends in mercury concentration or mercury deposition were not observed for the Clifty Falls data. A potential explanation for this geographic difference is that annual mercury emissions from sources in the vicinity of Clifty Falls were higher than those at the other stations. Other factors may help explain the differences in total mercury concentrations, such as the types of mercury emissions, mercury transport from stationary sources outside Indiana, and meteorological conditions. Additional data are needed to assign a localized or regional boundary to the area affected by high deposition of mercury near Clifty Falls.

Indiana↗

Distribution, persistence, and hydrologic characteristics of salmon spawning habitats in clearwater side channels of the Matanuska River, southcentral Alaska

Turbid, glacially influenced rivers are often considered to be poor salmon spawning and rearing habitats and, consequently, little is known about salmon habitats that do occur within rivers of this type. To better understand salmon spawning habitats in the Matanuska River of southcentral Alaska, the distribution and characteristics of clearwater side-channel spawning habitats were determined and compared to spawning habitats in tributaries. More than 100 kilometers of clearwater side channels within the braided mainstem of the Matanuska River were mapped for 2006 from aerial images and ground-based surveys. In reaches selected for historical analysis, side channel locations shifted appreciably between 1949 and 2006, but the relative abundance of clearwater side channels was fairly stable during the same period. Geospatial analysis of side channel distribution shows side channels typically positioned along abandoned bars at the braid plain margin rather than on bars between mainstem channels, and shows a strong correlation of channel abundance with braid plain width. Physical and geomorphic characteristics of the channel and chemical character of the water measured at 19 side channel sites, 6 tributary sites, 4 spring sites, and 5 mainstem channel sites showed conditions suitable for salmon spawning in side channels and tributaries, and a correlation of side channel characteristics with the respective tributary or groundwater source water. Autumn-through-spring monitoring of intergravel water temperatures adjacent to salmon redds (nests) in three side channels and two tributaries indicate adequate accumulated thermal units for incubation and emergence of salmon in side channels and relatively low accumulated thermal units in tributaries.

Alaska↗

Dam failure analysis for the Lago de Matrullas Dam, Orocovis, Puerto Rico

The U.S. Geological Survey, in cooperation with the Puerto Rico Electric Power Authority, completed a hydrologic and hydraulic study to assess the potential hazard to human life and property associated with the hypothetical failure of the Lago de Matrullas Dam, located within the headwaters of the Río Grande de Manatí. The hydrologic study yielded outflow hydrographs and peak discharges for Lago de Matrullas and other subbasins in the Río Grande de Manatí hydrographic basin for three extreme rainfall events: (1) a 6-hour probable maximum precipitation (PMP) event, (2) a 24-hour PMP event, and (3) a 100-year-recurrence, 24-hour rainfall event. The hydraulic study simulated the hypothetical dam failure of Lago de Matrullas using hypothetical flood hydrographs generated from the hydrologic study and selected dam breach parameters. The flood wave resulting from the failure was downstream-routed through the lower reaches of the Río Matrullas, the Río Toro Negro, and the Río Grande de Manatí for determination of water-surface profiles developed from the event-based hydrologic scenarios and “sunny day” (no precipitation) conditions. The Hydrologic Modeling System (HEC–HMS) and the River Analysis System (HEC–RAS) computer programs, developed by the Hydrologic Engineering Center (HEC) of the U.S. Army Corps of Engineers, were used for the hydrologic and hydraulic modeling, respectively. The flow routing in the hydraulic analyses was performed using the unsteady-state flow module available in the HEC–RAS model. At the Lago de Matrullas Dam, inflow peak discharges of about 1,104 and 1,032 cubic meters per second (m 3 /s) were estimated with HEC–HMS for the 6- and 24-hour PMP events, respectively. The 100-year recurrence, 24-hour rainfall event simulation resulted in a peak discharge of about 418 m 3 /s. For the hydrologic analysis, no dam failure conditions were considered with the model. The results of the hydrologic simulations indicated, however, that the dam would be overtopped by the simulated 6- and 24-hour PMP events. Unlike the 6- and 24-hour PMP events, no overtopping was observed during the simulated 100-year recurrence, 24-hour rainfall event. For the dam-breach hydraulic analysis, the hypothetical failures of the Lago de Matrullas Dam included two possible failure modes: overtopping and piping. Overtopping failure was evaluated in this study for the 6- and 24- hour probable-maximum-precipitation breach scenarios. Piping dam failure was simulated for sunny day conditions and for the 100-year-recurrence, 24-hour rainfall scenario. Results from the simulated dam failure of the Lago de Matrullas Dam using the HEC–RAS model for the 6- and 24-hour PMP events showed peak discharges at the dam of 3,149.33 and 3,604.70 m 3 /s, respectively. Dam failure during the 100-year-recurrence, 24-hour rainfall event resulted in a peak discharge of 2,103.12 m 3 /s directly downstream from the dam. Dam failure under sunny day conditions produced a peak discharge of 1,695.91 m 3 /s at the dam assuming the antecedent lake level was at the morning-glory spillway invert elevation. Flood-inundation maps prepared as part of the study depict the flood extent and provide valuable information for preparing an Emergency Action Plan. Results of the failure analysis indicate that a failure of the Lago de Matrullas Dam could cause flooding to many of the inhabited areas along stream banks from the Lago de Matrullas Dam to the mouth of the Río Grande de Manatí. Among the areas most affected are the low-lying regions in the vicinity of the towns of Ciales, Manatí, and Barceloneta. The delineation of the flood boundaries near the town of Barceloneta considered the effects of a levee constructed during 2000 at Barceloneta in the flood plain of the Río Grande de Manatí to provide protection against flooding to the near-by low-lying populated areas. The results showed overtopping can be expected in the aforementioned levee during 6- and 24-hour probable-maximum-precipitation dam failure scenarios. No overtopping of the levee was simulated, however, during dam failure scenarios under the 100-year recurrence, 24-hour rainfall event or sunny day conditions.

Orocovis↗

Assessment of contaminant trends in plumes and wells and monitoring network optimization at the Badger Army Ammunition Plant, Sauk County, Wisconsin

Soil and groundwater at the Badger Army Ammunition Plant (BAAP), Sauk County, Wisconsin, were affected by several contaminants as a result of production and waste disposal practices common during its operation from 1942 to 1975. Three distinct plumes of contaminated groundwater originate on BAAP property and extend off-site, as identified by previous studies. Routine sampling of groundwater quality from a network of monitoring wells and off-site private wells has been performed since 1990, although the number of wells monitored and the monitoring frequency have varied as the approved monitoring plan was modified. During the period of monitoring from 1990 to 2018, numerous site investigations and remedial actions were conducted to address the sources of contamination, contaminated soils, and groundwater. Concentrations of contaminants reportedly decreased between 2000 and 2012 within all three plumes. Five or six contaminants of concern (COCs) were identified for each of the three plumes. An independent assessment of the contaminant plumes and of the monitoring network was conducted using groundwater-quality data collected from more than 600 wells between 2000 and 2018. In a study conducted by the U.S. Geological Survey (USGS), in cooperation with the Army Environmental Command, a consistent data aggregation and interpolation scheme was applied to derive the likely maximum groundwater plume extents in four 3-year time periods between 2000 and 2018. The plume extent was defined by the Enforcement Standard for each COC and represents the maximum concentration observed in each 3-year time period. The plume boundary analysis shows that the spatial extent of groundwater contamination decreased for most COCs during the study period. Some plume boundaries are not well delineated by the existing monitoring network, particularly the downgradient edge of the Propellant Burning Ground plume. Maps identify the plume boundary in each time period, the sampling well network used to delineate the plume, and wells that were sampled in the 2010–12 period but not sampled in the 2015–18 period. A series of statistical analyses using the Monitoring and Remediation Optimization System, version 3.0, program were applied to the available COC concentration data for two distinct periods, 2000 to 2012 and 2013 to 2018, with the break between periods coinciding with changes to the monitoring network in 2013. Trends in the concentration of COCs in individual wells varied, although generally more wells had decreasing than had increasing concentrations for most COCs in both time periods. The exceptions were ethyl ether in the 2004–12 period and 2,6-dinitrotoluene in the 2013–18 period, for which more wells had an increasing trend. Spatial moment analysis of concentration data from the well network was used to assess the stability of each plume for the COCs. During the 2000–12 period, most of the contaminant plumes for which data were sufficient to complete the analysis were either decreasing or stable in mass and size. The exceptions were carbon tetrachloride (associated solely with the Propellant Burning Ground plume) and 2,4-dinitrotoluene and 2,6-dinitrotoluene (in the Deterrent Burning Ground plume), which showed an increasing trend in mass. No COCs showed an increasing trend in plume mass in the 2013–18 period. Some wells with increasing trends in concentration or with concentrations greater than the enforcement standard are near the tail of a plume, where increased monitoring may be of value to better define future plume boundaries. A spatial optimization analysis covering the 2013–18 period identified six wells that provided information redundant to that from other wells. A temporal optimization analysis identified optimal sampling frequencies for 125 wells. Remedial actions directed at the Propellant Burning Ground plume coincided with a general decrease in plume mass and size, although in specific areas and depths, the plume size for specific contaminants may still be increasing.

Wisconsin↗