USGS ScienceSearch

SEARCH · USGS Science

Results for “Soil Dynamics and Earthquake Engineering”

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.

1,696 records · Page 68Linked to original sources

Insights on gas hydrate formation and growth within an interbedded sand reservoir from well logging at the Qiongdongnan Basin, South China Sea

Although variable well log resolution and its control on saturation estimation has been studied, it has not been directly applied to a specific location to explore the nature of gas hydrate within a sand reservoir. We applied in-situ measurements of resistivities, neutron porosity, and gamma ray at two sites in the Qiongdongnan Basin, South China Sea (QDN-W05–2021 and QDN-W08–2021) to investigate the reservoir parameters of a hydrate-bearing sand reservoir. Our results show that gas hydrate is distributed in 5 zones with a total thickness of 10.7 m and an average saturation of 69% at the QDN-W05–2021 site, while they are distributed in 2 zones with a total thickness of 4.3 m and an average saturation of 49% at the QDN-W08–2021 site. We found that variances in saturations estimated from lateral-extra deep button (RX), phase shift (P40H-P40L), and attenuation (A40H-A40L) resistivities within the laterally mapped continuous sand body were affected by the nature of gas hydrate occurrences. Results indicate gas hydrate forms and accumulates at the center of the sand layer and tends to be less or not present toward the top and base. Integrated with seismic data, the in-situ measurements provide insights in the evolution of a mushroom-shaped, hydrate-gas reservoir system. In the system, free gas is likely horizontally transported from the top-center of the gas chimney to the surrounding areas in the early stage dominated by a warm-gas environment, whereas hydrate forms in the opposite pathway starting from the surrounding areas in the following stage with temperature reducing. Our study suggests that high-resolution in-situ measurements not only are a tool to identify the physical properties, but also can be used to help explain the physical process of hydrate growth and accumulation.

Marine Geology

Estimating the probability of export restrictions to inform mineral criticality

To assess risks associated with advanced technologies’ supply chain disruptions, governmental agencies and others have developed mineral “criticality” assessments, with criticality described using the economic impact and probability of supply chain disruptions. Previous work developed subjective supply risk indicators to approximate this probability, typically combining several factors such as supply diversity and trading partners’ political stability, where indicator weightings can substantially impact results. This work explicitly quantifies export barrier probability using an ensemble of machine learning classifiers, with probability estimates informed by exogenous variables, including prior barrier implementation and global export dominance. Major differences in high-probability countries and commodities are observed across models, but the ensemble method highlights Indonesia, China, Tanzania, and the United States as particularly high risk. The Supplementary Data File provides export barrier probability estimates for each analyzed country-commodity pair, enabling a direct, quantitative, objective contribution to assessing mineral criticality, enhancing risk identification and prioritization for policymakers.

Resources, Conservation, and Recycling

Flood-inundation maps for 14.8 miles of Little and Big Papillion Creeks in Omaha, Nebraska, 2023

Digital flood-inundation map libraries for two reaches that constitute 14.8 miles of Little and Big Papillion Creeks in Omaha, Nebraska, were created by the U.S. Geological Survey (USGS) in cooperation with the Papio-Missouri River Natural Resource District. The flood-inundation maps, which can be accessed through the USGS Flood Inundation Mapping Program website at https://www.usgs.gov/mission-areas/water-resources/science/flood-inundation-mapping-fim-program , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at Little Papillion Creek at Irvington, Nebr. (USGS station 06610750), Little Papillion Creek at Ak-Sar-Ben at Omaha, Nebr. (USGS station 06610765), and Big Papillion Creek at Q Street at Omaha, Nebr. (USGS station 06610770) streamgages. Near-real-time stages at these streamgages may be obtained from the USGS National Water Information System database at https://doi.org/10.5066/F7P55KJN or from the National Weather Service Advanced Hydrologic Prediction Service at https://water.weather.gov/ahps/ . Flood profiles were computed for two different reaches that constitute 14.8 miles of stream length in the study area by using hydraulic models. The models were calibrated by adjusting roughness coefficients to best represent the current (2022) stage-streamflow relation at the streamgages within the study reach. The hydraulic models were then used to compute water-surface profiles at 1-foot stage intervals for selected stage ranges to represent various flooding scenarios at the streamgages in each reach. The simulated water-surface profiles then were combined with a digital elevation model using a geographic information system, which had a 10-foot grid spacing to delineate the flooding extents and water depths for each stage. The availability of these flood-inundation maps, along with information regarding current stage from the USGS streamgages, can provide emergency management personnel and residents with information that is critical for flood response activities and post flood recovery efforts.

Nebraska

Prioritizing resource protection and understanding potential susceptibility of springs to surficial changes in a low-temperature geothermal system

