SEARCH · USGS Science

Results for “The Seismological Record”

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

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21Linked to original sources

USGS2024

Using crustal-scale refraction data of joint inversions of Rayleigh-wave dispersion curves and H/V spectral ratios for Atlantic Coastal Plain velocity structure, eastern U.S.

Shallow shear‐wave velocities ( ⁠⁠ V s ) sometimes are estimated from joint inversions of horizontal‐to‐vertical (H/V) spectral ratios and surface‐wave dispersion curves derived from ambient noise or small active sources. Here, we evaluate carrying out these inversions using Rayleigh‐wave dispersion curves computed from crustal‐scale P ‐wave seismic refraction data. We use data from the 2014–2015 Eastern North American Margin (ENAM) experiment in Virginia and North Carolina, but similar seismic refraction data sets have been acquired over sedimentary basins of interest for seismic hazard studies, including in major urban areas. The ENAM project deployed a pair of ∼215 km long, northwest–southeast linear arrays with ∼300 m receiver spacing to record 11 dynamite shots, and 80 continuously recording seismometers with 5–6 km spacing along the same arrays to record offshore airguns. The arrays crossed the onland portion of the Atlantic Coastal Plain sediments, which are a seaward‐thickening wedge of Cretaceous and younger sediments deposited mostly on crystalline bedrock. We compute Rayleigh‐wave dispersion curves from 3 to 9 km long portions of the receiver arrays on each side of the dynamite shots, and we compute ambient‐noise H/V ratios from the continuously recording seismometers. We use a genetic inversion algorithm in which forward velocity models in each “generation” are evaluated for misfits compared to the observed data, with subsequent generations constructed from the models with the smallest misfits. Velocities to depths of 500 m are defined well, as shown by a narrow range of velocities in the best‐fit models, by the consistency between multiple inversion runs at a site, and by forward modeling of site responses. The resulting velocity cross‐section of the Coastal Plain strata has seaward‐dipping contours in the thinner portions of the Coastal Plain but smaller dips in the deeper portions. We interpret these results as showing that velocity contours in the ACP strata are influenced by a combination of lithology and overburden pressure. Results demonstrate that existing seismic refraction data have the potential for determining detailed shallow shear‐wave velocity profiles.

eastern United States↗

USGS2015

On the reliability of Quake-Catcher Network earthquake detections

Over the past two decades, there have been several initiatives to create volunteer‐based seismic networks. The Personal Seismic Network, proposed around 1990, used a short‐period seismograph to record earthquake waveforms using existing phone lines ( Cranswick and Banfill, 1990 ; Cranswick et al. , 1993 ). NetQuakes ( Luetgert et al. , 2010 ) deploys triaxial Micro‐Electromechanical Systems (MEMS) sensors in private homes, businesses, and public buildings where there is an Internet connection. Other seismic networks using a dense array of low‐cost MEMS sensors are the Community Seismic Network ( Clayton et al. , 2012 ; Kohler et al. , 2013 ) and the Home Seismometer Network ( Horiuchi et al. , 2009 ). One main advantage of combining low‐cost MEMS sensors and existing Internet connection in public and private buildings over the traditional networks is the reduction in installation and maintenance costs ( Koide et al. , 2006 ). In doing so, it is possible to create a dense seismic network for a fraction of the cost of traditional seismic networks ( D’Alessandro and D’Anna, 2013 ; D’Alessandro, 2014 ; D’Alessandro et al. , 2014 ).

Seismological Research Letters↗

USGS2025

Seismic intensities, ground motions and basin amplification from the 2022 Mw 6.4 Ferndale, California intraslab earthquake

