An updated and extended earthquake catalogue for Greece and adjacent areas since 1900

Abstract. A homogeneous earthquake catalogue for Greece and adjacent areas covering the period 1900–2009 is presented, to be used for reliable seismic hazard studies. The catalogues of Makropoulos and Burton (1981) and Makropoulos et al. (1989), covering the time span 1900–1985, were updated for the period 1986–2009 using instrumentally determined focal coordinates, except for the magnitude from the bulletin of the ISC. For Ms, which is the magnitude scale included in the previous versions, the same procedure applied since 1964 was adopted, using the ISC body wave magnitude (mb) and a regression equation. In the present update, Mw is also calculated for the entire period 1900–2009 using a formula derived from all available moment magnitudes and directly determined by the moment tensor inversion method. Thus, a magnitude homogeneous catalogue concerning both Ms and Mw scales is presented. The extended catalogue contains 7352 events, 70% more than the 4310 events of the previous published (1989) version. The completeness test revealed that the catalogue is complete for magnitudes above 4 for the last 34 yr and that no earthquake with magnitude 6 or greater has been omitted in the whole instrumental era (1900–2009).

Download Full-text

The waveform inversion of mainshock and aftershock data of the 2006 M6.3 Yogyakarta earthquake

Geoscience Letters ◽

10.1186/s40562-021-00176-w ◽

2021 ◽

Vol 8 (1) ◽

Author(s):

Hijrah Saputra ◽

Wahyudi Wahyudi ◽

Iman Suardi ◽

Ade Anggraini ◽

Wiwit Suryanto

Keyword(s):

Joint Inversion ◽

Body Wave ◽

Moment Tensor ◽

Near Field ◽

Source Parameters ◽

Waveform Inversion ◽

Moment Tensor Inversion ◽

Slip Distribution ◽

Aftershock Distribution ◽

Velocity Models

AbstractThis study comprehensively investigates the source mechanisms associated with the mainshock and aftershocks of the Mw = 6.3 Yogyakarta earthquake which occurred on May 27, 2006. The process involved using moment tensor inversion to determine the fault plane parameters and joint inversion which were further applied to understand the spatial and temporal slip distributions during the earthquake. Moreover, coseismal slip distribution was overlaid with the relocated aftershock distribution to determine the stress field variations around the tectonic area. Meanwhile, the moment tensor inversion made use of near-field data and its Green’s function was calculated using the extended reflectivity method while the joint inversion used near-field and teleseismic body wave data which were computed using the Kikuchi and Kanamori methods. These data were filtered through a trial-and-error method using a bandpass filter with frequency pairs and velocity models from several previous studies. Furthermore, the Akaike Bayesian Information Criterion (ABIC) method was applied to obtain more stable inversion results and different fault types were discovered. Strike–slip and dip-normal were recorded for the mainshock and similar types were recorded for the 8th aftershock while the 9th and 16th June were strike slips. However, the fault slip distribution from the joint inversion showed two asperities. The maximum slip was 0.78 m with the first asperity observed at 10 km south/north of the mainshock hypocenter. The source parameters discovered include total seismic moment M0 = 0.4311E + 19 (Nm) or Mw = 6.4 with a depth of 12 km and a duration of 28 s. The slip distribution overlaid with the aftershock distribution showed the tendency of the aftershock to occur around the asperities zone while a normal oblique focus mechanism was found using the joint inversion.

Download Full-text

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

Bulletin of the Seismological Society of America ◽

10.1785/bssa0730020419 ◽

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.

Download Full-text

A New Uniform Moment Tensor Catalog for Yunnan, China, from January 2000 through December 2014

Seismological Research Letters ◽

10.1785/0220190242 ◽

2020 ◽

Vol 91 (2A) ◽

pp. 891-900

Author(s):

Yan Xu ◽

Keith D. Koper ◽

Relu Burlacu ◽

Robert B. Herrmann ◽

Dan-Ning Li

Keyword(s):

Stress Field ◽

Seismic Hazard ◽

Moment Tensor ◽

Horizontal Stress ◽

Strike Slip ◽

Moment Tensors ◽

Yunnan Region ◽

Seismic Waveforms ◽

The Moment ◽

Maximum Horizontal Stress

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.

Download Full-text