Geothermal systems are vulnerable to changes in water budget and composition, requiring science-based management. This study uses a dataset of spring water temperatures, time series of groundwater residence time tracers (tritium and carbon-14), and stable isotopes of water to understand geothermal flow in a low-temperature geothermal system in north west Colorado, United States (Steamboat Springs). The geothermal system is bisected by the Yampa River, necessitating a stream mass balance approach to quantify total discharge. Time series analysis of water temperature data provides a ranked list of features more susceptible to surficial changes, which is corroborated using time series of tritium which indicate spatially distinct patterns of mixing between modern and pre-modern groundwater. All springs contain a portion of pre-modern groundwater that is thousands to tens of thousands of years old, a period coinciding with melting of extensive Pleistocene glaciers that was likely one of the recharge sources to the geothermal system. Stream mass balance indicates that greater than 80% of the total geothermal discharge is derived from diffuse or small springs, highlighting the extensive nature of the geothermal outflow zone and the association with local geologic structures. This study provides baseline data to support management of the Steamboat Springs geothermal system and indicates the utility of these approaches in developing science-based geothermal management.

Colorado

Ductile and brittle Rio Grande Rift deformation in Oligocene granite records a two-stage rift history in southern Colorado

The timing and nature of early deformation in the Rio Grande Rift remains poorly constrained. We present evidence for the earliest structural signature of rift extension in the Sangre de Cristo Range, southern Colorado, based on new geologic mapping, structural analysis, rock magnetic data, and thermochronology. These analyses focus on the ~30.0 Ma granite of Chokecherry Canyon, which hosts discrete low-angle mylonitic shear zones and a distributed, gently SW-dipping protomylonitic fabric. Incremental stretching axes, stretching lineations, and Kmax magnetic lineations plunge gently WSW. Quartz microstructures and crystallographic orientations indicate dominantly coaxial strain in the protomylonite and general shear in the discrete shear zones. Quartz c-axis opening-angle thermometry suggests deformation at ~420–540°C. Thermal modeling of ⁴⁰Ar/³⁹Ar K-feldspar data indicates rapid post magmatic cooling below the brittle–plastic transition, supporting shear-zone formation immediately after emplacement. Slow cooling from ~20–13 Ma was followed by renewed rapid cooling at ~13 Ma, interpreted as the onset of extensional exhumation along the Sangre de Cristo Fault System. These results show that extension in the northern Rio Grande Rift was active by ~30 Ma, earlier than previously recognized. We propose a two-stage model for northern Rio Grande Rift evolution: Stage I (30–23 Ma) records ENE–WSW extension localized in low-angle mylonitic shear zones associated with mid-crustal intrusions; Stage II (≤18 Ma) reflects brittle high-angle normal faulting, focused exhumation, and rift narrowing. Stage I magmatism and deformation along the western range front likely established crustal weaknesses that guided later fault development.

Colorado

Chemical evolution and estimated flow velocity of water in the Trinity Aquifer, south-central Texas

Three permeable zones with varying lithology and water chemistry compose the Trinity aquifer, a principal source of water in the 5,500- square-mile study area in south-central Texas. The upper permeable zone locally yields small quantities of water to wells and was not included in this study. The middle permeable zone primarily is composed of limestone with minor amounts of dolostone. Terrigenous sand and marine limestone, with minor amounts of dolostone, are the principal lithologic units in the lower permeable zone. Dissolved solids concentrations range from 329 to 1,820 milligrams per liter in water samples from the middle permeable zone and from 518 to 3,030 milligrams per liter in water samples from the lower permeable zone. Principal hydrochemical facies in the middle permeable zone are calcium magnesium bicarbonate and calcium magnesium sulfate. Hydrochemical facies in ground-water samples from the lower permeable zone vary. Tritium concentrations as large as 5.3 tritium units in the southeastern part of the study area are indicative of relatively recent recharge. Results of a geochemical mass balance simulation along a flowpath in the middle permeable zone indicate a mass transfer of 4.25 millimoles per liter of dolomite dissolved, 5.74 millimoles per liter of gypsum dissolved, 0.46 millimole per liter of sodium chloride dissolved, 8.07 millimoles per liter of calcite precipitated, and 0.67 millimole per liter of calcium-for-sodium cation exchange between solid and aqueous phases. These results support dedolomitization as a principal chemical process in the middle permeable zone of the Trinity aquifer. Results of a simulation along a flowpath in the lower permeable zone indicate a mass transfer of 0.41 millimole per liter of dolomite dissolved, 0.001 millimole per liter of gypsum dissolved, 9.58 millimoles per liter of sodium chloride dissolved, 1.09 millimoles per liter of calcite precipitated, and 1.11 millimoles per liter of sodium-forcalcium cation exchange between solid and aqueous phases. Lower permeable zone processes indicate sodium chloride dissolution, dedolomitization, and cation exchange. Ground-water-flow velocities determined from adjusted carbon-14 ages, calculated using NETPATH, for selected flowpaths in the middle and lower permeable zones were about 1.7 feet per year and less than about 4.4 feet per year, respectively.

