USGS Science⌕ Search

SEARCH · USGS Science

Results for “Article”

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.

770 records · Page 7Linked to original sources
USGS2024

Automated, near real-time ground-motion processing at the U.S. Geological Survey

We describe automated ground‐motion processing software named gmprocess that has been developed at the U.S. Geological Survey (USGS) in support of near‐real‐time earthquake hazard products. Because of the open‐source development process, this software has benefitted from the involvement and contributions of a broad community and has been used for a wider range of applications than was initially envisioned. Here, we give an overview and introduction to the software, including how it has leveraged other open‐source libraries. We highlight some key features that gmprocess provides, compare response spectra calculated with the automated processing approach of gmprocess to the response spectra provided by the Next Generation Attenuation projects, and summarize projects that have utilized gmprocess. These use‐cases demonstrate that this software development effort has been successfully leveraged in earthquake research activities both within and outside the USGS.

Seismological Research Letters↗
USGS2025

Bølling-Allerød productivity in the subarctic Pacific driven by seasonal upwelling

The Bølling-Allerød deglacial event is marked by high diatom productivity and opal deposition throughout the subarctic Pacific. This opal could either constitute a strengthened biological pump and thus carbon sequestration, or a weakened biological pump and release of marine-sequestered CO 2 to the atmosphere. We quantify silicic acid supply at IODP Site U1340 in the Bering Sea using biogenic opal and δ 30 Si of Coscinodiscus , a diatom genus. These records, along with diatom environmental indicators, suggest the Bølling-Allerød had high silicic acid availability related to a shift from stratification to seasonal upwelling dynamics. We thus propose the primary cause of the high productivity event was increased macronutrient supply from vertical exchange that injected old, nutrient-rich, CO 2 -rich waters into the surface. Enhanced CO 2 release from the subarctic Pacific may help explain critical intervals of CO 2 rise that occur at the onsets of the Bølling-Allerød and PreBoreal.

Geophysical Research Letters↗
USGS2025

Ultralong, supershear rupture of the 2025 Mw 7.7 Mandalay earthquake reveals unaccounted risk

The 28 March 2025 moment magnitude ( M w ) 7.7 earthquake in Mandalay, Burma (Myanmar), ruptured 475 kilometers of the Sagaing Fault, which was more than twice the length predicted by magnitude scaling relationships. Kinematic slip models and observation of a Rayleigh Mach wave that passed through parts of Thailand confirmed that rupture occurred at supershear velocities of greater than 5 kilometers per second. The anomalous length exposed a vast population to violent near-fault shaking. The Mandalay earthquake is a modern analog for the M w 7.9 1906 San Francisco earthquake, another atypically long and fast rupture. Probabilistic seismic hazard analyses use scaling relations that do not account for such long ruptures at moderate magnitudes. This limitation, in conjunction with a likely increased population and infrastructure exposure for atypically long ruptures, contributes to a potential mischaracterization of seismic risk.

Mandalay↗
USGS2026

Generating geochemical and mineralogy distributions of soil in the conterminous United States using Bayesian hierarchical spatial models

Characterizing geochemical and mineralogical soil distributions across large spatial extents is essential for understanding mineral resources, ecosystem processes, and environmental risks. Rasters of soil geochemical distributions for the conterminous United States, however, are limited. We present a Bayesian modeling workflow and tool for generating predictive geochemical and mineralogy distribution maps for the conterminous United States using integrated nested Laplace approximation (INLA) with the stochastic partial differential equation approach. By modeling soil geostatistical data with environmental covariates (soil properties, topography, climate, and land cover), we generate predictive distributions of soil geochemistry that can be mapped or extracted for further analyses. As an example, we model the spatial distribution of trace elements in soil relevant to vertebrate health (cobalt, copper, iron, manganese, selenium, and zinc) and provide a workflow that can be used to generate and visualize predictive distributions of 39 other major and trace elements and 21 minerals of the soil survey, supporting a variety of ecological, environmental, and agricultural applications.

MethodsX↗
USGS2025

Ubiquity and causes of soil water preferential flow across 17 ecoregions

Preferential flow (PF) in soil causes the rapid transport of water, nutrients, and contaminants into the subsurface, influencing groundwater recharge and streamflow. Data scarcity has hindered the quantification of PF occurrence and the identification of its drivers across diverse ecoregions. We address this gap by analyzing high-frequency, multi-depth soil moisture data across 17 ecoregions in the USA, using ∼1,500 sensors at 40 sites. We discovered that PF is widespread, with sites experiencing PF in up to 60% of rainfall events ≥2 mm. Multiple approaches consistently show that PF is more likely to occur with increased peak rainfall intensity, finer textured material, low soil moisture variability, humid climate, and higher net primary productivity. This suggests that PF patterns could shift with projected climate changes, increasing uncertainty in predictions of groundwater recharge, water quality, and streamflow generation.

