Search USGS⌕ Search

SEARCH · Search USGS

Results for “Geochemical Perspectives Letters”

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.

235 records · Page 6Linked to original sources
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

Special issue “The 2024 M7.6 Noto Peninsula earthquake and seismic swarm”

The 2024 Mw7.5 (M JMA 7.6) Noto Peninsula earthquake struck the Noto Peninsula, central Japan on January 1, 2024, causing widespread damage. This Mw7.5 mainshock, together with the prolonged seismic swarm that persisted in the northeastern Noto Peninsula since the end of 2020, of which the largest event was the 2023 M JMA 6.5 earthquake, constitutes one of the most significant seismic sequences in Japan in recent decades. The sequence provides a rare opportunity to investigate the complex interplay among swarm activity, fluid migration, fault geometry and distribution, and the preparatory conditions required for a large inland crustal earthquake. Motivated by this scientific context, this special issue brings together studies examining the sequence from diverse disciplinary perspectives, incorporating seismic, geodetic, and tsunami-related observations and modelling.

Noto Peninsula↗
USGS2026

Dynamic drainage reorganization in Eastern Tibet: Insights from the Yangtze River first bend

The modern drainage network of eastern Tibet is widely believed to have developed through a series of river capture and flow reversal events; however, the timing and mechanisms driving this reorganization remain contentious. Among these events, the river capture that formed the First Bend of the Yangtze River (YFB) stands out as both iconic and particularly debated. Here we present sedimentary provenance data from the Late Miocene–Quaternary Dali Basin, located south of the YFB, which indicate that a southward-flowing Jinsha River (i.e., the present-day upper Yangtze River) sourced sediment to the Dali basin at ∼7.4–6.4 Ma in a drainage configuration different from that of today. Because this interval postdates the initial establishment of a near-modern Jinsha River system prior to the Miocene, our results imply at least two discrete fluvial reorganizations occurred at the YFB—one preceding ∼7.4 Ma and another following ∼6.4 Ma. By integrating these findings with landscape evolution modeling, we infer that the initiation of rapid uplift of the Yulong-Haba Mountains and the Diancang Shan may have been responsible for these drainage reorganizations. These results underscore that Cenozoic drainage systems on the eastern Tibetan Plateau have evolved dynamically on a short timescale of ∼10 5 –10 6 -year, rather than remaining in a long-term stationary configuration on ∼10 7 -year timescales.

eastern Tibetan Plateau, first bend of the Yangtze↗
USGS2026

Corundum discovered by SuperCam and the Perseverance rover at Jezero crater, Mars

Mars is primarily composed of mafic mineral assemblages and their alteration products, but small, scattered rocks strewn across the landscape offer clues to greater petrological diversity. While traversing the Jezero crater rim, the Perseverance rover encountered several plagioclase-rich light-toned float rocks. SuperCam identified the distinctive signature of corundum (α-Al 2 O 3 ) in these rocks using time-resolved luminescence spectroscopy. Two strong peaks (692.7 and 694.1 nm) with millisecond lifetimes, and additional supporting lines, are consistent with Cr 3+ substitution for Al 3+ in corundum. Corundum forms in Al-rich, Si-depleted environments through magmatic or metamorphic processes. Given the rocks' small size, association with plagioclase, and location on the crater rim, we interpret the most plausible formation scenario to be impact induced metamorphism at the interface of a felsic and a mafic/ultramafic member with the likely action of fluids at some stage, although other possibilities are not excluded.

Geophysical Research Letters↗
USGS2024

Onset of aftershocks: Constraints on the Rate-and-State model

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R (t)=K(c+t) −p ⁠ , with parameters K , c , and p . The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c ‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced stresses strongly vary in space, the c ‐value should vary accordingly. Short‐time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a ‐positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c ‐value. The onset of the power‐law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ <10 kPa and a power‐law decay of the background rate with distance to the fault, with exponent 2.7.

Seismological Research Letters↗
USGS2024

Testing rate‐and‐state predictions of aftershock decay with distance