Texas

Where will the cat cross the road? Comparing camera and GPS-based models for identifying wildlife corridors

Designing effective wildlife corridors is a critical conservation challenge in fragmented landscapes. GPS-based step selection functions strongly predict dispersal corridors and connectivity, but GPS collaring can be expensive and invasive. Camera-based occupancy models are widely used for connectivity analyses but may involve trade-offs in data resolution. Despite widespread use of both approaches, few studies have directly compared them using concurrent datasets. We developed a stacked single-species, single-season occupancy model and a Circuitscape connectivity surface for mountain lions (Puma concolor) on Washington’s Olympic Peninsula, USA, and compared them with a connectivity surface from an existing integrated step selection function. Both models predicted mountain lion GPS locations well, with binned Spearman rank correlations of 1 for Circuitscape and 0.96 for the step selection function, though step selection better identified habitat use by dispersers. Connectivity predictions were moderately correlated across the landscape ( r = 0.26), but agreement was strongest in human-dominated areas most critical for corridor planning. We conclude that GPS-based approaches are advantageous when data collection is feasible and the focus is on dispersal or fine-scale movement. However, camera-based approaches may be preferable for multi-species monitoring, large spatial and temporal scales, noninvasive sampling, when resources are limited, or when fine-scale or dispersal-specific inference is not required.

Washington

Light-dependent activity of deepwater sculpin (Myoxocephalus thompsonii) across substrates with and without predation risk

Light availability strongly influences predator–prey interactions in deepwater ecosystems, where visual constraints shape both foraging success and prey behavior. The behavioral response of deepwater sculpin ( Myoxocephalus thompsonii ) to siscowet lake trout ( Salvelinus namaycush siscowet ) was studied under ecologically relevant light intensities spanning several orders of magnitude typical of daytime downwelling light in the 20–100 m depth range of Lake Superior. Trials were conducted over varying substrates (gravel, sand, and black fabric). Deepwater sculpin showed a significant preference for gravel over sand and black fabric. In the absence of siscowet, sculpin movement frequency increased as light intensity decreased. Sculpin reaction distance to siscowet was influenced by both light and substrate. Reaction distance was shortest at low light intensities, peaked at intermediate light intensities (3.05 × 10⁹ to < 6.0 × 10⁹ photons m⁻2 s⁻1), and declined again at the highest light intensities tested. In the presence of siscowet, sculpin activity was suppressed at the upper end of the light range (≥ 6.0 × 10⁹ photons m⁻2 s⁻1). The greatest increase in movement occurred between 6.0 × 10⁹ and 3.05 × 10⁹ photons m⁻2 s⁻1, which corresponded to the range where siscowet prey capture declined, suggesting sculpin exploit this low-light window to move with reduced risk. Reduced activity in the presence of predators is common among cryptic species, and our findings suggest that sculpin restrict movement at higher light levels to avoid detection by siscowet.

Wisconsin

Simulated effects of future water availability and protected species habitat in a perennial wetland, Santa Barbara County, California

This study evaluates the potential water availability in Barka Slough and the effects of changing hydrological conditions on the aquatic habitat of five protected species. Barka Slough is a historically perennial wetland at the downstream western end of the San Antonio Creek Valley watershed (SACVW). A previously published hydrologic model of the SACVW for 1948–2018 was extended to include 2019–2021 and then modified to simulate the future years of 2022–2051. Two models simulating the future years of 2022–2051 were constructed, each with different climate inputs: (1) a repeated historical climate and (2) a 2070-centered Drier Extreme Warming climate (2070 DEW). The model with the 2070 DEW climate had warmer temperatures and an increase in average annual precipitation driven by larger, albeit more infrequent, precipitation events than the model with the historical climate. Simulated groundwater pumpage resulted in cumulative groundwater storage depletion and groundwater-level decline in Barka Slough in both future models. The simulations indicate that Barka Slough may transition from a perennial to an ephemeral wetland. Streamflow, stream disconnection, and depth to groundwater are key habitat metrics for federally listed species in Barka Slough. Future seasonal conditions for each metric are more likely to affect federally listed species’ habitats under 2070 DEW climatic conditions. Future seasonal streamflow volume may negatively impact unarmored threespine stickleback ( Gasterosteus aculeatus williamsoni ) and tidewater goby ( Eucyclogobis newberryi) habitats. Future seasonal stream disconnection may negatively impact the unarmored threespine stickleback habitat. Future groundwater-level decline may negatively impact Gambel’s watercress ( Nasturtium gambelii ) and La Graciosa thistle ( Cirsium scariosum var. loncholepis ) habitats and could influence the ability to use Barka Slough as a restoration or reintroduction site for these species. Results from this study can be used to inform water management decisions to sustain future groundwater availability in the SACVW.

