Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismologic observations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

70 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The geometry of signal space: a case study of direct mapping between seismic signals and event distribution

SUMMARY Under favourable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal-space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection and scaling, with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning and possible application to studies of repeating events in induced, volcanic and glacial seismicity. Our results place limits on the widely held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗
DOE OSTI2022

Global centroid moment tensor solutions in a heterogeneous earth: the CMT3D catalogue

SUMMARY For over 40 yr, the global centroid-moment tensor (GCMT) project has determined location and source parameters for globally recorded earthquakes larger than magnitude 5.0. The GCMT database remains a trusted staple for the geophysical community. Its point-source moment-tensor solutions are the result of inversions that model long-period observed seismic waveforms via normal-mode summation for a 1-D reference earth model, augmented by path corrections to capture 3-D variations in surface wave phase speeds, and to account for crustal structure. While this methodology remains essentially unchanged for the ongoing GCMT catalogue, source inversions based on waveform modelling in low-resolution 3-D earth models have revealed small but persistent biases in the standard modelling approach. Keeping pace with the increased capacity and demands of global tomography requires a revised catalogue of centroid-moment tensors (CMT), automatically and reproducibly computed using Green's functions from a state-of-the-art 3-D earth model. In this paper, we modify the current procedure for the full-waveform inversion of seismic traces for the six moment-tensor parameters, centroid latitude, longitude, depth and centroid time of global earthquakes. We take the GCMT solutions as a point of departure but update them to account for the effects of a heterogeneous earth, using the global 3-D wave speed model GLAD-M25. We generate synthetic seismograms from Green's functions computed by the spectral-element method in the 3-D model, select observed seismic data and remove their instrument response, process synthetic and observed data, select segments of observed and synthetic data based on similarity, and invert for new model parameters of the earthquake’s centroid location, time and moment tensor. The events in our new, preliminary database containing 9382 global event solutions, called CMT3D for ‘3-D centroid-moment tensors’, are on average 4 km shallower, about 1 s earlier, about 5 per cent larger in scalar moment, and more double-couple in nature than in the GCMT catalogue. We discuss in detail the geographical and statistical distributions of the updated solutions, and place them in the context of earlier work. We plan to disseminate our CMT3D solutions via the online ShakeMovie platform.

58 GEOSCIENCES↗
DOE OSTI2021

LL18-Signal Analysis-NDD2Ad (Final Report)

Under favorable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially-distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection, and scaling with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning, and possible application to studies of repeating events in induced, volcanic, and glacial seismicity. Our results place limits on the widely-held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗
DOE OSTI2020

A viscoelastic model for seismic attenuation using fractal mechanical networks

SUMMARY Seismic attenuation (quantified by the quality factor Q) has a significant impact on the seismic waveforms, especially in the fluid-saturated rocks. This dissipative process can be phenomenologically represented by viscoelastic models. Previous seismological studies show that the Q value of Earth media exhibits a nearly frequency-independent behaviour (often referred to as constant-Q in literature) in the seismic frequency range. Such attenuation can be described by the mathematical Kjartansson constant-Q model, which lacks of a physical representation in the viscoelastic sense. Inspired by the fractal nature of the pore fluid distribution in patchy-saturated rocks, here we propose two fractal mechanical network (FMN) models, that is, a fractal tree model and a quasi-fractal ladder model, to phenomenologically represent the frequency-independent Q behaviour. As with the classic viscoelastic models, the FMN models are composed of mechanical elements (spring and dashpots) arranged in different hierarchical patterns. A particular parametrization of each model can produce the same complex modulus as in the Kjartansson model, which leads to the constant-Q. Applying the theory to several typical rock samples, we find that the seismic attenuation signature of these rocks can be accurately represented by either one of the FMN models. Besides, we demonstrate that the ladder model in particular exhibits the realistic multiscale fractal structure of the saturated rocks. Therefore, the FMN models as a proxy could provide a new way to estimate the microscopic rock structure property from macroscopic seismic attenuation observation.