We analyze aftershocks of the 2019 M 7.1 Ridgecrest mainshock and isolated M 5–6 mainshocks in southern California to test predictions made by the rate‐and‐state friction model of Dieterich (1994) . Rate‐and‐state friction predicts that the seismicity rate after a stress step follows Omori decay, where the Omori c ‐value, which is the saturation in aftershock rate observed at small times, is larger for smaller stress steps. Put in the context of an aftershock sequence, this predicts that the Omori c ‐value will be systematically larger at greater distances from the mainshock. To our knowledge, this predicted effect has not been observed. In part this may be because the Omori c ‐value is difficult to measure because it often reflects short‐term catalog incompleteness rather than a true saturation in aftershock rate. We explore the dependence of the Omori c ‐value on the distance to the mainshock by applying the “ a ‐positive” method ( van der Elst and Page, 2023 ). This method is insensitive to short‐term aftershock incompleteness and allows resolution of the true aftershock rate deep into the mainshock coda. For aftershocks of the Ridgecrest mainshock and stacked M 5–6 mainshocks, we observe systematic differences in early aftershock rates, relative to mainshock distance, consistent with the predictions of rate‐and‐state friction. Furthermore, for the larger Ridgecrest dataset, we observe that aftershocks nearer to the mainshock start earlier, and we resolve a flattening of the Omori curve consistent with a larger Omori c ‐value for the farthest aftershocks, as predicted by Dieterich (1994) .

Seismological Research Letters↗
USGS2026

Moment magnitude for small earthquakes in the Delaware basin of west Texas and southeast New Mexico, USA

The Delaware Basin region of west Texas and southeast New Mexico has become one of the most prolific regions of seismic activity in the continental United States due to widespread hydraulic fracturing and wastewater disposal injection. In response to the increased number of earthquakes in this region, rapid and accurate characterization of earthquake sources is necessary to understand the evolution of seismic activity and level of seismic hazard associated with these earthquakes. This study re-evaluates earthquake magnitudes, estimating moment magnitude (MW) for small earthquakes in the Delaware Basin using 1) moment-rate spectra derived from S-wave coda envelopes, and 2) a relative magnitude method that relies exclusively on the ratio of waveform amplitudes between highly correlated waveform pairs. The coda-envelope method produces accurate M W estimates for small earthquakes ( M 1.5 – 3) that are consistent with independent, waveform modeled moment magnitudes for events with M W > 3 . Using the relative amplitudes method to extend these M W magnitudes to many other events, we successfully provide relative moment magnitude ( M W,rel ) values for 81% of the Texas Seismological Network catalog in the Delaware Basin region, and 45% of the USGS Induced Seismicity Project’s catalog of events in southeast New Mexico. The adoption and integration of the calibrated M W,rel method with current magnitude estimation methods offers valuable insights into the relationships between local and moment magnitude and will contribute to improved characterization of widespread induced seismicity.

New Mexico, Texas↗
USGS2025

Updating regional‐scale geospatial liquefaction models with locally available geotechnical data

We present a method to update the geospatial liquefaction model used by the U.S. Geological Survey’s near‐real‐time ground failure product with subsurface geotechnical data. The geospatial model estimates liquefaction probability from peak ground velocity (via ShakeMap) and geospatial susceptibility proxies. In many regions, additional information relevant to constraining liquefaction likelihood is also available, including surface geology maps and subsurface geotechnical measurements. There is currently no mechanism to use these data in the ground failure product liquefaction model, even though these data could provide more precise constraints on spatial variations in the lithologic character of the soil (surface geology) and direct measurements of the subsurface mechanical properties that affect liquefaction occurrence and severity (geotechnical measurements). In this study, we develop a method to integrate these data with the geospatial model and assess how these data can improve regional‐scale predictions. We develop a Bayesian updating framework and apply it to the 1989 magnitude 6.9 Loma Prieta, California, earthquake, for which mapped observations are available to evaluate performance. We constrain the Bayesian framework with 373 Northern California cone penetration tests and liquefaction susceptibility classes based on the mapped surface geology. This Bayesian model incorporates geotechnical information into the geospatial model and more accurately predicts liquefaction occurrences than the geospatial model, while sacrificing less accuracy in terms of predicting the absence of liquefaction than the geotechnical model. In future applications, this approach could be adapted to update other geospatial models using locally available subsurface data.

California↗
USGS2024

