Claim Missing Document
Check
Articles

Subsurface Structure in Japan Based on P and S waves Travel Time Analysis Using Genetic Algorithm in Japan Seismological Network A. M. Miftahul Huda; Bagus Jaya Santosa
International Journal of Science and Engineering Vol 6, No 1 (2014)
Publisher : Chemical Engineering Diponegoro University

Show Abstract | Download Original | Original Source | Check in Google Scholar | Full PDF (739.853 KB) | DOI: 10.12777/ijse.6.1.70-74

Abstract

Experiment to obtain the subsurface structure in Japan is conducted using seismograms analysis of earthquakes in Japan. All 101 data was used from events in 2012, selected by a maximum depth of 60 km and magnitude between 4.2 to 5.5 Mj. Determination of 1-D subsurface structure is done by utilizing the inversion method with genetic algorithm approach. P wave and S wave velocity structure are determined based on arrival times at receiver. The crustal thickness is known of 33,66 km. P wave velocity for the upper and lower crust, are 6,03 km/s and 6,92 km/s, respectively, and velocity in the upper mantle is 8,18 km/s. S wave velocity for the upper and lower crust are given 3,38 km/s and 3,89 km/s respectively, and the velocity in the upper mantle is 4,59 km/s. If the range integrated to the stable parameter of velocity structure, it shows stable result and the subsurface structure has sufficiently high compatibility.
S-Speed Structure : Seismogram Analysis and Fitting of Earthquake C122297A on RAR Observation Station Bagus Jaya Santosa
Jurnal Pendidikan Matematika dan Sains No 1 (2006): Jurnal Pendidikan Matematika dan Sains Tahun XI
Publisher : Faculty of Mathematics and Natural Sciences, Universitas Negeri Yogyakarta

Show Abstract | Download Original | Original Source | Check in Google Scholar | DOI: 10.21831/jpms.v11i1.12271

Abstract

Penelitian ini bertujuan menemukan struktur kecepatan S antara Pulau New Brittain, PNG dan stasiun observasi RAR di kepulauan Cook, Pacific Barat Daya.Analisis dilakukan dengan seismogram observasi akibat gempa C122297A, dalam domain waktu dan ketiga komponen Cartesian secara simultan. Seismogram sintetik dihitung dengan program GEMINI. Kedua ragam seismogram dikenakan filter low-pass order tujuh dengan frekuensi corner pada 20 mHz.Hasil analisis menunjukkan bahwa penyimpangan yang kuat dan tidak sistematis terjadi pada waktu tiba gelombang S, SS, ScS, ScS-ScS dan gelombang permukaan Love dan Rayleigh. Dengan metode ini terlihat bagaimana pekanya waveform terhadap struktur perlapisan bumi, dibandingkan dengan metode lain dalam analisis seismogram. Untuk menyelesaikan diskrepansi yang dijumpai diperlukan koreksi atas struktur bumi, meliputi ketebalan kulit bumi, gradient kecepatan h, dan besar koefisien-koefisien untuk h dan v di upper mantle, dan sedikit perubahan pada struktur kecepatan S di lapisan-lapisan bumi di bawah upper mantle hingga CMB. Fitting seismogram diperolah dengan baik pada waveform fase-fase gelombang, baik waktu tempuhnya maupun waveformnya. Struktur kecepatan pada model bumi PREMAN dan Ocean memiliki transversal isotrop dari kedalaman 25 hingga 220 km, setelah penelitian ini ternyata model transversal isotrop tersebut dapat dijumpai pada lapisan-lapisan yang lebih dalam, hingga CMB. Ini diperlukan untuk memperoleh fitting pada tiga jenis gelombang sekunder, yaitu S, SS dan ScS-ScS, dimana ketiganya merambat hingga kedalaman yang berbeda-beda. Sifat anisotropi di lapisan dangkal dapat dideskripsikan dengan lebih baik melalui fitting pada gelombang permukaan Love dan Rayleigh.Kata kunci: seismogram, model kecepatan S dari upper mantle - CMB, lautan
Seismic Anisotropy Analysis Beneath Sumatra Revealed by Shear-Wave Splitting Arya Dwi Candra; Bagus Jaya Santosa; Gazali Rachman
Indonesian Journal on Geoscience Vol 4, No 3 (2017)
Publisher : Geological Agency

Show Abstract | Download Original | Original Source | Check in Google Scholar | Full PDF (1955.14 KB) | DOI: 10.17014/ijog.4.3.169-179

