Source parameters of the 1949 magnitude 7.1 south Puget Sound, Washington, earthquake as determined from long-period body waves and strong ground motions

1987 ◽  
Vol 77 (5) ◽  
pp. 1530-1557
Author(s):  
Glenn Eli Baker ◽  
Charles A. Langston
Keyword(s):  
Focal Mechanism ◽  
Source Parameters ◽  
Puget Sound ◽  
Vertical Component ◽  
Test Laboratory ◽  
Body Waves ◽  
Polarization Angle ◽  
S Wave ◽  
Fault Surface ◽  
Strong Ground

Abstract Teleseismic P, SH, and SV first motions and SH to SV amplitude ratios recorded at eight teleseismic receivers from the 1949 magnitude 7.1 Olympia, Washington, earthquake in combination with data from three stations at regional distances were utilized in a grid testing routine to constrain focal mechanism. Identification of the pP phase places the event at 54 km depth. Distinct pulses, assumed to be source effects, are observed in the far-field waveforms. Analysis of these pulses for directivity made possible discrimination between the fault and auxiliary planes. The plane taken to represent the fault surface strikes east-west ± 15°, dips 45° ± 15° to the north, and has nearly pure left-lateral slip. The preferred source model has an eastward propagation of 40 km. Surface reflections of successive source pulses suggest an upward component of propagation of 5 km. Bounds on the earthquake location and rupture of the 13 April event were determined using depth and source mechanism constraints from the teleseismic study and characteristics of local strong ground motion recordings. The 9-sec S-instrument trigger time seen in the Seattle acceleration recordings places the event at least 60 km from Seattle. Strong motion velocity at the Olympia Highway Test Laboratory is characterized by an impulsive and rectilinear S wave. The low amplitude of the vertical component of initial S motion suggests that either the epicenter is within 5 km of the Olympia Highway Test Laboratory for a pure incident SV wave or located along an azimuth of N159° if the wave is SH. The combined constraint of minimum distance from Seattle and the S polarization angle implied by the teleseismic data focal mechanism places the initiation of rupture 5 to 10 km north to north-northwest of the Olympia Highway Test Laboratory at 47.13°N, 122.95°W. This is approximately 20 km west of previously determined epicenters. The T axis, gently dipping to the southeast, supports other evidence that the Juan de Fuca plate dips to the southeast in a zone between segments of the plate north and south of the event's location. The fault plane's slip is taken to indicate that subduction is still active beneath Washington and that motion of the two segments is probably independent.

Strain to Ground Motion Conversion of DAS Data for Earthquake Magnitude and Stress Drop Determination

2021 ◽  
Author(s):  
Itzhak Lior ◽  
Anthony Sladen ◽  
Diego Mercerat ◽  
Jean-Paul Ampuero ◽  
Diane Rivet ◽  
...  
Keyword(s):  
Early Warning ◽  
Stress Drop ◽  
Source Parameters ◽  
Simulated Data ◽  
Earthquake Early Warning ◽  
Body Waves ◽  
Ocean Bottom ◽  
Correlated Noise ◽  
S Wave ◽  
Earthquake Source Parameters

<p>The use of Distributed Acoustic Sensing (DAS) presents unique advantages for earthquake monitoring compared with standard seismic networks: spatially dense measurements adapted for harsh environments and designed for remote operation. However, the ability to determine earthquake source parameters using DAS is yet to be fully established. In particular, resolving the magnitude and stress drop, is a fundamental objective for seismic monitoring and earthquake early warning. To apply existing methods for source parameter estimation to DAS signals, they must first be converted from strain to ground motions. This conversion can be achieved using the waves’ apparent phase velocity, which varies for different seismic phases ranging from fast body-waves to slow surface- and scattered-waves. To facilitate this conversion and improve its reliability, an algorithm for slowness determination is presented, based on the local slant-stack transform. This approach yields a unique slowness value at each time instance of a DAS time-series. The ability to convert strain-rate signals to ground accelerations is validated using simulated data and applied to several earthquakes recorded by dark fibers of three ocean-bottom telecommunication cables in the Mediterranean Sea. The conversion emphasizes fast body-waves compared to slow scattered-waves and ambient noise, and is robust even in the presence of correlated noise and varying wave propagation directions. Good agreement is found between source parameters determined using converted DAS waveforms and on-land seismometers for both P- and S-wave records. The demonstrated ability to resolve source parameters using P-waves on horizontal ocean-bottom fibers is key for the implementation of DAS based earthquake early warning, which will significantly improve hazard mitigation capabilities for offshore and tsunami earthquakes.</p>