Are field observations of surface rupture useful? An example from the 2023 Mw 7.8 Pazarcık, Turkey (Türkiye), earthquake

Field investigations have long been an important component of the scientific response to surface‐faulting earthquakes. However, in light of advances in remote data and models, the question arises whether field‐based observations of surface rupture remain useful for understanding rupture processes and seismic hazards. We approach this question using a field‐based study of the central 2023 M w 7.8 Pazarcık, Turkey (Türkiye), earthquake rupture, at the intersection of the east Anatolian fault (EAF) and Narlı fault. Our field observations include the surface rupture expression and extent of the central EAF and northernmost Narlı fault in generally forested and steep terrain and 68 measurements of left‐lateral surface displacement. These data improve our understanding of the Pazarcık rupture complexity, resolve the surface geometry of the Narlı fault–EAF connection, and exhibit a clear (>2 m) change in surface displacement across this intersection zone that confirms remote‐based coseismic slip models. Our study shows that focusing field efforts in areas of obscured or low‐resolution remote data can yield essential data for refining rupture extent, documenting perishable on‐fault displacement, and improving postearthquake situational awareness. A comparison of similarly large‐magnitude continental surface‐rupturing earthquakes indicates that displacement uncertainties relate to a complex set of factors, including measurement methods, rupture complexity, and displacement magnitude. Our study validates the need for postearthquake field observations, which, when driven by clear motivating questions and knowledge of methodological strengths and limitations, provide high‐resolution rupture data that complement remote‐based models.

Pazarcık↗
USGS2024

Seismicity zoning at Coso geothermal field and stress changes from fluid production and migration

The Coso geothermal field is a major geothermal power production site in the western United States. It has been observed that low-magnitude seismic events ( M < 3.71) are unevenly distributed in three distinct zones, namely, nearfield (<3 km), midfield (3–6 km), and farfield (> 6 km) from the Coso geothermal plant. These zones exhibit distinct changes in earthquake location before and during geothermal production episodes that began in 1986. After 1986, the midfield region of the main flank experiences a significantly lower seismicity rate than the surrounding areas before production episodes. During 2014–2019, the farfield earthquakes cluster in the eastern and western parts of the greater Coso area, which is discernably different from how those pre-production earthquake events were distributed along the conjugate NW-SE and SW-NW trending structures across the main flank. Here, we analyze the stage of stress with finite-element-based poroelastic simulations to illustrate how the spatiotemporal evolution of the seismicity is associated with the pattern of stress perturbations caused by fluid migration amid the operations of geothermal power plants. Generally, ∼70% of co-production seismicity is found in zones of increased Coulomb stress between 2014 and 2019 at >99% confidence. Meanwhile, the midfield zone of seismic paucity overlaps with the zone of decreasing pore-fluid pressure. Overall, the results provide a physical explanation of how decadal geothermal operations at Coso have perturbed stress-field changes and contributed to the evolving characteristic seismic pattern, shedding insights into assessing the seismic hazard in other geothermal settings.

California↗
USGS2026

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
USGS2025

Searching for seismic precursors - The Barry Landslide hazard clean up

The Barry Landslide, located in Barry Arm of Prince William Sound, Alaska, poses a major hazard due to its steep, unstable slopes and the potential for a massive landslide-generated tsunami. With an estimated volume of 500–700 million cubic meters, the Barry Landslide could trigger highly destructive waves. In this study, we focus on seismic signals from the Barry Landslide, which are critical for providing timely tsunami warnings. Since the summer of 2020, the region has been instrumented to monitor the landslide, but the seismic record is complicated by the presence of nearby glaciers and frequent regional earthquakes. Among these signals, we analyze a specific class of short-duration, high-frequency seismic events that exhibit strong seasonal variability, increasing in rate from late summer to mid-winter before ceasing abruptly in late winter or early spring. Our analysis suggests that the source of these signals is likely near or beneath Cascade Glacier, adjacent to the landslide, rather than within the landslide mass itself. We apply detection algorithms to construct a time history for this signal type, which we then compare with environmental factors like precipitation, temperature, and slope displacement data from ground-based radar and remote sensing. Correlations indicate that these seismic events may be driven by seasonal hydrological changes, particularly the freeze-up of subglacial water pathways. While these events are not directly linked to landslide motion, they serve as indirect markers of subsurface hydrological conditions that influence slope stability. Our findings highlight the complex interplay between glaciers, groundwater, and landslide dynamics, emphasizing the need for multi-parameter monitoring to assess evolving geohazards in the region.

