EXPLORE THE ARCHIVE

Faults & earthquakes

Find USGS research on faults & earthquakes. Use publication dates, geographic coverage, abstracts, and original report links to locate studies relevant to a particular place.

215 records · Page 2Linked to original sources

USGS2026

Efficient physics‐informed ground‐motion simulations with reduced‐order models: CyberShake implications and high‐resolution site terms for southern San Andreas fault earthquakes

Recent advances in Probabilistic Seismic Hazard Analysis (PSHA) leverage physics‐based ground‐motion simulations to estimate seismic hazard, such as the CyberShake project. However, computational costs quickly escalate when performing PSHA for numerous faults or sites and can become prohibitively expensive. To reduce computational demands, CyberShake uses reciprocity and interpolates physics‐informed corrections from simulations conducted at fewer locations, but the accuracy of these interpolations remains poorly quantified. To quantify the interpolation accuracy, we derive high‐resolution, frequency‐dependent site terms for southern California and compare them with interpolated site terms using the CyberShake approach. We accomplish this by performing a set of earthquake point‐source simulations distributed along the nonplanar fault geometry for the southern San Andreas fault (SSAF) extending from Bombay Beach to Lake Hughes. Using SeisSol, we simulate three minutes of viscoelastic seismic wave propagation for these sources and store the horizontal‐component Green’s functions for 480,000 sites. We then use a scientific machine learning approach based on interpolated proper orthogonal decomposition to construct an accurate reduced‐order model of the Green’s functions to efficiently predict effective amplitude spectra (EAS) for finite‐source rupture models of SSAF earthquakes. Using minimum curvature interpolation with tension, as used in CyberShake, we compare the interpolated site terms against our high‐resolution site terms. We identify local discrepancies with EAS differing by up to a factor of approximately three. Furthermore, we identify locations where unexpectedly high or low ground motions are missed when using the interpolated dataset for these earthquakes. We estimate that our approach may be used within CyberShake to reduce the time‐to‐solution by a factor of 336 for the entire earthquake rupture forecast. Our analysis of physics‐based site terms provides more insight into the seismic hazard due to SSAF ruptures and guides future developments by combining high‐performance computing and reduced‐order modeling techniques for PSHA.

California↗

USGS2026

High‐resolution surface deformation and slip distribution observations for the 2023 Kahramanmaraş, Türkiye, earthquake sequence help constrain the rupture process

Splay, or branch, faults are a common geometric feature of earthquake surface ruptures and may provide constraints on the rupture behavior of an earthquake. The 2023 M w 7.8 Pazarcık and M w 7.5 Elbistan, Türkiye, earthquakes are examples of ruptures with multiple small splays, and the Pazarcık earthquake nucleated on a splay fault, the Narlı fault, before rupturing bilaterally on the East Anatolian fault (EAF). Here, we present 3‐m‐resolution surface displacement from subpixel correlation of Planet Dove optical images for the entirety of both ruptures with corresponding surface slip distributions. For a 30‐km‐long study region spanning the Narlı‐EAF intersection, we compare surface slip derived from five data sets with different resolutions (on‐the‐ground, WorldView, Planet Dove, Sentinel‐2, and Sentinel‐1) to elucidate complementary information. In addition, we integrate information from the surface expression of faulting with published dynamic rupture simulations and rupture process studies to constrain a rupture evolution for the Pazarcık earthquake that is consistent across data sets. This work highlights the complementary nature of disparate surface slip data sets and the role that high‐resolution surface displacement information, including from fault splays, can play in constraining nonunique rupture models and refining understanding of the earthquake rupture process.

Bulletin of the Seismological Society of America↗

USGS2025

Spatiotemporal variations in strain release and seismic rupture in multifault systems: An example from Panamint Valley, southeastern California

