SEARCH · Search USGS

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 523 records · Page 29Linked to original sources

USGS2022

Testing the ShakeAlert earthquake early warning system using synthesized earthquake sequences

We test the behavior of the United States (US) West Coast ShakeAlert earthquake early warning (EEW) system during temporally close earthquake pairs to understand current performance and limitations. We consider performance metrics based on source parameter and ground‐motion forecast accuracy, as well as on alerting timeliness. We generate ground‐motion times series for synthesized earthquake sequences from real data by combining the signals from pairs of well‐recorded earthquakes (⁠4.4≤M≤7.1⁠) using time shifts ranging from −60 to +180 s. We examine fore‐ and aftershock sequences, near‐simultaneous events in different source regions, and simulated out‐of‐network and offshore earthquakes. We find that the operational ShakeAlert algorithms Earthquake Point‐source Integrated Code (EPIC) and Finite‐Fault Rupture Detector (FinDer) and the Propagation of Local Undamped Motion (PLUM) method perform largely as expected: EPIC provides the best source location estimates and is often fastest but can underestimate magnitudes or, in extreme cases, miss large earthquakes; FinDer provides real‐time line‐source models and unsaturated magnitude estimates for large earthquakes but currently cannot process concurrent events and may mislocate offshore earthquakes; PLUM identifies pockets of strong ground motion, but can overestimate alert areas. Implications for system performance are: (1) spatially and temporally close events are difficult to identify separately; (2) challenging scenarios with foreshocks that are close in space and time can lead to missed alerts for large earthquakes; and (3) in these situations the algorithms can often estimate ground motion better than source parameters. To improve EEW, our work suggests revisiting the current algorithm weighting in ShakeAlert, to continue developments that focus on using ground‐motion data to aggregate alerts from multiple algorithms, and to investigate methods to optimally leverage algorithm ground‐motion estimates. For testing and certification of EEW performance in ShakeAlert and other EEW systems where applicable, we also suggest that 25 of our 73 scenarios become part of the baseline data set.

California↗

USGS1996

On the coherence of ground motion in the San Fernando valley

We present an analysis of the coherence of seismic ground motion recorded on alluvial sediments in the San Fernando Valley, California. Using aftershocks of the 17 January 1994 M w 6.7 earthquake recorded at a quasi-dense array of portable stations, we analyze the coherence of three well-recorded magnitude 3.7 to 4.0 events over the frequency range 0.5 to 15 Hz and a distance range of 0.5 to 5.3 km. All stations are located at sites with broadly similar near-site geology, characterized by medium to fine-grain Quaternary alluvial sediments. On average, relatively high values of coherence are observed for distances up to 3 to 4 km and frequencies up to 2 to 3 Hz; coherence drops sharply at frequencies near and above 3 Hz. Although average coherence functions are described reasonably well by a log-linear relationship with frequency, the curves at all distances exhibit a flattening at low frequencies that is not consistent with previous observations of coherence at hardrock sites. The distance decay of coherence is also markedly less strong, with high coherence values observed over station separations corresponding to multiple wavelengths. This may reflect fundamental differences in shallow-wave propagation in the two environments, with high-frequency scattering relatively more dominant in regions of hard-rock near-surface geology. Within a sedimentary basin or valley, the site response itself generally reflects a resonance phenomenon that may tend to give rise to more uniform ground motions. However, previous studies have demonstrated the existence of pathological focusing and amplification effects within complex sedimentary basin environments such as the greater Los Angeles region; our results undoubtedly do not quantify the full range of ground-motion variability at all sites, but rather represent the level of that variability that can be expected, and quantified, for typical source/receiver paths.

California↗

USGS2001

Comparisons of ground motions from the 1999 Chi-Chi, earthquake with empirical predictions largely based on data from California