Probabilistic Moment Tensor Inversion for Hydrocarbon-Induced Seismicity in the Groningen Gas Field, the Netherlands, Part 2: Application

Bulletin of the Seismological Society of America ◽

10.1785/0120200076 ◽

2020 ◽

Vol 110 (5) ◽

pp. 2112-2123 ◽

Cited By ~ 3

Author(s):

Bernard Dost ◽

Annemijn van Stiphout ◽

Daniela Kühn ◽

Marloes Kortekaas ◽

Elmer Ruigrok ◽

...

Keyword(s):

Induced Seismicity ◽

Moment Tensor ◽

Moment Tensor Inversion ◽

Gas Field ◽

Double Couple ◽

Recent Developments ◽

Inversion Procedure ◽

Seismic Sections ◽

The Moment ◽

Groningen Gas Field

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

Download Full-text

Uncertainty in FPSs from moment-tensor inversion

Geophysics ◽

10.1190/geo2011-0024.1 ◽

2011 ◽

Vol 76 (6) ◽

pp. WC65-WC75 ◽

Cited By ~ 18

Author(s):

Jing Du ◽

Norm R. Warpinski

Keyword(s):

Moment Tensor ◽

Failure Process ◽

Uncertainty Propagation ◽

Moment Tensor Inversion ◽

Hydraulic Fractures ◽

Mechanism Analysis ◽

Fault Plane Solutions ◽

Waveform Data ◽

Amplitude Data ◽

The Moment

Although microseismic monitoring of hydraulic fractures has primarily been concerned with the dimensions, complexity, and growth of fractures or fracture systems, there is an ever-increasing desire to extract more information about the hydraulic-fracturing and/or natural fractures from microseismic data. Source mechanism analysis, which is concerned with deducing details of the failure process from the microseismic waveform data, is, therefore, attracting more attention. However, most of the studies focus more on the moment-tensor inversion than on extracting fault-plane solutions (FPSs) from inverted moment tensors. The FPSs can be extracted from the inverted moment-tensor, but there remains a question regarding how errors associated with the inversion of the moment-tensor affect the accuracy of the FPSs. We examine the uncertainties of FPS, given the uncertainties of the amplitude data, by looking into the uncertainty propagation from amplitude data into the moment-tensor and then into the resultant FPS. The uncertainty propagation method will be demonstrated using two synthetic examples.

Download Full-text

Seismic Moment Tensor Solutions from GeoNet Data to Provide a Moment Magnitude Scale for New Zealand

10.26686/wgtn.16947268.v1 ◽

2021 ◽

Author(s):

◽

Elizabeth de Joux Robertson

Keyword(s):

New Zealand ◽

Seismic Moment ◽

Moment Tensor ◽

Velocity Model ◽

Moment Tensor Inversion ◽

Moment Magnitude ◽

Seismic Moment Tensor ◽

Waveform Data ◽

Moment Tensors ◽

The Moment

<p>The aim of this project is to enable accurate earthquake magnitudes (moment magnitude, MW) to be calculated routinely and in near real-time for New Zealand earthquakes. This would be done by inversion of waveform data to obtain seismic moment tensors. Seismic moment tensors also provide information on fault-type. I use a well-established seismic moment tensor inversion method, the Time-Domain [seismic] Moment Tensor Inversion algorithm (TDMT_INVC) and apply it to GeoNet broadband waveform data to generate moment tensor solutions for New Zealand earthquakes. Some modifications to this software were made. A velocity model can now be automatically used to calculate Green's functions without having a pseudolayer boundary at the source depth. Green's functions can be calculated for multiple depths in a single step, and data are detrended and a suitable data window is selected. The seismic moment tensor solution that has either the maximum variance reduction or the maximum double-couple component is automatically selected for each depth. Seismic moment tensors were calculated for 24 New Zealand earthquakes from 2000 to 2005. The Global CMT project has calculated CMT solutions for 22 of these, and the Global CMT project solutions are compared to the solutions obtained in this project to test the accuracy of the solutions obtained using the TDMT_INVC code. The moment magnitude values are close to the Global CMT values for all earthquakes. The focal mechanisms could only be determined for a few of the earthquakes studied. The value of the moment magnitude appears to be less sensitive to the velocity model and earthquake location (epicentre and depth) than the focal mechanism. Distinguishing legitimate seismic signal from background seismic noise is likely to be the biggest problem in routine inversions.</p>

