Seismic Source Processes of 25 Earthquakes (Mw>5) in the Gulf of California

2022 ◽  
Author(s):  
Eduardo Huesca-Pérez ◽  
Edahí Gutierrez-Reyes ◽  
Luis Quintanar
Keyword(s):  
Gulf Of California ◽  
Moment Tensor ◽  
Source Parameters ◽  
Seismic Source ◽  
Focal Mechanisms ◽  
Process Models ◽  
Moment Tensors ◽  
Time Space ◽  
Double Couple ◽  
Broadband Seismograms

ABSTRACT The Gulf of California (GoC) is a complex tectonic boundary that has been instrumented in the past several decades to record broadband seismograms. This volume of data has allowed us to study several source parameters systematically. Before, only a few source parameters of earthquakes greater than magnitude five had been studied in the GoC area. We re-examined the focal mechanisms of several earthquakes in the southern GoC that occurred over the last 20 yr using local–regional distance broadband seismograms. These focal mechanisms were then used as input data to retrieve the time–space history of the rupture for each earthquake. This work contributes to the study of 25 rupture-process models computed with the method proposed by Yagi et al. (1999). To investigate more about the nature of the seismicity in the GoC, we also calculated the non-double-couple component of moment tensors for 45 earthquakes. Previous studies (e.g., Ortega et al., 2013, 2016) have shown that non-double-couple components from moment tensors in this region are associated with complex faulting, suggesting that oblique faults or several parallel faults are interacting simultaneously. Our results show that, at least for moderate earthquakes (5 < M < 6), rupture processes in the GoC show a complex interaction between fault systems. It is revealed on the important contribution of non-double-couple component obtained in the full moment tensor analysis.

scholarly journals Source parameters estimation from broadband regional seismograms for earthquakes in the Aegean region and the Gorda plate

2018 ◽  
Vol 40 (3) ◽  
pp. 1032
Author(s):  
A. Agalos ◽  
P. Papadimitriou ◽  
N. Voulgaris ◽  
K. Makropoulos
Keyword(s):  
Moment Tensor ◽  
Source Parameters ◽  
Parameters Estimation ◽  
Aegean Area ◽  
Moment Tensors ◽  
Time Period ◽  
Broadband Seismograms ◽  
The Moment ◽  
Gorda Plate ◽  
Good Agreement

Seismic moment tensors are estimated for earthquakes offshore Northern California and Greece using inversion of regionally recorded broadband seismograms. This study includes inversion results for the strongest events that occurred inside the Gorda plate and near the Mendocino triple junction from 1991 to 2005. The regional results are in good agreement with obtained teleseismic results. We finally applied the moment tensor inversion methodology and validation mainly to moderate sized earthquakes, with magnitude greater than M~4.0, in the Aegean area. The focal mechanisms of HI earthquakes that occurred during the time period between June 2003 and March 2007 were estimated using this procedure.

scholarly journals Fully probabilistic seismic source inversion – Part 1: Efficient parameterisation

2013 ◽  
Vol 5 (2) ◽  
pp. 1125-1162 ◽  
Author(s):  
S. C. Stähler ◽  
K. Sigloch
Keyword(s):  
Moment Tensor ◽  
Source Parameters ◽  
Practical Importance ◽  
Principal Component ◽  
Seismic Source ◽  
Empirical Orthogonal Functions ◽  
Time Function ◽  
Source Inversion ◽  
Source Time Function ◽  
Earthquake Parameters

Abstract. Seismic source inversion is a non-linear problem in seismology where not just the earthquake parameters themselves, but also estimates of their uncertainties are of great practical importance. Probabilistic source inversion (Bayesian inference) is very adapted to this challenge, provided that the parameter space can be chosen small enough to make Bayesian sampling computationally feasible. We propose a framework for PRobabilistic Inference of Source Mechanisms (PRISM) that parameterises and samples earthquake depth, moment tensor, and source time function efficiently by using information from previous non-Bayesian inversions. The source time function is expressed as a weighted sum of a small number of empirical orthogonal functions, which were derived from a catalogue of >1000 STFs by a principal component analysis. We use a likelihood model based on the cross-correlation misfit between observed and predicted waveforms. The resulting ensemble of solutions provides full uncertainty and covariance information for the source parameters, and permits to propagate these source uncertainties into travel time estimates used for seismic tomography. The computational effort is such that routine, global estimation of earthquake mechanisms and source time functions from teleseismic broadband waveforms is feasible.