This article has the modest goal of comparing the ground motions recorded during the 1999 Chi-Chi, Taiwan, mainshock with predictions from four empirical-based equations commonly used for western North America; these empirical predictions are largely based on data from California. Comparisons are made for peak acceleration and 5%-damped response spectra at periods between 0.1 and 4 sec. The general finding is that the Chi-Chi ground motions are smaller than those predicted from the empirically based equations for periods less than about 1 sec by factors averaging about 0.4 but as small as 0.26 (depending on period, on which equation is used, and on whether the sites are assumed to be rock or soil). There is a trend for the observed motions to approach or even exceed the predicted motions for longer periods. Motions at similar distances (30-60 km) to the east and to the west of the fault differ dramatically at periods between about 2 and 20 sec: Long-duration wave trains are present on the motions to the west, and when normalized to similar amplitudes at short periods, the response spectra of the motions at the western stations are as much as five times larger than those of motions from eastern stations. The explanation for the difference is probably related to site and propagation effects; the western stations are on the Coastal Plain, whereas the eastern stations are at the foot of young and steep mountains, either in the relatively narrow Longitudinal Valley or along the eastern coast-the sediments underlying the eastern stations are probably shallower and have higher velocity than those under the western stations.

Bulletin of the Seismological Society of America↗

USGS1994

Magnetic field observations in the near-field the 28 June 1992 Mw 7.3 Landers, California, earthquake

Recent reports suggest that large magnetic field changes occur prior to, and during, large earthquakes. Two continuously operating proton magnetometers, LSBM and OCHM , at distances of 17.3 and 24.2 km, respectively, from the epicenter of the 28 June 1992 M w 7.3 Landers earthquake, recorded data through the earthquake and its aftershocks. These two stations are part of a differentially connected array of proton magnetometers that has been operated along the San Andreas fault since 1976. The instruments have a sensitivity of 0.25 nT or better and transmit data every 10 min through the GOES satellite to the USGS headquarters in Menlo Park, California. Seismomagnetic offsets of −1.2 ± 0.6 and −0.7 ± 0.7 nT were observed at these sites. In comparison, offsets of −0.3 ± 0.2 and −1.3 ± 0.2 nT were observed during the 8 July 1986 M L 5.9 North Palm Springs earthquake, which occurred directly beneath the OCHM magnetometer site. The observations are generally consistent with seismomagnetic models of the earthquake, in which fault geometry and slip have the same from as that determined by either inversion of the seismic data or inversion of geodetically determined ground displacements produced by the earthquake. In these models, right-lateral rupture occurs on connected fault segments in a homogeneous medium with average magnetization of 2 A/m. The fault-slip distribution has roughly the same form as the observed surface rupture, and the total moment release is 1.1 × 10 20 Nm. There is no indication of diffusion-like character to the magnetic field offsets that might indicate these effects result from fluid flow phenomena. It thus seems unlikely that these earthquake-generated offsets and those produced by the North Palm Springs earthquake were generated by electrokinetic effects. Also, there are no indications of enhanced low-frequency magnetic noise before the earthquake at frequencies below 0.001 Hz.

California↗

USGS2016

Comparisons between vs30 and spectral response for 30 sites in Newcastle, Australia from collocated seismic cone penetrometer, active- and passive-source vs data

Although the time‐averaged shear‐wave velocity down to 30 m depth ( V S 30 ) can be a proxy for estimating earthquake ground‐motion amplification, significant controversy exists about its limitations when used as a single parameter for the prediction of amplification. To examine this question in absence of relevant strong‐motion records, we use a range of different methods to measure the shear‐wave velocity profiles and the resulting theoretical site amplification factors (AFs) for 30 sites in the Newcastle area, Australia, in a series of blind comparison studies. The multimethod approach used here combines past seismic cone penetrometer and spectral analysis of surface‐wave data, with newly acquired horizontal‐to‐vertical spectral ratio, passive‐source surface‐wave spatial autocorrelation (SPAC), refraction microtremor (ReMi), and multichannel analysis of surface‐wave data. The various measurement techniques predicted a range of different AFs. The SPAC and ReMi techniques have the smallest overall deviation from the median AF for the majority of sites. We show that V S 30 can be related to spectral response above a period T of 0.5 s but not necessarily with the maximum amplification according to the modeling done based on the measured shear‐wave velocity profiles. Both V S 30 and AF values are influenced by the velocity ratio between bedrock and overlying sediments and the presence of surficial thin low‐velocity layers (<2 m thick and <150 m/s), but the velocity ratio is what mostly affects the AF. At 0.2< T <0.4 s, the AFs are largely controlled by the surficial geology of a particular site. AF maxima are the highest in the hard classes, which is the inverse of the findings used in the Australian Building Code. Only for T >0.5 s do the amplification curves consistently show higher values for soft site classes and lower for hard classes.