Download Full-text

Complementing regional moment magnitudes to GCMT: a perspective from the rebuilt International Seismological Centre Bulletin

Earth System Science Data ◽

10.5194/essd-13-1957-2021 ◽

2021 ◽

Vol 13 (5) ◽

pp. 1957-1985

Author(s):

Domenico Di Giacomo ◽

James Harris ◽

Dmitry A. Storchak

Keyword(s):

Body Wave ◽

Moment Tensor ◽

Specific Area ◽

Centroid Moment Tensor ◽

Earthquake Catalogues ◽

Origin Time ◽

Earthquake Parameters ◽

Short Period ◽

The Moment ◽

International Seismological Centre

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

Download Full-text

Indikasi Sesar Naik di Plampang, Pulau Sumbawa Berdasarkan Analisis Gempa Bumi 13 Juni 2020

EKSPLORIUM ◽

10.17146/eksplorium.2021.42.2.6273 ◽

2021 ◽

Vol 42 (2) ◽

pp. 111

Author(s):

Priyobudi Priyobudi ◽

Mohamad Ramdhan

Keyword(s):

Active Fault ◽

Fault Plane ◽

Moment Tensor ◽

Moment Tensor Inversion ◽

Reverse Fault ◽

Coulomb Stress Changes ◽

Stress Changes ◽

The Moment ◽

Splay Fault ◽

Lower Angle

ABSTRAK. Keberadaan sesar aktif dengan pola sesar naik di daerah Plampang berhasil diungkap dari sebaran hiposenter terelokasi, hasil inversi momen tensor, dan pemodelan perubahan tegangan Coulomb. Studi ini juga berhasil mengungkap sumber gempa pada sesar aktif tersebut dengan kedalaman relatif dangkal yang bisa menjadi ancaman di Pulau Sumbawa jika magnitudo maksimumnya rilis di masa yang akan datang. Hasil relokasi hiposenter menunjukkan sebaran episenter berarah barat daya–timur laut. Hal ini didukung juga oleh hasil inversi momen tensor yang menunjukkan bidang sesar berarah barat daya–timur laut (N2240E) dengan dip cukup curam (490). Penampang seismisitas vertikal pada arah dip menunjukkan adanya pola sesar naik yang semakin landai seiring bertambahnya kedalaman. Bidang sesar yang landai menunjukkan struktur decollement pada kedalaman 10–15 km dan berangsur menjadi curam sebagai struktur splay fault pada kedalaman 0–10 km. Hal tersebut konsisten dengan hasil inversi momen tensor yang menunjukkan mekanisme pergerakan sesar naik terjadi pada kedalaman 7 km. Pemodelan perubahan tegangan Coulomb menunjukkan adanya penambahan stress di luar area bidang sesar sehingga memicu terjadinya aftershocks. Sebaran gempa susulan menunjukkan adanya bidang sesar hipotetik dengan panjang 19 km dan lebar 12 km. Sesar sebesar ini berpotensi membangkitkan gempa dengan kekuatan Mw 6,4. Gempa Sumbawa 13 Juni 2020 dengan magnitudo M 5,3 disebabkan oleh sebagian kecil aktivitas dari bidang sesar tersebut.ABSTRACT. The existence of an active fault with a reverse fault mechanism in the Plampang area is successfully delineated from the distribution of the relocated hypocenter, the moment tensor inversion, and the Coulomb stress changes. This study also reveals the source of the earthquake in the active fault with a relatively shallow depth which can be a threat on Sumbawa Island if the maximum magnitude is released in the future. Seismicity from hypocenter relocation shows the distribution of the epicenter with a southwest–northeast direction. It is also supported by the moment tensor inversion result which shows the fault plane trending southwest–northeast (N2240E) with a steep dip (490). The vertical section of seismicity in the dip direction shows that the slope of the plane has a lower angle with increasing depth. The lower angle of a fault plane shows a decollement structure at a depth of 10–15 km and gradually becomes steep as a splay fault structure at a depth of 0–10 km. It is consistent with the result of moment tensor inversion which shows the mechanism of a reverse fault that occurred at a depth of 7 km. The Coulomb stress changes show the stress increasing outside the fault plane area, which triggers aftershocks. The distribution of aftershocks shows a hypothetical fault plane of 19 km long and 12 km wide. A fault of this size has the potential to generate an earthquake with a magnitude maximum of Mw 6.4. The Sumbawa earthquake on June 13, 2020, having M 5.3 was caused by a small part of the activity from the fault.