Alaska↗
USGS2026

Making many out of one: Synthetic geologic deformation model distributions for use in USGS NSHM25‐PRVI Puerto Rico-U.S. Virgin Island update

A key use‐case of geologic slip rates is within deformation models used in probabilistic seismic hazard analyses. Field‐derived geologic slip rates have formed the cornerstone of deformation models in such applications for decades. Recent advancements in seismic hazard analyses have expanded the use of faults for which geologic slip rates are not well constrained using categorical slip rate estimates. Because of these advancements, application of a geologic deformation model for use in 2025 U.S. Geological Survey National Seismic Hazard Model Puerto Rico‐U.S. Virgin Islands (NSHM25‐PRVI) proved challenging due to: (1) a lack of field‐based geologic slip rates, and (2) a lack of epistemic uncertainty distributions within a broad range of estimated slip rates. Preliminary versions of the NSHM25‐PRVI model sampled these slip rate bins in a coincident manner along preferred and extreme value branches, which yielded untenable correlations in mean hazard results. To minimize the influence of correlated uncertainties amid these challenges, we develop a synthetic epistemic uncertainty distribution for deformation rate on each crustal fault. Each fault has a weighting schema across four possible distribution shapes: uniform, normal, triangular favoring local minima, and triangular favoring local maxima. The synthetic distributions are then sampled several times for each logic tree branch. The results provide a more realistic distribution of rates across the study region as compared with using correlated extrema sampling. This exploration of our method in a small region like PRVI can pave the way for larger‐scale, more complicated applications (e.g., western United States).

Puerto Rico, US Virgin Islands↗
USGS2026

Recycling radiogenic osmium by crustal foundering in subduction zones: Evidence from pyroxenite xenoliths in the north Andean arc

Removal of mafic-ultramafic lower crust (e.g., via delamination) is fundamental to making andesitic continental crust, yet direct evidence of this process remains elusive. A unique suite of garnet clinopyroxenite and hornblendite (arclogite) xenoliths from the northern volcanic zone of the Andes, erupted to the surface in the mid-Pleistocene, have bulk-rock geochemistry and osmium (Os) isotopic compositions that demonstrate their crustal affinity, yet have equilibration pressures and temperatures below the arc Mohorovičić discontinuity (~53 km) and some as deep as within the sub-arc mantle wedge (~105 km). Garnet websterites from the same xenolith suite, sourced from depths approaching the Wadati-Benioff zone (~140 km), have elemental and isotopic compositions indicative of a mantle origin, likely formed as products of peridotite-liquid reactions above the subducting slab. Variability in bulk-rock 143 Nd/ 144 Nd values and garnet oxygen isotope ratios ( 𝛿 18 O) for these samples is attributed to minor assimilation of subducted components and/or older crustal material, but assimilation alone cannot account for the highly radiogenic Os isotopic compositions. These results provide direct petrologic evidence for modern arclogite formation and foundering in the Andean orogen, the archetypal active continental subduction system. Elevated 187 Os/ 188 Os values of Mercaderes arclogites imply that recycling of gravitationally unstable arc cumulates during continental crust formation introduces highly radiogenic Os into the convective mantle, which has implications for mass fluxes across the crust-mantle boundary as well as the sources and evolution of mantle heterogeneity as seen in orogenic peridotite massifs and in the sources of oceanic basalts

Earth and Planetary Science Letters↗
USGS2026

Earthquake ground-motion model adjustments for the San Francisco Bay area