California

Ground deformation and gravity for volcano monitoring

Introduction When magma accumulates or migrates, it can cause pressurization and related ground deformation. Characterization of surface deformation provides important constraints on the potential for future volcanic activity, especially in combination with seismic activity, gas emissions, and other indicators. A wide variety of techniques and instrument types have been applied to the study of ground deformation at volcanoes (sidebar, p. 2; Dzurisin, 2000, 2003, 2007). Geodetic instruments include continuously recording Global Navigation Satellite System (GNSS; of which the United States’ Global Positioning System is one example) stations (fig. D1), borehole tiltmeters, and interferometric synthetic aperture radar (InSAR) measurements (from satellites, occupied and unoccupied aircraft systems, and ground-based sensors). Additional geodetic measurements like continuous- and survey-mode gravity (fig. D2) can contribute substantially to interpreting these data. Borehole strainmeters (see chapter K , this volume, by Hurwitz and Lowenstern, 2024) also have outstanding utility for monitoring deformation, although because of cost and permitting challenges, we do not include them as part of standard volcano monitoring networks for U.S. volcanoes. Still other techniques like light detection and ranging (lidar), structure from motion, and optical satellite data can be used to derive gross topographic changes, which can be used to map volcanic deposits, infer eruption rates, and gain insights into the source processes associated with eruptive activity (see chapter G , this volume, on tracking surface changes caused by volcanic activity; Orr and others, 2024). Experience has shown that no single geodetic monitoring technique is adequate to detect and track the entire range of ground-motion patterns that occur at volcanoes, primarily because of the temporal and spatial diversity of volcano deformation (fig. D3). Similarly, the magnitude of surface deformation varies widely. Geodetic monitoring strategies should therefore include multiple techniques and instrument types to cover a wide range of spatial and temporal scales. In identifying recommendations for geodetic instrumentation for volcano monitoring networks, we attempted to maximize the diversity of instrument types to measure the full range of deformation signals and minimize their expense and number; thus, we do not include several well-known deformation-monitoring techniques in our recommendations. Extensometers, for example, measure strains over distances of a few meters and have an excellent record of success in detecting changes in preeruptive localized ground motion across existing cracks, including at Mount St. Helens, Washington (Iwatsubo and others, 1992), and Piton de la Fournaise, Réunion Island (Peltier and others, 2006). Despite being relatively inexpensive, extensometers are best used primarily when localized ground displacements (for example, ground cracks) need to be tracked, and are not necessary at all volcanoes. In considering volcano deformation monitoring strategies, two complicating factors are deserving of special attention. First, not all deformation is driven by subsurface magmatic activity—for example, at many large stratovolcanoes (for example, Mount Rainier), flank collapses and landslides are significant geologic hazards (Reid and others, 2001) that may occur even in the absence of magmatic activity. Monitoring the stability of volcanoes is thus another critical application of geodetic monitoring networks to inform hazard assessment. One of the most famous examples of edifice instability is the large flank collapse that initiated the May 18, 1980, eruption of Mount St. Helens. Deformation monitoring had detected a bulge on the north flank of the mountain in April 1980 that was expanding by several meters per day (Lipman and others, 1981). Given that flank collapses can happen at any time during a period of volcanic unrest (or even outside a period of unrest), the capability to assess edifice stability is critical. Second, although volcanoes are commonly treated as idealized structures that erupt from single points, like centralvent stratovolcanoes, many are characterized by long rift zones from which eruptions may originate, and distributed volcanic fields are characterized by broadly spaced vents. For example, linear dikes are common at Kīlauea, Mauna Loa, and between Mount Shasta and Medicine Lake in California. At Kīlauea, one of these linear dikes emerged more than 40 kilometers (km) away from the summit of the volcano during the lower East Rift Zone eruption in 2018. Other volcanic fields, like Lassen volcanic center, California, or the San Francisco Volcanic Field, Arizona, have many small vents spread over a wide area. Although the instrumentation guidelines presented in this chapter remain phrased for central-vent volcanoes, they should be modified as needed in the context of the eruptive characteristics of each individual volcanic system. Spatial analysis of geodetic network coverage could help to ensure adequate instrumentation in areas where volcanism can occur over a broad area as opposed to a central vent. As an example, consider the adjacent volcanoes Mount Shasta and Medicine Lake. If station locations are chosen based only on the distance from the centers of the volcanoes, then any geodetic anomalies between the two volcanoes—an area of potential volcanism as indicated by the presence of volcanic features—may remain undetected by ground-based instrumentation. The spatial analysis is accomplished via a grid of pressure point sources (Mogi, 1958) evenly distributed across the map area, at a depth of 5 km in this example (fig. D4). Each source is inflated until predicted deformations exceed the GNSS white noise uncertainty estimates at one site (Langbein, 2017; Murray and Svarc, 2017). This volume of detectable magma provides a measure of the quality of the coverage (fig. D4). The results indicate that, as of 2022, there is a large area between Mount Shasta and Medicine Lake volcano with existing mapped dikes in which a substantial amount of magma could intrude without being detected geodetically. Applying this style of analysis to individual volcanic systems can provide a guide for designing network geometry given the expected locations of future eruptions.