Bulletin of the Seismological Society of America↗

USGS1994

Regional propagation characteristics and source parameters of earthquakes in northeastern North America

The vertical components of the S wave trains recorded on the Eastern Canadian Telemetered Network ( ECTN ) from 1980 through 1990 have been spectrally analyzed for source, site, and propagation characteristics. The data set comprises some 1033 recordings of 97 earthquakes whose magnitudes range from M ≈ 3 to 6. The epicentral distances range from 15 to 1000 km, with most of the data set recorded at distances from 200 to 800 km. The recorded S wave trains contain the phases S , S m S , S n , and L g and are sampled using windows that increase with distance; the acceleration spectra were analyzed from 1.0 to 10 Hz. To separate the source, site, and propagation characteristics, an inversion for the earthquake corner frequencies, low-frequency levels, and average attenuation parameters is alternated with a regression of residuals onto the set of stations and a grid of 14 distances ranging from 25 to 1000 km. The iteration between these two parts of the inversion converges in about 60 steps. The average attenuation parameters obtained from the inversion were Q = 1997 ± 10 and γ = 0.998 ± 0.003. The most pronounced variation from this average attenuation is a marked deamplification of more than a factor of 2 at 63 km and 2 Hz, which shallows with increasing frequency and increasing distance out to 200 km. The site-response spectra obtained for the ECTN stations are generally flat. The source spectral shape assumed in this inversion provides an adequate spectral model for the smaller events ( M o < 3 × 10 21 dyne-cm) in the data set, whose Brune stress drops range from 5 to 150 bars. For the five events in the data set with M o ≧ 10 23 dyne-cm, however, the source spectra obtained by regressing the residuals suggest that an ω 2 spectrum is an inadequate model for the spectral shape. In particular, the corner frequencies for most of these large events appear to be split, so that the spectra exhibit an intermediate behavior (where | ü ( ω )| is roughly proportional to ω ).

Bulletin of the Seismological Society of America↗

USGS2006

Some Key Features of the Strong-Motion Data from the M 6.0 Parkfield, California, Earthquake of 28 September 2004

