Moment-tensor solutions for the 24 November 1987 Superstition Hills, California, earthquakes

1989 ◽  
Vol 79 (2) ◽  
pp. 493-499
Author(s):  
Stuart A. Sipkin

Abstract The teleseismic long-period waveforms recorded by the Global Digital Seismograph Network from the two largest Superstition Hills earthquakes are inverted using an algorithm based on optimal filter theory. These solutions differ slightly from those published in the Preliminary Determination of Epicenters Monthly Listing because a somewhat different, improved data set was used in the inversions and a time-dependent moment-tensor algorithm was used to investigate the complexity of the main shock. The foreshock (origin time 01:54:14.5, mb 5.7, Ms 6.2) had a scalar moment of 2.3 × 1025 dyne-cm, a depth of 8 km, and a mechanism of strike 217°, dip 79°, rake 4°. The main shock (origin time 13:15:56.4, mb 6.0, Ms 6.6) was a complex event, consisting of at least two subevents, with a combined scalar moment of 1.0 × 1026 dyne-cm, a depth of 10 km, and a mechanism of strike 303°, dip 89°, rake −180°.

1990 ◽  
Vol 80 (5) ◽  
pp. 1205-1231
Author(s):  
Jiajun Zhang ◽  
Thorne Lay

Abstract Determination of shallow earthquake source mechanisms by inversion of long-period (150 to 300 sec) Rayleigh waves requires epicentral locations with greater accuracy than that provided by routine source locations of the National Earthquake Information Center (NEIC) and International Seismological Centre (ISC). The effects of epicentral mislocation on such inversions are examined using synthetic calculations as well as actual data for three large Mexican earthquakes. For Rayleigh waves of 150-sec period, an epicentral mislocation of 30 km introduces observed source spectra phase errors of 0.6 radian for stations at opposing azimuths along the source mislocation vector. This is larger than the 0.5-radian azimuthal variation of the phase spectra at the same period for a thrust fault with 15° dip and 24-km depth. The typical landward mislocation of routinely determined epicenters of shallow subduction zone earthquakes causes source moment tensor inversions of long-period Rayleigh waves to predict larger fault dip than indicated by teleseismic P-wave first-motion data. For dip-slip earthquakes, inversions of long-period Rayleigh waves that use an erroneous source location in the down-dip or along-strike directions of a nodal plane, overestimate the strike, dip, and slip of that nodal plane. Inversions of strike-slip earthquakes that utilize an erroneous location along the strike of a nodal plane overestimate the slip of that nodal plane, causing the second nodal plane to dip incorrectly in the direction opposite to the mislocation vector. The effects of epicentral mislocation for earthquakes with 45° dip-slip fault mechanisms are more severe than for events with other fault mechanisms. Existing earth model propagation corrections do not appear to be sufficiently accurate to routinely determine the optimal surface-wave source location without constraints from body-wave information, unless extensive direct path (R1) data are available or empirical path calibrations are performed. However, independent surface-wave and body-wave solutions can be remarkably consistent when the effects of epicentral mislocation are accounted for. This will allow simultaneous unconstrained body-wave and surface-wave inversions to be performed despite the well known difficulties of extracting the complete moment tensor of shallow sources from fundamental modes.


1982 ◽  
Vol 72 (2) ◽  
pp. 439-456
Author(s):  
Thorne Lay ◽  
Jeffrey W. Given ◽  
Hiroo Kanamori

Abstract The seismic moment and source orientation of the 8 November 1980 Eureka, California, earthquake (Ms = 7.2) are determined using long-period surface and body wave data obtained from the SRO, ASRO, and IDA networks. The favorable azimuthal distribution of the recording stations allows a well-constrained mechanism to be determined by a simultaneous moment tensor inversion of the Love and Rayleigh wave observations. The shallow depth of the event precludes determination of the full moment tensor, but constraining Mzx = Mzy = 0 and using a point source at 16-km depth gives a major double couple for period T = 256 sec with scalar moment M0 = 1.1 · 1027 dyne-cm and a left-lateral vertical strike-slip orientation trending N48.2°E. The choice of fault planes is made on the basis of the aftershock distribution. This solution is insensitive to the depth of the point source for depths less than 33 km. Using the moment tensor solution as a starting model, the Rayleigh and Love wave amplitude data alone are inverted in order to fine-tune the solution. This results in a slightly larger scalar moment of 1.28 · 1027 dyne-cm, but insignificant (<5°) changes in strike and dip. The rake is not well enough resolved to indicate significant variation from the pure strike-slip solution. Additional amplitude inversions of the surface waves at periods ranging from 75 to 512 sec yield a moment estimate of 1.3 ± 0.2 · 1027 dyne-cm, and a similar strike-slip fault orientation. The long-period P and SH waves recorded at SRO and ASRO stations are utilized to determine the seismic moment for 15- to 30-sec periods. A deconvolution algorithm developed by Kikuchi and Kanamori (1982) is used to determine the time function for the first 180 sec of the P and SH signals. The SH data are more stable and indicate a complex bilateral rupture with at least four subevents. The dominant first subevent has a moment of 6.4 · 1026 dyne-cm. Summing the moment of this and the next three subevents, all of which occur in the first 80 sec of rupture, yields a moment of 1.3 · 1027 dyne-cm. Thus, when the multiple source character of the body waves is taken into account, the seismic moment for the Eureka event throughout the period range 15 to 500 sec is 1.3 ± 0.2 · 1027 dyne-cm.