Geometrically complex, multifault ruptures have been observed in recent, damaging earthquakes in southeastern California, sparking renewed efforts to identify physical conditions that promote or inhibit fault discontinuity-spanning coseismic ruptures. The likelihood of ruptures propagating across fault discontinuities is thought to be partly controlled by fault geometries, rupture direction, and the history of strain release. However, these parameters vary in space and time over multiple earthquake cycles, making it difficult to forecast the likelihood that an earthquake on one fault will trigger rupture on a nearby fault. Here we use tectono-geomorphic mapping of a geometrically complex fault zone in Panamint Valley, southeastern California, to assess spatiotemporal variations of paleo-rupture patterns and geometries of fault discontinuities over multiple earthquake cycles. First, we identify ten generations of late Pleistocene to Holocene alluvium using geomorphic parameters and luminescence dating to constrain ages of alluvium and bracket late Holocene earthquake timing. Then, we quantify slip kinematics using high-resolution structure from motion digital surface models. We find the Panamint Valley transtensional relay (PVTR) hosted four late Holocene earthquakes, bracketed to ~5.8–3.4 ka, ~3.8–2.2 ka, ~2.4–0.6 ka, and ~0.64–0.16 ka, with ~0.6–1.1 m of slip per event, correlative to Mw ≈ 6.7–6.9 earthquakes. Additionally, we find similarities in earthquake timing on the Ash Hill, PVTR, and Panamint Valley faults and similarities in the slip magnitude and slip kinematics between the Ash Hill and PVTR faults, implying that the PVTR may co-rupture with nearby faults. Paleo-rupture patterns indicate that seismogenic strain transfer may occur through the PVTR, along different combinations of fault segments and jump distances, over multiple earthquake cycles. These data highlight the utility of tectono-geomorphic mapping in evaluating paleo-rupture patterns and suggest that the PVTR may act to propagate and/or arrest rupture between the Ash Hill and Panamint Valley faults.

California↗

USGS2026

Low-frequency earthquakes track the motion of a captured slab fragment

Accurate tectonic models are essential for assessing seismic hazard and fault interactions. However, the plate configuration at the complex Mendocino triple junction, where the San Andreas Fault and the Cascadia subduction zone meet, remains uncertain. We analyzed fault slip associated with a recently identified zone of tectonic tremor and low-frequency earthquakes (LFEs) near the southern edge of the subducting Gorda slab. Based on tidal sensitivity and P-wave first motions, we show that the LFEs are generated by dipping, strike-slip motion. This suggests that a former Farallon slab fragment, now captured by the Pacific plate, is translating northward beneath westernmost North America. This geometry effectively extends the slab interface fault, challenging prevailing interpretations of slab window formation and creating a potential unaccounted earthquake hazard in this region.

Science↗

USGS2025

Fault displacement model for surface principal rupture of strike-slip faults

The probability distribution model for principal displacement accommodated on the surface main trace is a critical input to the fault displacement hazard analysis. This article presents a new model for strike-slip ruptures in the moment magnitude ( M ) range of 6 to 8.3. The new model is the outcome of a multi-year research effort to update the widely used model developed by Petersen and others in 2011. Updates include the adoption of the Fault Displacement Hazard Initiative database and enhancements to rupture and displacement data preparation. Statistical formulation and estimation have also been updated substantially. A three-parameter modified normal distribution that we refer to as the negative Exponentially Modified Gaussian distribution is adopted to model the probability distribution of the natural logarithm of principal displacement. Formulation for the mean parameter of the modified normal includes a random earthquake term, a nonlinear scaling relation with M , and an ellipse function for along-main-trace variation. The aleatory variability of the updated model now depends on M as well as site’s along-main-trace position. These updates not only significantly improve the fit to the distribution of the observed displacements but also yield reasonable 95th percentile predictions for M > 7.5 events. Alternative models representing the estimation uncertainty of the M -scaling relation are also developed. These new models are compared to the previous model in terms of percentile predictions and the calculated hazard curves. The steeper hazard curves from the new models yield a lower exceedance rate than the normal-distribution based model developed previously by Petersen and others.

Earthquake Spectra↗

USGS2025

Probing lower-crustal fault properties with frequency-dependent tidal tremor triggering

The way seismicity responds to periodic stress perturbations offers crucial insights into the processes that can trigger an earthquake. Laboratory and theoretical analyses have shown that the period of imposed forcing and source properties affect the sensitivity to triggering, but frequency-dependent triggering of tectonic faults is poorly understood. The rate of low-frequency earthquakes (LFEs) near Parkfield, California has been found to be strongly correlated with solid earth tides. Tidal forcing acts over multiple frequencies, and the sensitivity to tidal triggering of LFEs therefore provides a unique opportunity to probe the physics of earthquake triggering and underlying fault properties. Here, we find that the response of LFEs to solid earth tides at diurnal and semi-diurnal frequencies is highly variable but spatially coherent along the San Andreas Fault. Using rate-state friction modeling, we find that the variation of the amplitude of tidal modulation is mainly affected by the spatial variation of the background effective normal stress, whereas the spatially varied tidal modulation at diurnal and semi-diurnal frequencies is mainly affected by the fault frictional property. The spatial variations of the LFE response reveal a heterogeneous lower-crustal geologic structure and complex physical faulting processes below the rupture zone of eventual great earthquakes.