The 2004 Parkfield, California, earthquake was recorded by an extensive set of strong-motion instruments well positioned to record details of the motion in the near-fault region, where there has previously been very little recorded data. The strong-motion measurements obtained are highly varied, with significant variations occurring over only a few kilometers. The peak accelerations in the near fault region range from 0.13g to over 1.8g (one of the highest acceleration recorded to date, exceeding the capacity of the recording instrument The largest accelerations occurred near the northwest end of the inferred rupture zone. These motions are consistent with directivity for a fault rupturing from the hypocenter near Gold Hill toward the northwest. However, accelerations up to 0.8g were also observed in the opposite direction, at the south end of the Cholame Valley near Highway 41, consistent with bilateral rupture, with rupture southeast of the hypocenter. Several stations near and over the rupturing fault recorded relatively weak motions, consistent with seemingly paradoxical observations of low shaking damage near strike-slip faults. This event had more ground-motion observations within 10 km of the fault than many other earthquakes combined. At moderate distances peak horizontal ground acceleration (PGA) values dropped off more rapidly with distance than standard relationships. At close-in distance the wide variation of PGA suggests a distance-dependent sigma may be important to consider. The near-fault ground-motion variation is greater than that assumed in ShakeMap interpolations, based on the existing set of observed data. Higher density of stations near faults may be the only means in the near future to reduce uncertainty in the interpolations. Outside of the near-fault zone the variance is closer to that assumed. This set of data provides the first case where near-fault radiation has been observed at an adequate number of stations around the fault to allow detailed study of the fault-normal and fault-parallel motion and the near-field S-wave radiation. The fault-normal motions are significant, but they are not large at the central part of the fault, away from the ends. The fault-normal and fault-parallel motions drop off quite rapidly with distance from the fault. Analysis of directivity indicates increased values of peak velocity in the rupture direction. No such dependence is observed in the peak acceleration, except for stations close to the strike of the fault near and beyond the ends of the faulting.

Bulletin of the Seismological Society of America↗

USGS2021

Basin and site effects in the U.S. Pacific Northwest estimated from small‐magnitude earthquakes

Earthquake hazards in the U.S. Pacific Northwest (PNW) are increased by the presence of deep sedimentary basins that amplify and prolong ground shaking. To better understand basin and site effects on ground motions, we compile a database of recordings from crustal and intraslab earthquakes. We process 8028 records with magnitudes from 3.5 to 6.8 and hypocentral depths up to 62 km to compute Fourier amplitude spectra of ground acceleration for frequencies of 0–20 Hz. We compute residuals relative to the Bayless and Abrahamson (2019; hereafter, BA18) ground‐motion model and perform a series of linear, crossed, mixed‐effects regressions. In addition to estimating the bias, event, and site terms, we incorporate groupings for broad regionalized site response in three different regions (Seattle basin, Puget Lowland, non‐Puget Lowland), for effects from seismotectonic regime (crustal and intraslab sources), and for interactions between the regions and seismotectonic regimes. We find that the scaling of site response with respect to V S 30 (time‐averaged shear‐wave velocity from the surface to a depth of 30 m) and to basin depth indicators Z 1.0 and Z 2.5 (depths to the 1.0 and 2.5 km/s shear‐wave velocity horizons) is generally consistent with BA18; however, the region terms display strong spatial amplification patterns. For frequencies less than 5 Hz, the Seattle basin amplifies ground motions up to a factor of four, relative to the non‐Puget Lowland, with a maximum amplification around near 0.5 Hz. Sites in the Puget Lowland amplify low frequencies up to a factor of 2.5. At higher frequencies ( ⁠ f > 5 Hz ⁠ ), the Puget Lowland and Seattle basin show regional deamplification of ground motions, with the smallest average amplification factor of 0.65 occurring at 10.0 Hz. Although we observe slight differences in the seismotectonic regime terms, we find that the region terms are significantly more important for modeling earthquake hazard in the PNW.

Oregon, Washington↗

USGS1996

Use of liquefaction-induced features for paleoseismic analysis

Liquefaction features can be used in many field settings to estimate the recurrence interval and magnitude of strong earthquakes through much of the Holocene. These features include dikes, craters, vented sand, sills, and laterally spreading landslides. The relatively high seismic shaking level required for their formation makes them particularly valuable as records of strong paleo-earthquakes. This state-of-the-art summary for using liquefaction-induced features for paleoseismic interpretation and analysis takes into account both geological and geotechnical engineering perspectives. The driving mechanism for formation of the features is primarily the increased pore-water pressure associated with liquefaction of sand-rich sediment. The role of this mechanism is often supplemented greatly by the direct action of seismic shaking at the ground surface, which strains and breaks the clay-rich cap that lies immediately above the sediment that liquefied. Discussed in the text are the processes involved in formation of the features, as well as their morphology and characteristics in field settings. Whether liquefaction occurs is controlled mainly by sediment grain size, sediment packing, depth to the water table, and strength and duration of seismic shaking. Formation of recognizable features in the field generally requires a low-permeability cap above the sediment that liquefied. Field manifestations are controlled largely by the severity of liquefaction and the thickness and properties of the low-permeability cap. Criteria are presented for determining whether observed sediment deformation in the field originated by seismically induced liquefaction. These criteria have been developed mainly by observing historic effects of liquefaction in varied field settings. The most important criterion is that a seismic liquefaction origin requires widespread, regional development of features around a core area where the effects are most severe. In addition, the features must have a morphology that is consistent with a very sudden application of a large hydraulic force. This article discusses case studies in widely separated and different geological settings: coastal South Carolina, the New Madrid seismic zone, the Wabash Valley seismic zone, and coastal Washington State. These studies encompass most of the range of settings and the types of liquefaction-induced features likely to be encountered anywhere. The case studies describe the observed features and the logic for assigning a seismic liquefaction origin to them. Also discussed are some types of sediment deformations that can be misinterpreted as having a seismic origin. Two independent methods for estimating prehistoric magnitude are discussed briefly. One method is based on determination of the maximum distance from the epicenter over which liquefaction-induced effects have formed. The other method is based on use of geotechnical engineering techniques at sites of marginal liquefaction, in order to bracket the peak accelerations as a function of epicentral distance; these accelerations can then be compared with predictions from seismological models.

Engineering Geology↗

USGS2015

Earthquake geology of the Bulnay Fault (Mongolia)

The Bulnay earthquake of July 23, 1905 (M w 8.3-8.5), in north-central Mongolia, is one of the world's largest recorded intracontinental earthquakes and one of four great earthquakes that occurred in the region during the 20th century. The 375-km-long surface rupture of the left-lateral, strike-slip, N095°E trending Bulnay Fault associated with this earthquake is remarkable for its pronounced expression across the landscape and for the size of features produced by previous earthquakes. Our field observations suggest that in many areas the width and geometry of the rupture zone is the result of repeated earthquakes; however, in those areas where it is possible to determine that the geomorphic features are the result of the 1905 surface rupture alone, the size of the features produced by this single earthquake are singular in comparison to most other historical strike-slip surface ruptures worldwide. Along the 80 km stretch, between 97.18°E and 98.33°E, the fault zone is characterized by several meters width and the mean left-lateral 1905 offset is 8.9 ± 0.6 m with two measured cumulative offsets that are twice the 1905 slip. These observations suggest that the displacement produced during the penultimate event was similar to the 1905 slip. Morphotectonic analyses carried out at three sites along the eastern part of the Bulnay fault, allow us to estimate a mean horizontal slip rate of 3.1 ± 1.7 mm/yr over the Late Pleistocene-Holocene period. In parallel, paleoseismological investigations show evidence for two earthquakes prior to the 1905 event with recurrence intervals of ~2700-4000 years.

Bulnay Fault↗

USGS2017

Evaluating a kinematic method for generating broadband ground motions for great subduction zone earthquakes: Application to the 2003 Mw 8.3 Tokachi‐Oki earthquake

We compare broadband synthetic seismograms with recordings of the 2003 M w "> M w Mw 8.3 Tokachi‐Oki earthquake to evaluate a compound rupture model, in which slip on the fault consists of multiple high‐stress‐drop asperities superimposed on a background slip distribution with longer rise times. Low‐frequency synthetics ( < 1    Hz "> < 1 Hz <1 Hz ) are calculated using deterministic, 3D finite‐difference simulations and are combined with high‐frequency ( > 1    Hz "> > 1 Hz >1 Hz ) stochastic synthetics using a matched filter at 1 Hz. We show that this compound rupture model and overall approach accurately reproduces waveform envelopes and observed response spectral accelerations (SAs) from the Tokachi‐Oki event. We find that sufficiently short subfault rise times (i.e., < ∼ 1 – 2    s "> < ∼ 1 – 2 s <∼1–2 s ) are necessary to reproduce energy ∼ 1    Hz "> ∼ 1 Hz ∼1 Hz . This is achieved by either (1) including distinct subevents with short rise times, as may be suggested by the Tokachi‐Oki data, or (2) imposing a fast‐slip velocity over the entire rupture area. We also include a systematic study on the effects of varying several kinematic rupture parameters. We find that simulated strong ground motions are sensitive to the average rupture velocity and coherence of the rupture front, with more coherent ruptures yielding higher response SAs. We also assess the effects of varying the average slip velocity and the character (i.e., area, magnitude, and location) of high‐stress‐drop subevents. Even in the absence of precise constraints on these kinematic rupture parameters, our simulations still reproduce major features in the Tokachi‐Oki earthquake data, supporting its accuracy in modeling future large earthquakes.