Geochemistry & Geophysics↗
DOE OSTI2025

Observations and seismoacoustic simulations of earthquake-generated infrasound waves in non-epicentral regions

We analysed infrasound waves associated with the Gyeongju earthquake (ML 5.8) that occurred on 2016 September 12, in the southeastern Korean Peninsula. For infrasound wave detection, the Progressive Multichannel Correlation method was applied to the infrasound data set recorded at seven arrays operating in South Korea at epicentral distances ranging from 178 to 472 km. Based on the back-projection method constrained by array-dependent celerity and azimuth deviation models, the source regions were identified in both the epicentral and non-epicentral regions. Remarkably, the non-epicentral secondary sources of this earthquake were located in regions with shallow water depths: (i) the western coastal area in the Yellow Sea and (ii) the shallow ocean basin and bank in the East Sea. The location results obtained from the earthquake could be corroborated through its foreshock (ML 5.1), yielding location results consistent with those of the main shock. The generation of infrasound waves over shallow water depths was fortuitously validated by direct recordings of dominant single-frequency (~0.3 Hz) infrasound waves at close range via temporary sensors near the ocean basin and bank. We interpreted that low-frequency infrasound signals could be generated from interactions among the ocean floor, shallow seawater and atmosphere. We performed numerical simulations of seismoacoustic fields to predict ground motions on the seafloor and acoustic transmission efficiency between the water and air interface. The simulations quantified the energy transfer through different media and clarified our observational results. We found that because this solid Earth‒water‒atmosphere coupled air wave has a relatively low frequency (~0.3 Hz), it can survive propagation over long distances compared with high-frequency infrasound waves generated in inland and mountain regions. In this study, we extend our understanding of water‒atmosphere coupling and the monitoring framework for earthquake-associated non-epicentral infrasound waves, encompassing not only inland ground shaking but also shallow sea regions located far from the epicentre.

Earthquake monitoring and test-ban treaty verifica↗
DOE OSTI2024

Physics informed neural network can retrieve rate and state friction parameters from acoustic monitoring of laboratory stick-slip experiments

Various machine learning (ML) and deep learning (DL) techniques have been recently applied to the forecasting of laboratory earthquakes from friction experiments. The magnitude and timing of shear failures in stick-slip cycles are predicted using features extracted from the recorded ultrasonic or acoustic emission (AE) signals. In addition, the Rate and State Friction (RSF) constitutive laws are extensively used to model the frictional behavior of faults. In this work, we use data from shear experiments coupled with passive acoustic (variance, kurtosis, and AE rate) interleaved with active source ultrasonic monitoring (transmitted wave amplitude) to develop physics-informed neural network (PINN) models incorporating the RSF law and AE rate generation equation with wave amplitude serving as a proxy for friction state variable. This PINN framework allows learning RSF parameters from stick-slip experiments rather than measuring them through a series of velocity step experiments. We observe that when the stick-slip cycles are irregular, the PINN models outperform the data-driven DL models. Transfer learning (TL) PINN models are also developed by pre-training on data collected at one normal stress level followed by forecasting shear failures and retrieving RSF parameters at other stress levels (i.e., with different recurrence intervals) after retraining on a limited amount of new data. Our findings suggest that TL models perform better compared to standalone models. Both standalone and TL PINN-estimated RSF parameters and their ground truth values show excellent agreements thus demonstrating that RSF parameters can be retrieved from laboratory stick-slip experiments using the corresponding acoustic data and that the transmitted wave amplitude provides a good representation of the evolving frictional state during stick-slips.

58 GEOSCIENCES↗
DOE OSTI2022

Comparisons Between Array Derived Dynamic Strain Rate (ADDS) and Fiber‐Optic Distributed Acoustic Sensing (DAS) Strain Rate