A teleseismic body-wave analysis of the May 1980 Mammoth Lakes, California, earthquakes

1983 ◽  
Vol 73 (2) ◽  
pp. 419-434
Author(s):  
Jeffery S. Barker ◽  
Charles A. Langston
Keyword(s):  
Body Wave ◽  
Moment Tensor ◽  
Normal Fault ◽  
Strike Slip ◽  
Body Waves ◽  
P Wave ◽  
Moment Tensors ◽  
Double Couple ◽  
The Moment ◽  
Mammoth Lakes

abstract Teleseismic P-wave first motions for the M ≧ 6 earthquakes near Mammoth Lakes, California, are inconsistent with the vertical strike-slip mechanisms determined from local and regional P-wave first motions. Combining these data sets allows three possible mechanisms: a north-striking, east-dipping strike-slip fault; a NE-striking oblique fault; and a NNW-striking normal fault. Inversion of long-period teleseismic P and SH waves for the events of 25 May 1980 (1633 UTC) and 27 May 1980 (1450 UTC) yields moment tensors with large non-double-couple components. The moment tensor for the first event may be decomposed into a major double couple with strike = 18°, dip = 61°, and rake = −15°, and a minor double couple with strike = 303°, dip = 43°, and rake = 224°. A similar decomposition for the last event yields strike = 25°, dip = 65°, rake = −6°, and strike = 312°, dip = 37°, and rake = 232°. Although the inversions were performed on only a few teleseismic body waves, the radiation patterns of the moment tensors are consistent with most of the P-wave first motion polarities at local, regional, and teleseismic distances. The stress axes inferred from the moment tensors are consistent with N65°E extension determined by geodetic measurements by Savage et al. (1981). Seismic moments computed from the moment tensors are 1.87 × 1025 dyne-cm for the 25 May 1980 (1633 UTC) event and 1.03 × 1025 dyne-cm for the 27 May 1980 (1450 UTC) event. The non-double-couple aspect of the moment tensors and the inability to obtain a convergent solution for the 25 May 1980 (1944 UTC) event may indicate that the assumptions of a point source and plane-layered structure implicit in the moment tensor inversion are not entirely valid for the Mammoth Lakes earthquakes.

scholarly journals Moving from 1-D to 3-D velocity model: automated waveform-based earthquake moment tensor inversion in the Los Angeles region

2019 ◽  
Vol 220 (1) ◽  
pp. 218-234 ◽  
Author(s):  
Xin Wang ◽  
Zhongwen Zhan
Keyword(s):  
Los Angeles ◽  
Moment Tensor ◽  
Velocity Model ◽  
Primary Control ◽  
Lateral Velocity ◽  
Real Time System ◽  
Moment Tensors ◽  
Velocity Models ◽  
Double Couple ◽  
The Moment

SUMMARY Earthquake focal mechanisms put primary control on the distribution of ground motion, and also bear on the stress state of the crust. Most routine focal mechanism catalogues still use 1-D velocity models in inversions, which may introduce large uncertainties in regions with strong lateral velocity heterogeneities. In this study, we develop an automated waveform-based inversion approach to determine the moment tensors of small-to-medium-sized earthquakes using 3-D velocity models. We apply our approach in the Los Angeles region to produce a new moment tensor catalogue with a completeness of ML ≥ 3.5. The inversions using the Southern California Earthquake Center Community Velocity Model (3D CVM-S4.26) significantly reduces the moment tensor uncertainties, mainly owing to the accuracy of the 3-D velocity model in predicting both the phases and the amplitudes of the observed seismograms. By comparing the full moment tensor solutions obtained using 1-D and 3-D velocity models, we show that the percentages of non-double-couple components decrease dramatically with the usage of 3-D velocity model, suggesting that large fractions of non-double-couple components from 1-D inversions are artifacts caused by unmodelled 3-D velocity structures. The new catalogue also features more accurate focal depths and moment magnitudes. Our highly accurate, efficient and automatic inversion approach can be expanded in other regions, and can be easily implemented in near real-time system.

