USGS Science⌕ Search

Geology topics

Browse USGS reports

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

At least 883 records · Page 49Linked to original sources
USGS2026

Discrete element investigation of the influence of shallow soil density on the manifestations of strike-slip surface fault rupture

This study investigates the influence of soil relative density on strike-slip surface fault rupture manifestations using three-dimensional numerical simulations performed with the discrete element method (DEM). The simulations capture the formation of distinctive fault strands within complex flower structures using tens of millions of grains. The tendency for dense soils to localize shear manifests as multiple localized shear bands within wide zones of deformation in strike-slip fault rupture, whereas diffuse shear deformation develops within narrow zones in loose soils. The spatial extents of soil deformation are consistent between simulations having similar relative densities but with different quantities of grains. However, individual shears are more distinguishable in assemblages of finer grains than in assemblages of coarser grains. The simulations show the progressive development of new shears within the bounds of previously developed shears. Shear activity transitions inward as fault activity diminishes along the outermost shears and continues along newly developed shears until a vertically dipping throughgoing shear structure develops that accommodates most of the fault displacement thereon. The throughgoing fault develops at smaller fault displacements in looser soils because the first shear rupture propagates closer to the vertical direction and does not undergo as much inward translation of shear activity, as is observed in denser soils. In all simulations, ground surface uplift develops between nonintersecting active shears, and ground surface subsidence tends to develop where new shears intersect previous shears. The surface traces in these simulations are shown to be consistent with analog models and case histories of surface fault rupture occurring in different shallow subsurface materials. Although computationally costly, these modeling results are valuable for providing a strong numerical supplement to traditional analog models used to represent the mechanics of strike-slip zones in soil, and they provide quantifiable stresses and large-strain deformations throughout the model domain.

Journal of Geotechnical and Geoenvironmental Engin↗
USGS2026

Using organic compounds and salt mixtures to investigate potential mechanisms of major ion toxicity to the mayfly Neocloeon triangulifer

Field assessments and laboratory studies have demonstrated that mayflies (Ephemeroptera) are relatively sensitive to increased concentrations of major geochemical ions (Na + , Ca 2+ , Mg 2+ , K + , Cl − , SO 4 2− , HCO 3 − /CO 3 2- ), but much remains to be learned about potential mechanisms of toxicity. In the present study, we sought to determine whether solution osmolarity is sufficient to predict the toxic effects of major ions to a mayfly, Neocloeon triangulifer . We further investigated whether binary mixtures of major ion salts with common cations or common anions provide evidence of concentration addition or independent action, potentially illuminating toxicity drivers for this species. Finally, we sought to determine if there are mechanisms for major ion toxicity to mayflies that are comparable to those for other species. The toxicity of single salts showed a Ca-dependent correlation with the total component concentration (the sum of individual ion concentrations, a close correlate of osmolarity), but the effect concentrations were well below those causing toxicity from mannitol (presumed to have only an osmotic effect), suggesting that a general osmotic effect is not responsible for major ion toxicity in this species. The mixture experiments consistently indicated additive toxicity and consistency with the relationship to total component concentration, which supports a single toxic mechanism and thus no evidence for ion-specific mechanisms, as reported for other species. Given the lack of evidence of ion-specific toxicity, the total component metric is likely useful for a broad range of ionic compositions.

Environmental Toxicology and Chemistry↗
USGS2026

Monitoring reservoir processes at The Geysers geothermal field: Joint time-lapse acquisition and cooperative joint inversion of passive seismic and magnetotelluric data