Abstract

DOI: 10.17014/ijog.4.3.169-179A shear-wave splitting analysis was determined to observe anisotropic structures of an upper mantle layer beneath Sumatra. The data were collected from 35 BMKG stations with the magnitude of more than 6.25 Mw and the epicentre of 85o - 140o. A shear-wave splitting measurement was calculated by using Splitlab based on three methods simultaneously. The result of the shear-wave splitting measurement in the Sumatra Forearc and Fault Zone shows that there are two anisotropic layers. The first layer has a dominant-fast-polarization direction that is parallel with a trench, and has the delay time duration of 0.5 - 0.9 s-. It is presumed that it is caused by a shear-strain as a result of the existences of Mentawai and Sumatra Fault Zones. The second layer has a dominant-fast-polarization direction that is perpendicular to the trench with the delay time duration of about 1.1 - 1.9 s-. It is presumed that it is caused by a movement of a subduction plate on a mantle wedge. The measurement in the backarc shows that there is only one anisotropic layer that is a subduction plate. It is also found that there is a transition of an orientation change on the subduction plate between Sumatra and Java. The change of the polarization direction is probably related to the age difference and the direction velocity of the absolute plate movement (APM) from Sumatra to Java.
INTERPRETASI METODE MAGNETIK UNTUK PENENTUAN STRUKTUR BAWAH PERMUKAAN DI SEKITAR GUNUNG KELUD KABUPATEN KEDIRI Bagus Jaya Santosa; Mashuri Mashuri; Wahyu Tri Sutrisno; Abdurrahman Wafi; Riski Salim; Radhiyullah Armi
Jurnal Penelitian Fisika dan Aplikasinya (JPFA) Vol. 2 No. 1 (2012)
Publisher : Universitas Negeri Surabaya

Show Abstract | Download Original | Original Source | Check in Google Scholar | DOI: 10.26740/jpfa.v2n1.p7-14

Abstract

Telah dilakukan pengukuran dengan metode magnetik untuk mengetahui struktur bawah permukaan di sekitar gunung kelud. Pengambilan data dilakukan secara acak pada area seluas 0,6 km x 1 km dengan jumlah titik yang diperoleh 244 titik ukur. Proses akusisi dilakukan dengan menggunakan Magnetometer Proton ENVI SCINTREX. Pengolahan data diawali dengan koreksi IGRF dan koreksi variasi harian untuk mendapatkan anomali medan magnet total. Kemudian reduksi bidang datar, kontinuasi ke atas pada ketinggian 100 meter hingga 400 meter di atas sferoida referensi dan hasilnya digunakan untuk pemisahan anomali lokal dan regional. Hasil interpretasi kualitatif menunjukkan adanya anomali dipole magnetik di sebelah timur yang membentang dari arah barat laut ke tenggara sebesar -2125 nT hingga 1863 nT. Metode Talwani 2-D digunakan untuk interpretasi kuantitatif. Model geologi yang dihasilkan adalah patahan atau sesar. Nilai suseptibilitas magnetik di bawah kubah kawah gunung kelud sampai ke gunung lirang (k=0,0124 emu/gram) didominasi batuan basalt, dan di gunung sumbing (k=0,0234 emu/gram 0,0239 emu/gram) yang didominasi batuan andesit.
Computational physics Using Python: Implementing Maxwell Equation for Circle Polarization M Madlazim; Bagus Jaya Santosa
Jurnal Penelitian Fisika dan Aplikasinya (JPFA) Vol. 1 No. 1 (2011)
Publisher : Universitas Negeri Surabaya

Show Abstract | Download Original | Original Source | Check in Google Scholar | DOI: 10.26740/jpfa.v1n1.p1-7

Abstract