On December 20, 2022, a Mw 6.4 earthquake occurred at a depth of 18 km within the subducting Gorda plate in the Mendocino Triple Junction (MTJ), one of the most seismically active regions of the contiguous United States, causing widespread damage to local communities. Here we document the seismic intensities, ground motions, and basin amplification effects recorded by this earthquake across Humboldt County as part of an ongoing scientific effort to understand subduction zone earthquake hazards. Modified Mercalli Intensity (MMI) values from our post-earthquake field survey report shaking intensities as high as VIII (Severe). Strong ground motion data from 54 seismic stations were processed to calculate amplitude and frequency content parameters across Northern California. The maximum calculated geometric mean PGA and PGV are 1053 cm/s 2 and 52 cm/s (respectively), both recorded within the 3 km deep Eel River sedimentary basin. Comparisons with four published Ground Motion Prediction Equations indicate that PGA and PGV measurements align with expected attenuation-distance patterns for intraslab earthquakes of this nature. Within the Eel River Basin, ground motions for frequencies above 1 Hz are amplified, with respect to reference stations, by factors greater than 2. Our findings suggest that peak ground motions were mainly caused by sedimentary basin site-effects within the Eel River basin, although rupture directivity may have also increased ground motion amplitudes locally. While the Ferndale area is primarily impacted by shallow intraslab earthquakes, our results also raise questions about site-response and basin amplification hazards from a potential megathrust earthquake. The high seismicity rates of the southernmost Cascades call for stronger regional preparedness and improved strategies to mitigate the effects of such a large-scale disaster.

California↗

USGS2001

Observations of earthquake source parameters at 2 km depth in the Long Valley Caldera, eastern California

To investigate seismic source parameter scaling and seismic efficiency in the Long Valley caldera, California, we measured source parameters for 41 earthquakes (M 0.5 to M 5) recorded at 2 km depth in the Long Valley Exploratory Well. Borehole recordings provide a wide frequency bandwidth, typically 1 to 200–300 Hz, and greatly reduce seismic noise and path effects compared to surface recordings. We calculated source parameters in both the time and frequency domains for P and S waves. At frequencies above the corner frequency, spectra decay faster than ω 3 , indicating that attenuation plays an important role in shaping the spectra (path averaged Q p = 100–400, Q s = 200–800). Source parameters are corrected for attenuation and radiation pattern. Both static stress drops and apparent stresses range from approximately 0.01 to 30 MPa. Although static stress drops do not vary with seismic moment for these data, our analyses are consistent with apparent stress increasing with increasing moment. To estimate tectonic driving stress and seismic efficiencies in the region, we combined source parameter measurements with knowledge of the stress field and a Coulomb failure criterion to infer a driving stress of 40–70 MPa. Subsequent seismic efficiencies are consistent with McGarr's (1999) hypothesis of a maximum seismic efficiency of 6%.

California↗

USGS1997

Water-level changes in response to the 20 December 1994 earthquake near Parkfield, California

We analyze co-seismic changes of water level in nine wells near Parkfield, California, produced by an M D 4.7 earthquake on 20 December 1994 in order to test the hypothesis that co-seismic water-level changes are proportional to co-seismic volumetric strain. For each well, a quantitative relationship between water level and volumetric strain can be inferred from water-level fluctuations due to earth tides and barometric pressure. The observed co-seismic water-level changes, which ranged from −16 to +34 cm, can therefore be compared with volumetric strain recorded by borehole strainmeters or calculated using a dislocation model of the earthquake. We were able to find a dislocation model of the earthquake rupture that predicts volumetric expansion at five of the six wells where water level fell co-seismically, as well as volumetric contraction at one of the two sites where water level rose. Strain predicted by the dislocation model is in good quantitative agreement with the strain inferred from water-level changes observed at four of the well sites, as well as strain recorded by three borehole strainmeters. Water-level changes at two more well sites correspond to strain somewhat greater than predicted by the model but agree in sign with model-calculated strains. At three of the well sites, however, water-level changes took place that cannot be explained as responses to co-seismic volumetric strain for any plausible dislocation model of the earthquake rupture. At two of these sites, one in and one near the San Andreas fault, large water-level drops are probably influenced by co-seismic fault creep. The third site has a history of large water-level rises in response to earthquakes at distances up to several hundred kilometers. This data set shows that co-seismic water-level changes in many wells are proportional to volumetric strain but that other wells exist in which different mechanisms dominate co-seismic response.

California↗

USGS2001

Distribution of fine-scale mantle heterogeneity from observations of Pdiff coda

We present stacked record sections of Global Seismic Network data that image the average amplitude and polarization of the high-frequency Pdiff coda and investigate their implications on the depth extent of fine-scale (~10 km) mantle heterogeneity. The extended 1-Hz coda lasts for at least 150 sec and is observed to a distance of 130??. The coda's polarization angle is about the same as the main Pdiff arrival (4.4 sec/deg) and is nearly constant with time. Previous studies show that multiple scattering from heterogeneity restricted to the lowermost mantle generates an extended Pdiff coda with a constant polarization. Here we present an alternative model that satisfies our Pdiff observations. The model consists of single scattering from weak (~1%) fine-scale (~2 km) structures distributed throughout the mantle. Although this model is nonunique, it demonstrates that Pdiff coda observations do not preclude the existence of scattering contributions from the entire mantle.