We develop adjustments to ergodic ground‐motion models (GMMs) to improve their performance in the San Francisco Bay Area (SFBA). GMMs are widely used in hazard assessments to estimate characteristics of ground shaking based on known properties of the source, path, and site. Such models are often developed using datasets containing records from various regions, resulting in models that represent median ground‐motion behavior, which may not adequately represent ground motions within subregions. This is true for the SFBA, where ground motions attenuate more rapidly with distance than in many other parts of California that dominate GMM databases. To support improved seismic hazard estimates in the SFBA, we calculate regional constants and anelastic attenuation coefficient adjustments relative to two commonly used ergodic GMMs: BSSA14 ( Boore et al. , 2014 ) and ASK14 ( Abrahamson et al. , 2014 ). These adjustments are obtained for a suite of ground‐motion intensity measures (peak ground acceleration, peak ground velocity, and 5%‐damped pseudospectral acceleration at oscillator periods ranging from 0.075 to 10 s) using mixed‐effects regression. Use of the regionally adjusted models reduces the overall bias by up to 0.5 natural log units for BSSA14 and up to 0.6 natural log units for ASK14. We demonstrate one application of our attenuation adjustments and their implications in an earthquake early warning case study of the 2014 M 6.0 South Napa earthquake. The predicted extent of shaking using the adjusted models better matches observed shaking at large source‐to‐site distances, especially for lower shaking intensities, thus potentially reducing overalerting. We encourage the use of our model adjustments when ergodic models are considered for seismic hazard studies in the SFBA.

California↗
USGS2024

Agricultural tile drains increase the susceptibility of streams to longer and more intense streamflow droughts

Streamflow droughts are receiving increased attention worldwide due to their impact on the environment and economy. One region of concern is the Midwestern United States, whose agricultural productivity depends on subsurface pipes known as tile drains to improve trafficability and soil conditions for crop growth. Tile drains accomplish this by rapidly transporting surplus soil moisture and shallow groundwater from fields, resulting in reduced watershed storage. However, no work has previously examined the connection between tile drainage and streamflow drought. Here, we pose the question: does the extent of watershed-level tile drainage lead to an increased susceptibly and magnitude of streamflow droughts? To answer this, we use daily streamflow data for 122 watersheds throughout the Midwestern United States to quantify streamflow drought duration, frequency, and intensity. Using spatial multiple regression models, we find that agricultural tile drainage generates statistically significant ( p < 0.05) increases in streamflow drought duration and intensity while significantly reducing drought frequency. The magnitude of the effect of tile drainage on streamflow drought characteristics is similar to that of water table depth and precipitation seasonality, both of which are known to influence streamflow droughts. Furthermore, projected changes in regional precipitation characteristics will likely drive the installation of additional tile drainage. We find that for each 10% increase in tile-drained watershed area, streamflow drought duration and intensity increase by 0.03 d and 12%, respectively, while frequency decreases by 0.10 events/year. Such increases in tile drainage may lead to more severe streamflow droughts and have a detrimental effect on the socio-environmental usage of streams throughout the Midwest.

Illinois, Indiana, Iowa, Michigan, Minnesota, Ohio↗
USGS2025

Reflections on a trio of North American earthquakes in 1925

In 1925, three moderately large damaging earthquakes occurred in North America over four months: the 28 February (local time; LT) M 6.2 Charlevoix, 27 June (LT) M 6.6 Montana, and 29 June M 6.5 Santa Barbara earthquakes. The centennial anniversaries of these events motivated this retrospective consideration focused on the ground motions generated by the three events, including a reconsideration of early intensity assignments for the Montana earthquake. At the time, these three earthquakes appeared to support the arguments of some geologists who downplayed the severity of seismic hazard in southern California relative to other parts of the country. Some of the arguments advanced at that time, for example that Los Angeles “has the least to fear from ‘Acts of God’ of any city under the American flag,” ( Hill, 1928 ) sound naïve if not laughable now, but a comparison of well‐constrained shaking distributions for the three earthquakes reveals the dramatic difference in wave propagation efficiency in western versus eastern North America (ENAM), which leads to moderate ENAM events being felt to much larger distances. At M 6.2, the 1925 Charlevoix earthquake was a notably large event in ENAM. This earthquake was the largest event in eastern Canada since 1870 and caused damage in the epicentral region in addition to towns as far away as 200 km, with felt shaking extending over 1000 km. In contrast, felt shaking from the Santa Barbara earthquake barely extended beyond ∼200 km. Compiling published intensity distributions for larger ENAM earthquakes, we show that perceptible earthquake shaking is not uncommon in ENAM over century time scales, but experience with weakly felt shaking may incline people to downplay potential earthquake risk.

California, Montana, Quebec↗
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↗