Python is a relatively new computing language, created by Guido van Rossum [A.S. Tanenbaum, R. van Renesse, H. van Staveren, G.J. Sharp, S.J. Mullender, A.J. Jansen, G. van Rossum, Experiences with the Amoeba distributed operating system, Communications of the ACM 33 (1990) 4663; also on-line at http://www.cs.vu.nl/pub/amoeba/, which is particularly suitable for teaching a course in computational physics. There are two questions to be considered: (i) For whom is the course intended? (ii) What are the criteria for a suitable language, and why choose Python? The criteria include the nature of the application. High performance computing requires a compiled language, e.g., FORTRAN. For some applications a computer algebra, e.g., Maple, is appropriate. For teaching, and for program development, an interpreted language has considerable advantages: Python appears particularly suitable. PythonŸs attractions include (i) its system of modules which makes it easy to extend, (ii) its excellent graphics (VPython module), (iii) its excellent on line documentation, (iv) it is free and can be downloaded from the web. Python and VPython will be described briefly, and some programs demonstrated numerical and animation of some phenomenal physics. In this article, we gave solution of circle polarization by solving Maxwell equation.
Identification of the Grindulu Fault in Pacitan, East Java using Magnetic Method Latifatul Cholifah; Nurul Mufidah; Eden Lazuardi; Bagus Jaya Santosa; Sungkono Sungkono; Arif Haryono
Jurnal Penelitian Fisika dan Aplikasinya (JPFA) Vol. 10 No. 1 (2020)
Publisher : Universitas Negeri Surabaya

Show Abstract | Download Original | Original Source | Check in Google Scholar | DOI: 10.26740/jpfa.v10n1.p22-33

Abstract

Magnetic method in geophysical surveys is common for its non-destructive use of sub-surface structure delineation. In this study, ground-based measurements of magnetic intensity were performed using a set of instruments in some regions of Pacitan, a city in the southern area of East Java province. Based on these measurements, data acquisition was used to identify the Grindulu faulting zone in the region of interest, potentially vulnerable to geohazards. The data were first corrected using the IGRF and diurnal corrections. A filtering technique of upward continuation at a height of 900 m was then applied to separate local anomalies from regional ones as the targeted sources in the present case. These separate anomalies and their corresponding reductions to the poles as further filtering processes were analyzed for predicting the location and direction of the fault. The results, extracted from data analysis and interpretation, show that the main path of the Grindulu is directed along the NE-SW fault line or N60oE. The resulting anomalies also reflect that the Grindulu is a normal fault with surrounding minor faults lying across the Grindulu, calling for increased awareness of vulnerability in the city to seismic threats.
Estimation of the Source Parameters of the Flores Earthquake and Their Correlation to Aftershocks Emanuel Destianus Banggut; Bagus Jaya Santosa
IPTEK The Journal of Engineering Vol 5, No 1 (2019)
Publisher : Lembaga Penelitian dan Pengabdian kepada Masyarakat

Show Abstract | Download Original | Original Source | Check in Google Scholar | DOI: 10.12962/j23378557.v5i1.a5021

Abstract

This study uses data provided by BMKG Indonesia which is downloaded from the website www.webdc.eu. The data processing is done by using some software, one of them is Matlab to run Isola software. The software provides information about the fault direction. In addition, this study also uses Coulomb 3.3 software to provide information on subsequent earthquake forecasts. The results showed that the earthquake that occurred on 27th February 2015 gave impact to Coulomb's increase and decrease of stress to some areas in the northern mainland of Flores. The impact of this major earthquake is expected to trigger an earthquake that occurred in northern Flores. Based on the results of data processing obtained form focal mechanism is oblique reverse with the first nodal value strike 139, dip 73, and rake 50. While on the second nodal has a strike value 30, dip 43, and rake 154. Main earthquake impact Coulomb stress increase of 0.15-0.2 bar. Meanwhile, Coulomb stress result obtained shear value -18.890, normal -0.217, and Coulomb -19.976.
3D Tomographic Imaging of P Wave Velocity Structure Beneath Java Island using Fast Marching Tomography Method Uswatun Chasanah; Bagus Jaya Santosa
IPTEK Journal of Science Vol 2, No 1 (2017)
Publisher : Lembaga Penelitian dan Pengabdian kepada Masyarakat

Show Abstract | Download Original | Original Source | Check in Google Scholar | Full PDF (671.178 KB) | DOI: 10.12962/j23378530.v2i1.a2256

Abstract

Java Island one of the best locations for geophysics research, as it’s located near the edge of the junction of continental Eurasia Plate and the Indo-Australia Plate with the movement of the plates 6 cm/year. Its needed the accurate imaging of subsurface structure to understanding complex tectonic setting. The first arrival time from local earthquake of Java (6,360-9,150 S and 105,890-115,540 E) with magnitude greater than 4,5 Mw occurred from 2011 until 2013 recorded by local network seismograph has been inverted for three dimensional variation of the depth to the P wave velocity in the Java island. At the same time, earthquake hypocenter location has been corrected simultaneous. This research apply a new tomographic inversion scheme FMTOMO that has been developed by Rawlinson (2004).  The results of inversion show that there are three layer in 100 km from the surface at continental plate side (Eurasia). The continental upper crust has P wave velocity variation about 4-5,5 km/s in the depth until 30 km from the surface. There are low velocity anomaly inclined towards the slab which probably have associated with shallow earthquake along the slab
Coulomb Stress Change Analysis Center of Celebes on 29th May 2017 6.6 Mw Earthquake and Aftershocks Distribution Ibrahim Ibrahim; Gazali Rachman; Bagus Jaya Santosa
IPTEK Journal of Proceedings Series No 1 (2019): 4th International Seminar on Science and Technology 2018 (ISST 2018)
Publisher : Institut Teknologi Sepuluh Nopember

Show Abstract | Download Original | Original Source | Check in Google Scholar | Full PDF (545.723 KB) | DOI: 10.12962/j23546026.y2019i1.5114

Abstract

Mechanism of earthquakes associated with the distribution of stress static that occurs in rocks. When the rocks elastic limit is exceeded there will be a release of energy as an earthquake result as rocks no longer able to withstand the stress that will disturb the stress field in the neighborhood. In this study, analysis stress of changes was done by taken earthquake data Center of Celebes 6.6 Mw on May 29th, 2017 with hypocenter 13 km using four earthquake recording stations namely BKB, TOLI, PMSI, and LUWI through website GFZ (Geo Forschungs Zentrum) and Global CMT (Global Centroid Moment Tensor), and then calculated the earthquake source parameters so as to obtain model focal mechanism using discretization methods are iterative wave numbers and analyzed using Coulomb 3.3 to obtain the value of Coulomb stress change and its aftershocks distribution. Analysis results showed that orientation the focal mechanism model of earthquake fault plane has been a normal fault type, fault length 27.54 km, width fault 14.06 km with slip shift of 79.06 cm. Coulomb stress Changes are generated ranging from 0.05 to 0.2 bar trending southwest-northwest and northeast-southeast of the epicenter, Based on Coulomb stress plot that center of Celebes 6.6 Mw on May  29,  2017 earthquake triggering aftershocks on May 29th, 2017 at 14:53:44 UTC with latitude and longitude -1.12 ° and 120.17 ° northwest trending with a range of values from 0.1 to 0.05 bar, on May 31, 2017 at 04:42:06 with latitude and longitude -1.17 ° and 120.79 ° north-east trending with a range of values from 0.15 to 0.1 bar and on November 25, 2017 at 09:14:51 UTC with latitude and longitude -1.18 ° and 119.93 ° and 11:11:24 UTC with latitude and longitude -1.19 ° and 119.94 westbound with a range of values from 0.05 to 0.01 bar
Coulomb Stress Analysis of Halmahera Earthquake on June 7th 2016, Mw=6.3 and Its Correlation with Aftershocks Fathony Akbar Pratikno; Gazali Rachman; Bagus Jaya Santosa
IPTEK The Journal for Technology and Science Vol 30, No 2 (2019)
Publisher : IPTEK, LPPM, Institut Teknologi Sepuluh Nopember

Show Abstract | Download Original | Original Source | Check in Google Scholar | Full PDF (737.048 KB) | DOI: 10.12962/j20882033.v30i2.5006

Abstract

Main earthquake with magnitude Mw 6.3, has been occurred in the West of Halmahera Island, Northern Molucca on June 7th, 2016. The activity of the Pacific plate, Eurasia plate and Indo-Australia plates trigger some earthquake around Halmahera. The purpose of the current study was to determine the Coulomb stress change of the main earthquake and predict the aftershocks’ location around it. For understanding an interaction of the faults, this study used Coulomb stress change model. Three component seismic waveform data recorded by TNTI station within 132 km, SANI station within 381.2 km, LUWI station within 478.3 km, and TOLI2 station within 617.6 km of the epicenter. This study used software ISOLA-GUI to obtain the parameters of the earthquake source. The parameters of the earthquake source were then used to determine fracture orientation, length, width and slip displacement. It was also used to calculate the Coulomb stress changes around the main shock, using software Coulomb 33. The result shows that the orientation of the earthquake fault has a dip angle of 48° to the horizontal plane and has a strike of 210° against the North. The fault length is 19.49 km, the fault width is 11.59 km and slip displacement is 49.43 cm. Based on plotting result, the main earthquake in Northern Molucca on June 7th, 2016 has positive Coulomb stress change that spread across the Northeast to the Southwest. This change is likely triggered three aftershocks with a range of Mw 4.7 to 6.2.