Bulletin of the Seismological Society of America↗

USGS2020

Damping values derived from surface-source, downhole-receiver measurements at 22 sites in the San Francisco Bay Area of central California and the San Fernando Valley of southern California

A method discussed in Gibbs, Boore, et al. (1994) was applied to surface‐source, downhole‐receiver recordings at 22 boreholes, in the San Francisco Bay area in central California and the San Fernando Valley of southern California, to determine the average damping ratio of shear waves over depth intervals ranging from about 10 m to as much as 245 m (at one site), with most maximum depths being between 35 and 90 m. The average damping values range from somewhat less than 1% to almost 8%, with little dependence on grain size for sites in sediments. Surprisingly, the average damping values for sites with average velocities greater than about 450    m / s "> 450 m / s 450 m/s ⁠ , including, but not limited to rock sites, are generally larger than for sites with lower average velocities. The combined effect of the higher damping and shorter travel times through the rock columns, however, leads to an effective attenuation that is generally comparable or smaller than for soil sites.

California↗

USGS2001

Observations of basin ground motions from a dense seismic array in San Jose, California

We installed a dense array of 41 digital seismographs in San Jose, California, to evaluate in detail the effects of a deep sedimentary basin and shallow sedimentary deposits on earthquake ground motions. This urban array is located near the eastern edge of the Santa Clara Valley and spans the Evergreen sedimentary basin identified by gravity data. Average station spacing is 1 km, with three stations initially spaced 110 m apart. Despite the high-noise urban environment, the stations of the array successfully triggered on and recorded small local earthquakes (M 2.5-2.8 at 10-25 km distance) and larger regional events such as the M 5.0 Bolinas earthquake (90 km distance), M 4.6-5.6 earthquakes near Mammoth Lakes (270 km distance), M 4.9-5.6 events in western Nevada (420 km distance) and the M 7.1 Hector Mine earthquake (590 km distance). Maps of spectral ratios across the array show that the highest amplitudes in all frequency bands studied (0.125-8 Hz) are generally observed at stations farther from the eastern edge of the Santa Clara Valley. Larger spectral amplitudes are often observed above the western edge of the Evergreen Basin. Snapshots of the recorded wavefield crossing the array for regional events to the east reveal that large, low-frequency (0.125-0.5 Hz) arrivals after the S-wave travel from south to north across the array. A moving-window, cross-correlation analysis finds that these later arrivals are surface waves traveling from the south. The timing and propagation direction of these arrivals indicates that they were likely produced by scattering of incident S waves at the border of the Santa Clara Valley to the south of the array. It is remarkable that the largest low-frequency phases at many of the valley sites for regional events to the east are basin surface waves coming from a direction about 70 degrees different from that of the epicenters. Basin surface waves emanating from the eastern edge of the valley are also identified by the cross-correlation analysis.

Bulletin of the Seismological Society of America↗

USGS2019

Using component ratios to detect metadata and instrument problems of seismic stations: Examples from 18 years of GEOSCOPE data

Replacement or deterioration of seismic instruments and the evolution of the installation conditions and sites can alter the seismic signal in very subtle ways, so it is notoriously difficult to monitor the signal quality of permanent seismic stations. We present a simple tool, energy ratios between each pair of the three recorded components, aimed at characterizing and monitoring signal quality, as a complement to existing methods. To calculate stable daily energy ratios over a large frequency range (0.01 Hz – 5 Hz), we use the daily median energy ratio over all 5-minute windows within the day. The method is applied to all GEOSCOPE stations, for continuous BH channel data collected since 2001. We show applications to identify past gain problems (stations ROCAM and CRZF), to provide feedback after field interventions at remote sites (Antarctic station DRV), and to shed light on complex instrument problems (stations ECH and KIP). Our results show that component energy ratios have excellent time resolution and that they are visually simple for identification of problems. They can be used both for ongoing continuous monitoring of the signal quality, or as a tool to identify past problems. The associated python codes are available on GitHub.