A computer program for focal mechanism determination combining P and S wave data

1969 ◽  
Vol 59 (2) ◽  
pp. 503-519
Author(s):  
Agustin Udias ◽  
Dieter Baumann
Keyword(s):  
Computer Program ◽  
Focal Mechanism ◽  
Total Error ◽  
Source Model ◽  
Polarization Angle ◽  
S Wave ◽  
Period Range ◽  
S Waves ◽  
Double Couple ◽  
The Relationship

abstract A computer program has been developed to find the orientation of a double couple source model for the mechanism of an earthquake which best satisfies the data from P and S waves. The relationship between the two axes of the solution given by the equations for the polarization angle of S is used in order to rapidly find the orientation of the source model for which a total error value involving the error of S and P data is a minimum. The program gives best results for data from homogeneous instruments of similar period range. Solutions for three earthquakes, selected because of the orientation of the source, are presented and the reliability of their solutions under ideal conditions is discussed.

scholarly journals Strain to Ground Motion Conversion of DAS Data for Earthquake Magnitude and Stress Drop Determination

2021 ◽  
Author(s):  
Itzhak Lior ◽  
Anthony Sladen ◽  
Diego Mercerat ◽  
Jean-Paul Ampuero ◽  
Diane Rivet ◽  
...  
Keyword(s):  
Early Warning ◽  
Stress Drop ◽  
Source Parameters ◽  
Simulated Data ◽  
Earthquake Early Warning ◽  
Body Waves ◽  
Ocean Bottom ◽  
Correlated Noise ◽  
S Wave ◽  
Earthquake Source Parameters

Abstract. The use of Distributed Acoustic Sensing (DAS) presents unique advantages for earthquake monitoring compared with standard seismic networks: spatially dense measurements adapted for harsh environments and designed for remote operation. However, the ability to determine earthquake source parameters using DAS is yet to be fully established. In particular, resolving the magnitude and stress drop, is a fundamental objective for seismic monitoring and earthquake early warning. To apply existing methods for source parameter estimation to DAS signals, they must first be converted from strain to ground motions. This conversion can be achieved using the waves' apparent phase velocity, which varies for different seismic phases ranging from fast body-waves to slow surface- and scattered-waves. To facilitate this conversion and improve its reliability, an algorithm for slowness determination is presented, based on the local slant-stack transform. This approach yields a unique slowness value at each time instance of a DAS time-series. The ability to convert strain-rate signals to ground accelerations is validated using simulated data and applied to several earthquakes recorded by dark fibers of three ocean-bottom telecommunication cables in the Mediterranean Sea. The conversion emphasizes fast body-waves compared to slow scattered-waves and ambient noise, and is robust even in the presence of correlated noise and varying wave propagation directions. Good agreement is found between source parameters determined using converted DAS waveforms and on-land seismometers for both P- and S-wave records. The demonstrated ability to resolve source parameters using P-waves on horizontal ocean-bottom fibers is key for the implementation of DAS based earthquake early warning, which will significantly improve hazard mitigation capabilities for offshore and tsunami earthquakes.

The High-Frequency Decay Slope of Spectra (Kappa) for M≥3.5 Earthquakes on Rock Sites in Eastern and Western Canada

2020 ◽  
Vol 110 (2) ◽  
pp. 471-488 ◽  
Author(s):  
Samantha M. Palmer ◽  
Gail M. Atkinson
Keyword(s):  
Ground Motion ◽  
Classical Method ◽  
Vertical Component ◽  
Western Canada ◽  
S Wave ◽  
Motion Modeling ◽  
Eastern Canada ◽  
General Observation ◽  
Anelastic Attenuation ◽  
Two Parameters