2017 ◽  
Vol 2019 (16) ◽  
pp. 5087-5126 ◽  
Author(s):  
Yavar Kian ◽  
Lauri Oksanen

Abstract Given $(M,g)$, a compact connected Riemannian manifold of dimension $d \geq 2$, with boundary $\partial M$, we study the inverse boundary value problem of determining a time-dependent potential $q$, appearing in the wave equation $\partial_t^2u-\Delta_g u+q(t,x)u=0$ in ${\overline M}=(0,T)\times M$ with $T>0$. Under suitable geometric assumptions we prove global unique determination of $q\in L^\infty({\overline M})$ given the Cauchy data set on the whole boundary $\partial {\overline M}$, or on certain subsets of $\partial {\overline M}$. Our problem can be seen as an analogue of the Calderón problem on the Lorentzian manifold $({\overline M}, dt^2 - g)$.


2020 ◽  
Vol 224 (3) ◽  
pp. 2044-2059
Author(s):  
G M Petersen ◽  
P Niemz ◽  
S Cesca ◽  
V Mouslopoulou ◽  
G M Bocchini

SUMMARY Clusty is a new open source toolbox dedicated to earthquake clustering based on waveforms recorded across a network of seismic stations. Its main application is the study of active faults and the detection and characterization of faults and fault networks. By using a density-based clustering approach, earthquakes pertaining to a common fault can be recognized even over long fault segments, and the first-order geometry and extent of active faults can be inferred. Clusty implements multiple techniques to compute a waveform based network similarity from maximum cross-correlation coefficients at multiple stations. The clustering procedure is designed to be transparent and parameters can be easily tuned. It is supported by a number of analysis visualization tools which help to assess the homogeneity within each cluster and the differences among distinct clusters. The toolbox returns graphical representations of the results. A list of representative events and stacked waveforms facilitate further analyses like moment tensor inversion. Results obtained in various frequency bands can be combined to account for large magnitude ranges. Thanks to the simple configuration, the toolbox is easily adaptable to new data sets and to large magnitude ranges. To show the potential of our new toolbox, we apply Clusty to the aftershock sequence of the Mw 6.9 25 October 2018 Zakynthos (Greece) Earthquake. Thanks to the complex tectonic setting at the western termination of the Hellenic Subduction System where multiple faults and faulting styles operate simultaneously, the Zakynthos data set provides an ideal case-study for our clustering analysis toolbox. Our results support the activation of several faults and provide insight into the geometry of faults or fault segments. We identify two large thrust faulting clusters in the vicinity of the main shock and multiple strike-slip clusters to the east, west and south of these clusters. Despite its location within the largest thrust cluster, the main shock does not show a high waveform similarity to any of the clusters. This is consistent with the results of other studies suggesting a complex failure mechanism for the main shock. We propose the existence of conjugated strike-slip faults in the south of the study area. Our waveform similarity based clustering toolbox is able to reveal distinct event clusters which cannot be discriminated based on locations and/or timing only. Additionally, the clustering results allows distinction between fault and auxiliary planes of focal mechanisms and to associate them to known active faults.


Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. KS11-KS24 ◽  
Author(s):  
Xin Yu ◽  
Scott Leaney ◽  
Jim Rutledge ◽  
Chris Chapman

Moment-tensor inversion under single monitoring well geometries becomes unstable due to the singularity of the inversion matrix. But microseismic events observed during hydraulic fracturing commonly show clusters of events with similar source mechanisms despite differences in the origin time and the magnitude. If the events with similar source mechanisms can be grouped and inverted for a single common moment tensor, the singularity can be eliminated. We have developed a normalized multievent moment-tensor inversion (NME-MTI) method, which does the MTI simultaneously for multiple events, to test the feasibility of this multievent approach. First, the scalar moment for each event was estimated using the far-field low-frequency level at each receiver. Then, the displacements measured at the receivers were normalized by the scalar moment and used to invert for the common moment tensor simultaneously for all the events in the group. We introduced a gradient search method to minimize the overall misfit by adjusting the scalar moment for each event to reduce the errors introduced by the scalar moment estimation. The algorithm was tested with a synthetic data set with four monitoring wells and a field data set with dual monitoring wells. It was proved that the NME-MTI method can retrieve the moment tensors of the event groups with data from a single monitoring well. The effects of uncertainties on the inversion were examined with data noise, scalar moment uncertainty, and event location uncertainty. The results showed that the ME-MTI result is much less sensitive to the data noise and the scalar moment uncertainty than the single-event approach. The results also determined that although the bias to the solutions increases when the event location uncertainty increases, the bias can be controlled by reducing the event location uncertainties using a more accurate location algorithm.