To monitor time-lapse changes in reservoir water and steam concentrations, as well as subsurface flow in the geysers geothermal field in California, we implemented a multi-physics geophysical imaging approach. This approach utilized continuously monitored microseismic events coupled with repeat magnetotelluric (MT) observations. Travel time data from microseismic events recorded at surface stations were collected over an area of approximately 75 km² continuously from 2018 to 2023, alongside three repeat MT surveys conducted in 2021, 2022, and 2023. The resulting seismic and MT data were inverted both separately and jointly to generate time-lapse images of geophysical properties that correspond to temporal changes in water and steam saturation within the reservoir, as well as flow paths and barriers. By correlating the geophysical images with known water injection and steam production volumes, we were able to calibrate the data and gain confidence in the results, which can now be applied throughout the reservoir where borehole data are unavailable. The joint imaging results, combined with reservoir data derived from borehole observations, allow for the interpretation of temporal changes in the reservoir. This enables the operator to optimize reservoir management and their drilling program.

California↗
USGS2026

Bank-dwelling beavers contribute to the wood regime in a dryland river

Our findings suggest that the influence of dwelling beavers on large wood budgets may substantially determine the structure and function of medium to large rivers, particularly in dryland systems such as the White River. Recognizing the role of bank-dwelling beavers expands our understanding of biotic drivers of riverscape complexity and provides new opportunities to integrate ecology with fluvial geomorphology in managing and restoring dynamic river corridors.

Utah↗
USGS2026

Syn-emplacement crystallization of Mauna Loa 2022 lava flows, Hawaiʻi

Advances in near-real-time forecasts of lava flow advance rates, final travel distance, and areal coverage are transforming hazard response, particularly in locations such as Hawaiʻi that experience frequent lava flow activity. The most severe threats are posed by rapidly advancing channelized ʻaʻā flows, such as those that characterized the November–December 2022 eruption of Mauna Loa volcano. To constrain rheological inputs to flow forecasting models during the eruption, samples were collected from the two most persistent lava flows and rapidly assessed for downflow changes in crystallinity and vesicularity. As observed in other channelized lava flows, the lava lost bubbles and cooled during transport, with consequent increases in groundmass crystallinity and bulk viscosity. Lava rheology, however, is controlled by more than bulk bubble and crystal contents; also important are the melt viscosity and the size and shape distributions of bubbles and crystals. Here we further interrogate the sample suite to document a down flow increase of almost two orders of magnitude in melt viscosity alone (a function of cooling combined with compositional change) and assess the extreme anisotropy of the plagioclase crystals as a function of size and flow type (pāhoehoe or ʻaʻā). When combined with previous studies, these data constrain observed relations between surface morphology, flow temperature and crystal content and indicate that remote assessments of morphological transitions and core lava temperature could be used for channel-wide monitoring of rheological evolution for assimilation into flow models.

Hawaii↗
USGS2026

A 40,000-year packrat midden series from Cataviña, central Baja California, Mexico

Thirty packrat ( Neotoma spp.) middens collected from boulder fields near Cataviña, Baja California, Mexico, at 640–680 m elevation provide the first long chronology of macrofossils and pollen spanning the late Quaternary in the Central Desert of Baja California. Midden plant macrofossil and pollen assemblages document a rich chaparral/woodland assemblage during the last glacial and early Holocene dominated by Parry pinyon ( Pinus quadrifolia ) and California juniper ( Juniperus californica ) until 11,630 cal yr BP. This indicates chaparral/woodland had a much more extensive distribution in what are now desert elevations in northern and central Baja California. In contrast to late glacial and early Holocene midden records from northeastern Baja, the Cataviña middens of the same age lack plants adapted to warm season precipitation, suggesting that decreased temperatures and evapotranspiration during the growing season and enhanced winter precipitation, with little contribution from summer rains, supported the lowering of chaparral/woodland species distributions in central Baja California. Cataviña middens also record endemic desert plant taxa mixed in with chaparral/woodland species during the Pleistocene, persisting throughout the Holocene, followed by the quick arrival of other desert species after ∼11,000 cal yr BP. Baja California remains a high-potential yet poorly sampled area for packrat midden research in North America.

Baja California↗
USGS2026

Landscape connectivity and wildlife access to water across an international border: Barriers and opportunities for facilitating transboundary movement

