USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water-Resources Investigations Report”

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,206 records · Page 55Linked to original sources
USGS2023

An inventory of three-dimensional geologic models—U.S. Geological Survey, 2004–22

A database of spatial footprints and characteristics of three-dimensional geological models that were constructed by the U.S. Geological Survey between 2004 and 2022 was compiled as part of ongoing development of subsurface geologic information by the USGS National Cooperative Geologic Mapping Program. This initial inventory resulted in the compilation of 38 three-dimensional geological models that vary widely in their spatial extent, the type and purpose of the model, the number of subsurface units characterized by the model, and the software platforms used to create the model. This Data Report provides the scientific rationale and explanation of the contents of a companion USGS digital data release of spatial data and attributes associated with each three-dimensional model.

Data Report↗
USGS2025

Bathymetric and acoustic-backscatter mapping of Lake Sammamish, Washington, during USGS field activity 2021-656-FA

In 2021, scientists from the U.S. Geological Survey (USGS) acquired bathymetric and acoustic-backscatter data within Lake Sammamish, Washington. Mapping was completed to find evidence of past earthquakes such as underwater landslides and is part of a larger USGS project to understand the overall geologic hazards history of the Cascadia Margin region. The survey was conducted using the USGS research vessel Parke Snavely outfitted with an interferometric sidescan-sonar system for swath mapping. Three estimates of bathymetric uncertainty were calculated from the post-processed bathymetric data. This report provides a summary of the mapping mission.

Washington↗
USGS2026

Archean to Mesozoic–Cenozoic seismic crustal structure: Implications for geological and biological evolution

We use >4500 measurements of crustal structure to investigate the seismic structure of continental crust, Archean to Mesozoic–Cenozoic. The mean crustal thickness of continents, including their margins, is 36.5 km. We find that Archean, Paleoproterozoic, and Mesoproterozoic crust have similar mean crustal thickness (41 km), seismic velocities, Vp/Vs ratio, and density. Crusts of these ages span 3 Ga, from 4.0 to 1.0 Ga, and their similarity of physical properties (thickness, density, and seismic velocities) suggests that the process of crustal formation may have been similar during this time period, which covers 66% of Earth's history. Neoproterozoic crustal properties differ significantly from Mesoproterozoic and older crusts. The mean thickness of Neoproterozoic crust is 32 km, some 9 km thinner than the mean thickness of Archean, Paleoproterozoic, and Mesoproterozoic crust. A lithospheric root with a thickness of 150–200 km underlies Archean, Paleoproterozoic, and Mesoproterozoic crusts, and this root resists lithospheric rifting and crustal extension and thinning. Based on these observations, particularly the lithospheric thickness, we conclude that Archean, Paleoproterozoic, and Mesoproterozoic lithosphere are unique and together form the stable cratonic nuclei, defined as the thick (150–250 km), long-lived cores of continents. Higher mantle temperatures during the Archean, Paleoproterozoic, and Mesoproterozoic may have played a key role in the formation of the thick lithospheric roots. The second key finding is that the similarities of mean seismic properties indicate that the process of crustal formation operating in the Archean continued during the Paleoproterozoic and Mesoproterozoic. The thin (mean value < 135 km) lithospheric root beneath Neoproterozoic and younger crust may be related to the steady decrease in mantle temperature through time. Neoproterozoic and Paleozoic crust have similar physical properties, and these eras are characterized by pronounced biodiversification, including the renowned Garden of Ediacara, the Cambrian Explosion, and the Great Ordovician Biodiversification. Mesozoic–Cenozoic crust is the most diverse and reflects the current tectonic and magmatic processes of crustal formation.

Book chapter↗
USGS2024

2024 Crustal Deformation Modeling Workshop report

The 2024 Crustal Deformation Modeling Workshop was held June 10–14 at the Colorado School of Mines. The workshop included two days of tutorials on using PyLith for crustal deformation modeling, followed by three days of science talks and discussions. The workshop focused on four primary themes: ● Constraining long-term fault slip rates and their uncertainties using geodetic and geologic data; ● Earthquake cycle modeling with a focus on constraining models using seismic and geodetic data; ● Interaction of fluids and faulting; and ● Separating contributions of surface loading and tectonic loading in crustal deformation. The complete agenda is available on the CIG website.