ABSTRACT Spectral decay of ground-motion amplitudes at high frequencies is primarily influenced by two parameters: site-related kappa (κ0) and regional Q (quality factor, inversely proportional to anelastic attenuation). We examine kappa and apparent Q-values (Qa) for M≥3.5 earthquakes recorded at seismograph stations on rock sites in eastern and western Canada. Our database contains 20 earthquakes recorded on nine stations in eastern Canada and 404 earthquakes recorded on eight stations in western Canada, resulting in 105 and 865 Fourier amplitude spectra, respectively. We apply two different methods: (1) a modified version of the classical S-wave acceleration method; and (2) a new stacking method that is consistent with the use of kappa in ground-motion modeling. The results are robust with respect to the method used and also with respect to the frequency band selected, which ranges from 9 to 38 Hz depending on the region, event, and method. Kappa values obtained from the classical method are consistent with those of the stacked method, but the stacked method provides a lower uncertainty. A general observation is that kappa is usually larger, and apparent Q is smaller, for the horizontal component in comparison to the vertical component. We determine an average regional κ0=7  ms (horizontal) and 0 ms (vertical) for rock sites in eastern Canada; we obtain κ0=19  ms (horizontal) and 14 ms (vertical) for rock sites in western Canada. We note that kappa measurements are quite sensitive to details of data selection criteria and methodology, and may be significantly influenced by site effects, resulting in large site-to-site variability.

Inversion of regional surface-wave spectra for source parameters of aftershocks from the 1992 Petrolia earthquake sequence

1995 ◽  
Vol 85 (3) ◽  
pp. 705-715
Author(s):  
Mark Andrew Tinker ◽  
Susan L. Beck
Keyword(s):  
Surface Wave ◽  
Focal Mechanism ◽  
Source Parameters ◽  
Strike Slip ◽  
Earthquake Sequence ◽  
Fault System ◽  
Accretionary Wedge ◽  
Wave Spectra ◽  
The North ◽  
Gorda Plate

Abstract Regional distance surface waves are used to study the source parameters for moderate-size aftershocks of the 25 April 1992 Petrolia earthquake sequence. The Cascadia subduction zone had been relatively seismically inactive until the onset of the mainshock (Ms = 7.1). This underthrusting event establishes that the southern end of the North America-Gorda plate boundary is seismogenic. It was followed by two separate and distinct large aftershocks (Ms = 6.6 for both) occurring at 07:41 and 11:41 on 26 April, as well as thousands of other small aftershocks. Many of the aftershocks following the second large aftershock had magnitudes in the range of 4.0 to 5.5. Using intermediate-period surface-wave spectra, we estimate focal mechanisms and depths for one foreshock and six of the larger aftershocks (Md = 4.0 to 5.5). These seven events can be separated into two groups based on temporal, spatial, and principal stress orientation characteristics. Within two days of the mainshock, four aftershocks (Md = 4 to 5) occurred within 4 hr of each other that were located offshore and along the Mendocino fault. These four aftershocks comprise one group. They are shallow, thrust events with northeast-trending P axes. We interpret these aftershocks to represent internal compression within the North American accretionary prism as a result of Gorda plate subduction. The other three events compose the second group. The shallow, strike-slip mechanism determined for the 8 March foreshock (Md = 5.3) may reflect the right-lateral strike-slip motion associated with the interaction between the northern terminus of the San Andreas fault system and the eastern terminus of the Mendocino fault. The 10 May aftershock (Md = 4.1), located on the coast and north of the Mendocino triple junction, has a thrust fault focal mechanism. This event is shallow and probably occurred within the accretionary wedge on an imbricate thrust. A normal fault focal mechanism is obtained for the 5 June aftershock (Md = 4.8), located offshore and just north of the Mendocino fault. This event exhibits a large component of normal motion, representing internal failure within a rebounding accretionary wedge. These two aftershocks and the foreshock have dissimilar locations in space and time, but they do share a north-northwest oriented P axis.

Stress Drop Derived from Spectral Analysis Considering the Hypocentral Depth in the Attenuation Model: Application to the Ridgecrest Region, California

2021 ◽  
Author(s):  
Dino Bindi ◽  
Hoby N. T. Razafindrakoto ◽  
Matteo Picozzi ◽  
Adrien Oth
Keyword(s):  
Stress Drop ◽  
Spectral Decomposition ◽  
Epicentral Distance ◽  
Source Parameters ◽  
S Wave ◽  
Qualitative Comparison ◽  
Ground Shaking ◽  
Attenuation Model ◽  
Hypocentral Distance ◽  
The Impact

