1 / 27

RAPID SOURCE PARAMETER DETERMINATION AND EARTHQUAKE SOURCE PROCESS IN INDONESIA REGION

RAPID SOURCE PARAMETER DETERMINATION AND EARTHQUAKE SOURCE PROCESS IN INDONESIA REGION. Iman Suardi Seismology Course Indonesia. Supervisor: DR. Akio KATSUMATA (MRI-JMA) Prof. Yuji YAGI (Tsukuba University). Final Presentation of Master Report August 28, 2007.

tricia
Download Presentation

RAPID SOURCE PARAMETER DETERMINATION AND EARTHQUAKE SOURCE PROCESS IN INDONESIA REGION

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. RAPID SOURCE PARAMETER DETERMINATION AND EARTHQUAKE SOURCE PROCESS IN INDONESIA REGION Iman Suardi Seismology Course Indonesia Supervisor: DR. Akio KATSUMATA (MRI-JMA) Prof. Yuji YAGI (Tsukuba University) Final Presentation of Master Report August 28, 2007 International Institute of Seismology and Earthquake Engineering, Building Research Institute, Tsukuba, Japan, 2006-2007

  2. CONTENTS • INTRODUCTION • DATA • THEORY AND METHODOLOGY • RESULTS AND DISCUSSIONS • CONCLUSION • RECOMMENDATION

  3. INTRODUCTION Indonesia Earthquake Information and Tsunami Warning Center (InaTWS end 2008)

  4. Tectonic Setting of Indonesia

  5. Purpose of This Study • Propose a new rapid method to determine earthquake magnitude • release a new BMG magnitude tsunami early warning, • Determine focal mechanism and seismic moment, • Determine moment magnitude, • Compare moment magnitude with the above magnitude.

  6. DATA Distribution Of The New Accelerometer Network Of BMG data range from 4 March 2007 – 25 April 2007

  7. Data Preparation Calculated Magnitude of BMG Continuous waveform data from 4 March – 25 April 2007 Raw Data from Accelerometer TSA100S from Nanometrics Inc. (continuous waveform data). Convert to XFiles format Convert to Seed format Convert to SAC (by using rdseed4.6) using software Data Playback prepared by Nanometrics Inc. Moment Tensor Inversion Event data collection from Global CMT Retrieve waveform data from internet (IRIS network / FDSN-GSN stations) Remove instrument responses

  8. before The performance of waveform before converting and after converting into displacement record for event5 (2007/03/17, M6.2, depth 40.7 km, 1.30N; 126.37E). after

  9. Earthquake event list from the Global CMT (4 March 2007 – 25 April 2007). Data Retrieval from IRIS-DMC IRIS-DMC stations

  10. THEORY AND METHODOLOGY Magnitudes Scales • Body Wave Magnitude, mb • mb = Log10(A/T) + q(∆,h), where: A(µm), Amplitude (SP); T(s), Period; • ∆ (deg), Epicentral Distance; h(km), Focal Depth. • Based on amplitudes of P waves (T~1 sec)  not appropriate for large earthquake. • Surface Wave Magnitude, Ms, Gutenberg (1945) • Ms = Log10(A/T) + 1.66Log10∆ + 3.3, where: A(µm), Amplitude (LP); T(s), Period; ∆ (deg), Epicentral Distance. • Based on amplitudes of surface waves (T~20 sec)  easy to determine, underestimate size for great earthquakes. • Moment Magnitude, Mw, Kanamori (1977) • Mw = (Log10Mo – 9.1)/1.5, Mo (Nm) = µDS : Seismic Moment • where: µ rigidity; D Displacement; S Fault Area. • Based on moment tensor solutions such as USGS MT solutions and Global CMT • solution  No saturation, Difficult to determine. • Broadband Moment Magnitude, Tsuboi et al. (1995). • Mwp = Mw + 0.2

  11. Continued … Magnitude JMA, MJMA using accelerometer • For distance within 2000 km • , Tsuboi’s formula (1954). • Shallow earthquake < 60 km. • , Katsumata’s formula (2004). • R > 30 km, ∆ < 700 km • , Funasaki et al. (2004) • R > 5 km, ∆ < 700 km • where: AD is the maximum amplitude (µm), AZ isvertical component of velocity , , Attenuation Function depends on epicentral distance ∆ and on the focal depth, H, CVdepends on instrumental conditions of seismographs.

  12. Empirical Formula Find coefficients Least Square Method A New Magnitude for BMG Empirical Formula of Richter’s Magnitude

  13. unknown unknown station event minimum constant unknown constant minimum Least Square Method find coefficient a1, a2, and a3

  14. Green’s function Observed station Propagate wave Unit impulsive slip Moment Tensor Inversion • Determination of Source Parameter: • faulting type  Focal mechanism, • earthquake size  Seismic Moment, Moment Magnitude. Green’s Function The response function, effect of propagation, with unit impulsive slip and / or force.

  15. RESULTS AND DISCUSSIONS The Magnitude Determination of BMG Empirical Equation Type1: Type2:

  16. Type2: A Calculated Magnitude of BMG or Tsuboi’s formula

  17. Graph of Comparison of maximum amplitudes of three type of magnitude calculation plotted against epicentral distances for each events.

  18. Graph of comparison between the moment magnitude from Global CMT and the calculated magnitude of BMG (type2).

  19. Difference between the calculated magnitude of BMG and the moment magnitude of Global CMT plotted against the moment magnitude of Global CMT. Difference between the calculated magnitude of BMG and the moment magnitude from Global CMT plotted against the depth of earthquakes.

  20. Graph of difference between observed maximum amplitudes and estimated ones plotted against hypocentral distances for all the events.

  21. (a) (b) Moment Tensor Inversion The results of moment tensor inversion of 5 earthquake events. (a) Assumed source time function and focal mechanism. (b) Theoretical waveforms (red curves) from a point source and observed waveforms (black curves).

  22. Graph of comparison between the moment magnitude from Global CMT and the moment magnitude obtained from moment tensor inversion.

  23. Graph of comparison between the magnitude determination of BMG and the moment magnitude obtained from moment tensor inversion. Comparison of focal mechanism from Global CMT (blue color) and from moment tensor inversion (red color). Green stars denote epicenters; red triangles denote accelerometer station used for calculating the magnitude of BMG.

  24. CONCLUSION • The magnitude determination of BMG: • The consistency with the other magnitude determination, the magnitude determination of BMG is still sufficient and reliable for shallow earthquakes less than 150 kilometers and the hypocentral distance less than 1,000 kilometers, • The seismic moments and moment magnitude estimated by moment tensor inversion are generally consistent with those of Global CMT, • The magnitude determination of BMG is systematically consistent with the results from moment tensor inversion; • For the rapid earthquakes indicated by short source duration time, the magnitude determination of BMG will become overestimated, • The slow earthquakes indicated by long source duration time, the magnitude determination of BMG will become underestimated.

  25. RECOMMENDATION • The formula of the magnitude determination of BMG could be improved after seismic records are accumulated in BMG. Collecting and storing data should be well managed. The long time collected data is of an urgent necessity in order to make reliable coefficients of this formula. • The dense distribution of accelerometer station should be realized in order to obtain high quality of data. • Intensive study of magnitude determination of tsunami earthquakes will give advantage to improve the reliability of this formula.

  26. Thank you very much

More Related