Rapid global acceleration in the construction of physical barriers along international borders has greatly influenced biodiversity and animal movement. Physical barriers can fragment landscapes, hinder access to essential resources, impact long-distance migrations, and inhibit dispersal and gene flow. The effects of physical barriers on animal movement and landscape connectivity can be exacerbated in dryland environments where access to water is a limiting factor. In recent decades, the construction of border barrier infrastructure has accelerated along the international boundary between the United States and Mexico. Here, we used a landscape connectivity model to investigate the effects of barriers on wildlife access to the river in the Lower Rio Grande Valley. We used a modified omnidirectional connectivity model to compare access to the river for three large, terrestrial mammal species across three border barrier scenarios: (1) a landscape without border barriers; (2) a landscape with the existing barrier system; and (3) a potential future landscape with a continuous barrier system. The existing barrier system includes many discrete sections of barrier within tracts of the Lower Rio Grande Valley National Wildlife Refuge or on lands associated with the region's flood control system. Our results indicate that the existing border barriers can impede connectivity and wildlife access to the river in some areas, while some existing gaps between border barrier sections can serve as conduits for wildlife movement and river access. Our future scenario results show how a potential continuous border barrier system could further impede wildlife access to the river. We discuss management and landscape conservation options for enhancing wildlife access to water and riverine habitats. Collectively, our results illustrate the potential effects of border barriers on wildlife movement and access to water, providing information that can be used to better anticipate and lessen the ecological impacts of transboundary barriers.

Tamaulipas, Texas↗
USGS2026

Can surface treatments and climate matching enhance restoration success in the Sonoran Desert?

Introduction Restoring native plant communities in arid ecosystems through seeding is a critical yet often unsuccessful strategy due to severe environmental stressors, including degraded soils, low and variable rainfall, and seed predation. Objectives To improve restoration outcomes, we examined the influence of seed mix types and soil surface treatments on seeded and unseeded plant establishment across three Sonoran Desert sites, Arizona, United States. Methods Using a factorial design, we tested two seed mix types (cool versus warm adapted species) and four soil surface treatments (pits, mulch, connectivity modifiers or ConMods, and controls) on plant establishment. We examined if patterns of plant establishment could be predicted by whether the climate envelope of species in the seed mix types was similar to the site climate where they were seeded (climate matching). We monitored seeded and unseeded species establishment for 3 years and analyzed these effects using generalized linear mixed models and Tukey-adjusted multiple comparisons. Results Pit treatments significantly enhanced seeded species density by nearly 3× and cover by 2× compared with ConMods and increased unseeded plant density by 2× to 4× times during post-drought springs. During these same seasons, species suited to cooler temperatures (e.g. Salvia columbariae and Lupinus sparsiflorus ) had higher cover at the site receiving the highest precipitation, and species suited to warmer temperatures (e.g. Sphaeralcea ambigua and Senna covesii ) established better at the drier two sites. Conclusions In arid systems, surface treatments like pits that retain soil moisture were most promising for supporting seed-based restoration, and matching species that have high temperature tolerance with hotter, drier sites may enhance restoration success.

Arizona↗
USGS2026

An improved empirical model for predicting postfire debris-flow volume in the western United States

Reliable estimates of debris-flow volume can be used to help predict the magnitude of debris-flow hazards following wildfire in the western United States. In this study, we compiled and used a database of 227 postfire debris-flow volumes that were collected across the western United States to develop a multiple linear regression model for predicting postfire debris-flow volume. We explored 36 predictor variables related to rainfall, terrain, and fire characteristics, and selected the model with the combination of variables that yielded the most accurate predictions of debris-flow volume. We evaluated model performance against the entire volume database, as well as against four subsets of volume data from southern California, the Intermountain West, the Southwest, and regions with limited volume data, such as northern California and Washington. We also compared model performance against 3 existing postfire debris-flow volume models that were developed for use in southern California, the Intermountain West, and the Southwest. We demonstrate that the new volume model performs as well as the regional models in the regions for which they were developed and outperforms existing models when applied to volumes from data-limited regions in the western United States. These results indicate that the debris-flow volume model introduced in this study can be used to improve postfire hazard assessments across the western United States, especially outside of southern California.