Seismological Research Letters↗

USGS2009

Earthquake stress drops and inferred fault strength on the Hayward Fault, east San Francisco Bay, California

We study variations in earthquake stress drop with respect to depth, faulting regime, creeping versus locked fault behavior, and wall-rock geology. We use the P-wave displacement spectra from borehole seismic recordings of M 1.0-4.2 earthquakes in the east San Francisco Bay to estimate stress drop using a stack-and-invert empirical Green's function method. The median stress drop is 8.7 MPa, and most stress drops are in the range between 0.4 and 130 MPa. An apparent correlation between stress drop and magnitude is entirely an artifact of the limited frequency band of 4-55 Hz. There is a trend of increasing stress drop with depth, with a median stress drop of ~5 MPa for 1-7 km depth, ~10 MPa for 7-13 km depth, and ~50 MPa deeper than 13 km. We use S=P amplitude ratios measured from the borehole records to better constrain the first-motion focal mechanisms. High stress drops are observed for a deep cluster of thrust-faulting earthquakes. The correlation of stress drops with depth and faulting regime implies that stress drop is related to the applied shear stress. We compare the spatial distribution of stress drops on the Hayward fault to a model of creeping versus locked behavior of the fault and find that high stress drops are concentrated around the major locked patch near Oakland. This also suggests a connection between stress drop and applied shear stress, as the locked patch may experience higher applied shear stress as a result of the difference in cumulative slip or the presence of higher-strength material. The stress drops do not directly correlate with the strength of the proposed wall-rock geology at depth, suggesting that the relationship between fault strength and the strength of the wall rock is complex.

Bulletin of the Seismological Society of America↗

USGS2015

High‐resolution trench photomosaics from image‐based modeling: Workflow and error analysis

Photomosaics are commonly used to construct maps of paleoseismic trench exposures, but the conventional process of manually using image‐editing software is time consuming and produces undesirable artifacts and distortions. Herein, we document and evaluate the application of image‐based modeling (IBM) for creating photomosaics and 3D models of paleoseismic trench exposures, illustrated with a case‐study trench across the Wasatch fault in Alpine, Utah. Our results include a structure‐from‐motion workflow for the semiautomated creation of seamless, high‐resolution photomosaics designed for rapid implementation in a field setting. Compared with conventional manual methods, the IBM photomosaic method provides a more accurate, continuous, and detailed record of paleoseismic trench exposures in approximately half the processing time and 15%–20% of the user input time. Our error analysis quantifies the effect of the number and spatial distribution of control points on model accuracy. For this case study, an ∼87 m 2 exposure of a benched trench photographed at viewing distances of 1.5–7 m yields a model with <2 cm root mean square error (rmse) with as few as six control points. Rmse decreases as more control points are implemented, but the gains in accuracy are minimal beyond 12 control points. Spreading control points throughout the target area helps to minimize error. We propose that 3D digital models and corresponding photomosaics should be standard practice in paleoseismic exposure archiving. The error analysis serves as a guide for future investigations that seek balance between speed and accuracy during photomosaic and 3D model construction.

Bulletin of the Seismological Society of America↗

USGS2010

Temporal variations in Global Seismic Stations ambient noise power levels

Recent concerns about time-dependent response changes in broadband seismometers have motivated the need for methods to monitor sensor health at Global Seismographic Network (GSN) stations. We present two new methods for monitoring temporal changes in data quality and instrument response transfer functions that are independent of Earth seismic velocity and attenuation models by comparing power levels against different baseline values. Our methods can resolve changes in both horizontal and vertical components in a broad range of periods (∼0.05 to 1,000 seconds) in near real time. In this report, we compare our methods with existing techniques and demonstrate how to resolve instrument response changes in long-period data (>100 seconds) as well as in the microseism bands (5 to 20 seconds). High quality broadband data recorded by the GSN are fundamental to characterizing a wide range of Earth science issues including: the size and rupture of large earthquakes ( e.g. , Tsai et al. 2005 ); imaging the interior of the Earth ( e.g. , Van der Hilst et al. 1997 ); tracking global climate variation ( Aster et al. 2008 ); and monitoring calving glaciers (Ekström et al. 2003 , 2006a ). Recent studies based on theoretical Earth models ( Ekström et al. 2006b ; Davis and Berger 2007 ) suggest that broadband seismometer gain levels can vary with time. This has also been confirmed, for the STS-1 sensor, experimentally ( Yuki and Ishihara 2002 ). It therefore has become necessary to systematically check for temporal changes in amplitude at GSN stations. Many of these changes are frequency-dependent in nature and not a priori predictable ( Ekström et al. 2006b ). Robust methods that can be applied to a large number of stations in a broad range of frequency bands are necessary.