A Student’s Guide to and Review of Moment Tensors

1989 ◽  
Vol 60 (2) ◽  
pp. 37-57 ◽  
Author(s):  
M. L. Jost ◽  
R. B. Herrmann
Keyword(s):  
Point Source ◽  
Moment Tensor ◽  
Tensor Decomposition ◽  
Seismic Source ◽  
Order Moment ◽  
Review Of Literature ◽  
Moment Tensors ◽  
Numerical Examples ◽  
Testing Algorithms ◽  
The Moment

Abstract A review of a moment tensor for describing a general seismic point source is presented to show a second order moment tensor can be related to simpler seismic source descriptions such as centers of expansion and double couples. A review of literature is followed by detailed algebraic expansions of the moment tensor into isotropic and deviatoric components. Specific numerical examples are provided in the appendices for use in testing algorithms for moment tensor decomposition.

Seismic source inversion using Hamiltonian Monte Carlo and a 3-D Earth model in the Japanese Islands

2020 ◽  
Author(s):  
Saulė Simutė ◽  
Lion Krischer ◽  
Christian Boehm ◽  
Martin Vallée ◽  
Andreas Fichtner
Keyword(s):  
Monte Carlo ◽  
Moment Tensor ◽  
Seismic Source ◽  
Monte Carlo Algorithm ◽  
Earth Model ◽  
Model Parameters ◽  
Hamiltonian Monte Carlo ◽  
Double Couple ◽  
Full Moment Tensor ◽  
Japanese Islands

<p>We present a proof-of-concept catalogue of full-waveform seismic source solutions for the Japanese Islands area. Our method is based on the Bayesian inference of source parameters and a tomographically derived heterogeneous Earth model, used to compute Green’s strain tensors. We infer the full moment tensor, location and centroid time of the seismic events in the study area.</p><p>To compute spatial derivatives of Green’s functions, we use a previously derived regional Earth model (Simutė et al., 2016). The model is radially anisotropic, visco-elastic, and fully heterogeneous. It was constructed using full waveforms in the period band of 15–80 s.</p><p>Green’s strains are computed numerically with the spectral-element solver SES3D (Gokhberg & Fichtner, 2016). We exploit reciprocity, and by treating seismic stations as virtual sources we compute and store the wavefield across the domain. This gives us a strain database for all potential source-receiver pairs. We store the wavefield for more than 50 F-net broadband stations (www.fnet.bosai.go.jp). By assuming an impulse response as the source time function, the displacements are then promptly obtained by linear combination of the pre-computed strains scaled by the moment tensor elements.</p><p>With a feasible number of model parameters and the fast forward problem we infer the unknowns in a Bayesian framework. The fully probabilistic approach allows us to obtain uncertainty information as well as inter-parameter trade-offs. The sampling is performed with a variant of the Hamiltonian Monte Carlo algorithm, which we developed previously (Fichtner and Simutė, 2017). We apply an L2 misfit on waveform data, and we work in the period band of 15–80 s.</p><p>We jointly infer three location parameters, timing and moment tensor components. We present two sets of source solutions: 1) full moment tensor solutions, where the trace is free to vary away from zero, and 2) moment tensor solutions with the isotropic part constrained to be zero. In particular, we study events with significant non-double-couple component. Preliminary results of ~Mw 5 shallow to intermediate depth events indicate that proper incorporation of 3-D Earth structure results in solutions becoming more double-couple like. We also find that improving the Global CMT solutions in terms of waveform fit requires a very good 3-D Earth model and is not trivial.</p>

The Meckering earthquake of 14 October 1968: A possible downward propagating rupture

1987 ◽  
Vol 77 (5) ◽  
pp. 1558-1578
Author(s):  
Kristín S. Vogfjörd ◽  
Charles A. Langston
Keyword(s):  
Moment Tensor ◽  
Source Parameters ◽  
Time History ◽  
Vertical Displacement ◽  
Body Waves ◽  
Finite Fault ◽  
Short Period ◽  
Double Couple ◽  
The Moment ◽  
Fault Length