Arizona, California, Colorado, New Mexico, Utah, W↗
USGS2026

Genomic structural variation rescues a classic biological invader from a population bottleneck

Invasion genetics presents a classic paradox: how do species successfully spread despite severe population bottlenecks? The brown treesnake ( Boiga irregularis ) in Guam represents a striking example of this phenomenon, having been introduced with only a handful of individuals. We show that the population endured an extreme bottleneck, with roughly half of the genome exhibiting runs of homozygosity, comparable to species of conservation concern. Despite this, we uncovered extensive diversity in the form of nearly 19,000 genomic structural variants, which affect almost eight times more of the genome than single-nucleotide variants and provide material for ‘rescuing’ the population from inbreeding-driven declines. Structural variant density was highest in gene promoters, where recombination and DNA repair often occur, providing a mechanism for rapid evolution of gene-linked diversity. This diversity is enriched in genes vital for adaptive immunity and olfaction, suggesting genomic diversity in key chromosomal regions can rescue populations from inbreeding. This work has critical implications for invasion biology and conservation genetics practitioners.

BioRxiv↗
USGS2026

Cook Inlet beluga whale calling varies by group characteristics, behavior, and tidal state

Communication allows social species to exchange information among group members. In aquatic environments, acoustic signals are among the most effective forms of communication and are important for many species, including cetaceans. Beluga whales ( Delphinapterus leucas ) are highly social and vocal, yet little is known about the functionality of their social calls. To examine context-dependent vocal behavior in belugas, we collected passive acoustic data and fine-scale behavioral observations for the endangered Cook Inlet beluga population. The resulting dataset includes 1,720 annotated vocalizations collected over 21 behaviorally encoded encounters. We fit generalized linear mixed models to these data to investigate the effect of behavioral state, group size, calf presence, and tidal state on (1) calling rate (number of calls/minute) and (2) call category (whistles, pulsed calls, combined calls). Belugas were more likely to call when traveling and had higher calling rates during flood tides. Group-level calling rate increased sublinearly with group size, suggesting that individuals called less in larger groups, possibly reflecting increased listening, vocal coordination, or a strategy to avoid acoustic masking. Group calling rate increased before transitions between traveling and milling, suggesting a possible link between communication and behavioral transitions. Whistles were more prevalent when traveling, while pulsed calls were more prevalent when milling. Combined calls occurred only when calves were present, indicating the importance of these calls in communication with calves. Identifying these communication patterns and the contexts in which they occur can enhance our understanding of beluga whale ecology and aid in conservation efforts via passive acoustic monitoring.

Alaska↗
USGS2026

Consolidation and permeability of the B1 and D1 gas hydrate bearing sands and associated seal sediments of the extended-duration gas production test site on the Alaska North Slope

Gas hydrate, a solid combination of gas (mostly methane in nature) and water molecules stable at low temperatures and elevated pressures, occurs naturally in marine and permafrost-associated environments. Gas hydrate reservoirs, such as those in the Alaska North Slope, have been considered potential energy resources for gas production. To understand the petrophysical and geo-mechanical characteristics of the reservoir, core samples retrieved from the site of the JOGMEC-DOE-USGS collaborative gas hydrate R&D project have been analyzed in the laboratory for their hydraulic and mechanical properties. This paper focuses on both seal and reservoir samples associated with the B1 and D1 sands, which are evaluated for index properties (including porosity, grain size distribution, liquid and plastic limits, specific surface area, and specific gravity), consolidation, permeability, and water retention. Furthermore, the reservoir core samples were tested with pore-filling, laboratory-grown tetrahydrofuran hydrate, in order to assess reservoir behavior during gas production from hydrates. Under simulated in situ stress conditions, the seal and hydrate-free reservoir cores had a permeability anisotropy ratio of k h / k v = 3.0–5.0, and k h / k v = 2.4–3.0 for the reservoir tetrahydrofuran hydrate-bearing cores. The data suggest that depressurizing the reservoir to induce hydrate dissociation alters the reservoir effective permeability in three ways: permeabilities decrease due to porosity lost (e.g., the initial reservoir thickness can decrease by up to 5% upon 7 MPa depressurization), permeability increases due to the loss of solid hydrate in the pore space, and permeability anisotropy k h / k v decreases in response to the evolving pore-space geometry. We show that given the simulated in situ gas hydrate saturations (i.e., S h = 32% in core 7P-2E and S h = 21% in core 20P-4), gas production from the dissociation of tetrahydrofuran hydrate in the two tested cores results in a net increase in effective permeability and a decrease in k h / k v . This study highlights the importance of investigating seal and reservoir sediments and the impacts of depressurization on the porosity and permeability responses during production.