Geophysics ◽  
2020 ◽  
pp. 1-61
Author(s):  
Claudia Finger ◽  
Erik Saenger

An approach is presented to determine the time-dependent moment tensor and the origin time in addition to commonly derived locations of seismic events using time-reverse imaging (TRI). It is crucial to locate and characterize the occurring micro-seismicity without making a priori assumptions about the sources to fully understand the subsurface processes inducing seismicity. Low signal-to-noise ratios often force standard methods to make assumptions about sources or only characterize selected larger-magnitude events. In TRI, micro-earthquakes are located by back propagating the full recorded time-reversed wavefield through a velocity model until it ideally convergences on the source location. Therefore, it is less affected by low signal-to-noise ratios and potentially locates and characterizes most of the events. After distinguishing artificial convergence locations from source locations, the quality of the source location and the moment tensors are derived by recording the stress at the determined source locations during the back propagation of the time-reversed wavefield. A robust workflow is derived using synthetic test cases in a realistic scenario with velocity models that only approximate the true velocity model and/or noisy displacement traces. The influence of a rudimentary velocity model on the source-location accuracy and characterisation is significant. The proposed workflow handles these less-than optimal station distributions and velocity models. Finally, the derived workflow is successfully applied to field data recorded at the geothermal field of Los Humeros, Mexico. Although only a one-dimensional velocity model is currently available, source locations and (time-dependent) moment tensors could be determined for selected events.


2021 ◽  
Author(s):  
Claudia Finger ◽  
Erik H. Saenger

<p>In addition to stable and accurate hypocenters of seismic events, the characterisation of events is crucial for the investigation of seismicity in the context of geothermal reservoirs, CO2-sequestration and other geotechnical applications. Since the origin and nature of the seismicity in such cases is still under investigation, tools should rely on as few a priori assumptions about the sources as possible. Here, an approach is presented to determine the time-dependent moment tensor and origin time in addition to commonly derived hypocenter locations of seismic events using time-reverse imaging (TRI). The full six component moment tensor is derived and may be used to display for example focal mechanisms. The workflow consists of determining the location of potential sources, discriminating artificial and true source locations and obtaining the time-dependent moment tensors by recording the stress components at the derived source locations. Since TRI does not rely on the identification of seismic phases but on the simulation of the time-reversed wavefield through an adequate velocity model, no assumptions about the source location or the type of source mechanism is made. TRI is less affected by low signal-to-noise ratios and is thus promising for noisier sites and quasi-simultaneous events. However, a sufficient number of seismic stations are needed to accurately sample the wavefield spatially. The proposed workflow is demonstrated by locating and characterising microseismic events in the geothermal field of Los Humeros, Mexico. Although higher levels of noise are present and only a one-dimensional velocity model is available at this time, selected events could be located and characterised.</p>


1962 ◽  
Vol 52 (3) ◽  
pp. 673-682
Author(s):  
J. Cleary ◽  
H. Doyle

Abstract A description is given of the application of a close network of seismograph stations and of an IBM 650 computer program to the study of a sequence of near earthquakes in southeastern Australia. The epicenters, depths, and origin times of a moderate sized earthquake and its aftershocks were determined with unusual accuracy, revealing an interesting pattern of foci. The position and origin time of the main shock were then used to study arrivals at more distant stations; Pn arrivals were found to fit a linear travel-time equation closely to a distance of 15 degrees, with a surface velocity of 8.16 ± 0.03 km/sec. The Sn phase was not clear and gave an approximate velocity of 4.7 ± 0.2 km/sec. The possibility of a small velocity decrease in the upper mantle below eastern Australia is discussed.


2006 ◽  
Vol 13 (4) ◽  
pp. 393-400 ◽  
Author(s):  
E. De Lauro ◽  
S. De Martino ◽  
M. Falanga ◽  
M. Palo

Abstract. We analyze time series of Strombolian volcanic tremor, focusing our attention on the frequency band [0.1–0.5] Hz (very long period (VLP) tremor). Although this frequency band is largely affected by noise, we evidence two significant components by using Independent Component Analysis with the frequencies, respectively, of ~0.2 and ~0.4 Hz. We show that these components display wavefield features similar to those of the high frequency Strombolian signals (>0.5 Hz). In fact, they are radially polarised and located within the crater area. This characterization is lost when an enhancement of energy appears. In this case, the presence of microseismic noise becomes relevant. Investigating the entire large data set available, we determine how microseismic noise influences the signals. We ascribe the microseismic noise source to Scirocco wind. Moreover, our analysis allows one to evidence that the Strombolian conduit vibrates like the asymmetric cavity associated with musical instruments generating self-sustained tones.


Sign in / Sign up

Export Citation Format

Share Document