Determining the Double Couple Moment Tensor for Kamchatka Earthquakes from Regional Seismic Waveforms

2021 ◽  
Vol 57 (3) ◽  
pp. 332-347
Author(s):  
I. R. Abubakirov ◽  
V. M. Pavlov
1983 ◽  
Vol 73 (2) ◽  
pp. 419-434
Author(s):  
Jeffery S. Barker ◽  
Charles A. Langston

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.


Author(s):  
D. Chebrov ◽  
A. Chebrova ◽  
I. Abubakirov ◽  
E. Matveenko ◽  
S. Mityushkina ◽  
...  

The seismicity review of Kamchatka and surrounding territories for 2014 is given. In Kamchatka earthquake catalogue minimum local magnitude of completeness is MLmin=3.5, and for earthquakes under the Okhotsk sea with h≥350 kmMLmin=3.6. The Kamchatka earthquake catalogue for 2014 with ML3.5, published in the Appendix to this annual, includes 1114 events. 86 earthquakes of the catalogue with ML=3.35–6.2 were felt in Kamchatka and surrounding areas with seismic intensity I ranged from 2 to 5 according the MSK-64 scale. For all events with ML5.0 occurred in the area of responsibility of the KB GS RAS in 2014, an attempt to calculate the seismic moment tensor (SMT) was made. There are 40 such events in the regional catalogue. For 36 earthquakes, the SMT and depth h of the equivalent point source were calculated successfully. The calcu-lations were performed for the SMT double-couple model using a nonlinear algorithm. In 2014, a typical location of the earthquake epicenters was observed in the Kamchatka zone. In 2014, the seismicity level in all selected zones and in the region as a whole corresponded to the background one according to the “SESL’09” scale. The number of recorded events with ML3.6 and strong earthquakes with ML5.1 is close to the average annual value. Anomalous and outstanding events were not recorded.


2019 ◽  
Vol 220 (1) ◽  
pp. 218-234 ◽  
Author(s):  
Xin Wang ◽  
Zhongwen Zhan

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.


2020 ◽  
Vol 91 (2A) ◽  
pp. 891-900
Author(s):  
Yan Xu ◽  
Keith D. Koper ◽  
Relu Burlacu ◽  
Robert B. Herrmann ◽  
Dan-Ning Li

Abstract Because of the collision of the Indian and Eurasian tectonic plates, the Yunnan Province of southwestern China has some of the highest levels of seismic hazard in the world. In such a region, a catalog of moment tensors is important for estimating seismic hazard and helping understand the regional seismotectonics. Here, we present a new uniform catalog of moment tensor solutions for the Yunnan region. Using a grid-search technique to invert seismic waveforms recorded by the permanent regional network in Yunnan and the 2 yr ChinArray deployment, we present 1833 moment tensor solutions for small-to-moderate earthquakes that occurred between January 2000 and December 2014. Moment magnitudes in the new catalog vary from Mw 2.2 to 6.1, and the catalog is complete above Mw∼3.5–3.6. The moment tensors are constrained to be purely double-couple and show a variety of faulting mechanisms. Normal faulting events are mainly concentrated in northwest Yunnan, while farther south along the Sagaing fault the earthquakes are mostly thrust and strike slip. The remaining area includes all three styles of faulting but mostly strike slip. We invert the moment tensors for the regional stress field and find a strong correlation between spatially varying maximum horizontal stress and Global Positioning System observations of horizontal ground velocity. The stress field reveals clockwise rotation around the eastern Himalayan syntaxis, with northwest–southeast compression to the east of the Red River fault changing to northeast–southwest compression west of the fault. Almost 88% of the centroid depths are shallower than 16 km, consistent with a weak and ductile lower crust.


2016 ◽  
Vol 207 (2) ◽  
pp. 1137-1141 ◽  
Author(s):  
Zoë K. Mildon ◽  
David J. Pugh ◽  
Jon Tarasewicz ◽  
Robert S. White ◽  
Bryndís Brandsdóttir

Abstract Moment tensor analysis with a Bayesian approach was used to analyse a non-double-couple (non-DC) earthquake (Mw ∼ 1) with a high isotropic (implosive) component within the Krafla caldera, Iceland. We deduce that the earthquake was generated by a closing crack at depth. The event is well located, with high signal-to-noise ratio and shows dilatational P-wave first arrivals at all stations where the first arrival can be picked with confidence. Coverage of the focal sphere is comprehensive and the source mechanism stable across the full range of uncertainties. The non-DC event lies within a cluster of microseismic activity including many DC events. Hence, we conclude that it is a true non-DC closing crack earthquake as a result of geothermal utilization and observed magma chamber deflation in the region at present.