Scientific Investigations Report

Estimated hydrogeologic, spatial, and temporal distribution of self-supplied domestic groundwater withdrawals for aquifers of the Virginia Coastal Plain

Water use from private-domestic wells accounts for nearly 40 percent of total groundwater withdrawals in the Virginia Coastal Plain Physiographic Province (henceforth called the Virginia Coastal Plain). However, because self-supplied domestic water use generally falls below the Virginia Department of Environmental Quality (VDEQ) reporting and management threshold of 300,000 gallons per month, quantifying these withdrawals is challenging. This report builds upon the foundation of previous U.S. Geological Survey investigations by providing revised techniques to improve estimates of the aquifer source, spatial distribution, and monthly magnitude of these groundwater withdrawals. The aquifer sources of private-domestic wells in the Virginia Coastal Plain were estimated by cross-referencing 8,264 well records from the VDEQ and the Virginia Department of Health to a digital model of the Virginia Coastal Plain hydrogeologic framework. This analysis highlights the regional importance of the Yorktown-Eastover, Potomac, and surficial aquifers. Collectively, these three aquifers account for 80 percent of self-supplied domestic groundwater withdrawals. The population using self-supplied domestic water was estimated using census blocks, well-use ratios, building footprints, and land-use and land-cover data to produce a high-resolution, disaggregated, raster-based dataset. This approach improves upon previous models at the census-block or road-network scale by reducing the low-density spread of the self-supplied domestic population across undeveloped areas and concentrating the population and its corresponding water use in the areas where it is most likely to occur. Results show that an estimated 475,332 people comprise the 2020 self-supplied domestic population of the Virginia Coastal Plain, an increase of 5.7 percent since 2010, and the greatest concentrations of self-supplied domestic population surround large cities. Estimates could be further refined with the addition of current and complete spatial data on public water-system service areas. The quantity of water used by the self-supplied domestic population was estimated by modifying published state per-capita water-use coefficients with the corresponding monthly variability assessed from Virginia Coastal Plain public water-system withdrawal data. This analysis estimates an average increase of 12 percent from June through August and an average decrease of 8 percent from December through March from the baseline annual average of 80 gallons per day per capita, which generally matches similar studies in the eastern United States. The application of these revised methodologies for the estimation of private-domestic wells and the self-supplied domestic population improves understanding of domestic groundwater use in the Virginia Coastal Plain across hydrogeologic, spatial, and temporal scales. These revisions help better inform water-resource managers and decision makers and support higher resolution groundwater modeling. Furthermore, these methods are transferrable to other areas where self-supplied domestic water withdrawals are important to the overall water budget.

Virginia

ECCOE Landsat quarterly Calibration and Validation report—Quarter 1, 2024

Executive Summary The U.S. Geological Survey Earth Resources Observation and Science Calibration and Validation (Cal/Val) Center of Excellence (ECCOE) focuses on improving the accuracy, precision, calibration, and product quality of remote-sensing data, leveraging years of multiscale optical system geometric and radiometric calibration and characterization experience. The ECCOE Landsat Cal/Val Team continually monitors the geometric and radiometric performance of active Landsat missions and makes calibration adjustments, as needed, to maintain data quality at the highest level. This report provides observed geometric and radiometric analysis results for Landsats 8 and 9 for quarter 1 (January–March), 2024. All data used to compile the Cal/Val analysis results presented in this report are freely available from the U.S. Geological Survey EarthExplorer website: https://earthexplorer.usgs.gov . This quarterly report is the third to include analysis results for Landsat 9, which was launched in September 2021. The inclusion of Landsat 9 analysis results was dependent on two factors: a complete reprocessing of the Landsat 9 data archive and enough time elapsing to begin formulating lifetime trends. In April 2023, all Landsat 9 image data acquired since the satellite’s launch were reprocessed to take advantage of calibration updates identified by the ECCOE Landsat Cal/Val Team. Additional information about the Landsat 9 reprocessing effort is available at https://www.usgs.gov/landsat-missions/news/upcoming-reprocessing-all-landsat-9-data . Additional information about Landsat 9 prelaunch, commissioning, and early on-orbit imaging performance is available at https://www.mdpi.com/journal/remotesensing/special_issues/15B4V2K92K . This quarterly report is the first to not include analysis results for Landsat 7 because Enhanced Thematic Mapper Plus imaging was suspended on January 19, 2024, after the satellite transitioned into full sunlight. The satellite has been drifting since early 2022 after being lowered from the nominal orbit altitude, and the transition into full sunlight is a result of the satellite operating in its extended science mission. Additional information about the imaging suspension is available at https://www.usgs.gov/landsat-missions/news/landsat-7-imaging-suspended . Additional information about the Landsat 7 extended science mission is available at https://www.usgs.gov/landsat-missions/landsat-7-extended-science-mission .

