USGS Science⌕ Search

SEARCH · USGS Science

Results for “Interfaces”

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

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

At least 1,207 records · Page 67Linked to original sources
USGS2021

Relational database for horizontal‐to‐vertical spectral ratios

Frequency‐dependent horizontal‐to‐vertical spectral ratios (HVSRs) of Fourier amplitudes from three‐component recordings can provide useful information for site response modeling. However, such information is not incorporated into most ground‐motion models, including those from Next‐Generation Attenuation projects, which instead use the time‐averaged shear‐wave velocity ( ⁠ V S "> V S VS ⁠ ) in the upper 30 m of the site and sediment depth terms. To facilitate utilization of HVSR, we developed a publicly accessible relational database. This database is adapted from a similar repository for V S "> V S VS data and provides microtremor‐based HVSR data (mHVSR) and supporting metadata, but not parameters derived from the data. Users can interact with the data directly within a web portal that contains a graphical user interface (GUI) or through external tools that perform cloud‐based computations. Within the database GUI, the median horizontal‐component mHVSR can be plotted against frequency, with the mean and mean ± one standard deviation (representing variability across time windows) provided. Using external interactive tools (provided as a Jupyter Notebook and an R script), users can replot mHVSR (as in the database) or create polar plots. These tools can also derive parameters of potential interest for modeling purposes, including a binary variable indicating whether an mHVSR plot contains peaks, as well as the fitted properties of those peaks (frequencies, amplitudes, and widths). Metadata are also accessible, which includes site location, details about the instruments used to make the measurements, and data processing information related to windowing, antitrigger routines, and filtering.

Seismological Research Letters↗
USGS2022

Prospective and retrospective evaluation of the U.S. Geological Survey public aftershock forecast for the 2019-2021 Southwest Puerto Rico Earthquake and aftershocks

The M w "> M w 6.4 Southwest Puerto Rico Earthquake of 7 January 2020 was accompanied by a robust fore‐ and aftershock sequence. The U.S. Geological Survey (USGS) has issued regular aftershock forecasts for more than a year since the mainshock, available on a public webpage. Forecasts were accompanied by interpretive and informational material, published in English and Spanish. Informational products included narrative “scenarios” for how the aftershock sequence could play out, infographics, and a report on the potential duration of the aftershock sequence through the next decade. Forecasts are based on the epidemic‐type aftershock sequence (ETAS) model and generated using the USGS AftershockForecaster software—an interactive graphical user interface built on the OpenSHA platform ( Field et al. , 2003 ). The initial forecast is based on past sequences in similar tectonic environments; subsequent forecasts are tuned to the ongoing sequence via Bayesian model updating. Probabilistic aftershock forecasts for the next day, week, month, and year were publicly released and archived at a daily to monthly tempo, allowing for a truly prospective test of the forecast. Here, we evaluate the forecast over the first year of the recorded aftershocks. The ETAS‐based forecast performed well overall, successfully capturing both the chance of having at least one earthquake of a given magnitude in a forecast interval as well as the non‐Poissonian distribution of the total number of aftershocks within an interval. A retrospective analysis shows that the ETAS model is a substantial improvement over the existing Reasenberg and Jones (1989) forecast model. The exercise also reveals some limitations of the current model, in particular, with respect to nonstationarities in the aftershock magnitude distribution and model parameters throughout the evolving sequence.

Puerto Rico↗
USGS2023

Earthquake scenarios for Quito, Ecuador; Cali, Colombia; and Santiago De Los Caballeros, Dominican Republic

Earthquake risk associated with Quito, Ecuador; Cali, Colombia; and Santiago de los Caballeros, Dominican Republic is examined by generating a set of hypothetical earthquake scenarios considering seismic sources, recent seismicity, and major historical earthquakes recorded in the vicinity. In this study, particular focus is given to the development of earthquake scenarios for use in emergency planning in each city as well as stimulating discussion with respect to addressing the gaps in current understanding of built stock inventory and their vulnerability when subjected to strong shaking. Exposure and vulnerability models from the Global Earthquake Model foundation, hazard and local site information available for the cities, and the U.S. Geological Survey near‐real‐time products are utilized to estimate potential consequences for postearthquake response planning. Results showed that the historic city centers remain the most susceptible to experiencing severe damage resulting in widespread casualties. Similarly, the scenarios highlight areas susceptible to shaking induced ground failure hazards, which may pose additional challenges when responding to such earthquakes. Moderate earthquakes originating from nearby seismic sources, for example, Quito fault system for Quito or the Septentrional fault zone in the case of Santiago de Los Caballeros, could potentially be of greater consequence in terms of direct economic impact and disruption to the city when compared to very large distant subduction interface earthquakes.

Cali, Santiago De Los Caballeros, Quito↗
USGS2024

SKHASH: A python package for computing earthquake focal mechanisms

We introduce a Python package for computing focal mechanism solutions. This algorithm, which we refer to as SKHASH, is largely based on the HASH algorithm originally written in Fortran over 20 yr ago. HASH innovated the use of suites of solutions, spanning the expected errors in polarities and takeoff angles, to estimate focal mechanism uncertainty. SKHASH benefits from new features with flexible input formats and allows users to take advantage of recent advances in constraining focal mechanisms for small magnitude or poorly recorded earthquakes. The 3D locations of earthquakes and the velocity models used are varied when finding acceptable solutions. As a result, source–receiver azimuths are reflective of errors from the earthquake locations and velocity models, in addition to the takeoff angles. Users can consider weighted P ‐wave first‐motion polarities derived from traditional or machine‐learning picks, cross‐correlation consensus, and/or imputation techniques using SKHASH. Focal mechanism solutions can also be further constrained using traditional, machine learning, and/or cross‐correlation consensus S / P amplitude ratios. With improved reporting of individual and collective P polarity and S / P amplitude misfits, users can better evaluate the success of the solutions and the quality of the measurements. The reporting also makes it easier to identify potential issues with metadata, including incorrectly reported station polarity reversals. In addition, by leveraging vectorized operations, taking advantage of an efficient backend Python C Application Programming Interface, and the use of a parallel environment, the Python SKHASH routine may compute mechanisms quicker than the HASH routine.

Seismological Research Letters↗
USGS2025

High-pass corner frequency selection and review tool for use in ground-motion processing

Raw seismological waveform data contain noise from the instrument’s surroundings and the instrument itself that can dominate recordings at low and high frequencies. To use these data in ground‐motion modeling, the effects of noise on the signals must be reduced and the signals’ usable frequency range identified. We present automated procedures to efficiently reduce low‐frequency noise that are implemented in the software package gmprocess. These procedures check for, and as needed remove, low‐frequency artifacts in the displacement record using polynomial fits, which can be used in combination with existing signal‐to‐noise ratio (SNR)‐based corner‐frequency selection procedures. The automated selections are then efficiently verified and refined using a graphical user interface (GUI) that plots relevant ground‐motion time series and spectra and tracks modifications to signal processing parameters. We demonstrate these procedures using recordings from the 2020 M 5.1 Sparta, North Carolina, and the 2013 M 4.7 southern Ontario earthquakes. Data processed with the SNR‐only and polynomial criteria for these events contain displacement artifacts in 37% and 23% of processed traces, respectively. Records with remaining artifacts are corrected manually using the GUI. These processing steps illustrate the workflow for efficient data processing with quality control.

Seismological Research Letters↗
USGS2025

Using gridded seismicity to forecast the long-term spatial distribution of earthquakes for the 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

Gridded (or background) seismicity models are a critical component of probabilistic seismic hazard assessments, accounting for off‐fault and smaller‐magnitude earthquakes. They are typically developed by declustering and spatially smoothing an earthquake catalog to estimate a long‐term seismicity rate that can be used to forecast future earthquakes. Here, we present new gridded seismicity models for use in the 2025 National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The previous PRVI NSHM was released in 2003, and our new models incorporate updates to both data and methodology. We utilize an updated earthquake catalog based on improved Puerto Rico Seismic Network data with newly characterized completeness epochs. The catalog is divided into crustal, subduction interface, and intraslab seismicity using new methods and Slab2 subduction zone geometries. To forecast the long‐term spatial distribution of earthquakes, we use an updated methodology developed for the 2023 U.S. 50‐state NSHM, considering three declustering methods and two spatial smoothing methods based on 2D Gaussian kernels. To adapt it for the complex seismotectonics of the region, we also adopt probabilistic methods to account for events with unknown depths and uncertainties in tectonic classification, and develop a new method for spatial scaling to counteract the effects of spatial variability in network coverage while maintaining the use of smaller events. Finally, we test the performance of these spatial models in forecasting the location of M w ≥ 5earthquakes in the region. Our updated methodology improves the representation of epistemic uncertainty relative to the 2003 model, and our results demonstrate the effectiveness of the new measures we have introduced to address heterogeneities in network detection and systematically evaluate forecast performance.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

Wavelet Inversion for SliP (WISP): Open-source earthquake slip modeling software

Models of the spatiotemporal evolution of earthquake slip, termed finite-fault models, are a critical component of rapid earthquake and tsunami response, earthquake forecasting, seismic ground-motion estimates, and studies of earthquake kinematics. Here, we detail a newly released finite-fault modeling software, Wavelet Inversion for SliP (WISP), in use at the U.S. Geological Survey’s National Earthquake Information Center (NEIC) and available to the public. WISP version 1.1.0 allows inversion of teleseismic body and surface waves, as well as local strong-motion, static and dynamic Global Navigation Satellite System, and satellite imagery (e.g., Interferometric Synthetic Aperture Radar) observations on single or multiple planar fault segments. The software is used in NEIC rapid response of earthquakes M w ≥ 7, generally resulting in a published model within the first few hours after the event origin time. The rupture location and dimensions are then used as inputs to downstream products to estimate earthquake shaking, predict loss, and model the likelihood of secondary hazards, namely landslides and liquefaction. WISP is also used in research studies to evaluate the characteristics of complex ruptures including multifault ruptures and earthquake doublets, among others. The WISP version 1.1.0 software release is composed of Python-wrapped FORTRAN code to accomplish the inversion procedure. A simple command line interface facilitates ease of use even for those with only a cursory knowledge of Python scripting. WISP version 1.1.0 includes a Jupyter Notebook tutorial demonstrating use of the software for modeling the 2015 M w 8.3 Illapel, Chile, earthquake. In parallel with the tutorial, we demonstrate the typical usage of the WISP software using the M w 8.3 Illapel earthquake example here.

Seismological Research Letters↗
USGS2025

Subduction zone earthquake catalog separation tool: Implementation in the USGS 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

The U.S. Geological Survey (USGS) periodically releases updates to National Seismic Hazard Model (NSHM) for the United States and its territories leveraging current scientific knowledge and methodologies to guide public policy, building codes, and risk assessments regarding potential ground shaking due to earthquakes that may result in infrastructure damage. In subduction zones, there is a need to separate the earthquake catalog into tectonic regimes to create specific seismicity models for which the most appropriate ground‐motion models are then applied. Here, we describe newly developed methods and software, called CatSep, that classifies subduction zone events into three primary tectonic regimes: crustal, interface, and intraslab. This method incorporates information about the location of the earthquake relative to the subducting slab, the depth of the Mohorovičić discontinuity, and the earthquake’s moment tensor. Applying this method is a first step in the NSHM workflow for regions covering U.S. subduction zones. Results using this subduction zone earthquake catalog separation tool for the 2025 Puerto Rico and U.S. Virgin Islands NSHM earthquake catalog are presented and analyzed.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

The digital archivist: Automating legacy macroseismic data processing using large language models

Macroseismic data are a key resource to investigate shaking and damage from preinstrumental and early instrumental eras. However, data are often stored as inconsistently formatted reports describing observed shaking and damage, making manually parsing and interpreting accounts labor‐intensive. We introduce a novel workflow using Google’s Gemini 2.5 Pro large language model (LLM) to automate the extraction and structuring of macroseismic observations from summary reports. We apply this workflow to the 22 March 1957 M 5.3 Daly City, California, earthquake as a case study. We used Gemini to extract addresses, originally assigned modified Mercalli intensity values, and descriptions from each report. To address coordinate precision limits, addresses were geocoded via Google’s Geocoding application programming interface. This workflow yielded over 2300 geocoded intensity reports for the Daly City earthquake. We use the geocoded accounts, with the original report intensity assignments, to develop a shaking intensity map that in some respects rivals modern Did You Feel It? Maps. We also extract and present data for the 9 February 1971 M L 6.7 Sylmar, California, earthquake. Our results demonstrate the potential of LLMs for reliably extracting and analyzing large, unstructured macroseismic datasets. LLMs offer a scalable solution for rapidly digitizing macroseismic archives, enabling their broader use to constrain ground‐motion models in modern seismic hazard analysis and to improve our understanding of site effects in urban areas. The concepts explored here may also be applied to the handling of other legacy seismological and earth science data.

Seismological Research Letters↗
USGS2026

Amploc: A complete Python workflow for amplitude-based locations

Amplitude‐based locations are well suited for characterizing long‐duration seismic sources that may migrate with time, such as debris flows, avalanches, rockfall, and volcanic tremor. Although amplitude‐based location methods have existed for decades, their adoption in real‐time monitoring systems remains limited. We present a complete Python‐based workflow that includes station corrections, amplitude‐based location algorithms, and interactive review tools that are appropriate for real‐time use. The algorithms themselves are not new; however, Amploc integrates several previously published methods into a unified operational framework. The workflow also introduces a browser‐based graphical interface that allows users to review real‐time location results using formats familiar to many observatories. By integrating calibration, location algorithms, and visualization tools, Amploc lowers the barrier to operational use of amplitude‐based locations and improves characterization of phenomena that inform hazard assessments such as long‐period earthquakes, volcanic tremor, and debris flows.

British Columbia, Oregon, Washington↗
USGS2024

Distinct yet adjacent earthquake sequences near the Mendocino Triple Junction: 20 December 2021 Mw 6.1 and 6.0 Petrolia, and 20 December 2022 Mw 6.4 Ferndale

Two earthquake sequences occurred a year apart at the Mendocino Triple Junction in northern California: first the 20 December 2021 �w 6.1 and 6.0 Petrolia sequence, then the 20 December 2022 �w 6.4 Ferndale sequence. To delineate active faults and understand the relationship between these sequences, we applied an automated deep‐learning workflow to create enhanced and relocated earthquake catalogs for both the sequences. The enhanced catalog newly identified more than 14,000 M 0–2 earthquakes and also found 852 of 860 already cataloged events. We found that deep‐learning and template‐matching approaches complement each other to improve catalog completeness because deep learning finds more M 0–2 background seismicity, whereas template‐matching finds the smallest M < 0 events near already known events. The enhanced catalog revealed that the 2021 Petrolia and 2022 Ferndale sequences were distinct in space and time, but adjacent in space. Though both the sequences happened in the downgoing Gorda slab, the shallower Ferndale sequence ruptured within the uppermost slab near the subduction interface, while the onshore Petrolia sequence occurred deeper in the mantle. Deep‐learning‐enhanced earthquake catalogs could help monitor evolving earthquake sequences, identify detailed seismogenic fault structures, and understand space–time variations in earthquake rupture and sequence behavior in a complex tectonic setting.

California↗
USGS2026

Rapid characterization of the 2025 Mw 8.8 Kamchatka, Russia earthquake

The 29 July 2025 M w 8.8 Kamchatka, Russia, earthquake was the sixth largest instrumentally recorded earthquake. This event was seismically well observed at regional and teleseismic distances, but publicly available near‐source data were sparse at the time of the event, presenting unique challenges for rapid source and impact characterization. The U.S. Geological Survey (USGS) National Earthquake Information Center provides global real‐time monitoring for earthquakes, including rapid response information products that estimate source characteristics, shaking, and the resulting impacts. We describe the USGS rapid response earthquake information products following the Kamchatka event and discuss their implications for ongoing hazards in the region. We describe potential improvements to our response workflows motivated by this event, including more rapid constraints on source geometries and the automated selection of fault geometries for finite‐fault inversions. The rapid response products together support the interpretation of a unilateral southwestward rupture with significant slip on the southwestern end of the rupture extent. The M w 8.8–9.0 event in 1952, which ruptured a comparable extent of the Kuril–Kamchatka subduction interface, has many similarities to the 2025 rupture. This illustrates that slip deficits may remain following great earthquakes and highlights the usefulness of comparative studies between historic and modern events.

Kuril-Kamchatka subduction zone↗
USGS1977

Reflection and refraction of type-II S waves in elastic and anelastic media

The general theory of viscoelasticity, which accounts for elastic as well as anelastic linear behavior of materials, predicts that two types of S waves propagate in anelastic earth materials. The particle motion for an inhomogeneous plane S wave of type I is elliptical in the plane defined by the directions of propagation and attenuation, while the particle motion for an inhomogeneous plane S wave of type II is linear perpendicular to this plane. The general theory predicts that an S -wave incident upon a plane boundary perpendicular to the plane defined by the directions of propagation and attenuation generates S waves only of the same type. General characteristics of the type-II S waves reflected and refracted at plane anelastic boundaries are: (a) velocities and maximum attenuations which depend on the angle of incidence and frequency, (b) maximum energy flow at a different velocity and in a different direction than phase propagation, (c) energy flow across the boundary due to interaction of the incident and reflected waves. The general theory predicts these characteristics for the waves whenever a plane type-II S wave interacts with a plane anelastic boundary such as a soil-bedrock, crust-mantle, or core-mantle interface. None of these characteristics are predicted for the plane SH waves described by elasticity theory.

Bulletin of Seismological Society of America↗
USGS1988

Crustal velocities near Coalinga, California, modeled from a combined earthquake/explosion refraction profile

Crustal velocity structure for the region near Coalinga, California, has been derived from both earthquake and explosion seismic phase data recorded along a NW-SE seismic-refraction profile on the western flank of the Great Valley east of the Diablo Range. Comparison of the two data sets reveals P-wave phases in common which can be correlated with changes in the velocity structure below the earthquake hypocenters. In addition, the earthquake records reveal secondary phases at station ranges of less than 20 km that could be the result of S- to P-wave conversions at velocity interfaces above the earthquake hypocenters. Two-dimensional ray-trace modeling of the P-wave travel times resulted in a P-wave velocity model for the western flank of the Great Valley comprised of: (a) a 7- to 9-km-thick section of sedimentary strata with velocities similar to those found elsewhere in the Great Valley (1.6 to 5.2 km/sec); (b) a middle crust extending to about 14 km depth with velocities comparable to those reported for the Franciscan assemblage in the Diablo Range (5.6 to 5.9 km/sec); and (c) a 13- to 14-km-thick lower crust with velocities similar to those reported beneath the Diablo Range and the Great Valley (6.5 to 7.30 km/sec). This lower crust may have been derived from subducted oceanic crust that was thickened by accretionary underplating or crustal shortening.

California↗
USGS1994

Wide-band analysis of the 3 March 1985 central Chile earthquake: Overall source process and rupture history

We apply a linear, finite-fault waveform inversion scheme to the near-source strong-motion records, the teleseismic body waves, and the long-period Rayleigh waves recorded for the 3 March 1985 Chile earthquake to recover the mainshock rupture history. The data contain periods between about 2 and 350 sec and are inverted by allowing a variable dislocation rise time at each point on the fault. The results indicate that the mainshock had a seismic moment of 1.5 × 10 28 dyne-cm ( M w 8.0) and ruptured mainly updip and to the south of the hypocenter for a distance of about 150 km along the Nazca-South America plate boundary. A smaller northward component of propagation is also evident, giving a total rupture length of about 200 km. The total source duration of the mainshock is 70 sec, with the majority of the slip occurring within the first 40 sec in a broad 100-km-wide zone in the northern half of the rupture area. Slip in this region extends from a depth of 55 km to within about 10 km of the surface and contains two areas of maximum slip (2.3 and 2.9 m) with rise times of approximately 14 sec. Slip in the southern portion of the fault reaches lower peak values (1.8 m) and extends downdip to depths no greater than 30 km. An independent variable rise-time inversion of the teleseismic body waves alone yields similar results, indicating that a significant component of slow fault motion is not required for this earthquake. The mainshock was preceded by several smaller precursors, the largest of which is an M w ∼ 6.6 thrust earthquake occurring at a depth of 22 km in the shallow 15° dipping portion of the plate interface.

Bulletin of the Seismological Society of America↗
USGS1997

Estimating the diminution of shear-wave amplitude with distance: Application to the Los Angeles, California, urban area

The rate of decay with distance of shear-wave amplitude, computed from 20-sec S -wave spectra, is determined from TERRA scope records of small earthquakes in the greater Los Angeles area. Piecewise log-linear interpolation functions and traditional diminution functions are used to fit spectral decay to a maximum distance of 150 km. Simultaneously, isotropic source and receiver terms are determined. Separate branches of the spectral decay function are found for two categories of source depth: greater than 10 km and less than 10 km. In the hypocentral distance range of 20 to 150 km and in the frequency range of 0.5 to 8.0 Hz, an important result of the investigation is that the horizontal-component decay rate associated with deeper-crustal sources is generally greater than that associated with shallower sources and is greater than that which is estimated using more traditional models of spectral decay with distance. The same behavior generally holds for vertical-component spectra. The variation in apparent attenuation rate with source depth should affect seismic-hazard estimates associated with the rupture of blind thrust faults in the Los Angeles basin and vicinity. The results of the inversions suggest that interpolation function representations of spectral decay are sensitive to perturbations of S -wave amplitude due to crustal reflectors, such as post-critical S -wave arrivals from mid-crustal to deep-crustal velocity interfaces.

California↗
USGS1999

Fault-slip distribution of the 1995 Colima-Jalisco, Mexico, earthquake

Broadband teleseismic P waves have been analyzed to recover the rupture history of the large ( M S 7.4) Colima-Jalisco, Mexico, shallow interplate thrust earthquake of 9 October 1995. Ground-displacement records in the period range of 1-60 sec are inverted using a linear, finite-fault waveform inversion procedure that allows a variable dislocation duration on a prescribed fault. The method is applied using both a narrow fault that simulates a line source with a dislocation window of 50 sec and a wide fault with a possible rise time of up to 20 sec that additionally allows slip updip and downdip from the hypocenter. The line-source analysis provides a spatio-temporal image of the slip distribution consisting of several large sources located northwest of the hypocenter and spanning a range of rupture velocities. The two-dimensional finite-fault inversion allows slip over this rupture-velocity range and indicates that the greatest coseismic displacement (3-4 m) is located between 70 and 130 km from the hypocenter at depths shallower than about 15 km. Slip in this shallow region consists of two major sources, one of which is delayed by about 10 sec relative to a coherent propagation of rupture along the plate interface. These two slip sources account for about one-third of the total P -wave seismic moment of 8.3 × 10 27 dyne-cm ( M w 7.9) and may have been responsible for the local tsunami observed along the coast following the earthquake.

Colima-Jalisco↗
USGS2005

Hazard calculations for the WGCEP-2002 earthquake forecast using OpenSHA and distributed object technologies

We present seismic-hazard calculations for what is arguably the most sophisticated earthquake forecast ever developed—the model by the 2002 Working Group on California Earthquake Probabilities (2003), or WGCEP-2002 as referred to hereafter. These calculations have been made possible by developments in both OpenSHA ( Field et al. , 2003 ) and the Information Technology Research (ITR) Collaboration of the Southern California Earthquake Center (SCEC) ( Jordan et al. , 2003 ). In particular, we demonstrate the use of a freely available, platform-independent, and graphical-user-interface-based application for computing hazard curves. This application utilizes distributed-object technologies, meaning the code running on the client's computer communicates over the Internet with code hosted elsewhere on server computers. This interoperability, which is invisible to the user, makes the application fast and executable without heavy installation requirements. Instructions on how to reproduce the calculations presented here are available under “Publications” at the Web site http://www.OpenSHA.org/ .

Seismological Research Letters↗