conterminous United States and Puerto Rico↗
USGS2025

A comprehensive geologic framework of the National Crustal Model for seismic hazard studies in the conterminous United States

A three-dimensional (3D) geologic framework has been developed for the conterminous United States (U.S.) as part of the U.S. Geological Survey National Crustal Model to enhance seismic hazard modeling. The geologic framework is created from geologic maps and multiple subsurface geologic unit boundaries including the base of the Miocene, Cenozoic, Phanerozoic, and the Mohorovičić discontinuity. Modifications are made to surficial geologic maps to remove discontinuities across state and country borders. The subsurface distribution of rock type and age is extrapolated from the surface, seeded with subsurface geologic information, and constrained by a map of basement geology. The framework provides the basis for estimates of subsurface seismic velocity and density that is needed to improve estimates of earthquake ground shaking and seismic hazard. The present framework greatly expands and updates a previously published 3D geologic framework of the western part of the U.S. that was itself a first-of-its-kind digital 3D portrayal of the nation.

conterminous United States↗
USGS2025

Long‐period ground motions from dynamic rupture simulations of large earthquakes on the creeping Hayward–Calaveras–Rodgers Creek fault system

he Hayward, Calaveras, and Rodgers Creek faults in the San Francisco Bay region of California have a high probability of producing a large earthquake in the next decades. Although these faults creep, the creep is insufficient to keep up with their relatively rapid slip rates on their deepest sections, so they have been storing tectonic strain since their last large earthquakes, with the Hayward’s and Rodgers Creek’s more than 150 yr ago. We do not know what the next large Hayward–Calaveras–Rodgers Creek earthquakes will look like or how strongly they will shake the San Francisco Bay region. Harris et al. (2021) used the 3D dynamic (spontaneous) rupture method to simulate large earthquakes on these creeping faults. In this article, we examine the resulting simulated long‐period ( T > 1 s) ground shaking from 0 to 50 km distance, for earthquakes nucleating on the Hayward fault and earthquakes nucleating on the Rodgers Creek fault. We compare these simulated long‐period ground motions with the Boore et al. (2014) well‐established empirically based ground‐motion model suitable for the slowest material velocity in our 3D velocity structure. We find that the simulated long‐period ground motions from the creeping‐fault earthquake scenarios produce a reasonable agreement with the empirical expectations if frictional cohesion is included only where it is appropriate.

California↗
USGS2025

Prospectivity modeling of the NASA VIPER landing site at Mons Mouton near the Lunar South Pole

We use a high-resolution digital elevation model and a numerical thermal model to produce a variety of inputs for a water-ice prospectivity model for the Volatiles Investigating Polar Exploration Rover (VIPER) landing site. These input data are maps of topography, surface slope, surface aspect, surface curvature, maximum temperature, depth to ice stability, permanently shadowed regions (PSRs), distance to PSRs, and PSR density. This model predicts where water ice is most likely within the top meter of regolith, assuming plausible relationships between ice concentration and the various inputs. The model is designed to be adjusted in near-real time as data are collected during the VIPER mission. As such, it is a tool for both analyzing data from the mission as well as planning operations. Since the current model, at this point, relies only on orbital remote sensing, the final version will also be a tool to extrapolate the VIPER mission results across the lunar poles.

Planetary Science Journal↗
USGS2025

Zircon as a pathfinder to REE mineralization

Carbonatites and alkaline silicate rocks are major primary sources of the rare earth elements (REE) and other critical metals, such as Nb. Despite the economic significance of these rocks, their formation and the processes of REE enrichment are poorly understood. Here, statistical analysis of a global dataset demonstrates that zircon geochemistry is a powerful recorder of REE metallogenesis and a potential pathfinder for REE deposits. Zircons from REE and Nb fertile intrusions lack Eu anomalies and have elevated Gd/Yb and Th/Yb, indicating they crystallised from magmas that originated from deep, oxidised and enriched mantle sources. Complexes with Nb enrichment have low U/Nb, reflecting an enriched mantle source, whereas high U/Nb in REE-only fertile intrusions suggest a subduction-metasomatised mantle source. Machine learning models demonstrate high accuracy in classifying zircon from barren and fertile deposits. Classification of detrital zircons shows that REE-enriched deposits correlate with supercontinent assembly, whereas Nb fertile complexes are associated with supercontinent breakup. This approach offers a new, mineral to global scale, petrologic and exploration tool that enhances understanding of REE metallogenesis.