Abstract Distributed acoustic sensing (DAS) strain rate and particle velocity can be compared through approximate scaling with medium velocity. We instead performed a direct comparison between array derived dynamic strain (ADDS) rate and DAS strain rate for six frequency bands. The PoroTomo project at Brady's Hot Springs, Nevada, deployed a 240‐geophone 3C array co‐located with fiber‐optic DAS system and 8.7 km of buried cable. We selected subsets of the geophone array to create four smaller arrays and computed ADDS. The horizontal components of the ADDS were rotated into the direction of the fiber‐optic cable and then compared with the observed DAS strain rates. From three example regional earthquakes of local magnitudes 2.9, 4.1, and 4.3, the ADDS are found to be coherent with DAS for frequencies ≤1 Hz. For frequencies >1‐Hz, this correlation decays quickly. Small differences between linear and areal dynamic strains at 1‐Hz suggest poor signal‐to‐noise or localized strain that is perturbed by shallow heterogeneities compare to the average strain propagating across the geophone array. The implication is that around 1‐Hz, straight fiber DAS is measuring axial strain along the fiber and can provide good approximations to translational particle motions. However, above 1‐Hz, DAS becomes more sensitive to shallow velocity gradients that can be beneficial for geophysical imaging yet becomes a limitation for traditional seismic analysis methods depending on absolute amplitude and phase from translational particle motions.

58 GEOSCIENCES↗
DOE OSTI2025

Bayesian event categorization matrix approach for explosion monitoring

Current efforts to correctly categorize natural events from suspected explosion sources with data that is collected by ground- or space-based sensors presents historical challenges that remain unaddressed by the Event Categorization Matrix (ECM) model. Smaller historical events (lower yield explosions) may have data available from fewer measurement techniques than are available today, and therefore, a historical event record can lack a complete set of discriminants. The covariance structures can also differ between such observations of event (source-type) categories. Both obstacles are problematic for the classic ECM model. Our work addresses this gap and presents a Bayesian update to the previous ECM model, termed the Bayesian Event Categorization Matrix model, which can be trained on partial observations and does not rely on a pooled covariance structure. We further augment the ECM model with Bayesian Decision Theory so that false negative or false positive rates of an event categorization can be reduced in an intuitive manner. To demonstrate improved categorization rates for the Bayesian Event Categorization Matrix model, we compare an array of Bayesian and classic models with multiple performance metrics using Monte Carlo experiments. We use both synthetic and real data. Our Bayesian models show consistent gains in overall accuracy and lower false negative rates relative to the classic ECM model. Here, we propose future avenues to improve Bayesian Event Categorization Matrix models’ decision making and predictive capability.

58 GEOSCIENCES↗
DOE OSTI2025

Subsurface fluvial sediments beneath InSight on Mars from geophysical constraints

Subsurface structure investigation on Mars is crucial for understanding its geological evolution and past hydrological conditions. Elysium Planitia (EP), located near the hypothesized ancient ocean shorelines, could contain clues for past water activity and paleoclimate. Here we present better-constrained subsurface models beneath InSight extending to ~800 m depth, obtained from joint inversion of seismic and seismoacoustic coupling data, and use the well-resolved subsurface structure to explore the lithological profile through rock physics models. The derived subsurface lithology agrees well with local geological context and exhibits a shallow 60-m-thick low-rigidity layer consistent with hydrated sedimentary materials. Despite possible contributions of aeolian and volcanic deposits, we favor the interpretation that the low-rigidity layer originated from fluvial activity in EP during the Hesperian or Hesperian-to-Amazonian epoch, as supported by adjacent paleo-shoreline morphology observations. These results hint at a period of warmer paleoclimate at low latitudes, possibly during high-obliquity phases of Mars’ rotational axis.

58 GEOSCIENCES↗
DOE OSTI2022

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗
DOE OSTI2022

Evaluating the location capabilities of a regional infrasonic network in Utah, US, using both ray tracing-derived and empirical-derived celerity-range and backazimuth models