Seismological Research Letters↗

USGS2020

Near-field ground motions from the July, 2019 Ridgecrest, California, earthquake sequence

The 2019 Ridgecrest, California, earthquake sequence, including an Mw 6.4 event on 4 July and an Mw 7.1 approximately 34 hr later, was recorded by 15 instruments within 55 km nearest‐fault distance. To characterize and explore near‐field ground motions from the Mw 6.4 foreshock and Mw 7.1 mainshock, we augment these records with available macroseismic information, including conventional intensities and displaced rocks. We conclude that near‐field shaking intensities were generally below modified Mercalli intensity 9, with concentrations of locally high values toward the northern and southern termini of the mainshock rupture. We further show that, relative to near‐field ground motions at hard‐rock sites, instrumental ground motions at alluvial near‐field sites for both the Mw 6.4 foreshock and Mw 7.1 mainshock were depleted in energy at frequencies higher than 2–3 Hz, as expected from ground‐motion models. Both the macroseismic and instrumental observations suggest that sediments in the Indian Wells Valley experienced a pervasively nonlinear response, which helps explain why shaking intensities and damage in the closest population center, Ridgecrest, were relatively modest given its proximity to the earthquakes.

California↗

USGS2010

Seismology of the continental crust and upper mantle

More seismological studies of the continental crust and sub‐crustal lithosphere of the United States have been completed in the past four years than at any other similar period, and a continued growth in activity is likely to continue for years to come. Several trends account for this phenomenon. First, the interest in seismic reflection studies generated initially by COCORP results in this country [ Brown and others , 1986], and later by the British BIRPS results [ Matthews and Gheadie , 1986], has led to the development of several other seismic reflection programs. Among the most active of these research programs are those of the University of Wyoming [ Smithson and others , 1986], Virginia Polytechnic Institute, CALCRUST [ Henyey , 1986], and the U.S. Geological Survey (USGS) [ Hamilton , 1986]. In Canada, the Lithoprobe program has achieved remarkable results in a variety of geographic locations [ Green et al. , 1986]. A second trend is the resurgence of interest in seismic‐refraction/wide‐angle reflection profiling. The year 1978 marked the beginning of this increased activity when several large projects were conducted in the western U.S., such as the Yellowstone‐Snake River Plain experiment organized by the University of Utah. Since 1979, a large amount of refraction/wide‐angle reflection data has been collected by the USGS and by university groups utilizing large numbers of state‐of‐the‐art industry seismographs in cooperative experiments. A third trend is the increased sophistication of other seismic methods such as teleseismic delay‐time methods, tomography, and receiver‐transfer functions. In these studies, permanent or temporary recording arrays have been used to determine local and regional crustal and upper‐mantle structure with impressive resolution.

Reviews of Geophysics↗

USGS2009

The puzzle of the 1996 Bárdarbunga, Iceland, earthquake: no volumetric component in the source mechanism

A volcanic earthquake with M w 5.6 occurred beneath the Bárdarbunga caldera in Iceland on 29 September 1996. This earthquake is one of a decade-long sequence of events at Bárdarbunga with non-double-couple mechanisms in the Global Centroid Moment Tensor catalog. Fortunately, it was recorded well by the regional-scale Iceland Hotspot Project seismic experiment. We investigated the event with a complete moment tensor inversion method using regional long-period seismic waveforms and a composite structural model. The moment tensor inversion using data from stations of the Iceland Hotspot Project yields a non-double-couple solution with a 67% vertically oriented compensated linear vector dipole component, a 32% double-couple component, and a statistically insignificant (2%) volumetric (isotropic) contraction. This indicates the absence of a net volumetric component, which is puzzling in the case of a large volcanic earthquake that apparently is not explained by shear slip on a planar fault. A possible volcanic mechanism that can produce an earthquake without a volumetric component involves two offset sources with similar but opposite volume changes. We show that although such a model cannot be ruled out, the circumstances under which it could happen are rare.