Geochemical Perspectives Letters↗
USGS2025

Potential for continental scientific drilling to inform fault mechanics and earthquake science

Our understanding of fault mechanics and earthquake processes remains limited, largely due to minimal direct observations near active faults at seismogenic depths. This lack of data restricts our ability to accurately assess and mitigate both natural and human-induced seismic hazards. However, recent advancements in drilling capabilities and downhole sensing technologies offer an opportunity: the ability to observe the physical conditions within a volume near active fault zones. In this contribution, we highlight how scientific drilling can provide access to the near-fault environment, enabling measurements of the stress, temperature, fluid pressure, and rock properties at depths where ruptures initiate, propagate, and arrest. These observations are essential to refine models of earthquake nucleation and dynamic rupture, bridging gaps between laboratory experiments, numerical simulations, and surface observations. These insights can advance fundamental understanding in earthquake science but also support the development of more effective seismic hazard assessments and risk mitigation strategies.

Seismica↗
USGS2025

Dating a medieval tsunami with uranium-series techniques on Caribbean corals

Uranium-series dates from coral boulders constrain the timing of a medieval tsunami from the Puerto Rico Trench. Previously reported evidence for this tsunami includes hundreds of coral boulders that came to rest hundreds of meters inland on Anegada, British Virgin Islands. New U-series dates on these coral boulders provide limiting dates for the tsunami. The narrowest limits were by dating interior bands of a coral that retains the hemispherical form of a living coral colony, and which include adjustments for the number of annual density band couplets between the dated samples and the boulder exteriors. By those limits, the tsunami dates between 1381 and 1391 CE, and likely occurred during summer or fall. The tsunami is important as the only reported sign that the eastern Puerto Rico Trench has produced a great earthquake. The dating may aid in defining the earthquake source and in communicating tsunami hazards.

British Virgin Islands, Puerto Rico↗
USGS2025

Rapid fault healing from cementation controls the dynamics of deep slow slip and tremor

Despite its status as one of the most important discoveries in geophysics, the physical mechanism(s) responsible for slow slip events (SSEs) are not well understood. Here, we synthesize observations of deep SSEs in the Cascadia Subduction Zone and argue that rapid, cohesive fault strengthening may control the dynamics of deep SSEs. Cohesive strength is frequently ignored in constitutive laws used to describe fault rheology in numerical simulations of earthquakes and SSEs alike. To demonstrate its importance, we perform and analyze a suite of petrological experiments that simulate fault healing under representative pressure and temperature conditions. We show that significant cohesive strength recovery caused by dissolution-precipitation processes occurs on timescales of just a few hours. Together, our experimental and observational results support the idea that cohesion is a key component of fault strength under SSE conditions and highlight the need for its inclusion in both future experiments and numerical models of fault slip.

British Columbia, Washington↗
USGS2025

Long-term dynamics of earthquake swarms in the Yellowstone caldera

The factors controlling the spatial distribution and temporal evolution of earthquake swarms in volcanic systems remain unclear. We leverage leading-edge deep learning algorithms and a detailed three-dimensional velocity model to construct a 15-year high-resolution earthquake catalog of the Yellowstone caldera region. More than half of the region’s earthquakes are clustered into swarm-like families characterized by episodes of hypocenter expansion and migration. Adjacent earthquake swarms, separated by long quiescent periods, are found to be a dominant feature. We suggest that these swarms are controlled by the interplay between slowly diffusing aqueous fluids and rapid episodic fluid injections, which may result from the breaking of permeability seals. Our analyses also indicate that clustered seismicity beneath the caldera occurs on relatively immature, rougher fault structures, compared to more planar faults outside. Our results provide additional context for understanding seismicity in hydrothermal systems, highlighting the key role played by long-term fluid diffusion processes in driving the occurrence of earthquake swarms.

Idaho, Montana, Wyoming↗
USGS2025

Unveiling coseismic deformation from differenced legacy aerial photography and modern lidar topography: The 1983 M6.9 Borah Peak earthquake, Idaho, USA