Open-File Report

Hydrogeology, karst, and groundwater availability of Monroe County, West Virginia

Monroe County is in southeastern West Virginia, encompassing an area of 474 square miles. The area consists of karst and siliciclastic aquifers of Ordovician, Silurian, Devonian, and Mississippian age and is in parts of two physiographic provinces: the Valley and Ridge Province to the east of Peters Mountain, and the Appalachian Plateau Province to the west of Peters Mountain. This study was developed in response to inquiries from the Monroe County Commission requesting assessment of the water resources of the county to better understand the quantity of the county’s groundwater resources, for both current [2023] and future demand, and to provide information to support protection and management of the county’s valuable groundwater resources. Various products were developed for this study that provide knowledge with respect to water availability and contamination susceptibility of the karst aquifers within the county. U.S. Geological Survey (USGS) geologists conducted extensive geologic mapping in support of the project, producing (1) a countywide bedrock geologic map, (2) a countywide hydrogeologic map, and (3) a light detection and ranging (lidar)-derived countywide digital elevation model and associated sinkhole map. A significant part of this work was to map in detail the Greenbrier Group at the formation level, which prior to this study had only partially been completed. The report also includes (4) a description of the lithologic units identified as part of the geologic mapping process. U.S. Geological Survey hydrologists completed several additional products for the hydrology part of the effort, including development of (1) a countywide potentiometric surface (water-table) map, (2) a countywide base-flow stream assessment, (3) countywide water-budget estimates, (4) well log surveys for 15 wells to better understand subsurface controls on groundwater flow within the study area, (5) two groundwater tracer tests to better refine the groundwater divide from the northern and southern parts of the karst aquifer in Monroe County; and finally, based on all available data collected for the study including the potentiometric surface map, geologic map, current [2023] and legacy fluorometric groundwater tracer tests, and base-flow stream assessments, (6) groundwater-basin delineations were reassessed for principal groundwater basins within the Greenbrier aquifer. In Monroe County, four principal hydrogeologic settings produce large yields of water for residential, agricultural, and other uses. The most relied upon water-bearing zone with respect to current [2023] public water supply is from springs along Peters Mountain. These springs are derived from intervals of fractured sandstone and resultant alluvial deposits. Groundwater flows downslope through these permeable alluvial deposits and discharges at the contact with less permeable strata, such as the Reedsville Shale. The second most relied upon water-bearing zone in Monroe County is within the karstic Greenbrier Group aquifer, in which the basal Hillsdale Limestone overlies the less permeable Maccrady Shale. This geologic contact between the Hillsdale Limestone and Maccrady Shale is not only targeted as a source of water for agricultural supply but also is targeted as a source of water for residential supply. The third most relied upon water-bearing zone is composed of shallow perched aquifers within the Greenbrier Group. The discontinuous nature of these perched aquifers makes mapping their extent impossible, but they are related to permeable geologic strata, such as karstified limestones with solutionally enhanced permeability that overlies less permeable shale or chert bedrock. During geologic mapping of the county, several of these perched aquifers were documented in the Pickaway, Union, and Alderson Limestones. A fourth zone consists of springs from Ordovician carbonates at the base of Peters Mountain, which are influenced by sinking streams as well as upwelling along faults. In terms of water quantity, the most sustainable springs are those having deeper-sourced flows. Public supplies are a principal source of water used for residential and commercial supply in the region, accounting for 0.49 million gallons per day (Mgal/d) of fresh-water withdrawals (0.14 Mgal/d of groundwater and 0.35 Mgal/d of surface water) for residential and commercial use and serving 6,645 individuals (49.2 percent of the population). An estimated 6,861 people, (50.8 percent of the population) primarily rely on private wells or other unregulated sources, such as springs, and withdraw 0.55 Mgal/d of groundwater for their residential use. Public water supply in the region is primarily (71.4 percent) derived from springs and augmented by stream withdrawals (backup sources mainly during low-flow periods), with the remaining portion (28.6 percent) derived from groundwater withdrawals from wells. For rural residents, however, 100 percent of their withdrawals are derived from groundwater (wells or springs).