Bulletin of the Seismological Society of America↗

USGS2025

Magnitude, depth and methodological variations of spectral stress drop within the SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest Earthquake Sequence

We present the first ensemble analysis of the 56 different sets of results submitted to the ongoing Community Stress Drop Validation Study using the 2019 Ridgecrest, California, earthquake sequence. Different assumptions and methods result in different estimation of the source contribution to recorded seismograms, and hence to the source parameters (principally corner frequency, f c ⁠ , spectral stress drop, Δσ, and seismic moment, M 0 ⁠ ) obtained from modeling calculated source spectra. For earthquakes smaller than magnitude (M) 2.5 there is negligible correlation between the f c values obtained by different studies, implying that no present method is reliable using available data. For larger magnitude events, correlation between f c measurements of different studies, within even a small M range is always higher than spectral ⁠Δσ , because the f c measurements simply reflect the underlying physical decrease in f c with increasing M. We model the observed trends of submitted f c with both magnitude and depth. Most methods report an increase in spectral Δσ with M, although a magnitude‐invariant spectral Δσ is within the confidence limits. The depth dependence is smaller and depends on whether a study allows attenuation to vary with source depth; a combination of depth‐dependent attenuation correction, and depth‐dependent shear‐wave velocity can compensate for reported depth trends. We model the submitted values to remove differing M and depth variation to investigate the relative interevent variability. We find consistent relative variation between individual events, and also lower relative spectral Δσ in the northwest of the aftershock sequence, and higher on the cross fault and in the region of main fault intersection. This large‐scale comparison implies that absolute spectral Δσ estimates are dependent on the methods used; studies of different regions or using different methods should not be directly compared and improved constraints on path and site corrections are needed to resolve these absolute spectral Δσ differences.