Alaska↗
USGS2026

Riverine pesticide trends in the United States: Assessing a decade of national-scale monitoring

Pesticides in freshwater systems can compromise water availability by degrading water quality, with implications for human health and aquatic life. Despite recognition of the need for national-scale monitoring and analysis, few studies have documented long-term trends in surface water pesticide contamination across the US. This study addresses that need by analyzing temporal trends and acute and chronic benchmark exceedances for aquatic life and human health from 81 river sites sampled from 2013 to 2022 using an analytical method targeting 80 pesticides. The majority (79%) of single site and pesticide combinations had too few pesticide detections to estimate trends. When detections were more frequent, increasing trends in concentration were twice as common as decreasing trends. Increasing pesticide concentrations were common in primary drainages of the Mississippi River Basin. Aquatic life benchmarks were exceeded by 19 pesticides, and exceedances were geographically widespread, with both acute and chronic aquatic life benchmark exceedances at 62% of sites. The herbicides atrazine and metolachlor and the insecticide imidacloprid were identified as the greatest threats to surface water availability based on their trends and aquatic life benchmark exceedances. These findings demonstrate the need for continued monitoring and trend analysis, driver investigation, and management strategies to protect freshwater resources.

conterminous United States↗
USGS2026

Two hundred years of historical spawning and nursery data for coregonine fishes in the Laurentian Great Lakes

Historical data can provide critical ecological information for species across the globe, many of which are facing unprecedented rates of ecosystem change. Yet, historical information related to freshwater species, especially fishes, remains scattered, often in original formats, and underutilized for informing conservation and restoration activities. Here, we present a Data Descriptor called Coregonine Spawning History (CORHIST), a database designed to house diverse data related to past spawning and nursery areas for fishes in the family Salmonidae, subfamily Coregoninae (ciscoes and whitefishes), in the Laurentian Great Lakes and their tributaries. Data for 11 species of coregonines historically occurring in the Great Lakes are included in CORHIST. Over 3,400 occurrence records at the coordinate scale have been entered, over 2,200 of which are for Cisco ( Coregonus artedi ) and Lake Whitefish ( C. clupeaformis )—two focal species for which there is either multinational conservation interest or restoration efforts underway in the Laurentian Great Lakes. CORHIST is already proving useful for several studies developing habitat suitability models and delineating spatial units for conservation or restoration planning.

Laurentian Great Lakes↗
USGS2026

Life history traits and population dynamics of Freshwater Drum across large river gradients