West Virginia

Evaluation of best management practices at an edge-of-field site in the Eagle Creek watershed, Ohio, 2012–20

In 2010, the U.S. Geological Survey worked in partnership with the Great Lakes Restoration Initiative and the Natural Resources Conservation Service to identify farm fields in priority watersheds—watersheds critically important to Great Lake health. Intensive best management practices were implemented in five States to test the efficacy of best management practices, also referred to as agricultural conservation practices, to reduce sediment and nutrient runoff on agricultural fields. These fields were chosen as representatives of the priority watersheds because their farming practices and geographical conditions were common among farms in those watersheds. A farm site located near Findlay, Ohio, in the Eagle Creek watershed used two conservation practices, cover crop (Natural Resources Conservation Service practice standard 340) and variable rate technology (VRT; Natural Resources Conservation Service practice standard 590), which represent common practices in the Maumee priority watershed. This study monitored surface runoff and subsurface tile runoff at the Eagle Creek watershed site. The effects of cover crop and VRT nutrient application on sediment and nutrient runoff were assessed from October 2012 to September 2020. Cover crops were applied in fall of 2016, 2017, and 2018, coinciding with VRT applications. The following parameters were analyzed as part of the study: runoff discharge; peak discharge; and concentrations and loads of total phosphorus, particulate phosphorus, dissolved reactive phosphorus, and suspended sediment. High peak runoff events disproportionally affected nutrient concentrations and loads in runoff. A threshold of the 85th percentile of the peak discharge was chosen to assess the effect of cover crop and VRT nutrient application on nutrient mitigation at each gage. Runoff events below the 85th percentile were considered small, and runoff events over the 85th percentile were considered large. Concentrations from collected water-quality samples were used for the analysis; estimated concentrations were excluded from the statistical evaluation. All surface runoff parameters during small runoff events, except surface runoff volume, were significantly lower in the period with cover crops and VRT implementation. Contrarily, subsurface runoff parameters during small events were not statistically different between the periods with and without cover crops and VRT implementation. Surface and subsurface parameters during large events were also not statistically different between the periods with and without cover crops and VRT implementation. These results indicate cover crops and VRT may improve water quality during small runoff events, but additional best management practices that mitigate large runoff events may lead to greater water-quality improvement given the contribution of large events to overall losses.

Ohio

ECCOE Landsat quarterly Calibration and Validation report—Quarter 2, 2024

Executive Summary The U.S. Geological Survey Earth Resources Observation and Science Calibration and Validation (Cal/Val) Center of Excellence (ECCOE) focuses on improving the accuracy, precision, calibration, and product quality of remote-sensing data, leveraging years of multiscale optical system geometric and radiometric calibration and characterization experience. The ECCOE Landsat Cal/Val Team continually monitors the geometric and radiometric performance of active Landsat missions and makes calibration adjustments, as needed, to maintain data quality at the highest level. This report provides observed geometric and radiometric analysis results for Landsats 8 and 9 for quarter 2 (April–June) of 2024. All data used to compile the Cal/Val analysis results presented in this report are freely available from the U.S. Geological Survey EarthExplorer website: https://earthexplorer.usgs.gov . This is the fourth quarterly report to include analysis results for Landsat 9, which was launched in September 2021. The inclusion of Landsat 9 analysis results was dependent on two factors: a complete reprocessing of the Landsat 9 data archive and enough time elapsing to begin formulating lifetime trends. In April 2023, all Landsat 9 image data acquired since the satellite’s launch were reprocessed to take advantage of calibration updates identified by the ECCOE Landsat Cal/Val Team. Additional information about the Landsat 9 reprocessing effort is available at https://www.usgs.gov/landsat-missions/news/upcoming-reprocessing-all-landsat-9-data . Additional information about Landsat 9 prelaunch, commissioning, and early on-orbit imaging performance is available at https://www.mdpi.com/journal/remotesensing/special_issues/15B4V2K92K . This is the second quarterly report that does not include analysis results for Landsat 7 because Enhanced Thematic Mapper Plus imaging was suspended on January 19, 2024, after the satellite transitioned into full sunlight. The satellite has been drifting since early 2022 when it was lowered from the nominal orbit altitude, and the transition into full sunlight is a result of the satellite operating in its extended science mission. Additional information about the imaging suspension is available at https://www.usgs.gov/landsat-missions/news/landsat-7-imaging-suspended . Additional information about the Landsat 7 extended science mission is available at https://www.usgs.gov/landsat-missions/landsat-7-extended-science-mission .

Open-File Report

Airborne radiometric data map alteration of porphyry copper systems in the Elkhorn district, MT