California↗

USGS1994

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

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

Bulletin of the Seismological Society of America↗

USGS2017

Influence of lithostatic stress on earthquake stress drops in North America

We estimate stress drops for earthquakes in and near the continental United States using the method of spectral ratios. The ratio of acceleration spectra between collocated earthquakes recorded at a given station removes the effects of path and recording site and yields source parameters including corner frequency for, and the ratio of seismic moment between, the two earthquakes. We determine stress drop from these parameters for 1121 earthquakes greater than M ∼3 in 60 earthquake clusters. We find that the average Brune stress drop for the few eastern United States (EUS) tectonic mainshocks studied (2.6–36 MPa) is about three times greater than that of tectonic mainshocks in the western United States (WUS, 1.0–7.9 MPa) and five times greater than mainshocks potentially induced by wastewater injection in the central United States (CUS, 0.6–5.6 MPa). EUS events tend to be deeper thrusting events, whereas WUS events tend to be shallower but have a wide range of focal mechanisms. CUS events tend to be shallow with strike‐slip to normal‐faulting mechanisms. With the possible exception of CUS aftershocks, we find that differences in stress drop among all events can be taken into account, within one standard deviation of significance, by differences in the shear failure stress as outlined by Mohr–Coulomb theory. The shear failure stress is a function of vertical stress (or depth), the fault style (normal, strike slip, or reverse), and coefficient of friction (estimated here to be, on average, 0.64). After accounting for faulting style and depth dependence, we find that the average Brune stress drop is about 3% of the failure stress. These results suggest that high‐frequency shaking hazard (>∼1 Hz) from shallow induced events and aftershocks is reduced to some extent by lower stress drop. However, the shallow hypocenters will increase hazard within several kilometers of the source.

Bulletin of the Seismological Society of America↗

USGS2002

Triggered surface slips in the Salton Trough associated with the 1999 Hector Mine, California, earthquake