2020 ◽  
Vol 110 (5) ◽  
pp. 2112-2123 ◽  
Author(s):  
Bernard Dost ◽  
Annemijn van Stiphout ◽  
Daniela Kühn ◽  
Marloes Kortekaas ◽  
Elmer Ruigrok ◽  
...  

ABSTRACT Recent developments in the densification of the seismic network covering the Groningen gas field allow a more detailed study of the connection between induced seismicity and reactivated faults around the gas reservoir at 3 km depth. With the reduction of the average station distance from 20 km to 4–5 km, a probabilistic full-waveform moment tensor inversion procedure could be applied, resulting in both improved hypocenter location accuracy and full moment tensor solutions for events of M≥2.0 recorded in the period 2016–2019. Hypocenter locations as output from the moment tensor inversion are compared to locations from the application of other methods and are found similar within 250 m distance. Moment tensor results show that the double-couple (DC) solutions are in accordance with the known structure, namely normal faulting along 50°–70° dipping faults. Comparison with reprocessed 3D seismic sections, extended to a depth of 6–7 km, demonstrate that (a) most events occur along faults with a small throw and (b) reactivated faults in the reservoir often continue downward in the Carboniferous underburden. From non-DC contributions, the isotropic (ISO) component is dominant and shows consistent negative values, which is expected in a compacting medium. There is some indication that events connected to faults with a large throw (>70  m) exhibit the largest ISO component (40%–50%).


2020 ◽  
Author(s):  
Saulė Simutė ◽  
Lion Krischer ◽  
Christian Boehm ◽  
Martin Vallée ◽  
Andreas Fichtner

<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>


Author(s):  
D. Chebrov ◽  
V. Saltikov ◽  
E. Matveenko ◽  
S. Droznina ◽  
E. Romasheva ◽  
...  

The seismicity review of Kamchatka and surrounding territories for 2015 is given. In the Kamchatka earthquake catalogue, the minimum local magnitude of completeness is MLmin=3.5, and for earthquakes with h≥350 km under the Okhotsk sea MLmin=3.6. The Kamchatka earthquake catalogue for 2015 with ML3.5, published in the Appendix to this issue, includes 1213 events. 92 earthquakes of the catalogue with ML=3.0–6.5 were felt in Kamchatka and surrounding areas with seismic intensity I=2–6 according to the MSK-64 scale. For all events with ML5.0 that occurred in 2015 in the KB GS RAS area of responsibility, an attempt to calculate the seismic moment tensor (SMT) was made. There are 32 such events in the regional catalogue. For 28 earthquakes the SMT and depth h of the equivalent point source were calculated successfully. The calculations were performed for the SMT double-couple model using a nonlinear algorithm. In 2015, a typical location of the earthquake epicenters was observed in the Kamchatka zone. In 2015, the seismicity level in all selected zones and in the region as a whole correspond to the background one according to the “SESL’09” scale. The number of recorded events with ML3.5 and strong earthquakes with ML5.0 is close to the average annual value. Anomalous and significant events were not recorded.


2020 ◽  
Vol 110 (3) ◽  
pp. 1125-1133
Author(s):  
William Menke ◽  
Joshua B. Russell

ABSTRACT The non-double-couple (non-DC) components of the moment tensor provide insight into the earthquake processes and anisotropy of the near-source region. We investigate the behavior of the isotropic (ISO) and compensated linear vector dipole (CLVD) components of the moment tensor for shear faulting in a transversely ISO medium with an arbitrarily oriented symmetry axis. Analytic formulas for ISO and CLVD depend on the orientation of the fault relative to the anisotropy symmetry axis as well as three anisotropic parameters, which describe deviations of the medium from isotropy. Numerical experiments are presented for the preliminary reference Earth model. Both ISO and CLVD components are zero when the axis of symmetry is within the fault plane or the auxiliary plane. For any orientation in which the ISO component is zero, the CLVD component is also zero, but the opposite is not always true (e.g., for strong anisotropy). The relative signs of the non-DC components of neighboring earthquakes may help distinguish source processes from source-region anisotropy. We prove that an inversion of ISO and CLVD components of a set of earthquakes with different focal mechanisms can uniquely determine the orientation and strength of anisotropy. This study highlights the importance of the ISO component for constraining deep slab anisotropy and demonstrates that it cannot be neglected.


Sign in / Sign up

Export Citation Format

Share Document