California↗

USGS2026

Dynamic coupling between faulting, rifting and magmatism during 2021-2025 unrest on Reykjanes Peninsula, Iceland

Interactions among faulting, earthquakes, and eruptions are fundamental to plate tectonics and hazard forecasting yet rarely observed along mid-ocean ridges. On Iceland's Reykjanes Peninsula, seismotectonic–volcanic unrest resumed after nearly 800-year hiatus, providing an opportunity to observe these interactions during 2021–2025 activity. By integrating high-resolution seismicity, focal mechanisms, satellite geodesy, surface deformation, and eruption data, we document ∼4 m of total extension accommodated through 14 rifting episodes. The largest, in 2023, involved graben reactivation and diking, with seismic swarms and earthquake faulting that matched the surface ruptures, where strike-slip faulting preceded normal-faulting earthquakes and extension. The accrued extension was released by extension fractures triggered by magma accumulation. Long-term observations show no correlation between erupted magma volume, seismicity, and crustal extension. This highlights dynamic relation between rifting, faulting, and magmatism in transtensional settings and their implications for hazard assessment.

Reykjanes Peninsula↗

USGS2025

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth↗

USGS2026

The U.S. Geological Survey 2025 Puerto Rico and U.S. Virgin Islands time-independent earthquake rupture forecast

We present the 2025 U.S. Geological Survey Puerto Rico and U.S. Virgin Islands (PRVI) time‐independent earthquake rupture forecast (ERF), developed for the 2025 update to the National Seismic Hazard Model (NSHM) for PRVI. The updated ERF improves upon a prior model from 2003, including an expanded fault inventory with slip‐rate estimates, updated seismicity catalogs, and refined subduction zone geometries and deformation models. It applies the fault‐system inversion methodology to solve for rates of ruptures on modeled faults, adapted from the 2023 NSHM (NSHM23) for the western United States, including the first application of the inversion to model rates on a U.S. subduction interface. Off‐fault and intraslab seismicity are constrained by observed seismicity and use updated methods developed for NSHM23. Uncertainties in model components are substantial, and the ERF represents epistemic uncertainties through a comprehensive logic tree consisting of 1.7 billion logic‐tree branches combined across all sources.

Puerto Rico, U.S. Virgin Islands↗

USGS2026

Ground motions from finite faults with unknown geometry

Earthquake emergency response and loss modeling require accurate models of ground‐motion intensity metrics (IMs), such as those produced by ShakeMap. In the critical minutes and hours after an earthquake, however, knowledge of the earthquake source geometry limits the accuracy of these models, especially in areas where seismic instrumentation is insufficient to constrain the shaking intensity pattern. Previous attempts to overcome this limitation have centered on computing the mean distance metrics from a suite of assumed source geometries. That approach is inadequate with modern ground‐motion models (GMMs) because they typically include terms for the hanging wall effect, depth to top of rupture, and other nonlinear terms. Such approaches also propagate the distance uncertainty through the GMMs and may not properly allocate the additional uncertainty into the between‐ and within‐event components. Here we present Finite‐Fault SIMulation with Mixed Effects Regression (FFSIMMER), a method to address these limitations by computing IM distributions from a suite of rupture realizations that fully accounts for all geometry‐related GMM terms and allows for the separation of the between‐ and within‐event components of the total uncertainty. The accuracy of this method is limited by the available constraints on the source parameters. To address this, FFSIMMER can incorporate constraints based on knowledge of the source, such as the approximate strike and dip based on knowledge of the plate tectonics and focal mechanism of the earthquake. We demonstrate the improvement in accuracy that can be achieved with the 2025 M 7.7 Mandalay, Burma (Myanmar), earthquake by incorporating constraints on the range of plausible source parameters.

Bulletin of the Seismological Society of America↗

USGS2025

Evidence for fluid pressurization of fault zones and persistent sensitivity to injection rate beneath the Raton Basin

Subsurface wastewater injection has increased the seismicity rate within the Raton Basin over more than two decades, with the basin-wide injection rate peaked between 2009-2015. To understand the evolution of injection-induced earthquakes, we systematically analyzed 2016-2024 broadband recordings with a machine-learning-based phase picker and constructed a catalog with 95,993 earthquakes (-1≤ M L ≤4.3). We then inverted for full centroid moment tensors (CMT) for 90 M L ≥ 2 events, with a special interest in constraining the non-double-couple components via probabilistic metrics. Both relocations and CMT solutions support basement-rooted normal faults, including graben and half-graben structures. Furthermore, we observe the non-double-couple components that imply elevated pore pressure in the fault zones. An earthquake cluster emerged in the north-central basin in 2023, preceded by ~1-yr of increased injection volume from wells within 15km. Despite a basin-wide decrease in the injection volume, we highlights the persistence of seismicity that remains to sensitive to injection rates within the Raton Basin.