Surface fracturing occurred along the southern San Andreas, Superstition Hills, and Imperial faults in association with the 16 October 1999 (Mw 7.1) Hector Mine earthquake, making this at least the eighth time in the past 31 years that a regional earthquake has triggered slip along faults in the Salton Trough. Fractures associated with the event formed discontinuous breaks over a 39-km-long stretch of the San Andreas fault, from the Mecca Hills southeastward to Salt Creek and Durmid Hill, a distance from the epicenter of 107 to 139 km. Sense of slip was right lateral; only locally was there a minor (~1 mm) vertical component of slip. Dextral slip ranged from 1 to 13 mm. Maximum slip values in 1999 and earlier triggered slips are most common in the central Mecca Hills. Field evidence indicates a transient opening as the Hector Mine seismic waves passed the southern San Andreas fault. Comparison of nearby strong-motion records indicates several periods of relative opening with passage of the Hector Mine seismic wave-a similar process may have contributed to the field evidence of a transient opening. Slip on the Superstition Hills fault extended at least 9 km, at a distance from the Hector Mine epicenter of about 188 to 196 km. This length of slip is a minimum value, because we saw fresh surface breakage extending farther northwest than our measurement sites. Sense of slip was right lateral; locally there was a minor (~1 mm) vertical component of slip. Dextral slip ranged from 1 to 18 mm, with the largest amounts found distributed (or skewed) away from the Hector Mine earthquake source. Slip triggered on the Superstition Hills fault commonly is skewed away from the earthquake source, most notably in 1968, 1979, and 1999. Surface slip on the Imperial fault and within the Imperial Valley extended about 22 km, representing a distance from the Hector Mine epicenter of about 204 to 226 km. Sense of slip dominantly was right lateral; the right-lateral component of slip ranged from 1 to 19 mm. Locally there was a minor (~1-2 mm) vertical component of slip; larger proportions of vertical slip (up to 10 mm) occurred in Mesquite basin, where scarps indicate long-term oblique-slip motion for this part of the Imperial fault. Slip triggered on the Imperial fault appears randomly distributed relative to location along the fault and source direction. Multiple surface slips, both primary and triggered slip, indicate that slip repeatedly is small at locations of structural complexity.

Bulletin of the Seismological Society of America↗

USGS2010

Global significance of a sub-Moho boundary layer (SMBL) deduced from high-resolution seismic observations

We infer the fine structure of a sub-Moho boundary layer (SMBL) at the top of the lithospheric mantle from high-resolution seismic observations of Peaceful Nuclear Explosions (PNE) on superlong-range profiles in Russia. Densely recorded seismograms permit recognition of previously unknown features of teleseismic propagation of the well known Pn and Sn phases, such as a band of incoherent, scattered, high-frequency seismic energy, developing consistently from station to station, apparent velocities of sub-Moho material, and high-frequency energy to distances of more than 3000 km with a coda band, incoherent at 10 km spacing and yet consistently observed to the end of the profiles. Estimates of the other key elements of the SMBL were obtained by finite difference calculations of wave propagation in elastic 2D models from a systematic grid search through parameter space. The SMBL consists of randomly distributed, mild velocity fluctuations of 2% or schlieren of high aspect ratios (???40) with long horizontal extent (???20 km) and therefore as thin as 0.5 km only; SMBL thickness is 60-100 km. It is suggested that the SMBL is of global significance as the physical base of the platewide observed high-frequency phases Pn and Sn. It is shown that wave propagation in the SMBL waveguide is insensitive to the background velocity distribution on which its schlieren are superimposed. This explains why the Pn and Sn phases traverse geological provinces of various age, heat flow, crustal thickness, and tectonic regimes. Their propagation appears to be independent of age. temperature, pressure, and stress. Dynamic stretching of mantle material during subduction or flow, possibly combined with chemical differentiation have to be considered as scale-forming processes in the upper mantle. However, it is difficult to distinguish with the present sets of Pn/Sn array data whether (and also where) the boundary layer is a frozen-in feature of paleo-processes or whether it is a response to an on-going processes; nevertheless, the derived quantitative estimates of the SMBL properties provide important constraints for any hypothesis on scale-forming processes. Models to be tested by future numerical and field experiments are, for example, repeated subduction-convection stretching of oceanic lithosphere (marble-cake model) and schlieren formation at mid-ocean ridges. It is also proposed that the modeling of the observed blocking of Sn and Pn propagation at active plate margins offers a new tool to study the depth range of tectonics below the crust-mantle boundary. Finally, the deduced schlieren structure of the SMBL closes an important scale gap of three to four orders of magnitude between structural dimensions studied in petrological analysis of mantle samples (xenoliths or outcrop of oceanic lithosphere) and those imaged in classical seismological studies of the lithosphere.