Download Full-text

Penentuan Karakteristik Mekanisme Gempa Tahun 2018-2019 Di Nusa Tenggara Menggunakan Metode Inversi Momen Tensor

Kappa Journal ◽

10.29408/kpj.v5i1.3315 ◽

2021 ◽

Vol 5 (1) ◽

pp. 31-39

Author(s):

Sri Rizki Eka Putri ◽

◽

Hiden Hiden ◽

Suhayat Minardi ◽

◽

...

Keyword(s):

Moment Tensor ◽

Inversion Method ◽

Eurasian Plate ◽

Waveform Data ◽

Fault Parameters ◽

The Indian Ocean ◽

The Moment ◽

High Level ◽

Earthquake Characteristics ◽

Orientation Parameters

Nusa Tenggara is one of the areas with a high level of seismic activity in Indonesia because this area is located between the Indian Ocean plate, which moves northward and pushes the Eurasian plate. One method that is often used to determine an earthquake's epicentre is the Tensor Moment Inversion method. This study aims to determine the moment tensor magnitude of each earthquake event and determine earthquake characteristics based on the earthquake focus mechanism in Nusa Tenggara from 2018 to 2019. The earthquake with a magnitude of ≥ 5.7 SR and to find out the fault parameters, namely strike, dip, and rake using waveform data. One method that is often used to determine an earthquake's epicentre is the Tensor Moment Inversion method. The results showed that the fault planes formed were reverse faults and oblique faults. It has been calculated the moment tensor for each of the six components, namely Mxx, Myy, Mzz, Mxy, Myz and Mxz. From the results of the focal analysis of the 2018-2019 Nusa Tenggara earthquake mechanism, the values of the fault plane orientation parameters such as strike, dip and rake are obtained. For strikes in Nusa Tenggara on area 1, namely: 73° to 122°, Dip: 20° to 72° and Rake: 53° to 139°. While in field 2 for a strike, it is 232° to 280°, Dip 28° to 75°, Rake: 52° to 102°.

Download Full-text

Adaptive moment-tensor joint inversion of clustered microseismic events for monitoring geological carbon storage

Geophysical Journal International ◽

10.1093/gji/ggz293 ◽

2019 ◽

Vol 219 (1) ◽

pp. 80-93

Author(s):

Yu Chen ◽

Lianjie Huang

Keyword(s):

Carbon Storage ◽

Oil Recovery ◽

Joint Inversion ◽

Moment Tensor ◽

Crack Opening ◽

Cost Effective ◽

Microseismic Monitoring ◽

Moment Tensor Inversion ◽

Inversion Method ◽

Microseismic Events

SUMMARY Moment-tensor inversion of induced microseismic events can provide valuable information for tracking CO2 plumes at geological carbon storage sites, and study the physical mechanism of induced microseismicity. Accurate moment-tensor inversion requires a wide-azimuthal coverage of geophones. Cost-effective microseismic monitoring for geological carbon storage often uses only one geophone array within a borehole, leading to a large uncertainty in moment-tensor inversion. We develop a new adaptive moment-tensor joint inversion method to reduce the inversion uncertainty, when using limited but typical geophone receiver geometries. We first jointly invert a number of clustered microseismic events using a uniform focal mechanism to minimize the waveform misfit between observed and predicted P and S waveforms. We then invert the moment tensor for each event within a limited searching range around the joint inversion result. We apply our adaptive joint inversion method to microseismic data acquired using a single borehole geophone array at the CO2-Enhanced Oil Recovery field at Aneth, Utah. We demonstrate that our inversion method is capable of reducing the inversion uncertainty caused by the limited azimuthal coverage of geophones. Our inverted strikes of focal mechanisms of microseismic events are consistent with the event spatial distribution in subparallel pre-existing fractures or geological imperfections. The large values up to 40 per cent of the CLVD components might indicate crack opening induced by CO2/wastewater injection or rupture complexity.

Download Full-text