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

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.

1981 ◽  
Vol 71 (5) ◽  
pp. 1423-1444
Author(s):  
Jeffrey S. Barker ◽  
Charles A. Langston

abstract Seismograms from WWSSN and Canadian network stations were modeled to determine the source parameters of the 20 June 1978 Thessaloniki, Greece, earthquake (Ms = 6.4). The depth of the initial rupture was constrained to 11 ± 1 km by comparison of the arrival times of surface reflections with synthetic short-period seismograms. A focal sphere plot of first motion polarities provided little constraint on other focal parameters, except to indicate that predominantly normal faulting was involved. A generalized inverse technique utilizing the moment tensor formalism was applied to teleseismic P and SH waves for six increments of depth. The moment tensor obtained indicated a nearly horizontal, N-trending tension axis and a nearly vertical compression axis, and yielded the following double-couple source parameters: strike 280° ± 7°; dip 55° ± 3°; rake −65° ± 5°; seismic moment 5.7 × 1025 dyne-cm; and a skewed triangular source time function with a rise time of about 1 sec and duration of 6 to 8 sec. Due to indications of multiple or finite source effects for this event, and the assumption in the moment tensor formalism of a point source, a low-pass filter was applied to the data and the inversions were repeated. The results were nearly identical with those of the original inversion, suggesting that any individual sources had similar mechanisms, or that the point source model is sufficient for this earthquake.


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.


By recording several components of tilt, strain and acceleration at one location, one can determine the focal mechanism, or moment tensor, of an earthquake. Alternatively, recordings made at several locations can be used. The moment tensor can be decomposed into its isotropic part and its deviatoric part. When the eigerrvalues of the deviator are in the sequence (— 1, 0, 1) the equivalent double couple can be found.


2020 ◽  
Vol 92 (1) ◽  
pp. 212-226
Author(s):  
Sam Wimpenny ◽  
C. Scott Watson

Abstract We present a compilation of 2131 high-fidelity mechanisms and centroid depths of moderate-magnitude earthquakes derived using synthetic body-waveform modeling (the Global Waveform-Modelled Earthquake Catalog v1.0—gWFM), which can be visualized and downloaded online (see Data and Resources). In this article, we describe the methods used to construct the gWFM and present a comparison between the earthquake depths and focal mechanisms in the gWFM with those derived by the International Seismological Centre, Global Centroid Moment Tensor (Global CMT) project, and the U.S. Geological Survey (USGS) W-phase, as well as 60 events studied using geodesy. We find that 20%–30% of the earthquakes in routine global catalogs have depths that differ by more than 10 km from those in the gWFM. Shallow-crustal earthquakes of Mw 5–6 are typically the worst located in depth by routine catalogs. Over 90% of the earthquakes in the gWFM are within ±15° in strike, ±5° in dip, and ±15° in rake of the Global CMT and USGS W-phase best double-couple moment tensor solutions. However, the mechanisms of shallow Mw 5–6 earthquakes in the routine catalogs can be inaccurate, due to the well-known insensitivity of long-period surface waves to the vertical dip-slip components of the moment tensor. The gWFM is an archive of well-constrained earthquake source parameters, though it will continue to update as new earthquake mechanisms and depths are published, thereby remaining an up-to-date research tool.


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.


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 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%).


2021 ◽  
Vol 13 (5) ◽  
pp. 1957-1985
Author(s):  
Domenico Di Giacomo ◽  
James Harris ◽  
Dmitry A. Storchak

Abstract. Seismologists and geoscientists often need earthquake catalogues for various types of research. This input usually contains basic earthquake parameters such as location (longitude, latitude, depth, and origin time), as well as magnitude information. For the latter, the moment magnitude Mw has become the most sought after magnitude scale in the seismological community to characterize the size of an earthquake. In this contribution we provide an informative account of the Mw content for the newly rebuilt Bulletin of the International Seismological Centre (ISC, http://www.isc.ac.uk, last access: May 2021), which is regarded as the most comprehensive record of the Earth's seismicity. From this data, we extracted a list of hypocentres with Mw from a multitude of agencies reporting data to the ISC. We first summarize the main temporal and spatial features of the Mw provided by global (i.e. providing results for moderate to great earthquakes worldwide) and regional agencies (i.e. also providing results for small earthquakes in a specific area). Following this, we discuss their comparisons, by considering not only Mw but also the surface wave magnitude MS and short-period body wave magnitude mb. By using the Global Centroid Moment Tensor solutions as an authoritative global agency, we identify regional agencies that best complement it and show examples of frequency–magnitude distributions in different areas obtained both from the Global Centroid Moment Tensor alone and complemented by Mw from regional agencies. The work done by the regional agencies in terms of Mw is fundamental to improve our understanding of the seismicity of an area, and we call for the implementation of procedures to compute Mw in a systematic way in areas currently not well covered in this respect, such as vast parts of continental Asia and Africa. In addition, more studies are needed to clarify the causes of the apparent overestimation of global Mw estimations compared to regional Mw. Such difference is also observed in the comparisons of Mw with MS and mb. The results presented here are obtained from the dataset (Di Giacomo and Harris, 2020, https://doi.org/10.31905/J2W2M64S) stored at the ISC Dataset Repository (http://www.isc.ac.uk/dataset_repository/, last access: May 2021).


2019 ◽  
Vol 219 (2) ◽  
pp. 1148-1162
Author(s):  
Jiun-Ting Lin ◽  
Wu-Lung Chang ◽  
Diego Melgar ◽  
Amanda Thomas ◽  
Chi-Yu Chiu

SUMMARY We test the feasibility of GPS-based rapid centroid moment tensor (GPS CMT) methods for Taiwan, one of the most earthquake prone areas in the world. In recent years, Taiwan has become a leading developer of seismometer-based earthquake early warning systems, which have successfully been applied to several large events. The rapid determination of earthquake magnitude and focal mechanism, important for a number of rapid response applications, including tsunami warning, is still challenging because of the limitations of near-field inertial recordings. This instrumental issue can be solved by an entirely different observation system: a GPS network. Taiwan is well posed to take advantage of GPS because in the last decade it has developed a very dense network. Thus, in this research, we explore the suitability of the GPS CMT inversion for Taiwan. We retrospectively investigate six moderate to large (Mw6.0 ∼ 7.0) earthquakes and propose a resolution test for our model, we find that the minimum resolvable earthquake magnitude of this system is ∼Mw5.5 (at 5 km depth). Our tests also suggest that the finite fault complexity, often challenging for the near-field methodology, can be ignored under such good station coverage and thus, can provide a fast and robust solution for large earthquake directly from the near field. Our findings help to understand and quantify how the proposed methodology could be implemented in real time and what its contributions could be to the overall earthquake monitoring system.


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