Objective Monitoring and assessment of nongame native fishes is limited, but conservation interest in these species is growing. Freshwater Drum Aplodinotus grunniens are a wide-ranging species that serve important functional roles and could serve as an indicator for similar but less common species. Our overall objectives were to quantify and compare population dynamic rates and life history of Freshwater Drum among study reaches in the upper Mississippi and Illinois rivers and relate these metrics to hypothesized environmental and anthropogenic factors. Methods We integrated recently collected age data with monitoring data to estimate age and size distributions, growth curves, maturation schedules, mortality rates, and young-to-adult ratios of Freshwater Drum in six study reaches spanning 1,500 km of river. Principal component analyses and linear regression were used to relate environmental and anthropogenic gradients (latitude, commercial harvest, hydrologic dynamics, primary productivity) to life history traits and population dynamic rates. Results We found latitudinal gradients in life history traits and population dynamic rates whereby Freshwater Drum in upstream, higher-latitude study reaches generally exhibited later maturity, slower growth, smaller maximum size, and lower mortality rates compared with those in lower-latitude study reaches. Further, young-to-adult ratios positively corresponded with chlorophyll- a concentration. No clear relationships were apparent between population dynamic rates and hydrologic variation or commercial harvest. Conclusions Latitude is an important structuring component of life history traits and population dynamics of Freshwater Drum in the upper Mississippi and Illinois rivers likely due to both temperature seasonality and disturbance regimes. The presence of demographic structure in a widespread, common species such as Freshwater Drum suggests similar patterns likely exist in other long-lived native fishes.

Illinois, Iowa, Minnesota, Missouri, Wisconsin↗
USGS2026

Acoustic measurements and modeling of sub-bottom layers on the summit of the Atlantis II Seamount

Sub-bottom profiler images taken from the summit plateau of the Atlantis II Seamount reveal distinct seabed stratigraphy, including marine sediment, limestone, and basalt layers. Acoustic reflection data also show arrivals reflecting from this sub-bottom structure. A wavenumber integration model with elastic geoacoustic properties is able to reproduce the arrival pattern of seabed reflections and particularly the phase inversion of the sub-bottom return. The reflection model suggests that the limestone layer is eroded with high porosity and possesses a lower compressional velocity than a well-cemented layer. The model results also highlight the necessity of incorporating elastic effects for realistic geoacoustic characterization.

JASA Express Letters↗
USGS2026

Reconstructing ancient sedimentary source-to-sink systems – Examples from southern Laurentia’s Proterozoic accretionary orogens

Provenance analysis is a powerful tool for investigating sediment delivery networks, constraining magmatic histories, and reconstructing the tectonic evolution of orogenic belts and basins. Basin analysis studies increasingly use detrital zircon (DZ) U-Pb forward mixture modeling to enhance provenance interpretations by quantifying the relative contributions of different sources. Forward mixture modeling requires significant a priori knowledge that limits deep-time applications. This challenge is overcome with an inverse mixture modeling approach non-negative matrix factorization to reconstruct the number and age distributions of paleo-source regions of Proterozoic metasedimentary rocks in the southwestern United States. This analysis indicates eight reconstructed end-member distributions representing unique sediment sources: two multi-modal end members characterized by ages older than ca. 1.8 Ga from cratonic Laurentia, five unimodal age distributions between ca. 1.80 Ga and 1.65 Ga consistent with Paleoproterozoic arc magmatic sources, and a ca. 1.6−1.5 Ga end member likely derived from exotic cratons in supercontinent Nuna (Columbia). Sediments deposited between ca. 1.80 Ga and 1.73 Ga yield heterogeneous age distributions suggesting multiple arc-backarc systems and several phases of slab roll back, contraction, and accretionary orogenesis, including input from pre−1.8 Ga Laurentian cratons. Homogenization of DZ signatures during the Yavapai orogeny (ca. 1.72−1.68 Ga) reflect crustal assembly as well as the uplift of Paleoproterozoic arcs in the orogenic hinterland. Detrital zircon age distributions from strata deposited during the Mazatzal orogeny (ca. 1.65−1.60 Ga) suggest the Mazatzal Province is a continental arc constructed on older crust. Mesoproterozoic samples are consistent with multiple basins derived from local recycling and long-distance sediment transport. Collectively, these data record the tectonic transition from the episodic accretion of disparate crustal domains to an increasingly integrated continental margin. These results provide new insights into the Proterozoic tectonic and paleogeographic evolution of the southwestern United States at basin to orogen scales and highlight the power of inverse DZ modeling to extract geologically meaningful quantitative mixture models from sedimentary records alone, offering a powerful tool for deep-time tectonic and basin analysis.

Arizona, California, Colorado, New Mexico, Nevada,↗