Conference Paper↗
USGS2024

Getting to the surface of the problem: A dynamic rupture benchmark for shallowly dipping faults near Earth’s surface

The Statewide California Earthquake Center (SCEC) Dynamic Rupture Group workshop “Getting to the Surface of the Problem: A Dynamic Rupture Benchmark for Shallowly-Dipping Faults Near Earth’s Surface” was convened virtually, in an online meeting on November 4, 2024. 52 people participated, including two who sent pre-recorded lightning talks earlier, but were unable to join us that day due to the international time-zone difference. One-half of our participants were students or postdocs, and the group included scientists from 11 countries (USA, Canada, China, Czech Republic, France, Germany, Iceland, Japan, New Zealand, Saudi Arabia, Switzerland).

Final Report↗
USGS2026

Preliminary geology of the North Meadow Creek Area, Tobacco Root Mountains, southwest Montana: The North Meadow Creek Fault

This report documents a previously unmapped fault, informally referred to herein as the North Meadow Creek fault, on the east flank of the Tobacco Mountains of southwestern Montana. This fault has an apparently long and complex history, including Quaternary offset of an older alluvial terrace, offset of a debris flow of presumed Pleistocene age, offset of rhyolite lava flows of presumed Tertiary age from their presumed source several kilometers east, on the opposite side of the fault, and offset and possible rotation of blocks of Archean gneiss on northeast side of the fault. The western part of the fault is covered by glacial moraines of two apparent ages, both Pleistocene.

Montana↗
USGS2024

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↗
USGS2026

Ground-motion simulations for the 2024 Mw 4.8 Tewksbury, New Jersey, earthquake

Ground-motion simulations of notable earthquakes in the central and eastern United States are limited and typically assume one-dimensional (1D) Earth structure. In this study, we use a three-dimensional (3D) seismic velocity model to better constrain the depth and focal mechanism of the April 5th, 2024, moment magnitude 4.8 Tewksbury earthquake and investigate the spatial variability of earthquake ground motions and the effects of nearby sedimentary basins. We perform earthquake ground-motion simulations up to 0.5 Hz using the 3D spectral-element wave-propagation solver SPECFEM3D over a region 280-km wide by 260-km long by 77-km deep. Topography and subsurface geophysical structure are assigned using the U.S. Geological Survey National Crustal Model with a minimum shear-wave velocity of 200 m/s. We use earthquake time series from 13 broadband seismic stations in the region that have a uniform azimuthal distribution and epicentral distances ranging from 76 to 131 km to compare with synthetics and explore the effects of 1D versus 3D seismic structure on focal mechanism and depth solutions. Ground-motion intensity metrics are also presented relative to the NGA-East ground-motion models (GMMs) currently used in seismic hazard assessments for the region. We find that the 3D model, which reveals a wide spatial variability of period-dependent ground motions, yields better predictions of earthquake ground motions relative to the 1D model and the NGA-East ergodic ground-motion model, with 76 percent reduction of residual variance in observed ground motions averaged over 3-, 5-, 7-, and 10-second periods. Use of the 3D model to solve for a focal mechanism yields a shallower focal depth at 4 km and a shallower east-dipping focal plane relative to the U.S. Geological Survey regional moment tensor and Global Centroid Moment Tensor. Our study demonstrates that use of 3D seismic velocity models can improve estimates of earthquake focal mechanisms, ground motions, and seismic hazard.

New Jersey↗
USGS2026

Cruise summary—Samoa Basin abyssal mapping—Box coring leg, Pago Pago, Territory of American Samoa to Pago Pago, American Samoa, April 11– May 1, 2026

Expedition Summary U.S. Geological Survey scientists led a box coring effort to the Samoa Basin to characterize minerals and the surrounding abyssal sediments and fauna. Thirty-eight box cores were deployed between April 13, 2026, and April 28, 2026. Thirty-six box cores recovered sufficient material to determine nodule density, and 35 recovered sufficient material for subcores to be collected. The purpose of this Data Report is to provide a summary of samples collected, initial results that were obtained shipboard, and briefly mention planned future analyses from this expedition.