The 1983 M6.9 Borah Peak, Idaho, earthquake is one of the largest historical normal fault earthquakes in the western United States. We quantified meter-scale vertical change along the 35 km-long rupture using topographic differencing of 1966 aerial imagery and 2019 lidar-derived data. The initial differencing results are largely obscured by horizontal and vertical georeferencing errors and flight-line stripes. Our error corrections are designed to be insensitive to the coseismic deformation and reduced error by 50%. We calculated vertical separation and resolved a maximum of 2.02 ± 0.46 m at Doublespring Pass. Our vertical separation measurements are generally consistent with those from prior studies using field data and post-earthquake topographic data. However, the differencing measurements are a few decimeters lower than these prior measurements, indicating that differencing can isolate historical from prehistoric earthquake deformation. Our study demonstrates that revisiting historical earthquakes can provide new insights into the magnitude and patterns of coseismic deformation.

Idaho↗
USGS2025

Mapping a Carrington storm

A map is presented of median 1-min-resolution peak geoelectric-field strength across the United States as would be induced by magnetic storms as intense as the 2 September 1859 Carrington storm. The map is constructed from two data sets: Magnetometer time series from 22 ground-based observatories recording 40 magnetic storms, and surface impedance tensors derived from magnetotelluric measurements acquired at 1616 survey sites across the contiguous United States. Carrington-class storm geoelectric fields are likely to be very strong in the United States East and Midwest; > 5.00 V/km at many places. In Virginia, strengths would likely range from 30.30 V/km, with a 68% confidence interval of [19.44,47.20] V/km, to as low as 0.05 [0.03,0.07] V/km. Comparison of model geopotentials with those measured on 30 long lines, indicates errors of about 18%. A Carrington-class storm would likely induce geoelectric fields with strengths 55% greater than for the 13–14 March 1989 storm.

contiguous United States↗
USGS2026

The saline groundwater legacy of a large buried coastal paleo-estuary

Elevated groundwater salinity in coastal regions threatens the beneficial use of fresh groundwater. Coastal groundwater management typically focuses on preventing intrusion from modern sources of seawater; however, past geological processes can also leave a legacy of saline groundwater now hidden in the subsurface. Here, multiple extensive airborne electromagnetic surveys provide detailed evidence of residual salinity from a paleo-estuary filling a late Pleistocene incised valley impacting more than 10,000 km 2 that is now hidden beneath coastal Louisiana’s deltaic plain. Our results show that the three-dimensional pattern of saline groundwater beneath Louisiana mimics that of near-surface aquifers surrounding the modern Delaware Bay estuary, fingerprinting the signature of the past drowning of a large, incised valley of the Mississippi River following post-glacial sea-level rise. These findings demonstrate a new framework for understanding legacy sources of saltwater critical for managing stressed water resources along global coastlines.

Louisiana↗
USGS2024

Evaluation of an open earthquake early warning system in Mexico, and laboratory tests of their sensors

Microelectromechanical system (MEMS) accelerometers are useful for seismological and engineering applications because of their ability to record unsaturated large seismic signals. Recent advances in MEMS technologies enable the design of instruments with improved capabilities that also allow the recording of small signals. As a result, MEMS can be useful across a broad dynamic range and for both major earthquakes and smaller magnitude events. Here, we analyze the performance of a MEMS‐based accelerometer with the capability of real‐time data transmission. We discuss the self‐noise level, dynamic range, and sensitivity, along with the timing precision and data transmission latency, and compare these parameters to other MEMS‐based instruments. We also summarize observations from a regional network deployed in southern Mexico over a period of 3+ yr for the purpose of earthquake early warning. In addition to the sensor evaluation, we present a large, openly available data set of strong motion data from the Mexican network that comprises continuous ground‐motion records from 24 accelerometers in the period between 2017 and 2022.

Seismological Research Letters↗
USGS2026

Subsurface vertical connectivity shapes solute transport to montane streams: Insights from watershed-scale geophysics-informed modeling

Headwater catchments provide essential water and nutrients to downstream ecosystems. The timing and evolution of their transport is shaped in part by the “invisible” subsurface structure but have been largely unexplored at the watershed scale. This raises a fundamental question: how and to what extent does subsurface structure influence streamflow and solute transport in montane watersheds? Here, we test controls of three-dimensionally resolved subsurface structure on solute transport, which reflects source waters and water residence time. We integrate the geophysically mapped architecture of Coal Creek, a headwater catchment in the Colorado River Basin, into numerical flow-transport models. Our results reveal greater vertical connectivity elongates flow paths and enhances deep groundwater contributions to streams, therefore substantially influencing the timing and magnitude of solute transport. In contrast, watershed-integrated outlet discharge can be predicted without high-resolution subsurface information. These findings demonstrate incorporating subsurface architecture is critical for predicting ecosystem health and contaminant transport.

Colorado↗