Bardarbunga↗

USGS2017

Detection and characterization of pulses in broadband seismometers

Pulsing - caused either by mechanical or electrical glitches, or by microtilt local to a seismometer - can significantly compromise the long‐period noise performance of broadband seismometers. High‐fidelity long‐period recordings are needed for accurate calculation of quantities such as moment tensors, fault‐slip models, and normal‐mode measurements. Such pulses have long been recognized in accelerometers, and methods have been developed to correct these acceleration steps, but considerable work remains to be done in order to detect and correct similar pulses in broadband seismic data. We present a method for detecting and characterizing the pulses using data from a range of broadband sensor types installed in the Global Seismographic Network. The technique relies on accurate instrument response removal and employs a moving‐window approach looking for acceleration baseline shifts. We find that pulses are present at varying levels in all sensor types studied. Pulse‐detection results compared with average daily station noise values are consistent with predicted noise levels of acceleration steps. This indicates that we can calculate maximum pulse amplitude allowed per time window that would be acceptable without compromising long‐period data analysis.

Bulletin of the Seismological Society of America↗

USGS1988

Seismicity in South Carolina

The largest historical earthquake in South Carolina, and in the southeastern US, occurred in the Coastal Plain province, probably northwest of Charleston, in 1886. Locations for aftershocks associated with this earthquake, estimated using intensities based on newspaper accounts, defined a northwest trending zone about 250 km long that was at least 100 km wide in the Coastal Plain but widened to a northeast trending zone in the Piedmont. The subsequent historical and instrumentally recorded seismicity in South Carolina images the 1886 aftershock zone. Except for a few scattered earthquakes and a swarm of shallow (≤ 4 km deep), small ( M L ≤ 2.5), primarily reverse faulting earthquakes that occurred along the flanks of a granite pluton about 60 km northwest of Columbia, the seismicity in the Piedmont province has been associated with water level changes in reservoirs. Reservoir induced seismicity (RIS) is shallow (≤ 6 km deep), primarily strike-slip or thrust faulting corresponding to an inferred maximum horizontal compressive stress oriented approximately N 60° E. Instrumentally recorded seismicity in the Coastal Plain province occurs in 3 seismic zones or clusters: Middleton Place-Summerville (MPSSZ), Adams Run (ARC), and Bowman (BSZ). Approximately 68% of the Coastal Plain earthquakes occur in the MPSSZ, a north trending zone about 22 km long and 12 km wide, lying about 20 km northwest of Charleston. The hypocenters of MPSSZ earthquakes range in depth from near the surface to almost 12 km. Thrust, strike-slip, and some normal faulting are indicated by the fault plane solutions for Coastal Plain earthquakes. The maximum horizontal compressive stress, inferred from the P-axes of the fault plane solutions, is oriented NE-SW in the shallow crust (< 9 km deep) but appears to be diffusely E-W between 9 to 12 km deep. Although there is localized variability, the current seismicity and associated faulting in South Carolina probably represent a regional response to the NE-SW maximum horizontal compressive stress prevalent throughout eastern North America.

South Carolina↗

USGS1996

Rupture directivity and slip distribution of the M 4.3 foreshock to the 1992 Joshua Tree earthquake, Southern California

Details of the M 4.3 foreshock to the Joshua Tree earthquake were studied using P waves recorded on the Southern California Seismic Network and the Anza network. Deconvolution, using an M 2.4 event as an empirical Green's function, corrected for complicated path and site effects in the seismograms and produced simple far-field displacement pulses that were inverted for a slip distribution. Both possible fault planes, north-south and east-west, for the focal mechanism were tested by a least-squares inversion procedure with a range of rupture velocities. The results showed that the foreshock ruptured the north-south plane, similar to the mainshock. The foreshock initiated a few hundred meters south of the mainshock and ruptured to the north, toward the mainshock hypocenter. The mainshock (M 6.1) initiated near the northern edge of the foreshock rupture 2 hr later. The foreshock had a high stress drop (320 to 800 bars) and broke a small portion of the fault adjacent to the mainshock but was not able to immediately initiate the mainshock rupture.

Bulletin of the Seismological Society of America↗