SUMMARY More realistic models for infrasound signal propagation across a region can be used to improve the precision and accuracy of spatial and temporal source localization estimates. Motivated by incomplete infrasound event bulletins in the Western US, the location capabilities of a regional infrasonic network of stations located between 84–458 km from the Utah Test and Training Range, Utah, USA, is assessed using a series of near-surface explosive events with complementary ground truth (GT) information. Signal arrival times and backazimuth estimates are determined with an automatic F-statistic based signal detector and manually refined by an analyst. This study represents the first application of three distinct celerity-range and backazimuth models to an extensive suite of realistic signal detections for event location purposes. A singular celerity and backazimuth deviation model was previously constructed using ray tracing analysis based on an extensive archive of historical atmospheric specifications and is applied within this study to test location capabilities. Similarly, a set of multivariate, season and location specific models for celerity and backazimuth are compared to an empirical model that depends on the observations across the infrasound network and the GT events, which accounts for atmospheric propagation variations from source to receiver. Discrepancies between observed and predicted signal celerities result in locations with poor accuracy. Application of the empirical model improves both spatial localization precision and accuracy; all but one location estimates retain the true GT location within the 90 per cent confidence bounds. Average mislocation of the events is 15.49 km and average 90 per cent error ellipse areas are 4141 km2. The empirical model additionally reduces origin time residuals; origin time residuals from the other location models are in excess of 160 s while residuals produced with the empirical model are within 30 s of the true origin time. We demonstrate that event location accuracy is driven by a combination of signal propagation model and the azimuthal gap of detecting stations. A direct relationship between mislocation, error ellipse area and increased station azimuthal gaps indicate that for sparse networks, detection backazimuths may drive location biases over traveltime estimates.

58 GEOSCIENCES↗
DOE OSTI2022

Hydromechanical impact of basement rock on injection-induced seismicity in Illinois Basin

Abstract The common explanation of observed injection-induced microseismicity is based on the measured stress state at the injection interval and the assumption that it remains the same in the vicinity. We argue here that representing the stress state in different geologic formations over the injection site with the single Mohr’s circle is insufficient due to local stratigraphic features and contrast in compressibilities of the involved formations. The role of hydromechanical coupling in the microseismic response is also crucial for the proper assessment of the problem. Thoroughly monitored Illinois Basin Decatur Project revealed the majority of CO 2 injection-associated microseismic events being originated in the crystalline basement. Even though basement faults can serve as the conduits for fluid flow—the predicted pressure increase seems to be insufficient to trigger seismicity. To address this issue, accurate laboratory measurements of rock properties from the involved formations are conducted. The pre-injection stress state and its evolution are evaluated with the hydromechanically coupled numerical model. It appears that the presence of an offset in a stiff competent layer affects the stress state in its vicinity. Therefore, both the pre-injection stress state and its evolution during the fluid injection should be addressed during the induced seismicity assessment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluation of Station Performance of the Idaho National Laboratory Seismic Monitoring Network Using Network Detection Thresholds

The Idaho National Laboratory (INL) Seismic Monitoring Network is located in eastern Idaho and monitors a portion of the intermountain seismic belt. It has been in place for 50 yr and has undergone several major changes, the most recent of which has been the transition to the Antelope real‐time acquisition system and the implementation of automatic phase picking algorithms to aid in analysis. This study discusses the efforts to evaluate the performance of the INL seismic monitoring network (and other surrounding stations) using the new real‐time acquisition system. The method outlined by Wilson et al. (2021) is used to develop an empirical relationship between the observability of local earthquakes as a function of magnitude and distance. This relationship is used to produce detection thresholds for Pwaves for all stations of interest. The INL seismic network has two main goals: monitor tectonic‐and volcanic‐related events and measure ground motions for input into seismic hazard analysis. Because of these two overall objectives, several seismic stations have been installed near critical facilities and, therefore, are not as quiet as stations that are used primarily for earthquake detection. This is reflected in their detection thresholds, which are much smaller for stations away from facilities. This study shows that the INL Seismic Monitoring Network is able to detect earthquakes near INL facilities with M L > 1.2, with redundancies built in to ensure this sensitivity even if data became unavailable from some stations. This study also shows “holes” in the monitoring network where the detection of smaller earthquakes is highly dependent on sparsely placed seismic stations. In conclusion, the results of this study will be used to govern plans for expansion of earthquake monitoring in Idaho and the surrounding region and to fine‐tune the detection thresholds for individual stations.