Data Report↗
USGS2026

Creation of aerial camera system boresight installations in Bangor, Maine; Monroe, Louisiana; and Bassville Park, Florida

Aerial camera systems rely on boresights for defining internal and external metrics. The U.S. Fish and Wildlife Service’s Migratory Bird Surveys Branch required boresights to calibrate their camera systems for creating georeferenced imagery that can then be used in geographic information systems. The U.S. Geological Survey created high-accuracy boresight installations in Bangor, Maine, in July 2024; Monroe, Louisiana, in September 2024; and Bassville Park, Florida, in January 2025 to provide this metric-defining capacity to the U.S. Fish and Wildlife Service. This U.S. Geological Survey Data Report details the methods and results of the boresight installation process and provides camera operators with the data needed to use each of the installations.

Florida, Louisiana, Maine↗
USGS2025

3D Converted wave reverse time migration imaging

We describe a newly developed method for recovering high-resolution images of seismic discontinuities, such as subducting slabs, in 3D. Our method makes use of converted S → P or P → S waves observed by dense arrays of seismometers to infer the locations and relative strengths of seismic discontinuities at depth in a target region. Observed direct and converted waves are backpropagated to their times of origin. The time-reversed wavefield is then separated into its constituent P and S components via the Helmholtz decomposition, and those separated wavefields are used to compute imaging functions that characterize the locations and relative strengths of seismic discontinuities. Imaging functions may be designed to use either S → P or P → S waves, so that users can target those arrivals expected to be most dominant in a given dataset. We have previously demonstrated the efficacy of our method in two dimensions, and we now present a 3D implementation of our technique which addresses the significant computational challenges posed by the size of volumetric wavefield data in three dimensions. Through a series of synthetic examples, we demonstrate that our method is capable of recovering the fine scale structure of a subducting slab given realistic station coverage and earthquake sources. We investigate optimal seismic station geometries for our technique and explore image interpretability in regions with poor data coverage. We find that linear station geometries yield more optimal, interpretable imaging functions than collections of small arrays can. We also show that our method can successfully recover bothS → P or P → S images when realistic shear earthquake sources are used, and we explore the additional computational challenges presented by the high frequency content of S waves. Our results demonstrate the potential for our technique to recover high-resolution information about subducting slabs in real-world regions, given that relatively sparse seismic arrays with only approximately 100 stations are capable of recovering interpretable imaging functions from just a few realistic earthquake sources for multiple discontinuities at significant depth in an area of approximately 400~sq~km.

Seismica↗
USGS2026

Overstorey mortality promotes juvenile piñon pine growth during favourable weather at cooler, wetter sites

Hotter droughts have resulted in widespread tree die-off events globally, frequently leading to regeneration failure. Dry forest recovery often depends on the growth and survival of extant juvenile trees. However, it is unclear how microenvironmental changes following overstorey tree die-off affect juvenile trees, particularly in dryland systems where tree recruitment is typically limited by water availability and heat stress. We simulated an overstorey tree die-off event by girdling trees in piñon-juniper woodlands across the south-western United States. We sampled juvenile piñon pine growth from live and dead overstorey treatments across six study sites spanning a regional latitudinal gradient and local elevational gradients. We examined how juvenile branch and needle growth differed between live and dead overstorey treatments, and whether responses varied with weather conditions and juvenile tree size following overstorey mortality. We found greater juvenile branch and needle growth under dead compared with live overstorey trees for 2 years following overstorey mortality at mid- and high-elevation sites which are typically cooler and wetter than the other sites. These observed growth releases were contingent on favourable post-mortality weather conditions. Higher growth under dead overstorey occurred at sites experiencing near-average climatic water deficits compared with sites experiencing above-average climatic water deficits. Growth at all sites increased from the first to second year after overstorey mortality. Across sites, growth was unrelated to juvenile tree size. Synthesis . Our results underscore differentiation in juvenile responses to overstorey tree die-off driven by local site conditions and weather across the range of Pinus edulis . Overstorey mortality resulted in consistently higher juvenile growth only at climatically favourable sites and during favourable weather, while unmeasured microsite differences could help account for variation observed at the hottest and driest site. Results from less climatically favourable sites suggest that overstorey trees neither directly limit nor facilitate juvenile growth, though further study over longer timeframes is needed to resolve the pace and magnitude of potential recovery or decline. Overall, juvenile vigour may be promoted following overstorey mortality only in a narrow spatial (site) and temporal (weather) environmental context, suggesting additional vulnerabilities for piñon populations under more arid conditions.