ABSTRACT We investigate the impact of considering a depth-dependent attenuation model on source parameters assessed through a spectral decomposition. In particular, we evaluate the effect of considering the hypocentral depth as an additional variable for the attenuation model, using as the target the tendency of the average stress drop to increase with depth, as observed in recent studies. We analyze the Fourier spectra of S-wave windows for about 1900 earthquakes with a magnitude above 2.5 recorded in the Ridgecrest region, southern California. Two different parameterizations of the attenuation term are implemented in the spectral decomposition, either as a function of the hypocentral distance alone or as a function of both epicentral distance and depth. The comparison of the spectral attenuation curves shows that, although the hypocentral model describes, on average, the range of values spanned by the attenuation curve for different depths, systematic differences with distance, depth, and frequency are observed. These differences are transferred to the source spectra and, in turn, to the source parameters extracted from the best-fitting ω−2 models. In particular, stress drops for events deeper than 7 km are, on average, almost double even when depth is introduced explicitly in the attenuation model. The increase of stress drop with depth is confirmed also after accounting for the increase of the shear velocity with depth, which absorbs about 30%–40% of the total increase. Moreover, a qualitative comparison with a model for the gradient of the effective normal stress confirms the reliability of the observed trend. Finally, the coherent spatial patterns shown by a simplified 2D tomographic representation of the spectral residuals highlights the impact on ground-shaking variability of the lateral variability of the crustal attenuation properties in the region.

Nonlinear behavior of soil sediments identified by using borehole records observed at the Ashigra Valley, Japan

1995 ◽  
Vol 85 (6) ◽  
pp. 1821-1834
Author(s):  
Toshimi Satoh ◽  
Toshiaki Sato ◽  
Hiroshi Kawase
Keyword(s):  
Shear Strain ◽  
Main Part ◽  
Nonlinear Behavior ◽  
Strong Motion ◽  
S Wave ◽  
Ground Shaking ◽  
Wave Velocities ◽  
The One ◽  
Soil Sediments ◽  
Strong Ground

Abstract We evaluate the nonlinear behavior of soil sediments during strong ground shaking based on the identification of their S-wave velocities and damping factors for both the weak and strong motions observed on the surface and in a borehole at Kuno in the Ashigara Valley, Japan. First we calculate spectral ratios between the surface station KS2 and the borehole station KD2 at 97.6 m below the surface for the main part of weak and strong motions. The predominant period for the strong motion is apparently longer than those for the weak motions. This fact suggests the nonlinearity of soil during the strong ground shaking. To quantify the nonlinear behavior of soil sediments, we identify their S-wave velocities and damping factors by minimizing the residual between the observed spectral ratio and the theoretical amplification factor calculated from the one-dimensional wave propagation theory. The S-wave velocity and the damping factor h (≈(2Q)−1) of the surface alluvial layer identified from the main part of the strong motion are about 10% smaller and 50% greater, respectively, than those identified from weak motions. The relationships between the effective shear strain (=65% of the maximum shear strain) calculated from the one-dimensional wave propagation theory and the shear modulus reduction ratios or the damping factors estimated by the identification method agree well with the laboratory test results. We also confirm that the soil model identified from a weak motion overestimates the observed strong motion at KS2, while that identified from the strong motion reproduces the observed. Thus, we conclude that the main part of the strong motion, whose maximum acceleration at KS2 is 220 cm/sec2 and whose duration is 3 sec, has the potential of making the surface soil nonlinear at an effective shear strain on the order of 0.1%. The S-wave velocity in the surface alluvial layer identified from the part just after the main part of the strong motion is close to that identified from weak motions. This result suggests that the shear modulus recovers quickly as the shear strain level decreases.

scholarly journals Triplicated P-wave measurements for waveform tomography of the mantle transition zone

Solid Earth ◽  
2012 ◽  
Vol 3 (2) ◽  
pp. 339-354 ◽  
Author(s):  
S. C. Stähler ◽  
K. Sigloch ◽  
T. Nissen-Meyer
Keyword(s):  
Transition Zone ◽  
Epicentral Distance ◽  
Body Wave ◽  
Source Parameters ◽  
Body Waves ◽  
P Wave ◽  
Theoretical Modelling ◽  
Finite Frequency ◽  
Mantle Transition Zone ◽  
Band Limited

Abstract. Triplicated body waves sample the mantle transition zone more extensively than any other wave type, and interact strongly with the discontinuities at 410 km and 660 km. Since the seismograms bear a strong imprint of these geodynamically interesting features, it is highly desirable to invert them for structure of the transition zone. This has rarely been attempted, due to a mismatch between the complex and band-limited data and the (ray-theoretical) modelling methods. Here we present a data processing and modelling strategy to harness such broadband seismograms for finite-frequency tomography. We include triplicated P-waves (epicentral distance range between 14 and 30°) across their entire broadband frequency range, for both deep and shallow sources. We show that is it possible to predict the complex sequence of arrivals in these seismograms, but only after a careful effort to estimate source time functions and other source parameters from data, variables that strongly influence the waveforms. Modelled and observed waveforms then yield decent cross-correlation fits, from which we measure finite-frequency traveltime anomalies. We discuss two such data sets, for North America and Europe, and conclude that their signal quality and azimuthal coverage should be adequate for tomographic inversion. In order to compute sensitivity kernels at the pertinent high body wave frequencies, we use fully numerical forward modelling of the seismic wavefield through a spherically symmetric Earth.