58 - GEOSCIENCES↗
DOE OSTI2021

Regional Moment Tensor Inversion Using Rotational Observations

There are benefits from the addition of rotational motions to translational displacements for moment tensor (MT) inversions. The rotational radiation pattern is orthogonal to the shear radiation pattern, thus incorporating rotations is equivalent to gaining another observation point on the focal sphere. We demonstrated this by simulating curl and displacement wavefields for a regional distance station. Thus, one 6-C station (3-Component translational + 3-Component rotational) gathers the same information on radiation pattern as two 3-C stations at 90° azimuth from one another along the focal plane axis, which is sometimes difficult to obtain when restricted to surface sensors. We added rotational Green's functions to a regional MT inversion scheme (long-period, time-domain, and linear inversion) by computing spatial gradients from f-κ reflectivity synthetics. For rotational data, we used Array Derived Rotations (ADR) from Piñon Flats Observatory Array and Golay array deployed during IRIS Community Wavefield Demonstration Experiment. The hope is to ultimately use compact and field deployable broadband rotational seismometers instead of ADRs. Rotational motions were predicted from well-constrained deviatoric MT solutions of nine earthquakes recorded by the arrays and other seismic networks. Additionally, we formed three station sparse datasets with (3-C and 6-C) and without rotational ground motions (3-C only) for eight earthquakes testing the benefits of including rotational waveforms in MT inversions when they are equally weighted with translational waveforms. Adding rotational motions improved the double-couple components and reduced compensated-linear-vector-dipole and isotropic components in full-MT solutions using sparse datasets with poor coverage.

58 GEOSCIENCES↗
DOE OSTI2023

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗
DOE OSTI2023

3-D Simulations of earthquakes rupture jumps: 1. Homogeneous pre-stress conditions

SUMMARY Observational and modelling studies indicate that earthquake ruptures can jump between fault sections as large as ∼3 and ∼5 km for compressional and extensional offsets, respectively. Here, we compare characteristics of the rupture jump process on parallel but offset fault sections from traditional 3-D dynamic rupture simulations governed by slip weakening friction using the finite element code, FaultMod, to those from quasi-dynamic simulations governed by rate- and state-dependent friction (rate-state friction) using the code RSQSim. These simulations use spatially uniform initial stresses. For a variety of measures the rupture renucleation position on the offset fault, the rate-state friction and slip weakening friction models produce very similar results. The principal difference is the additional occurrence of delayed rupture jumps that arise from the time- and stress-dependent nucleation that is characteristic of rate-state friction. For immediate rupture jumps, models with slip weakening friction span greater offsets than those with rate-state friction. However the jump distances are nearly identical when delayed rupture jumps are included in the comparisons. We propose that delayed rupture jumps are the likely mechanism for adjacent large-earthquake pairs and clusters. Based on the similarity of renucleation positions with both dynamic and quasi-dynamic models, we conclude that the renucleation positions for rupture initiation on the receiver fault (separated by less than ∼3 km from the source fault) are primarily controlled by static stress changes induced by slip on the initiating fault. However, in light of the slightly greater maximum jump distances (>3 km) seen with the dynamic slip weakening friction model, dynamic stress changes from seismic waves play an increasingly important role as offset distances increase.

, RSQSim↗