Abstract Average source parameters of the 1968 Meckering, Australia earthquake are obtained by the inversion of body waves. The objectives of the inversion are the elements of the moment tensor and the source-time history. An optimum source depth of 3 km is determined, but because of source complexity the point source assumption fails and the moment tensor obtained at that depth has a large nondouble-couple term, compensated linear vector dipole = 34 per cent. The source parameters of the major double-couple are: strike = 341°; dip = 37°; rake = 61°; and seismic moment = 8.2 ×1025 dyne-cm. The source-time function is of approximately 4 sec duration, with a long rise time and a sharp fall-off. The fault length is constrained on the surface by the observed surface break, and results from vertical displacement modeling suggest a width of approximately 10 km in the middle, assuming a dip of 37°. That restricts the entire faulted area to lie above 6 km depth. Two finite fault models for the earthquake are presented, with rupture initiating at a point (1) near the top of the fault and (2) at the bottom of the fault. Both models produce similar long-period synthetics, but based on the short-period waveforms, model 1 is favored. It is argued that such a rupture process is the most reasonable in this cold shield region.

Bayesian multiple rupture plane inversion to assess rupture complexity: application to the 2018 Mw 7.9 Alaska earthquake

2021 ◽  
Author(s):  
Angela Carrillo Ponce ◽  
Torsten Dahm ◽  
Simone Cesca ◽  
Frederik Tilmann ◽  
Andrey Babeyko ◽  
...  
Keyword(s):  
Moment Tensor ◽  
Real Data ◽  
Strike Slip ◽  
Focal Mechanisms ◽  
Strike Slip Fault ◽  
Body Waves ◽  
Source Inversion ◽  
Lateral Strike Slip ◽  
Double Couple ◽  
Alaska Earthquake

<p>When the earthquake rupture is complex and ruptures of multiple fault segments contribute to the total energy release, the produced wavefield is the superposition of individual signals produced by single subevents. Resolving source complexity for large, shallow earthquakes can be used to improve ground shaking and surface slip estimations, and thus tsunami models. The 2018 Mw 7.9 Alaska earthquake showed such complexity: according to previous studies, the rupture initiated as a right-lateral strike-slip fault on a N-S oriented fault plane, but then jumped onto a left-lateral strike-slip fault oriented westward. Rupture complexity and presence of multiple subevents may characterize a number of other earthquakes. However, even when individual subevents are spatially and/or temporally separated, it is very difficult to identify them from far field recordings. In order to model complex earthquakes we have implemented a multiple double couple inversion scheme within Grond, a tool devoted to the robust characterization of earthquake source parameters included in the Pyrocko software. Given the large magnitude of the target earthquake, we perform our source inversions using broadband body waves data (P and S phases) at teleseismic distances. Our approach starts with a standard moment tensor inversion, which allows to get more insights about the centroid location and overall moment release. These values can then be used to constrain the double source inversion. We discuss the performance of the inversion for the Alaska earthquake, using synthetic and real data. First, we generated realistic synthetic waveforms for a two-subevents source, assuming double couple sources with the strike-slip mechanisms proposed for the Alaska earthquake. We model the synthetic dataset both using a moment tensor and a double double couple source, and demonstrate the stability of the double double couple inversion, which is able to reconstruct the two focal mechanisms, the moment ratio and the relative centroid locations of the two subevents. Synthetic tests show that the inversion accuracy can be in some cases reduced, in presence of noisy data and when the interevent time between subevents is short. A larger noise addition affects the retrieval of the focal mechanism orientations only in some cases, but in general all the parameters were well retrieved. Then, we test our tool using real data for the earthquake. The single source inversion shows that the centroid is shifted 27 s in time and 40 km towards NE with respect to the original assumed location retrieved from the gCMT catalogue. The following double double couple source inversion resolves two subevents with right-lateral and left-lateral strike-slip focal mechanisms and Mw 7.6 and 7.8 respectively. The subevent centroids are separated by less than 40 km in space and less than 20 s in time.</p>

Focal mechanisms for small to intermediate earthquakes in the northern part of the Alps and their seismotectonic interpretation

2020 ◽  
Author(s):  
Thomas Plenefisch ◽  
Laura Barth ◽  
Keyword(s):  
Stress Field ◽  
Focal Mechanisms ◽  
Small Magnitude ◽  
Moment Tensors ◽  
Time Period ◽  
Double Couple ◽  
Geodynamic Processes ◽  
The Alps ◽  
Seismic Networks ◽  
Seismotectonic Interpretation