Colorado, New Mexico↗

USGS2025

Induced earthquakes are generally not tidally triggered in Oklahoma and Kansas

Human-induced earthquakes occur along critically stressed faults as injected wastewater simultaneously heightens fluid pressure and pushes faults to failure. We investigate the possibility that small stresses imposed by Earth tides could trigger earthquakes in the induced seismicity region of Oklahoma and Kansas from 2011 to 2018. We decluster a catalog consisting of ∼110,000 earthquakes using three methods (Reasenberg, nearest-neighbor distance, and phase-bin). We find no significant tidal earthquake triggering using Schuster's p -value test for the declustered catalogs as a whole. We search for localized triggering using discretized space-time cells and find ∼0–6% of cells have significant tidal triggering which is close to what is randomly expected (5%) and indicates there is an insignificant amount of tidal triggering for the full study region. One area that has significant p -values across multiple time windows, ∼2014–2016 is ∼15 km from a region of large wastewater injection volume. It is possible that localized tidal triggering occurs for this time and area because faults remain critically stressed and are particularly susceptible to slip under the small stress load from the semidiurnal tide. Possible explanations for the lack of tidal triggering in our broader study are that the pre-seismic stressing rate in the earthquake nucleation area is faster than the tidal stressing rate (∼3 kPa/day), faults are not close enough to critically stressed to be affected by tidal forcing, and that nucleation occurs over longer periods than the tides considered in this study (∼1, ∼14 days). Fluid injection could be the source of a higher pre-seismic stress rate.

Kansas, Oklahoma↗

USGS2025

Revisiting an enigma on California's north coast: The Mw6.5 Fickle Hill earthquake of 21 December 1954

Many earthquakes occur along the North Coast of California in the vicinity of the Mendocino Triple Junction (MTJ), where the Pacific, Gorda, and North American (NA) plates meet, and on the adjacent plate boundaries. The MTJ marks the nexus of the Mendocino and San Andreas faults with the Cascadia subduction zone (CSZ). Historically, most large earthquakes around the MTJ have been within the offshore Gorda plate and its subducted portion beneath the NA plate. North of the MTJ, active faults mapped in the NA plate are part of the CSZ fold‐and‐thrust belt. Although some events have been detected in the NA plate, no large historic events have been associated with mapped surface faults. The 21 December 1954 M w 6.5 earthquake in Humboldt County is one possible exception. Using published data from catalogs and articles, unpublished data from Berkeley’s archives, and S ‐ P times interpreted from two U.S. Coast and Geodetic Survey (USCGS) accelerometers, we determine a probability cloud for the earthquake’s hypocenter using NonLinLoc. The highest probability location lies beneath Fickle Hill just east of the city of Arcata, California, at 40.87° N, 124.03° W, and ∼11 km depth. Using P ‐wave polarities from Berkeley stations and the digitized waveforms from the accelerometers, we find that the focal mechanism most consistent with the data indicates thrust movement with strike, dip, and rake of 350°, 10°, and 90°, respectively, at a depth of 14 km. Given the depth uncertainties of both this event and the megathrust, this implies that the earthquake most likely took place on the subduction interface rather than on the mapped faults in the Mad River fault zone that trend 322° and dip to the northeast. The revisited intensity in the epicentral region also supports a location beneath Fickle Hill to the east of the city of Arcata, California.

California, Idaho, Nevada, Oregon, Washington↗

USGS2025

To heal or not to heal?: 2. The moment-recurrence time behavior of repeating earthquakes in the 2011 Prague, Oklahoma aftershock sequence is consistent with laboratory healing rates