The Elkhorn district in southwestern Montana is known for past production of high-grade base and precious metals from skarn and carbonate replacement deposits. The district also hosts a porphyry Cu-Mo deposit. A radiometric survey flown over the district and surrounds provides maps of surficial concentrations of potassium, thorium, and uranium. Handheld gamma ray spectrometer measurements on altered and unaltered rocks in the district and surrounding region provide a means to ground truth the airborne data. We demonstrate that the airborne and ground-based measurements show similar map patterns. We calculate ratios of the radioelement concentrations and present a potassium enhancement map that combines potassium and ratios of potassium to thorium and uranium. The results highlight the rocks previously mapped as having K-feldspar-biotite and quartz-sericite alteration assemblages in the Elkhorn district and map additional hydrothermal systems in the region, including known alteration in the Radersburg district. The data enhancement techniques can be used as a screening tool for mapping additional porphyry copper systems.

Montana

A joint Gaussian process model of geochemistry, geophysics, and temperature for groundwater TDS in the San Ardo Oil Field, California, USA

Decline in availability of fresh groundwater has expanded interest in brackish groundwater resources; however, the distribution of brackish groundwater is poorly understood. Water resources in sedimentary basins across the United States often overlie oil and gas development. Mapping of groundwater total dissolved solids (TDS) using data from oil well geophysical logs has become an important technique for identifying fresh and brackish groundwater. Existing geophysical log analysis methods use porosity and temperature to relate formation resistivity to TDS. Typically, natural geothermal gradients are used to estimate temperature at the location of collected resistivity. However, in thermally enhanced oil fields, steam is injected into the subsurface to mobilize high viscosity oil, creating variable temperature distributions. Furthermore, TDS derived from resistivity also depends on the fractions of dominant ions. Typically, chloride and bicarbonate fractions must be determined. It is also necessary to model TDS across many geologic units with heterogenous porosity distributions. Collectively, each quantity used to estimate TDS (resistivity, porosity, temperature, bicarbonate fraction) varies in space and time, and available data points are rarely collocated. Here, we present a new method of mapping groundwater TDS that continuously models each quantity together with a joint Gaussian process. This method enables mapping fresh and brackish water with practically available data. We apply this method to the San Ardo Oil Field in Monterey County, California, where steam injection occurs. In some areas of the aquifer system overlying the oil zone, the temperature is ∼75 °C, roughly twice the natural background value. Groundwater TDS is typically <1,500 mg/L in the aquifer and increases with depth to ∼9,000 mg/L in the oil-producing zone. A low-permeability clay layer delineates the fresh and brackish water, likely by inhibiting surface recharge from penetrating the deeper zones, allowing higher-TDS connate water to remain in place. Weaker lateral TDS trends may be controlled by recharge patterns associated with the Salinas River. Our model reveals with high certainty that groundwater has freshened in one localized part of the oil-producing zone and suggests with less certainty that more widespread freshening has also occurred. The lowering of TDS was possibly from decades of low-TDS steam injection and the associated fluid production and disposal operations.

California

Characterization of hydrogeologic units using matrix properties, Yucca Mountain, Nevada

Determination of the suitability of Yucca Mountain, in southern Nevada, as a geologic repository for high-level radioactive waste requires the use of numerical flow and transport models. Input for these models includes parameters that describe hydrologic properties and the initial and boundary conditions for all rock materials within the unsaturated zone, as well as some of the upper rocks in the saturated zone. There are 30 hydrogeologic units in the unsaturated zone, and each unit is defined by limited ranges where a discrete volume of rock contains similar hydrogeologic properties. These hydrogeologic units can be easily located in space by using three-dimensional lithostratigraphic models based on relation- ships of the properties with the lithostratigraphy. Physical properties of bulk density, porosity, and particle density; flow properties of saturated hydraulic conductivity and moisture-retention characteristics; and the state variables (variables describing the current state of field conditions) of saturation and water potential were determined for each unit. Units were defined using (1) a data base developed from 4,892 rock samples collected from the coring of 23 shallow and 8 deep boreholes, (2) described lithostratigraphic boundaries and corresponding relations to porosity, (3) recognition of transition zones with pronounced changes in properties over short vertical distances, (4) characterization of the influence of mineral alteration on hydrologic properties such as permeability and moisture-retention characteristics, and (5) a statistical analysis to evaluate where boundaries should be adjusted to minimize the variance within layers. This study describes the correlation of hydrologic properties to porosity, a property that is well related to the lithostratigraphy and depositional and cooling history of the volcanic deposits and can, therefore, be modeled to be distributed laterally. Parameters of the hydrogeologic units developed in this study and the relation of flow properties to porosity that are described can be used to produce detailed and accurate representations of the core-scale hydrologic processes ongoing at Yucca Mountain.

Nevada