International Geology Review↗

USGS2022

Empirical map-based nonergodic models of site response in the greater Los Angeles area

We develop empirical estimates of site response at seismic stations in the Los Angeles area using recorded ground motions from 414 M 3–7.3 earthquakes in southern California. The data are from a combination of the Next Generation Attenuation‐West2 project, the 2019 Ridgecrest earthquakes, and about 10,000 newly processed records. We estimate site response using an iterative mixed‐effects residuals partitioning approach, accounting for azimuthal variations in anelastic attenuation and potential bias due to spatial clusters of colocated earthquakes. This process yields site response for peak ground acceleration, peak ground velocity, and pseudospectral acceleration relative to a 760 m/s shear‐wave velocity ( ⁠ V S "> V V S "> s ) reference condition. We employ regression kriging to generate a spatially continuous site response model, using the linear site and basin terms from Boore et al. (2014) as the background model, which depend on V S "> V V S "> s V S "> 30 ⁠ and depth to the 1 km/s V S "> V V S "> s isosurface. This is different from past approaches to nonergodic models, in which spatially varying coefficients are regressed. We validate the model using stations in the Community Seismic Network (CSN) that are in the middle of our model spatial domain but were not considered in model development, finding strong agreement between the interpolated model and CSN data for long periods. Our model could be implemented in regional seismic hazard analyses, which would lead to improvements especially at long return periods. Our site response model also has potential to improve both ground‐motion accuracy and warning times for the U.S. Geological Survey ShakeAlert earthquake early warning (EEW) system. For a point‐source EEW simulation of the 1994 M 6.7 Northridge earthquake, our model produces ground motions more consistent with the ground‐truth ShakeMap and would alert areas with high population density such as downtown Los Angeles at lower estimated magnitudes (i.e., sooner) than an ergodic model for a modified Mercalli intensity 4.5 alerting threshold.

California↗

USGS2020

Capturing, preserving and digitizing legacy seismic data from the Montserrat Volcano Observatory analog seismic network, July 1995 – December 2004

An eruption of the Soufrière Hills Volcano (SHV) on the eastern Caribbean island of Montserrat began on 18 July 1995 and continued until February 2010. Within nine days of the eruption onset, an existing four‐station analog seismic network (ASN) was expanded to 10 sites. Telemetered data from this network were recorded, processed, and archived locally using a system developed by scientists from the U.S. Geological Survey (USGS) Volcano Disaster Assistance Program (VDAP). In October 1996, a digital seismic network (DSN) was deployed with the ability to capture larger amplitude signals across a broader frequency range. These two networks operated in parallel until December 2004, with separate telemetry and acquisition systems (analysis systems were merged in March 2001). Although the DSN provided better quality data for research, the ASN featured superior real‐time monitoring tools and captured valuable data including the only seismic data from the first 15 months of the eruption. These successes of the ASN have been rather overlooked. This article documents the evolution of the ASN, the VDAP system, the original data captured, and the recovery and conversion of more than 230,000 seismic events from legacy SUDS, Hypo71, and Seislog formats into Seisan database with waveform data in miniSEED format. No digital catalog existed for these events, but students at the University of South Florida have classified two-thirds of the 40,000 events that were captured between July 1995 and October 1996. Locations and magnitudes were recovered for ~10,000 of these events. Real-time seismic amplitude measurement, seismic spectral amplitude measurement, and tiltmeter data were also captured. The result is that the ASN seismic dataset is now more discoverable, accessible, and reusable, in accordance with FAIR data principles. These efforts could catalyze new research on the 1995–2010 SHV eruption. Furthermore, many observatories have data in these same legacy data formats and might benefit from procedures and codes documented here.

Soufrière Hills Volcano↗