The timing and failure conditions of an earthquake are governed by the interplay between fault reloading and restrengthening. The moment-recurrence time behavior of repeating earthquakes can give observational estimates of fault healing rates; however, it is difficult to link these observed healing rates to laboratory studies of frictional healing in part because of uncertainty in lithology. Here, we study the 2011 Prague earthquake sequence, which includes repeating earthquakes in the Arbuckle group and the granitic basement, and compare them to laboratory experiments on samples of the Arbuckle and Troy granite (representative of the basement rock) (Okamoto et al., 2025, https://doi.org/10.1029/2024JB030573 ). We find three spatially distinct groups of repeating earthquakes with different moment-recurrence behavior: (a) constant moment-recurrence time in the Arbuckle group, (b) scattered moment-recurrence time at the intersection of the foreshock-mainshock fault in the granitic basement, and (c) moment-predictable behavior outside of the foreshock-mainshock fault intersection also in the granitic basement. Our observation of stagnant healing for repeating sequences in the Arbuckle group is consistent with laboratory observations of low healing rates for moderately high pore fluid pressures in Arbuckle samples. For the moment-predictable group, the source radius that is required in order to match healing rates is consistent with source radius estimations when taking into account reasonable attenuation of the P- pulse width. Overall, we observe diverse healing behaviors in the seismic families that are consistent with laboratory healing rates, providing seismic evidence that contact-scale frictional mechanisms are relevant to large-scale earthquake dynamics.

Oklahoma↗

USGS2026

Temporal and spatial changes in seismic attenuation associated with inferred fluid migration in the 2016 central Apennines earthquake sequence

Prior work suggests that high‐frequency seismic attenuation acts as a highly sensitive proxy for crustal permeability and fluid mobility in fractured media. We test the hypothesis that the fault system responsible for the 2016–2017 Amatrice–Visso–Norcia–Capitignano sequence acted as an impermeable seal, compartmentalizing pressurized fluids until dynamic rupture triggered widespread fluid diffusion. By tracking across the sequence the spatiotemporal evolution of the S ‐wave anelastic attenuation parameter, we identify large, positive low‐frequency attenuation anomalies emerging within the hanging wall following the Amatrice mainshock and strictly preceding subsequent large ruptures. Conversely, we observe weaker, negative anomalies in the footwall, anticorrelated in time with those of the hanging wall, revealing a massive asymmetry in fluid redistribution and permeability evolution across the fault system. Furthermore, aftershock migration rates reveal distinct linear alignments in a distance‐reduced time space, allowing us to explicitly track and quantify episodes of lateral and upward fluid migration. These physically consistent patterns suggest that stress‐driven fluid diffusion directly weakens adjacent fault patches, dictating the spatiotemporal migration of seismicity. We conclude that near‐real‐time monitoring of seismic attenuation may help detect fluid redistribution in active fault systems and may provide useful information for time‐dependent seismic hazard assessment.

central Apennines↗

USGS2026

Surface rupture and slip distribution of the 2025 Mw7.7 Mandalay earthquake and updated length scaling of supershear earthquakes

The 2025 M w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes ( M w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log 10 (surface rupture length) = 0.89 M w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.

Mandalay↗

USGS2026

Distributed faulting of the northern West Napa Fault Zone in Napa Valley, California

Mapped surface ruptures from the 24 August 2014 M w 6.0 South Napa earthquake in the Napa Valley, California, show a 2‐km‐wide zone of distributed faulting in the southern and central West Napa fault zone (WNFZ). In the northern WNFZ at Hendry Winery (HW), however, the mapped 2014 surface ruptures encompass an ∼100‐m‐wide zone, implying significant narrowing of the near‐surface fault zone to the north. We present a tomographic shear‐wave velocity ( ⁠⁠ V S ) model and guided‐wave data that indicate the northern WNFZ is at least 400‐m wide, with multiple near‐surface fault traces. Our V S model shows that the 2014 surface ruptures are underlain by discrete low‐velocity zones (LVZs), and coincident guided‐wave data show that the LVZs carry fault‐zone guided waves. If nearby (<500 m) mapped faults to the east of HW are part of the WNFZ, the entire WNFZ is more than 1 km wide in the northern Napa Valley. WNFZ guided waves travel up to 38% slower than S body waves, and low‐strain guided‐wave shaking is up to five times stronger than the associated body‐wave shaking. Our data suggest that guided waves, traveling along distributed faults, may result in an increased shaking hazard over a 1‐km‐wide area of the northern Napa Valley during future significant earthquakes. In places, the 2014 surface ruptures were difficult to find one year after the earthquake, and paleoseismic trenching showed only weak evidence for faulting, which may not have been identified in trenches if the locations of the 2014 surface ruptures had not been previously mapped ( Prentice et al. , 2015 ). Guided‐wave and V S tomography data, however, show strong evidence for faulting beneath the 2014 surface ruptures and at locations to the east. Although paleoseismic trenching is the gold standard for identifying near‐surface faulting, methods such as peak ground velocities of guided waves may better identify immature near‐surface fault traces.

California↗

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↗