scholarly journals Metode Seismik dalam Pemodelan Fisis Letusan Gunung Lumpur Bledug Kuwu

2019 ◽  
Vol 2 (2) ◽  
pp. 61-66
Author(s):  
Ahmad Fauzi Pohan ◽  
Rusnoviandi Rusnoviandi
Keyword(s):  
Wave Velocity ◽  
Physical Modeling ◽  
Vertical Component ◽  
Dominant Frequency ◽  
P Wave ◽  
Physical Parameters ◽  
S Wave ◽  
Interesting Phenomenon ◽  
Central Java ◽  
P Wave Velocity

Aktivitas gunung lumpur Bledug Kuwu di Jawa  Tengah merupakan fenomena yang menarik dikaji menggunakan pemodelan fisis. Tujuan penelitian ini adalah mengetahui parameter dari medium gunung lumpur Bledug Kuwu. Adapun pemodelan fisis yang dilakukan dengan menggunakan media fisis akuarium berukuran 59 × 59 × 37,3 cm yang diisi material dari lumpur Bledug Kuwu. Sumber letusan dihasilkan dari tekanan kompresor yang dapat diatur kedalaman (10.5, 13, dan 15.5 cm) dan sudut (30o, 45o dan 60o) sumbernya. Sensor yang digunakan geophone komponen vertikal sebanyak 3 buah dengan durasi perekaman selama 5 dan 2,5 detik. Data diambil dengan frekuensi sampel 2 dan 4 kHz untuk masing-masing durasi perekaman. Konfigurasi sumber dan geophone dibuat sesuai dengan pemodelan fisisnya. Pengukuran desnsitas lumpur menunjukkan angka sebesar 1200 kg/m3. Berdasarkan hasil analisis seismogram model fisis diperoleh kecepatan perambatan gelombang-P pada medium lumpur Bledug Kuwu adalah sebesar 48,74 m/s,dan gelombang-S sebesar 28,14 m/s dengan frekuensi dominan antara 20 sampai 25 Hz.   Bledug Kuwu mud volcano activity in Central Java is an interesting phenomenon to be studied using both physical  modeling. The objective of this study was to determine the physical parameters of the medium of Bledug Kuwu. The Physical model was an aquarium with a dimension of 59 × 59 × 37.3 cm filled with Bledug Kuwu’s mud. The eruption source is generated by a compressor pressure that can be controled both the depth(10.5, 13, and 15.5 cm) and the angel of the source (30o, 45o and 60o). The resulting seismic signals were recorded by using 3 vertical component geophones for 10 and 5 seconds durations at a frequency of 2 and 4 kHz respectivel, mud density 1200 kg/m3 . The physical modeling shows that the P-wave velocity of the Bledug Kuwu’s medium is 48.7 m/s, S-wave velocity of Bledug Kuwu’s is 28,14 m/s  with a dominant frequency of 20 to 25 Hz.

scholarly journals The use of parallel computing for the high-resolution determination of earthquake source parameters

2019 ◽  
pp. 68-75
Author(s):  
A. S. Fomochkina ◽  
V. G. Bukchin
Keyword(s):  
Parallel Computing ◽  
Focal Mechanism ◽  
Source Parameters ◽  
Nodal Plane ◽  
Earthquake Parameters ◽  
Earthquake Source Parameters ◽  
Detailed Search ◽  
High Degree ◽  
Special Case

Alongside the determination of the focal mechanism and source depth of an earthquake by direct examination of their probable values on a grid in the parameter space, also the resolution of these determinations can be estimated. However, this approach requires considerable time in the case of a detailed search. A special case of a shallow earthquake whose one nodal plane is subhorizontal is an example of the sources that require the use of a detailed grid. For studying these events based on the records of the long-period surface waves, the grids with high degree of detail in the angles of the focal mechanism are required. We discuss the application of the methods of parallel computing for speeding up the calculations of earthquake parameters and present the results of studying the strongest aftershock of the Tohoku, Japan, earthquake by this approach.

Sign in / Sign up

Export Citation Format

Share Document