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

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.

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

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

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

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

Joint microseismic moment tensor inversion and location using P- and S-waves diffraction stacking

The microseismic location methods based on diffraction stacking which does not require arrival picking can yield accurate and reliable source location for data with a low signal-to-noise ratio. However, due to the complex radiation pattern from a rupturing source, variation in the waveform polarities brings challenges to the diffraction-stacking based methods. The current implementations of joint source mechanism inversion and location methods which only use P-wave amplitudes have limitations in noise resistance and location accuracy. To mitigate those issues, we develop a new method for joint microseismic moment tensor inversion and event location using diffraction stacking with P- and S-waves amplitudes, both of which are used to invert for the moment tensor of a microseismic event, and then the inverted moment tensor is used to correct the waveform polarity changes before stacking. In addition, to expedite the large amount of calculations required for moment tensor inversion at each potential source position and origin time, we develop an optimized grid search scheme and implement the algorithm with GPUs. The proposed location method does not require manual picking of the first arrivals, and can automatically detect and locate microseismic events from continuous data. We first validated the method with two synthetic examples, and then applied it to a surface monitoring dataset for hydraulic fracturing at a shale gas well pad in the southern Sichuan Basin, China, where billions of cubic meters of shale gas are being produced annually. The locations of the microseismic events are nicely correlated with the fracturing stages and the determined source mechanisms are also consistent with the expected fracture growth. The proposed method is feasible for microseismic surface monitoring with dense nodal arrays and can provide important information for fracture growth and regional stress characterization.

How to Assess the Moment Tensor Inversion Resolution for Mining Induced Seismicity: A Case Study for the Rudna Mine, Poland

Underground exploitation of georesources can be highly correlated with induced seismic activity. In order to reduce the risk and improve the mining operations safety, the mining activity is monitored by a dedicated seismic network. Moment tensor inversion is a powerful method to investigate the rupture process of earthquakes in mines, providing information on the geometry of the earthquake source and the moment release. Different approaches have been proposed to estimate the source mechanisms, with some advantages and limitations. One of the simplest and most used methods rely on the fit of the polarity and amplitude of first P wave onsets. More advanced techniques fit the full waveforms and their spectra. Here, we test and compare moment tensor and focal mechanism estimations for both inversion techniques. In order to assess the inversion resolution, we built realistic synthetic data, accounting for real seismic noise conditions and network geometry for the Rudna copper mine, SW Poland. The Rudna mine pertains to the Legnica-Glógow Copper District, where thousands of mining induced earthquakes are detected yearly, representing a serious hazard for miners and mining infrastructures. We simulate a range of different processes and locations, considering pure double couple, deviatoric and full moment tensors with different magnitudes and located in different mining panels. Results show that the P-wave first onset inversion is very sensitive to the geometry of the seismic network, which is limited by the existing underground infrastructure. On the other hand, the quality of the moment tensor solutions for the full waveform inversion is mainly determined by the strength of mining tremor and the signal-to-noize ratio. We discuss the performance of both inversion techniques and provide recommendations toward a reliable moment tensor analysis in mines.