Arizona, Colorado↗
USGS1976

A comparison of electronic particle counting and pipette techniques in routine mud analysis

Electronic particle counting and pipette size-analysis techniques are compared for their relative effectiveness in large-scale routine mud analysis. Optimum laboratory procedures for conducting electronic analyses with a Model TA Coulter Counter are presented. Comparative electronic and pipette analyses conducted on the same samples resulted in equivalent precision for both techniques. Results of electronic analyses consistently showed size distributions significantly coarser than those obtained by pipette analyses. The laboratory efficiency of the relatively rapid electronic technique is far superior to that of the pipette technique, although more extensive precautionary measures are required to maintain a high level of analytical accuracy and precision. The quality of electronic analyses largely depends upon the effectiveness of standardized operational procedures established within the laboratory.

Open-File Report↗
USGS2026

Increasing artificial light at night enhances salmonid predator foraging efficiency in an urbanized lake

Artificial light at night (ALAN) poses a threat to ecosystems globally. It includes both direct and indirect light, or skyglow, which occurs when ALAN scatters in the atmosphere, extending beyond its original source. We analyzed ALAN trends in Lake Washington, WA, from 2014 to 2023 using Visible Infrared Imaging Radiometer Suite (VIIRS) nighttime light measurements, evaluated the relationship between in situ and satellite measurements, and modeled predator search volumes for a juvenile salmon predator, cutthroat trout ( Oncorhynchus clarki ), as a function of ambient light and turbidity conditions. Open water regions experienced significant increases in ALAN, while nearshore areas primarily showed no or negative trends, revealing the role of skyglow in shaping open water light environments. Using a visual foraging model for cutthroat trout, we found that juvenile salmon at the shallow southern pelagic site experienced light 28 times brighter, resulting in a 168% greater predation vulnerability than those at the northern site. In the Ship Canal, a narrow corridor for outmigrating salmon, predator search volumes were 249% higher than at the southern site. These contrasts in predation vulnerability demonstrate how local conditions influence predator–prey dynamics and provide critical insight for targeting mitigation of both nearshore and distant light sources.

Washington↗
USGS2026

Estimating paleotemperature using stable isotopes of soil-formed phyllosilicates from paleosols: A review

Fossilized soils, or paleosols, contain soil-formed phyllosilicates whose stable isotopic compositions may be used to calculate paleotemperature and thus reconstruct ancient terrestrial environments. Though paleosols are common in the geologic record, the use of phyllosilicates as paleotemperature proxies is limited in the literature owing to difficulties with selecting optimal paleosols, isolation from non-clay minerals and organic materials, mixtures of phyllosilicates in natural samples, wide variations of chemical compositions for phyllosilicates, and limited to undefined equilibrium fractionation factors between phyllosilicates-water. Here, we address these challenges by examining and comparing methods used for sample selection, mineral isolation, pretreatments, mineral identification, conventional and developing methods for oxygen and hydrogen isotopic analyses, and determination of phyllosilicate-water equilibrium fractionation factors, concluding with recommendations for best approaches for paleotemperature estimation. Additionally, we discuss how to identify and avoid detrital phyllosilicates, the impacts of diagenesis, comparison of stable isotope and non-isotope paleosol paleotemperature proxies, and challenges and opportunities for broadly using paleosols as paleoclimate archives. With ongoing efforts to refine this multi-faceted paleotemperature approach, the stable isotope geochemistry of soil-formed phyllosilicates continues to be an invaluable proxy system, enhancing our understanding of terrestrial paleoenvironments and paleoclimate.

Kansas↗