<p>In the framework of the AlpArray project more than 600 broadband stations have been installed and operated in the Alps and the surroundings. Together with the permanent stations in the area it is one of the most densely spaced seismic networks worldwide. Thereby, it offers an excellent opportunity to investigate the seismicity and seismotectonics of the Alpine chain. Due to the huge number of stations focal mechanisms can be calculated even for small magnitude earthquakes with high accuracy. The focal mechanisms are one important key to reveal the contemporary stress field and thus contribute to a better understanding of the geodynamic processes of the Alps.</p><p>In our study we focus on small to intermediate earthquakes in the Northern Alps, namely on four distinct sub-regions. These are from West to East the Lake Constance, the Arlberg region, the area of Garmisch-Partenkirchen and the broader region of Innsbruck. In order to calculate the focal mechanisms, we apply the FOCMEC program (Snoke, 2003), which inverts for a pure double-couple source. P-polarities as well as amplitude ratios of SH to P are used as input parameters for the inversion. Thanks to the dense network a good coverage of the focal sphere is achieved in most cases.</p><p>Altogether, we calculated focal mechanisms for 25 earthquakes in the magnitude range between 2.5 and 3.5 from the time period 2016 to 2019. Most of the focal mechanisms represent reverse or strike-slip faulting, normal faulting events are rather rare. The mechanisms are analysed with respect to lateral changes along the Northern Alpine. On one hand we compare the mechanisms with mechanisms of older studies as well as with moment tensors of events of slightly larger magnitudes. Those events are the scope of another subproject in the framework of the AlpArray (Petersen et al., 2019). On the other hand, we compare our mechanisms with geological indicators, namely orientation of faults. Finally, the focal mechanisms are used as input to invert for the stress field.</p>

A user-friendly probabilistic earthquake source inversion framework for joint inversion of seismic, geodetic, and gravitational signals - The Grond toolkit

2020 ◽  
Author(s):  
Sebastian Heimann ◽  
Marius Isken ◽  
Daniela Kühn ◽  
Hannes Vasyura-Bathke ◽  
Henriette Sudhaus ◽  
...  
Keyword(s):  
Open Source ◽  
Joint Inversion ◽  
Moment Tensor ◽  
Source Parameters ◽  
Seismic Source ◽  
Magma Ascent ◽  
Caldera Collapse ◽  
Waveform Data ◽  
Trade Offs ◽  
Finite Fault

<p>Seismic source and moment tensor waveform inversion is often ill-posed or non-unique if station coverage is poor or signals are weak. Three key ingredients can help in these situations: (1) probabilistic inference and global search of the full model space, (2) joint optimisation with datasets yielding complementary information, and (3) robust source parameterisation or additional source constraints. These demands lead to vast technical challenges, on the performance of forward modelling, on the optimisation algorithms, as well as on visualisation, optimisation configuration, and management of the datasets. Implementing a high amount of automation is inevitable.</p><p>To tackle all these challenges, we are developing a sophisticated new seismic source optimisation framework, Grond. With its innovative Bayesian bootstrap optimiser, it is able to efficiently explore large model spaces, the trade-offs and the uncertainties of source parameters. The program is highly flexible with respect to the adaption to specific source problems, the design of objective functions, and the diversity of empirical datasets.</p><p>It uses an integrated, robust waveform data processing, and allows for interactive visual inspection of many aspects of the optimisation problem, including visualisation of the result uncertainties. Grond has been applied to CMT moment tensor and finite-fault optimisations at all scales, to nuclear explosions, to a meteorite atmospheric explosion, and to volcano-tectonic processes during caldera collapse and magma ascent. Hundreds of seismic events can be handled in parallel given a single optimisation setup.</p><p>Grond can be used to optimise simultaneously seismic waveforms, amplitude spectra, waveform features, phase picks, static displacements from InSAR and GNSS, and gravitational signals.</p><p>Grond is developed as an open-source package and community effort. It builds on and integrates with other established open-source packages, like Kite (for InSAR) and Pyrocko (for seismology).</p>

Sign in / Sign up

Export Citation Format

Share Document