Untitled

SEACG 2022 Program Book
Welcome to SEACG 2022 Head of Geophysical Engineering

It is our pleasure to welcome you Magister & Doctorate Degree
to the 2022 4th Southeast Asian Program’s Welcome
Conference on Geophysics. Our Ladies and gentlemen,
current theme is "Great It’s my pleasure to welcome all of
Challenges and Opportunities of participants, presenters, and in
Geophysics Today and Future". particular, keynote and invited
The main objective of this speakers from all over the
conference is to lead the country and to introduce the
development of new ideas and two-day 4th Southeast Asian
technologies in geophysical Conference on Geophysics
methods and related science. This (SEACG) in my capacity as head
conference will be the meeting of Geophysical Engineering
point between young and senior generations of Department. As a matter of fact, we still must be aware with
geophysicists so that the sharing experiences can be social distancing norm due to the COVID-19 pandemic,
conducted effectively and efficiently. therefore the conference is conducted by hybrid. We hope
The SEACG 2022 will be held in a hybrid conference that all of you will share the knowledge and experience with
(physical and virtual), live from Bandung, Indonesia. The others. Along with a wonderful organizing committee, we
program will last for two days, from August 9, 2022, hope the conference will be remembered for many years to
consisting of the keynote and invited sessions with come.
preeminent speakers and parallel sessions covering a wide The theme of two-day conference is “Great challenges and
range of topics, including geosciences research and opportunities of geophysics today and future”. We are
development. There are 11 invited international speakers looking forward to talking in plenary seasons on issues
and 88 contributed speakers who will present their ideas. about recent development in geophysical technology and
With five panel rooms, 35 contributed speakers will deliver related sciences. The conference is created to provide an
their presentations offline and 53 virtually. Please visit our opportunity for all of us to participate in the discourse
website (https://seacg.fttm.itb.ac.id) for further through research, experiment, and field studies.
information.
We certainly encourage all of you to attend and engage with
We would thank the organization's team members, the other participants and presenters. I do hope this will be an
reviewers, and the faculty staff. They have worked very interesting experience. We hope that the use of the hybrid
hard to make this event possible. We also would like to conference technology will pave the way for different types
thank the sponsors of this conference: The Institute of of future conferences that can integrate physical and virtual
Research and Community Services – ITB and The World conference, giving wider opportunities to those who would
Class University Program – ITB. find it difficult to attend in person.
We conclude by reflecting those better times will come and
We wish all attendees of SEACG 2022 have an exciting
we sincerely hope that in the not-too-distant future you will
experience and enjoy the moment in this scientific forum.
come and visit us in Bandung. We look forward to seeing
We look forward to seeing you in the SEACG 2024.
you at Bandung’s SEACG 2022, which I am sure that it will
be not only exciting technically but also good fun. Finally, I
Dr. Indra Gunawan hope all of you will enjoy the two-day conference. Thank
Chairman SEACG 2022 you very much.
Dr. Ir. Fatkhan, M.T.

Head of Geophysical Engineering
Magister & Doctorate Degree Program
Faculty of Mining and Petroleum Engineering
Institut Teknologi Bandung
2
SEACG 2022 Program Book
Organizing Committee
Steering Committee Prof. Djoko Santoso
Prof. Nanang T. Puspito
Time Zone
Prof. Sri Widiyantoro Please mind that the time zone for this hybrid conference is
Prof. Wawan G. A. Kadir Western Indonesia Time (UTC+7)
Prof. Sigit Sukmono
Prof. Andri D. Nugraha
Dr. Fatkhan Conferences Session
Dr. Warsa
Dr. Moh. Rachmat Sule Keynote or Invited Session will be held at the beginning of
each day conference (9 and 10 August 2022). Please find
Editor Committees Prof. Hendra Grandis the schedule and details on the Highlights Page.
Prof. Satria Bijaksana
Ign. Sonny Winardhi, Ph.D A parallel scientific session will be held on two days (9 and
Dr. Tedi Yudistira 10 August 2022). Three physical rooms and five virtual
rooms (hosted by Zoom links) will be used simultaneously
Chairman Dr. Indra Gunawan for the presentation. Please mind your zoom room number,
General Affair Manager date, time, and session.
Dr. Zulfakriza
Rizka, MT
Secretariat Officer
Aditya Lesmana, MT
Oral Presentation Guidelines
Dr. David P. Sahara Pre-online oral presentation
Treasurer
Dian Kusumawati, M.Si The presenter will conduct the oral presentation offline or
Public Relation & Ekkal Dinanto, MT online using the sharing screen feature in the Zoom
Information Dr. Fernando Lawrens application. Therefore, to avoid connection problems that
cause delayed schedule, the presenter may send the
Event Silvia Jannatul Fajar, MT
presentation material or recorded video to the committee
Dr. Shindy Rosalia
if they want the committee to play it during the oral
Riskiray Ryannugroho, MT
Logistic presentation session.
Vigo Denis Radestian
Rifki The slide presentation is in Microsoft PowerPoint (*.pptx)
at least Office 2010.
Documentation Faridz Nizar Ahmady, MT
Dr. Dadi Abdurrahman Online oral presentation (Live)
Official Team
Dona Sita Ambarsari, ST, MT 1. The presenter must be ready in the online meeting
Dr. Endra Gunawan room at least 15 minutes before the presentation
Madaniya Oktariena MT begins. The presenter should check the oral schedule
Firda Haifa Fadhillah at: https://seacg.fttm.itb.ac.id/schedule-seacg-2022/
Sindi Hajah Pattimah 2. Presenters use English in their oral presentations.
Maulidia Ain Bening 3. The committees determine the order of presentation
Ilona Jocelynn Untung 4. Each presenter must ensure that their presentations
Velly Amalia Hoja P. are not more than 15 minutes, and it will be followed
Ahmad Toifur by 5 minutes of discussion.
Ziyad Fakhri Kunadi 5. The committee will give time remainder in five, three
Ibra Kumara and one minute before the presentation time is over.
Widjra Cyiena Christi N. 6. All participants are prohibited to unmute the mic
Sari Ratna Dewi during oral presentation.
Syaiful Apri Kurniawan
During the online sessions, all participants and presenters
are expected to maintain peace and order.
Important Dates
[Extended to]
Abstract Submission
20 Jun 2022
Abstract Accepted
30 Jun 2022
Notification
Extended Abstract
01 Jul - 30 Aug 2022
Submission
Registration 01 Jul - 06 Aug 2022
SEACG 2022 Conference 9 - 10 Aug 2022
3
SEACG 2022 Schedule
Highlights
Day-1
Tuesday, 9 August 2022 SCIENTIFIC PARALLEL SESSION I
SEACG 2022 OPENING TIME 10.30 – 11.50 (UTC+7)
Room Topics
Opening and Welcome Speech
ZOOM 1 Fault and Deformation Study (1)
ZOOM 2 Passive Seismic
ZOOM 3 Environmental Geophysics (1)
ZOOM 4 Geophysical Approaches in
Hydrocarbon (1)
ZOOM 5 Geophysical Approaches in
Geothermal
<< 11.50 - 12.45 SESSION BREAK >>

Dr. Indra Gunawan as SEACG 2022 Chairman
Dr. Fatkhan as Head of Geophysical Engineering INVITED SESSION II
Magister & Doctorate Degree Program
TIME 08.00 – 08.30 (UTC+7)
ROOM ZOOM 1
KEYNOTE SESSION
TIME 12.45 – 13.45 (UTC+7)

ROOM ZOOM 1
MODERATOR Rexha verdhora Ry, S.T., M.T.
SPEAKERS Prof. Phil Cummins
TIME 08.30 – 09.15 (UTC+7) Australian National University
ROOM ZOOM 1 “A New Seismotectonic Model and
MODERATOR Andri Hendriyana, S.T, M.T. Earthquake Hazard Assessment for Papua
SPEAKERS Dr. Rer. Nat. R. Mohammad Rachmat New Guinea”
Sule Prof. Sri Widiyantoro
Institut Teknologi Bandung, Indonesia Institut Teknologi Bandung / Maranatha
Christian University
“Fault Locking south of Jakarta inferred
INVITED SESSION I from an investigation of seismic activity
along the Baribis fault, northwestern Java,
Indonesia.”
<< 13.45 - 14.00 COFFEE BREAK >>
SCIENTIFIC PARALLEL SESSION II

TIME 14.00 – 15.40 (UTC+7)
Room Topics
TIME 09.15 – 10.15 (UTC+7) ZOOM 1 Fault and Deformation Study (2)
ROOM ZOOM 1 ZOOM 2 Advance geophysics (1)
MODERATOR Dr. Zulfakriza ZOOM 3 Passive seismic (2)
SPEAKERS Dr. Tedi Yudistira ZOOM 4 Geophysical approaches in hydrocarbon (2)
Institut Teknologi Bandung, Indonesia ZOOM 5 Volcanology (1)
“Seismic Velocity Structures Beneath Merapi
Volcano Complex from Ambient Seismic
Noise Tomography”
Dr. Erdinc Saygin
CSIRO, Australia
“4D Subsurface Monitoring & Imaging
with Novel Sensor Platforms and
Correlation Wavefield”
<< 10.15 - 10.30 COFFEE BREAK >>
4
SEACG 2022 Schedule
Day-2r
Wednesday, 10 August INVITED SESSION IV

INVITED SESSION III
TIME 12.45 – 14.15 (UTC+7)

TIME 08.00 – 09.30 (UTC+7) ROOM ZOOM 1
ROOM ZOOM 1 MODERATOR Dr. Endra Gunawan
MODERATOR Dr. Fernando Lawrens SPEAKERS Prof. Wawan Gunawan
SPEAKERS Prof. Takeshi Tsuji Institut Teknologi Bandung, Indonesia
University of Tokyo “An opportunity and contributions of
“Continuous monitoring based on Portable geophysical approach in developing
Active Seismic Source (PASS).” Carbon Capture, Utilization, and Storage
Prof. Wahyu Srigutomo (CCUS)”
Institut Teknologi Bandung, Indonesia Dr. Zhijun Du
“Electromagnetic modeling and its applications Ocean University of China
for volcanological and geothermal studies in “Hydrocarbon Mapping and CO2
Indonesia” Monitoring from Integrated Analysis of
Multi-Method Geophysical Data”
Moeto Fujisawa, M.Sc.
JGI, Inc. Dr. Jean-Philippe Metaxian
“Diffraction Imaging using Microseismic Institut de Physique du Globe de Paris
Record Obtained by DAS.” “Strombolian seismic activity
characterization using fibre-optic cable and
distributed acoustic sensing.”
<< 09.30 - 09.45 COFFEE BREAK >>
SCIENTIFIC PARALLEL SESSION III << 14.15 - 14.30 COFFEE BREAK >>
TIME 09.45 – 11.45 (UTC+7)
Room Topics SCIENTIFIC PARALLEL SESSION IV
ZOOM 1 Tomography (1) TIME 14.30 – 15.30 (UTC+7)
ZOOM 2 Advance geophysics (2) Room Topics
ZOOM 3 Environmental Geophysics (2) ZOOM 1 Seismology
ZOOM 4 Seismic hazards (1) ZOOM 2 Tomography (2)
ZOOM 5 Geomagnetic ZOOM 3 Environmental Geophysics (3)
ZOOM 4 Volcanology (2)
<< 11.45 - 12.45 SESSION BREAK >> ZOOM 5 Volcanology (3)
CLOSING
TIME 15.30 – 16.00 (UTC+7)
The 4th Southeast Asian Conference on Geophysics

Hybrid Conference
09-10 August 2022
https://seacg.fttm.itb.ac.id
5
SEACG 2022 Room
SEACG Room Layout

Room Number Room Link Meeting ID Passcode
Zoom 1 https://bit.ly/SEACG2022_ZOOM-01 987 0782 7440 FTTM-ITB
Zoom 2 https://bit.ly/SEACG2022_ZOOM-02 937 4595 3793 SEACG2022
Zoom 3 https://bit.ly/SEACG2022_ZOOM-03 959 7792 0220 814255
Zoom 4 https://bit.ly/SEACG2022_ZOOM-04 933 0805 6724 762844
Zoom 5 https://bit.ly/SEACG2022_ZOOM-05 926 4900 8551 SEACG2022
*If the links don’t work, please visit seacg.fttm.itb.ac.id for newest updates
6
SEACG 2022 Parallel Scientific Session Day-1 | Session 1
TUE, 09 AUG 2022 | Session 1
Zoom 1 Zoom 2 Zoom 3 Zoom 4 Zoom 5
DAY 1 – TUE
Fault and Deformation Study (1) Passive Seismic (1) Environmental Geophysics (1) Geophysical approaches in hydrocarbon (1) Geophysical approaches in geothermal
(UTC+7) Chairperson: Dr. David Prambudi Sahara Chairperson: Dr. Rexha Verdhora Chairperson: Dr. Dadi Abdurrahman Chairperson: Ekkal Dinanto, MT Chairperson: Dr. Fernando Lawrens
SEACG-004 SEACG-030 SEACG-054 SEACG-017 SEACG-052
Investigation Of Postseismic Deformation Of The 2018 Data Processing of Surface Wave Tomography in Detection of Old Tram Tracks Near South Part of Saturated Fluid Variation Effect on Seismic Anisotropy: Characteristics of Subsurface Conditions of The
10.30 - Lombok Earthquake Sequence Using Gps Data Australia using Automated Surface Wave Phase- Surabaya Zoo Using Georadar Method A Case Study of Deep-Water Reservoir Selogajah Hotsprings Area Using Very-Low Frequency
1 Risya Fauziyyah, Endra Gunawan, Sri Widiyantoro, Velocity Measuring System (ASWMS) Rayhan Farisi Ramadhan, Mariyanto Mariyanto, JPGN Madaniya Oktariena, Wahyu Triyoso, Fatkhan Electromagnetic (VLF-EM) Method
10.50
Irwan Meilano, Syamsuddin Amira Saqina, Nicolas Silaen, Riyan Yusami, Yusrina Rochman Fatkhan, Sigit Sukmono, Dona Sita Ambarsari, Rahmi Azizah, Mariyanto Mariyanto, M. Singgih
Alfiani Haqi, Sandy K. Suhardja Erlangga Septama, Pongga D. Wardaya, Befriko S. Purwanto
Murdianto
Study On Estimation Of Fault Plane And Rupture Processing Receiver Function Data in Mentawai Steam-Vent Identification using Frequency Fault Delineation and Geometrical Modeling of the 3D Subsurface Characterization using Gravity Data of
10.50 - Duration For Potential Assessment Of Tsunami Region, Sumatra Characteristic Method in Ground Penetrating Radar Geological Carbon Storage (GCS) Candidate for the Mataloko Geothermal Field, Flores, Nusa Tenggara
2 Hazards In Southern Part Of Sumatera And Southern M.Aldi, Nicolas Silaen, Dayyan Dhaifullah, Sandy K (GPR) CCS and CCUS Facilities Based on Gravity Dataset. A Timur, Indonesia
11.10 Part Of West Java Case Study at the Northtern Part of East Java,
Suhardja Radithya A. Pambudi, Warsa, Wahyudi W. Parnadi Fadhil Farras, Susanti Alawiyah, Setianingsih, Eko
Ajeng Marina Utamie, Nanang T. Puspito, Pepen Indonesia Januari Wahyudi
Supendi Accep Handyarso
Moment tensor inversion implementation in Preliminary Results on Low Frequency Passive Detecting Subsurface Utilities with Multi-frequency and The Effects of Porosity on the Anisotropy Parameters Identification of Shallow Fault with Seismic Refraction
determining focal mechanism solution of Palu-Koro Seismic (LFPS) Study in Lake Toba Multi-antenna Configuration GPR of the Slope Fan Facies Sand Reservoirs in the Modeling in Gondang Area, Bojonegoro
and Matano fault events: processing strategy and Riyan Yusami, Sandy K Suhardja, Dzikri N Rochim, Wahyudi W Parnadi, Indra Gunawan, Warsa, Agus Deepwater Kutai Basin, East Kalimantan, Indonesia Nahari Rasif, Aditya Aulia Rachman, Mariyanto
preliminary result Nicolas Silaen, Amira Saqina Laesanpura, Djoko Santoso Rahma Nurhazlita, Teuku Abdullah Sanny, Dona Sita Mariyanto
11.10 -
3 Dian Kusumawati, David P. Sahara, Nanang T. Ambarsari
11.30 Puspito, Mudrik R. Daryono, Andri D. Nugraha,
Mohamad Ramdhan, Andri Kurniawan, Wahyu
Tanihaha, Leonardo Junior Johan Solihin, Muhamad
Diva Pratama, Sindi Hajah Patimah, Annisa Trisnia
Sasmi, Hilmy Muhammad, Hammam Muhammad,
Widjra Cyiena Christi Natafrisca, Maulidia Ain Bening
SEACG-058 SEACG-015 SEACG-079 SEACG-011
Geoelectric Analysis for Potential Slip Plane Sediment thickness in the eastern part of Kalimantan GPR Impedance Inversion for Imaging Groundwater Tomographic Seismic Model Using Ray Tracing
11.30 - Interpretation in Labuan Bulan Hill, Salem District, from high-frequency receiver functions Level: a Case Study at PPSDM Test Site, Cipatat, Algorithm Dijkstra to Visualize Structure In The
4 Brebes Regency, Central Java Syuhada Syuhada, Bayu Pranata, Titi Anggono, West Java Province Geological Model Of The Wayang Windu Geothermal
11.50 Field, West Java
Indra A. Dinata, M. Khalifah B.Y. Praja, Imam A. Mohamad Ramdhan, Febty Febriani, Aditya D Wahyu W. Parnadi
Sadisun, Agus Kuswanto Prasetio, Cinantya N Dewi, Mohammad Hasib Muhammad Rizqi Fadhlurrahman, Indriati Retno
Palupi
SESSION BREAK
7
SEACG 2022 Parallel Scientific Session Day-1| Session 2
TUE, 09 AUG 2022 | Session 2
DAY 1 – TUE Fault and deformation study (2) Advance geophysics (1) Passive seismic (2) Geophysical approaches in hydrocarbon (2) Volcanology (1)
(UTC+7) Chairperson: Dr. Endra Gunawan Chairperson: Dr. Indra Gunawan Chairperson: Dr. Shindy Rosalia Chairperson: Dona Sita Ambarsari, S.T., M.T. Chairperson: Dr. Zulfakriza
SEACG-022 SEACG-001 SEACG-059 SEACG-036
Late Quaternary Deformation of the Ambon fault, Early Results of Modeling Subsurface Study in Lake Application Of Seismic Multiattribute for Predicting The Geophysical Overview over the Paleo-Volcanic-
14.00 – Indonesia: Preliminary Results Toba, North Sumatra Using Receiver Function Distribution Of Reservoar Property In “A” Field, Kutai Caldera-Rims Anomaly, from Synthetically
6
14.20 Adi Patria, Danny Hilman Natawidjaja, Muhammad Lio Julio, Nicolas Silaen, Wulan S Ratna, Indra R Basin Modelling to the Field Cases
Hanif, Hiroyuki Tsutsumi Siregar, Sandy K Suhardja P D Wardaya, A Novianto, A K Pratiknyo, H Hamdalah, Accep Handyarso
Y Rizkianto, A M Handini
Neotectonics of the eastern Matano fault in Sulawesi, A Freeware Tool for Masw Data Processing In Seira Early Forward RF Model In South Central Java Delineation Of Gas Sand Reservoir Prospect in The “A” Analyses of Explosion Earthquakes at Sinabung
14.20 - Indonesia Software To Determine Shear Wave Velocity Field, Deep Water Area Of Kutai Basin Volcano: Characteristics of Waveform, Spectra, and
7 R. P. C. Tambunan, S. K. Suharja
14.40 Adi Patria, Danny Hilman Natawidjaja, Mudrik Putu Pradnya Andika, Tedi Yudistira, Ruhul Firdaus, P D Wardaya, A Novianto, A K Pratiknyo, H Hamdalah, Energy
Rahmawan Daryono, Muhammad Hanif, Anggraini Imam Catur Priambodo Y Rizkianto, A M Handini Aulia Dhita Pritaningtyas, Mohammad Hasib, Estu
Rizkita Puji, Hiroyuki Tsutsumi Kriswati, Tati Zera
Mataram Fault Section – New Active Fault Section One Dimensional Magnetotelluric Data Inversion Extracting Rayleigh wave phase velocities from Low Frequency Seismic Extrapolation in Full Waveform Geological Characteristics Identification of The
14.40 – Crosses East-West in The Centre of Yogyakarta City Modeling Using Convolutional Neural Network ambient noise cross-correlation: application to Inversion (FWI) using U-Net Western Tangkuban Parahu Volcanic Area Using
8
15.00 Mudrik R. Daryono, Astyka Pamumpuni, Dadan D. Muhammad Iqbal Khatami, Hendra Grandis Jakarta seismic network Asido Saputra Sigalingging, Sonny Winardhi, Ekkal H/V Analysis
Wardhana, Danny H. Natawidjaja Rexha V. Ry, Phil R. Cummins, Babak Hejrani, Sri Dinanto Vigo Denis Radestian, Zulfakriza, Andri Dian
Widiyantoro Nugraha, Aditya Lesmana
Geological slip rate study on Renun Segment of Numerical Example to Estimate Groundwater Level Preliminary Result of Shear-wave Velocity Estimation Data Conditioning for Seismic Reservoir The Interpretation of Subsurface Structure of Banda
15.00 – Sumatran Fault in Dolok Sanggul area, North Changes using Time-lapse Microgravity Data in Central Part of Sulawesi Using Microtremor Data Characterization of Thin Bedded Fluvio-Deltaic Api Volcano Based on Gravity Modelling
9 Sumatra using UAV photogrammetry and DEMNAS Eko Januari Wahyudi, Setianingsih, Dasapta Erwin Sari Ratna Dewi, Zulfakriza, Sri Widiyantoro Reservoirs Faridz Nizar Ahmady, Djoko Santoso, Fauzan
15.20
Mudrik R. Daryono, Muhammad Hanif, Dadan Dany Irawan, R. M. Rachmat Sule Andri Hendriyana, M. Noor Alamsyah, Dona Sita Ardana
Wardhana, Danny Hilman Natawidjaja Ambarsari, Istifani Ferucha, Diah Ayu Widianingsih,
Sigit Sukmono
Correlation of Stress reconstruction, Correlation Preliminary Results and Interpretation of Gravity Modeling Crustal Structure in West Sumatera Using 3D Check shot Modeling for Time-Depth Conversion Digitization of Analog Seismograms using Machine
15.20 – Dimension (Dc), and Seismic Moment Rate Loading Survey Over Pre-Tertiary Rock in Northern Kebumen Teleseismic Receiver Function Learning Methods
10 Jefri Bambang Irawan, Wahyu Triyoso
15.40 along the Sumatra subduction margin Faridz Nizar Ahmady, Djoko Santoso, Susanti Wulan Sri Ratna, Lio Julio, Indra Rivaldi S., M. Aldi, Martanto Martanto, Andri Dian Nugraha, David P.
Wahyu Triyoso, Dina A. Sarsito, Sigit Sukmono Alawiyah Mohamad Ramdhan, Sandy Kurniawan Suhardja, Sahara, Zulfakriza Zulfakriza, Devy Kamil
Iktri Madrinovella, Ida Herawati, Bayu Pranata Syahbana, Corentin Caudron
8
SEACG 2022 Parallel Scientific Session Day-2| Session 3
WED, 10 AUG 2022 | Session 3

DAY 2 – WED
Tomography (1) Advance geophysics (2) Environmental Geophysics (1) Seismic hazards (1) Geomagnetic
(UTC+7) Chairperson: Dr. David Prambudi Sahara Chairperson: Dr. Andri Hendriyana Chairperson: Dr. Dadi Abdurrahman Chairperson: Dr. Zulfakriza Chairperson: Faridz Nizar Ahmady, M.T
The early model of tomography in Eastern Indonesia Applying Map Filters and Inversion Model for Gravity Identification Of Aquifer Layer Using Geoelectrical Early Analysis of Vulnerability Study in Lake Toba Geomagnetic Data Analysis to Identify Subsurface
09.45 – using FMTOMO Data on Matlab-based Application Resistivity Method with Schlumberger Array In Bumi Dzikri N Rochim, Nicolas Silaen, Riyan Yusami, Condition of Gondang Area, Bojonegoro
1
10.05 Elta Purnama Wulan, Sinta Nur Asyidah, Yusrina M Tsaqif Ghifari, Fawwaz Byru Fitrianto, Indriati Harapan Village, New Capital City Of Indonesia Sandy K Suhardja, Soni Satiawan, Dicky A Zaky, M. Yolanda Mustika Bohal, Mariyanto Mariyanto,
Alfiani Haqi Retno Palupi, Wiji Raharjo Fahmi Saputra, Febrian Dedi Sastrawan, Meidi Ramdhan Ihsan Fermantono
Arisalwadi
2D And 3D Modelling Of Molucca Sea Double Tomographic Inversion Code Analysis using Characterization Of Agricultural Land Using Development of a Land Subsidence Fragility Curve in Characterization of Caldera Lakes’ Sediments
10.05 – Subduction Zone Based On Primary Waves Travel Regularization Iteratively Reweighted Norm (IRN) Microtremor Approach Bandung Basin: a Preliminary Result Through Rock Magnetic and Geochemical
2 Time Seismic Tomography Inversion Belista Ajeng Natasya, Ryandika Syahputra, Elta W, Edlyn Yoadan Nathania, T. A. Sanny Irwan Gumilar, Riantini Virtiana, Teguh P. Sidiq, Brian Analyses. Case Study: Lakes Buyan and
10.25 Gendis Ashella Nareswari Pratita, Indriati Retno Kamaludin Nursal, Aulia Rizki, Nicolas Silaen, Sandy Bramanto, Ramdan Sudarsono, Titis Idam Haq Tamblingan in The Island of Bali
Palupi, Prasetyadi K Suhardja Ni Komang Tri Suandayani, Satria Bijaksana,
Darharta Dahrin, Irwan
Early Results of Surface Tomography using An Attempt to Invert Tomography Delay Time Water Pipeline Leak Identification using Ground How does A Fast-Cascading Earthquake in a Multi- The Use of Magnetic Method to Identify of Rare
10.25 – Automated Surface-Wave Phase-Velocity Measuring Problem with Hybrid Lsqr Penetrating Radar Method Scale Fracture Network Affect Ground Motions? Earth Elements in Lateritic Soils from Southeast
3 System (ASWMS) in Hawaii Ryandika Syahputra, Elta W, Kamaludin Nursal, Audi Dian Bakti, Mariyanto, Eko Minarto, Tigar Ade Kadek Hendrawan Palgunadi, Alice-Agnes Gabriel, Sulawesi
10.45 Kamaludin Nursal, S.K. Suhardja Nicolas Silaen, Belista Ajeng N, Sandy K Suhardja Atma Wijaya Dmitry Garagash, Thomas Ulrich, Paul Martin Mai Khalil Ibrahim, Satria Bijaksana, La Ode
Ngkoimani, Putu Billy Suryanata, Silvia Jannatul
Fajar
Early Result of Imaging 3D Seismic Velocity Structure Ultrasonic Tomography of Core Plug Using Fresnel Application of Crosshole Seismic Tomography to Calculation of Logarithmic Ratio of Energy-Seismic REEs Content and Magnetic Properties in Placer
10.45 – In Central Java Using Double-Difference Tomography Volume Approach Inferred Cavities Condition Moment, Rupture Duration, and Mwpd for the Deposits from Central Kalimantan
4 Ryandika Saputra, Sandy Kurniawan Suhardja, Iktri Hafidz Naufal Aryan, Teuku Abdullah Sanny, Dona Fatkhan, Riskiray Ryannugroho, Fernando Lawrens, Tsunami Early Warning System in Bali and Nusa Gusfrimanuel, Satria Bijaksana, Putu Billy
11.05 Madrinovella, Mohammad Ramdhan Sita Ambarsari Andri Dian Nugraha, Mohammad Rachmat Sule Tenggara Suryanata, Silvia Jannatul Fajar
Muh. Soekarno Saputra Rahman, Nanang T. Puspito,
Pepen Supendi
Preliminary Result of Local Earthquake Tomography Constrained Global Optimization of 1D Horizontal Reconstruction of Retrogressive Landslides Using Seismic Hazard Study Of The 2018 Situbondo Inverse modeling using geomagnetic data of field
Around Palu Koro Fault using BMKG Seismic Electric Dipole Time Domain Electromagnetic Method Integrated Technique of Geoelectric Resistivity and Earthquake For Situbondo And Its Surrounding Cities camp geophysics in Karangsambung (2005-2019)
11.05 – Network (HED TDEM) using Multi-Objective Particle Swarm Borehole Logging in Sumampir, Purbalingga Sindi Hajah Patimah, Wahyu Triyoso, Endra Eko Januari Wahyudi
5 Optimization (MOPSO) Regency, Central Java
11.25 H Jayadi, B. J. Santosa, A. D. Nugraha, D. D. Gunawan
Warnana, R. Jamroni, P. Supendi, Z. Zulfakriza, F. I. G. P. Fadjar Soerya Djaja, Enjang J. Mustopa, Imam Achmad Sadisun, Indra A. Dinata, Gunawan
Muttaqy, G. Rachman, A. S. Putra, I. U. Meidji Wahyu Srigutomo Wibisono
SEACG-077 SEACG-086
Identification of Andesite Distribution in Pasrepan Magnetic Basement Modelling around Misool
11.25 – Sub-District, Pasuruan District Using Schlumberger Island Raja Ampat Regency With MagB_Inv
6 Configuration Resistivity Method Syaiful Apri Kurniawan, Satria Bijaksana, Silvia
11.45 Flavio Da Conceição Martins, Almer Harahap, Jannatul Fajar, Putu Billy Suryanata
Mariyanto, Nita Ariyanti, Mochammad Fauzan
Dwiharto, Ahmad Wafi Irsyad
SESSION BREAK
9
SEACG 2022 Parallel Scientific Session Day-2 | Session 4
WED, 10 AUG 2022 | Session 4

DAY 2 – WED
Seismology Tomography (2) Environmental Geophysics (3) Volcanology (2) Volcanology (3)
(UTC+7) Chairperson: Dr. Shindy Rosalia Chairperson: Dr. Rexha Verdhora Chairperson: Dr. Eko Januari Wahyudi Chairperson: Dr.Eng. Fernando Lawrens Chairperson: Dr. Endra Gunawan
14.30 – Preliminary Study Of Hyposenter Relocation Early Analisis on Teleseismic Shear Wave Splitting in The Correlation of Land Survey and Satellite Gravity Early Model on Teleseismic Study in Merapi Merbabu Semeru Volcano Eruption Activity in 2021 Based
6 Distribution Using Double-Difference Method In Central Java Data, Study Case: Timor Leste Area Sirilus C. Sibagariang, Sandy Kurniawan Suhardja, on Seismic and Deformation Monitoring Data
14.50 Bayah Dome Area, Lebak, Banten Eka W. A. R. P., Cicilia, Rita Saputri, Sandy K. Ahmad Toifur, Ziyad Fakhri Kunadi, Susanti Alawiyah, Sinta N. Asyda, Cecilia E.P, M. Kholil Kristianto, Ahmad Basuki, Heruningtyas Desi P.,
Agung Sabtaji, Hartanto, Fitri Afiadi Suhardja Indra Gunawan Devy Kamil Syahbana
Preliminary Results: Determination of Microseismic Early Model of Tomography in Timor using FMTOMO Predictive Modeling of Land Subsidence Due to Preliminary Results on Receiver Function Study in Mt. Observing Near-Real Time Volcanoes Deformation
14.50 - Event Locations on Anisotropic Medium Using Time Nabiilah Rahmawati, Sandy K Suhardja Groundwater Level Decline in Gedebage District, Merapi, Central Java, Indonesia Using InSAR – Time Series of LiCSBAS in
7 Reverse Modeling Bandung, Indonesia Ramadanti Fitria, Suhardja Sandy K. Indonesia
15.10 Syamsurizal Rizal, Awali Priyono, Andri Dian Dwi Sarah, Eko Soebowo, Arifan Jaya Syahbana, Erlangga Ibrahim Fattah, Asep Saepuloh, Andri
Nugraha, Mochamad Apri, Mochamad Agoes Nugroho Aji Satriyo, Zulfahmi, Wahyudin Dian Nugraha, Kristianto, Ahmad Basuki
Moelyadi, Jamhir Safani
Hypocenter relocation of eastern Manggarai Early Modelling of Seismic Velocity Structure P, S, Application of Vertical Electrical Sounding Method Volcanostratigraphy and Identification of Geological Ambient Seismic Noise Cross-correlation of the
15.10 – earthquake 2022 and its aftershocks based on BMKG and Vp/Vs Ratio in Bengkulu Region Using Damped Least-Square Inversion for Structure of the Karaha-Sadakeling Colcano Complex Tangkuban Parahu Volcano. West Java Indonesia:
8 seismic network Nicolas Silaen, Yusrina A Haqi, Ryandika Syahputra, Interpretation of Leachate in the Bakung Landfill, Tasikmalaya Regency, West Java Preliminary Result
15.30 Mohamad Ramdhan, Priyobudi, Adityo Mursityanto, Amira Saqina, Dzikri N Rochim, M F Firmansyah, Lampung Wilfridus FS Banggur, MN Kartadinata, Sofyan Aditya Lesmana, Awali Priyono, Andri Dian
Kadek Hendrawan Palgunadi, Andersen Lasbayakin Belista A Natasya, Riyan Yusami, Aulia Rizkiyana, Rizka, Djoko Santoso, Warsa Primulyana, Estu Kriswati, Azrie Ezziat Putera Nugraha, Zulfakriza, Shindy Rosalia, Ahmad
Panjaitan, Jajat Jatnika Sandy K Suhardja Muhammad, Gradita Trihadi Basuki, Kristianto
10
Zoom 1 – Fault and Deformation Study (1) Session 1
TUE-01 Zoom 1 / SEACG-004
Investigation of Postseismic Deformation of The 2018 Lombok Earthquake Sequence

Using GPS Data
Risya Fauziyyah1, Endra Gunawan2, Sri Widiyantoro2,3, Irwan Meilano4, Syamsuddin5
1
Geophysical Engineering, Faculty of Mining and Petroleum Engineering, Bandung Institute of Technology, Indonesia
2
Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Bandung Institute of Technology, Indonesia
3
Faculty of Engineering, Maranatha Christian University, Indonesia
4
Geodesy and Geomatics Engineering, Faculty of Earth Science and Technology, Bandung Institute of Technology, Indonesia
5
Physics Department, Faculty of Mathematic and Natural Sciences, University of Mataram, Indonesia
*corresponding author’s email :
ABSTRACT
A series of earthquakes struck the northern region of Lombok from July to August 2018 due to the tectonic processes of an east-
west trending fault beneath the northern part of Lombok Island with thrust fault mechanism. The responsible fault of the 2018
Lombok earthquake sequence is a parallel fault structure with the Flores thrust located further south towards Lombok. A
continuous GPS network was installed soon after the 2018 earthquake sequence with GPS stations located in northern Lombok.
This study aims to investigate the postseismic deformation of the 2018 Lombok earthquake sequence using one month of Global
Positioning System (GPS) data. A postseismic deformation strategy analysis was performed to simultaneously quantify the
multiple mechanisms of afterslip and viscoelastic relaxation. In our model, we construct viscoelastic deformation driven by
coseismic stress changes from the latest major Lombok earthquake using Maxwell viscoelastic response. In addition, we obtained
a viscoelastic model of an elastic layer thickness, i.e. 50 km and the Maxwell viscosity is 9.0 × 1017 Pa s. The afterslip model
shows the displacement is relatively close to the afterslip estimated, and relative to the north towards the earthquake source.
Finally, we conclude that afterslip dominates the postseismic deformation mechanism than the viscoelastic relaxation for one
month period after the 2018 Lombok earthquake sequence.
Keywords: 2018 Lombok earthquake sequence, afterslip, viscoelastic relaxation, rheology model, GPS
.....
Study On Estimation Of Fault Plane And Rupture Duration For Potential Assessment
Of Tsunami Hazards In Southern Part Of Sumatera And Southern Part Of West Java
Ajeng Marina Utamie1,2*, Nanang T. Puspito3, Pepen Supendi2
1
Master Program in Earth Science, Faculty of Earth Sciences and Technology, Institut Teknologi Bandung, Indonesia
2
Agency for Meteorology, Climatology, and Geophysics (BMKG), Indonesia
3
Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Indonesia
*corresponding author’s email : ajengmarina7@gmail.com
ABSTRACT
The southern part of Sumatra and the southern part of West Java have a high potential for earthquakes and tsunamis. Several
earthquakes and tsunamis in this area have caused many casualties. We use the earthquake catalog data sourced from the
Indonesian Agency for Meteorology, Climatology, and Geophysics (BMKG) between 2009 and 2019 to estimate the fault plane
and duration of fault rupture as an effort to understand the fault mechanism that causes tsunami-generating earthquakes in these
areas. Our results show that there are differences in fault plane patterns and mechanisms in earthquakes in the southern part of
Sumatra and the southern part of West Java. The October 25, 2010 (Mw 7.1) earthquake in Mentawai at 01:11 UTC showed a
reverse fault plane pattern and had the largest estimated duration of fault rupture at 182.9 s, generating a tsunami that caused
fatalities while on April 3, 2011 (Mw 6.6) earthquake in southwest Cilacap at 20:06 UTC showed a normal fault plane pattern,
had estimated duration of fault rupture at 61.2 s and did not triggered tsunami. The results of this study are expected to be
considered in efforts to mitigate earthquakes and tsunamis in the southern part of Sumatra and the southern part of West Java.
11
Zoom 1 – Fault and Deformation Study (1) Session 1
Keywords:
Moment Tensor Inversion Implementation In Determining Focal Mechanism Solution

Of Palu-Koro And Matano Fault Events: Processing Strategy And Preliminary Result
Dian Kusumawati1,2*, David P. Sahara1, Nanang T. Puspito1, Mudrik R. Daryono3, Andri D.
Nugraha1, Mohamad Ramdhan3, Andri Kurniawan2, Wahyu Tanihaha2, Leonardo Junior Johan
Solihin2, Muhamad Diva Pratama2, Sindi Hajah Patimah2, Annisa Trisnia Sasmi2, Hilmy
Muhammad2, Hammam Muhammad2, Widjra Cyiena Christi Natafrisca2, Maulidia Ain Bening2
1
Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Bandung, Indonesia
2
Geophysical Engineering, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Bandung, Indonesia
3
Research Center for Geological Disaster, National Research and Innovation Agency (BRIN), Bandung, Indonesia
*corresponding author’s email : diann.kusumawati03@gmail.com
ABSTRACT
Earthquake focal mechanisms are useful in analysing seismotectonic features in certain area, as they can depict subsurface
structure. Earthquake focal mechanism solution can be resulted from first polarity or moment tensor inversion method. The
moment tensor inversion can produce more complete information of focal mechanism solution than the first polarity method.
However, the method is slightly more complicated for it is based on waveform inversion. Several problems related to signal
processing might arise. Strategies should be applied to overcome the difficulties and obtain reliable focal mechanism solution. In
this study we applied moment tensor inversion to produce focal mechanisms data of Palu-Koro and Matano Fault events. We
processed focal mechanism of events with depth less than 60 km and with minimum magnitude of Mw 3.8. The moment tensor
inversion is conducted using Isolated Asperities (ISOLA) software. Full waveform inversion of stations located in the vicinity
of Palu-Koro and Matano Fault is applied. The cause of waveform fitting difficulties using the data in this study area and the
solution are discussed. For instance, the existence of unseen signal noise and clear-yet-disturbing signal are observed. We
noticed that the waveform fitting is easier to be applied in the Matano Fault than in the Palu-Koro Fault. Moreover, the overall
resulting focal mechanisms in Palu-Koro and Matano fault is presented.
Keywords: focal mechanism, Matano Fault, moment tensor inversion, Palu-Koro Fault
.....
Geoelectric Analysis for Potential Slip Plane Interpretation in Labuan Bulan Hill,
Salem District, Brebes Regency, Central Java
Indra A. Dinata1, M. Khalifah B.Y. Praja2, Imam A. Sadisun1,3*, Agus Kuswanto4
1
Applied Geology Research Group, Faculty of Earth Sciences and Technology, Institut Teknologi Bandung, Indonesia
2
Undergraduate Student of Geological Engineering Study Program, Faculty of Earth Sciences and Technology, Institut Teknologi
Bandung, Indonesia
3
Research Centre for Disaster Mitigation, Institut Teknologi Bandung, Indonesia
4
Geoscan Ltd.
*corresponding author’s email : indraandradinata@gmail.com
ABSTRACT
Landslide is a natural disaster which often occurs in Indonesia. Based on National Disaster Management Agency (BNPB),
there were 1,321 landslide incidents which occurred in 2021. Labuan Bulan Hill located on high susceptibility landslide seen
from Landslide Susceptibility Zone of Brebes Regency, furthermore, soil tensile cracks presence with 5 - 50 cm width and 15
- 125 cm height at upper slope. In this study, the identification of subsurface lithologies and potential slip plane was carried
out using geoelectric survey method with Wenner – Schlumberger configuration. The Wenner – Schlumberger configuration
has deeper penetration than other geoelectric survey configurations. Geoelectric surveys are carried out on two lines with 240
m length and 5 m electrode space. Geoelectric line is made parallel and perpendicular with slope direction to estimate
subsurface conditions easier with the help of tensile cracks presence in the surface. Results of measurements identified that the
upper layer consists of poorly graded silty sand with 285 Ωm – 2224 Ωm resistivity, the second layer is a moderate - highly
weathered volcanic breccia with 58.7 Ωm – 285 Ωm resistivity, and the third layer is fresh – slightly weathered volcanic breccia
with 1 Ωm – 1388 Ωm resistivity. A fresh – slightly weathered volcanic breccia layer was identified as potential slip plane with
20 – 30 m depth.
Keywords: landslide, geoelectric survey, potential slip plane
12
Zoom 2 – Passive Seismic (1) Session 1

Data Processing of Surface Wave Tomography in Australia using Automated Surface
Wave Phase-Velocity Measuring System (ASWMS)
Amira Saqina1*, Nicolas Silaen1, Riyan Yusami1, Yusrina Alfiani Haqi1, Sandy K. Suhardja1, 2
1
Faculty of Exploration and Production Technology, Pertamina University, Jl. Teuku Nyak Arief, Simprug, Kebayoran Lama,
Jakarta 12220, Indonesia
2
Agency for Meteorology, Climatology and Geophysics, Jalan Angkasa I, No. 2, Kemayoran, Indonesia
*corresponding author’s email : amirasaqina1407@gmail.com
ABSTRACT
The Australian continent is divided into tectonic blocks, including Archaean and Proterozoic units in the western and central
part, and Phanerozoic units in the eastern part. Surface wave tomography is an effective method of getting information about
velocity variations in a particular region. Surface wave tomography was used in this project to map the subsurface of the
Australian continent down to a depth of 100 km with the ASWMS (Automated Surface-Wave Phase-Velocity Measuring
System). From 2010 to 2015, we gathered waveform data from 767 events and 219 permanent and temporal seismic stations
dispersed across Australia. We selected a set of events that has a minimum magnitude of 6. Rayleigh wave extraction and
cross-correlation between two nearby stations are among the first processing processes. Following that, the apparent phase-
velocity value is calculated using the Eikonal and tomographic inversion. To generate a tomographic map with a true structural
phase velocity, each phase velocity data is stacked with amplitude correction. According to current tomographic imaging
studies, low-velocity values under Australia tend to be dispersed easterly, with Rayleigh velocity ranging from 3.7 to 3.9 km/s.
The tectonic history and evolution of eastern Australia are related to a sequence of orogenic events that are pushing toward the
continent's eastern boundary.
Keywords: Cross-Correlation, Phase-velocity, Surface-wave tomography
.....

Processing Receiver Function Data in Mentawai Region, Sumatra
M.Aldi1*, Nicolas Silaen1, Dayyan Dhaifullah1, Sandy K Suhardja1,2

1
Faculty of Exploration and Production Technology, University of Pertamina, Jl. Teuku Nyak Arief, Simprug, Kebayoran
Lama, Jakarta 12220, Indonesia
2
*corresponding author’s email : aldi.mhd07@gmail.com
ABSTRACT
The Mentawai Islands are located in the western part of Sumatra, in the Indo-Australian and Eurasian plate subduction zone.
Because of these tectonic settings, the Mentawai Region has a high level of seismicity, making it particularly appealing for
subsurface structure study. We used the receiver function method in this study to determine the depth of the oceanic crust and
identify the low-velocity zone (LVZ) under the Mentawai Islands. This method uses the conversion of P-to-S waves, which
have a larger amplitude in the radial component than vertical waves. The deconvolution technique is used to extract the two
components of the receiver function signal. Water level deconvolution and iterative deconvolution were employed in this study.
In practice, any deconvolution method cannot produce a suitable receiver function due to various methodologies that cannot
be used in all conditions. To get the best results, evaluate many deconvolution methods. This study uses two teleseismic stations
situated in the north and south of the Mentawai, which was obtained through the IRIS website. The Butterworth filter is used
before the data is rotated to get ordered and noise-free data. The first results show that the water level is stronger than the
iterative, which is more sensitive to the Gaussian filter width parameter. However, they produced almost identical results at
both stations, precisely the arrival time of Ps, which was in 3-4 seconds and was followed by a low velocity.
Keywords: Deconvolution, Receiver Function, Parameter Receiver Function, Iterative Deconvolution, Water Lever Deconvolution
.....
13
Zoom 2 – Passive Seismic (1) Session 1

Preliminary Results on Low Frequency Passive Seismic (LFPS) Study in Lake Toba
Riyan Yusami1*, Sandy K Suhardja1, Dzikri N Rochim1, Nicolas Silaen1, Amira Saqina1
1 Faculty of Exploration and Production Technology, University of Pertamina, Jl. Teuku Nyak Arief, Simprug, Kebayoran
*corresponding author’s email : riyanyusami03@gmail.com
ABSTRACT
Sumatra Island is a subduction zone island located between the Indo-Australian plate and the Eurasian plate. The Sumatran
Fault, also known as the Semangko Fault, is a geological structure. The presence of the Bukit Barisan Mountains, which are
characterized by active volcanic activity, makes Sumatra Island vulnerable to volcanic eruptions. This study was carried out
in Lake Toba, North Sumatra, using earthquake data from the IRIS network of recording stations. The goal of this study was
to use the Low Frequency Passive Seismic (LFPS) method to investigate the amplitude of knowing the presence of fluid
reservoirs and volcanic analysis in the Lake Toba region.The Vertical to Horizontal Spectral Ratio (VHSR) and Power
Spectral Density (PSD) methods were used to locate and connect reservoirs, as well as to perform temporal analysis to
observe changes in amplitude in the time domain. Because of the anomalous spectrum of the maximum value of V/H and a
relatively high maximum value of the PSD amplitude, the results obtained indicate that the research locations that have the
potential for the existence of trapped subsurface fluid reservoirs through the VHSR and PSD methods are along the
Southwest–Northeast region, Samosir Island, and several areas around Lake Toba. As a result of the VHSR and PSD
methods, it can be assumed that volcanic activity in the Lake Toba region has been relatively active and dynamic up to this
point, which supports previous research.
Keywords: North Sumatra, Lake Toba, Fluid Reservoir, Vertical to Horizontal Spectral Ratio, Power Spectral Density, Low Frequency
Passive Seismic
.....
Sediment thickness in the eastern part of Kalimantan from high-frequency receiver

functions
Syuhada Syuhada1, Bayu Pranata2, Titi Anggono1, Mohamad Ramdhan1, Febty Febriani1, Aditya D
Prasetio1, Cinantya N Dewi1, Mohammad Hasib1
1
Research Center for Geological Disaster, National Research and Innovation Agency (BRIN), Cisitu, Bandung, Indonesia
2
Meteorology, Climatology and Geophysical Agency (BMKG), Kemayoran, Jakarta, Indonesia
*corresponding author’s email : hadda9@gmail.com
ABSTRACT
We analyze high-frequency receiver functions from two seismic stations to characterize the sediment properties around the
eastern coast of Kalimantan. First, we compute high-frequency receiver functions from teleseismic events using the time-domain
deconvolution method. The obtained receiver functions are then stacked using the H-κ stacking approach to determine the
sediment properties. The computed receiver functions show the complicated waveforms for the seismic station located near
Kutai Basin. These waveforms may indicate the presence of multi sedimentary layers. The H-κ stacking results suggest that the
sediment is thinner in the northern part of the area (~2.40 km depth) than in the southern part of the region (~7.38 km depth).
These results are consistent with the previous geological and geophysical studies conducted in this region. The previous works
suggested that the sediment thickness around the study area can reach up to 9 km. Furthermore, the information about sediment
properties around this area may be useful for the accurate characterization of crustal properties using receiver function analysis.
Keywords: Receiver function, Kalimantan, Sediment, H-κ stacking
14
Zoom 3 – Environmental Geophysics (1) Session 1
Detection of Old Tram Tracks Near South Part of Surabaya Zoo Using Georadar
Method
Rayhan Farisi Ramadhan1*, Mariyanto Mariyanto1, JPGN Rochman1
1
Department of Geophysical Engineering, Faculty of Civil, Planning and Geo Engineering, Institut Teknologi Sepuluh Nopember,
Surabaya 60111, Indonesia
*corresponding author’s email : rayhanfr1@gmail.com
ABSTRACT
The development of public facilities and infrastructure in the city of Surabaya is carried out to improve public services and city
management. One aspect that is currently being focused on is the development of public transportation system with the
government's plan to reactivate the mass transportation system in the Dutch era in the form of trams. In order to reactivate the
series of tram tracks, detection of old tram tracks that have been buried using the Ground Penetrating Radar or Georadar method
was carried out. This very near surface exploration method utilizes electromagnetic wave propagation which is sensitive to the
high dielectric permittivity contrast of materials such as metals. In this study, GPR measurement using GPR GSSI Antenna 270
MHz around Terminal Joyoboyo Surabaya. GPR data processing techniques such as Adjust Signal Position, Remove DC,
Dewow, Inverse Amplitude Decay, Remove Global Background, Karhunen-Loeve Filter and Time to Depth Conversion are
conducted to facilitate the identification and interpretation process. In addition, validation was also carried out with a map of
the old Dutch tram tracks belonging to PT. KAI. From the analysis, it was found that the presence of old tram tracks was detected
1 meter beneath the JB07 and JB14 and also corresponding to the map of the old tram tracks belonging to PT. KAI.
Keywords: Georadar, Old Tram Tracks, Dielectric Permittivity

.....
Steam-Vent Identification Using Frequency Charasteristic Method In Ground

Penetrating Radar (GPR)
Radithya A. Pambudi1*, Warsa1, Wahyudi W. Parnadi1
1
Department of Geophysical Engineering, Institut Teknologi Bandung, Indonesia
*corresponding author’s email : rarifpambudi@gmail.com
ABSTRACT
Ground Penetrating Radar (GPR) is a geophysical method that widely used and developed in a lot of fields because of the non-
destructiveness and relatively high resolution of the method. Despite all that this method can produce a biased result because
of factors that could affect the anomaly. To prevent those advanced analytical methods are used, one of them is frequency
characteristic. Frequency characteristic method used GPR’s trace to then analyze it in frequency and time domain. The result
is quantified in various characters of frequency. This method is used to help identify steam-vents in geothermal environment
in Darajat, Garut. Finite difference time domain method (FDTD) is used to generate a synthetic model to then analyzed and
compared with the field data to identify the property of steam-vent. From the result we know the applicability of the method
to identify anomalies.
Keywords: Ground Penetrating Radar, frequency characteristic, steam vent, FDTD
15
Detecting Subsurface Utilities with Multi-frequency and Multi-antenna Configuration

GPR
Wahyudi W Parnadi1*, Indra Gunawan1*, Warsa1, Agus Laesanpura1, Djoko Santoso1
1
Geophysical Engineering Department, Institut Teknologi Bandung, Jalan Ganesha 10, Bandung 40132, Indonesia
*corresponding author’s email: wahyudi@gf.itb.ac.id
ABSTRACT
Ground-Penetrating Radar (GPR) technique is a powerful method to image subsurface utilities using its high-frequency
electromagnetic waves. Higher frequency signals deliver high-resolution detection tools, whereas lower frequency signals offer
deeper investigation depth. Combining radargram from both signals into one composite radargram will provide an effective
radargram for interpretation. In addition, applying acquisition techniques using different antenna configurations will provide
more complete data for interpretation. In this research, we conducted experiments with different antenna orientations:
perpendicular broadside, parallel end-fire, and cross-polarization using two different working frequency antennas. We
implemented these experiments using the Ramac Mala GPR system with 100 MHz and 200 MHz unshielded antennae in the
PPSDM Geominerba's Test Site, in which we buried various pipes of different sizes, different depths, and material types in the
ground. Our study showed that acquisition with two different antenna configurations in the same lines is needed to yield a better
subsurface utility interpretation. Combining two different radargrams acquired with 100 MHz and 200 MHz antennae enhances
its interpretability.
Keywords: GPR, multi-frequency, multi-configuration antenna, PPSDM Geominerba’s Test Site
....
GPR Impedance Inversion for Imaging Groundwater Level: a Case Study at PPSDM
Test Site, Cipatat, West Java Province
Wahyudi W Parnadi1*
1
Geophysical Engineering Department, Institut Teknologi Bandung, Jalan Ganesha 10, Bandung 40132, Indonesia
*corresponding author’s email : wahyudi@gf.itb.ac.id
ABSTRACT
Ground-Penetrating Radar (GPR) method is an appropriate tool for imaging the groundwater level. A conventional technique to
interpret processed GPR data is interpreting its radargram in wiggle format. That technique provides only geometry information
like groundwater level. In this research, we carried out GPR data analysis using impedance inversion techniques to know the
distribution of physical properties, leading to a better subsurface interpretation. We performed some GPR survey lines at the
PPSDM test site using the GPR Mala system with a 100 antenna. Then we processed the acquired data with conventional steps
like dewow, bandpass filtering, and deconvolution. The next step was to conduct GPR data inversion with the aid of secondary
data. The results provide clear images of the subsurface electromagnetic properties, including the groundwater level.
Keywords: Impedance inversion, Ground-Penetrating Radar, GPR, groundwater level
16
Zoom 4 – Geophysical approaches in hydrocarbon (1) Session 1
Saturated Fluid Variation Effect on Seismic Anisotropy: A Case Study of Deep-Water

Reservoir
Madaniya Oktariena1*, Wahyu Triyoso2, Fatkhan Fatkhan2, Sigit Sukmono2, Dona Sita Ambarsari2,
Erlangga Septama3, Pongga D. Wardaya3, Befriko S. Murdianto4
1
Graduate Program of Geophysical Engineering, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Jalan
Ganesha 10, Bandung 40132, Indonesia
2
Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Jalan Ganesha 10, Bandung 40132, Indonesia
3
PT Pertamina (Persero) Upstream Research & Technology Innovation, SopoDel Office Tower 50th Floor, Jalan Mega Kuningan Barat
III Jakarta 12950, Indonesia
4
Pertamina Hulu Kalimantan Timur, CIBIS NINE Tower 9th Floor, Jalan TB Simatupang 2, Jakarta Selatan 12560, Indonesia
*corresponding author’s email : madaniyaoktariena@outlook.com
ABSTRACT
Many factors can drive seismic anisotropy, whether it is the intrinsic condition of rock, external force disturbing the initial
elastic condition of the rock, or longwave anisotropy. Saturated fluid, as one of the reservoir's intrinsic elements, is known to
trigger an anisotropic effect in far offset seismic with Vertical Transverse Isotropic character. However, seismic anisotropy
is not a popular tool to be used as an attribute in reservoir characterization, and it is with good reasoning. More often than
not, observed anisotropy (η parameter) at a certain depth of reservoir level is an accumulation of various geological and
geophysical conditions. Anisotropy on seismic data presents as hockey stick curving that a combination of two-parameter
can correct during move-out: VNMO and the elliptical anisotropy η parameter. This research means observing how the
saturated fluid within the reservoir contributes to the observed anisotropy in seismic data. Rock physics modeling (RPM) and
fluid substitution (FRM) are conducted to investigate the saturated fluid effect at a particular temperature and pressure that
reflects the deep-water environment. First, a synthetic seismogram is produced by utilizing VTI Rüger Reflectivity on the
modified log and the mother wavelet from the existing 3D seismic survey. Then, the anisotropy parameter is extracted by
applying Tsvankin 4th Order NMO to obtain the elliptical anisotropy η parameter. Based on synthetic seismic modeling with
a modified log from the FRM process, saturated fluid has an effect on far offset data as a manifestation of seismic anisotropy
η. In this case of a deep-water environment, gas saturated sand needs the most considerable η correction, followed by oil and
water-saturated. Although the seismic anisotropy η value difference is relatively small and easily dismissed as insignificant,
saturated fluid as part of intrinsic rock element indeed does contribute to the observed seismic anisotropy.
Keywords: Fluid Replacement Modelling, Eta, Vertical Transverse Isotropy, Seismic Anisotropy, Deep-Water Reservoir
....
Fault Delineation and Geometrical Modeling of the Geological Carbon

Storage (GCS) Candidate for the CCS and CCUS Facilities Based On
Gravity Dataset. A Case Study at the Northtern Part of East Java, Indonesia
Accep Handyarso1*
1
Research Center for Geological Resources, National Research and Innovation Agency Jl. Sangkuriang, Bandung 40135, Indonesia
*corresponding author’s email: accep.handyarso@brin.go.id
ABSTRACT
Global warming due to the effects of greenhouse gas emissions has triggered climate change and caused various geological
and hydro-meteorological disasters in the world. Most countries have carried out decarbonization efforts following the
Paris Agreement, which aims to reduce carbon emissions from 29% up to 41% by 2030. Sustainable Development Goals
no. 13 which states combating the climate change is globally accepted framework. One element of the successful CCS and
CCUS facilities is the geological carbon storage (GCS) readiness. The GCS locates at a specific place underneath, it has
limited capacity, and the existing fault around the area must be considered to avoid leakage and uncontrollable behavior
after injection phase of the liquefied (high-pressure) carbon. This research exposes the existing fault and estimates the
geometry of the GCS candidate in the northern part of East Java, Indonesia. A synthetically strike-slip fault modeling was
made to help gravity data interpretation and the 3D gravity data inversion is conducted to estimate the storage capacity
which is constrained with the thickness of the targeted formation in the area as the prior information. According to the
research, the proposed GCS area has large-scale anticline closure without any strike-slip fault detected. The subsurface
3D modeling has been made to visualize the existing fault system which belongs to the Baribis-Kendeng Fault Zone and
calculating the potential capacity estimation of the GCS area. The study area proposed to be a potential GCS area with
low-leakage potential besides Central Processing Plant Gundih and Sukowati in East Java, Indonesia.
Keywords: Climate-Change, Decarbonization, Geological-Carbon-Storage, Fault, Capacity, Gravity
17
The Effects of Porosity on The Anisotropy Parameters of The Slope Fan Facies Sand
Reservoirs in The Deepwater Kutai Basin, East Kalimantan, Indonesia
Rahma Nurhazlita1*, Teuku Abdullah Sanny1, Dona Sita Ambarsari1
1
Geophysical Engineering, Faculty of Mining and Petroleum Engineering, Institute of Technology Bandung, Indonesia
*corresponding author’s email : rahmanurhazlita00@gmail.com
ABSTRACT
Knowledge about the effects of anisotropy and its relationship with reservoir properties is needed for good interpretation and
modeling. Related to this issue, this study evaluated the effects of porosity on the anisotropy parameters for sand reservoirs
deposited in the upper slope fan facies. The study objects are obtained from core plug samples of sand reservoirs in the deepwater
Kutai Basin. Core plug samples were collected from the depth of 13602.5 ft to 13704.5 MD. Thomsen's ε, γ, δ, and η anisotropy
parameters data were obtained by ultrasonic measurements and porosity data were obtained by laboratory measurement with
Coreval 700 apparatus and Boyle’s law . The analysis results show that the relationship of the anisotropy parameters with
porosity appears when high porosity sandstone and low porosity sandstone are separated. The plots of anisotropy ε, γ, and δ,
show trends that for greywacke, with increasing anisotropy value the porosity increases. The effect of porosity on the high
porosity (29%-37%) sandstone shows a steeper change than on the lower porosity (12%-13%) sandstone. The analysis also
shows that the higher composition of lithic mineral grain reduces the effect of anisotropy on porosity.
Keywords: Anisotropy, Core plug, Porosity, Reservoir, Sandstone
....
18
Zoom 5 – Geophysical approaches in geothermal Session 1

Characteristics of Subsurface Conditions of The Selogajah Hotsprings Area Using Very-
Low Frequency Electromagnetic (VLF-EM) Method
Rahmi Azizah1*, Mariyanto Mariyanto1, M. Singgih Purwanto1
1
Department of Geophysical Engineering, Faculty of Civil, Planning, and Geo Engineering
Institut Teknologi Sepuluh Nopember, Surabaya 60111, Indonesia
*corresponding author’s email : rahmiazizah05@gmail.com
ABSTRACT
Selogajah Hill is located in the North Slope area of Mount Pandan which has geothermal manifestations in the form of hot
springs. The area is located in the Kendeng zone which has several active faults and other geothermal manifestations. Based on
these manifestations, VLF-EM was used to identify the characteristics of subsurface conditions in the Selogajah Hill area,
Gondang District, Bojonegoro Regency. So far, there has been no intensive study conducted to map the characteristics of the
subsurface conditions around the Selogajah hot spring so that the VLF-EM method is expected to be a reference and assist local
communities in planning development and regional development in the area. Data acquisition was carried out on 6 lines in the
Selogajah hot spring area and VLF data processing was carried out with several filters including the NA-MEMD filter to
decompose data and reduce noise. From the qualitative processing results from the Karous-Hjelt filter, the conductivity
distribution with apparent depth is obtained. Where from the 6 parallel lines, there is a conductive anomaly which is
concentrated in the eastern part of the study area, in the middle area of the track and slightly appears in the western part. The
interpretation of the conductive anomaly is that the rock layer is saturated with water in the form of hydrothermal fluid and the
resistive anomaly in the area indicates the presence of rock fractures that gave rise to the Selogajah 1 and Selogajah 2 hot
springs.
Keywords: VLF-EM, NA-MEMD, Selogajah
.....

3D Subsurface Characterization using Gravity Data of Mataloko Geothermal Field,
Flores, Nusa Tenggara Timur, Indonesia
Fadhil Farras1*, Susanti Alawiyah1, Setianingsih1, Eko Januari Wahyudi1
1
Department of Geophysical Engineering, Faculty of Mining and Petroleum Engineering (FTTM), Bandung Institute of Technology,
Bandung, Indonesia
*corresponding author’s email : farrasfadhil5@gmail.com
ABSTRACT
The Mataloko geothermal field is located in the central part of Flores, there are many alteration zones including Nage centered
around the research area.There is a normal Wae Luja fault which causes the formation of a reservoir zone in the claycap zone.
This fault causes the formation of fractures in this zone and as a channel for fluid entry from the reservoir to the shallow
reservoir.The conceptual model begins with interpreting geoscientific data (geology, geochemistry, geophysics). The final
result of the modeling is expected to be 3D with the aim of clarifying the visualization of the subsurface model of the Mataloko
geothermal field.
Keywords: geothermal, Mataloko, fault, conceptual model
19
Zoom 5 – Geophysical approaches in geothermal Session 1

Identification of Shallow Fault with Seismic Refraction Modeling in Gondang Area,
Bojonegoro
Nahari Rasif1*, Aditya Aulia Rachman1, Mariyanto Mariyanto1
1
*corresponding author’s email : naharirasif18@gmail.com
ABSTRACT
There are various considerations in determining whether a field is said to have geothermal potential or not. One of them is the
appearance of shallow faults in the area. This shallow fault can function as a permeable layer that is used as a fluid flow path in
the geothermal system. Due to the emergence of manifestations around Mount Pandan in the Gondang area, Bojonegoro, seismic
refraction measurements are also needed that can support and validate the geothermal potential in the area with the hope of
shallow faults that are clearly visible through the interpretation of seismic wave propagation velocity. This measurement
produces an anomaly which is thought to be a shallow fault on the path that passes through the main fault, namely on paths 1,
2 and 5, not far from the manifestation of Banyukuning and Selogajah hot springs. This is also reinforced by direct observation
of the surface conditions as a result of the appearance of shallow faults.
Keywords: Seismic Refraction, Shallow Fault, Modelling
.....

Tomographic Seismic Model Using Ray Tracing Algorithm Dijkstra To Visualize
Structure In The Geological Model Of The Wayang Windu Geothermal Field, West
Java
Muhammad Rizqi Fadhlurrahman1*, Indriati Retno Palupi1
1
Jurusan Teknik Geofisika, Fakultas Teknologi Mineral, Universitas Pembangunan Nasional “Veteran” Yogyakarta
*corresponding author’s email : rizqi.fadhlurrahman@gmail.com
ABSTRACT
The seismic tomography method is one of the techniques used to image the subsurface conditions of the earth, in this case the
structure using earthquake waves that propagate. The tomographic seismic method in the Wayang Windu geothermal field uses
the parameter of S wave propagation (Vs). The processing carried out to display the image of the structure here uses the ray
tracing Dijkstra’s algorithm where this algorithm is used to find the shortest path, namely with the minimum weight in this case
from the hypocenter to the station. The methodology for data processing consists of several stages, namely, parameterization of
the model that will be made into an initial velocity model where the initial velocity model will determine the placement and
number of hypocenters and stations, for the placement of hypocenters placed in structural zones, followed by making synthetic
travel time. where at this stage the weighting is carried out using the Dijkstra’s algorithm, followed by performing Inverse
Modeling to display seismic tomography at the Wayang Windu geothermal field. In this initial speed model, it has a grid
measuring 30x50 by entering the initial speed based on the AK135 Table Time Travel with a depth of 1-30 km and has a speed
value of 2.5-8.05 km/s. By using the Dijkstra’s algorithm, it will quickly solve the ray tracing problem in the model created
where this algorithm will determine the blocks that will be passed by the wave based on the determination of weights either
from the addition of noise or other influences. This weighting will vary from one hypocenter to another due to the different
dimensions of each hypocenter to the station. Based on the processing carried out using the MatLab and Python programs, the
results of seismic tomography are shown where the structural pattern is quite clear. Where is this structural pattern related to
the number of hypocenters and by placing the hypocenter in the structural zone which will clarify the structure image on this
seismic tomography.
Keywords: tomography seismic, Dijkstra’s algorithm, hypocenter
20
Zoom 1 – Fault and deformation study (2) Session 2
Late Quaternary Deformation of the Ambon fault, Indonesia: Preliminary Results
Adi Patria1,2*, Danny Hilman Natawidjaja1, Muhammad Hanif1, Hiroyuki Tsutsumi3

1
2
Graduate School of Science and Engineering, Doshisha University, Kyotanabe, Japan
3
Department of Environmental Systems Science, Doshisha University, Kyotanabe, Japan
*corresponding author’s email: adip006@brin.go.id
ABSTRACT
The Banda Arc in eastern Indonesia is a tectonically complex region situated at the convergence zone between the Australian,
Eurasian, and Pacific plates. Our previous active fault investigation in the northern Banda Arc has revealed that the Ambon
fault can generate Mw 6.9 earthquakes. This capability poses a significant risk to Ambon Island because it is located near
populous areas. Recent seismicity in Ambon Island and surrounding areas indicates that this region is actively deforming.
Ambon Island has experienced many damaging earthquakes since the 17th century. However, the potential seismic hazard on
this island remains enigmatic due to a lack of detailed geological information, such as fault traces, slip rate, and paleoseismic
history. This study reports the preliminary result of tectonic geomorphic investigations and ground-penetrating radar (GPR)
surveys on the Ambon fault. Tectonic geomorphic features, such as fault scarps, warped surfaces, and sag ponds pronounce
the late Quaternary activity of the Ambon fault. The Ambon fault predominantly indicates a normal sense of motion with a dip
of ~60°. The Ambon fault is interpreted as an antithetic fault of the Banda Detachment in Weber Deep, accommodating the
extensional deformation in the Banda Sea.
Keywords: Active fault, Ambon fault, Ground-penetrating radar (GPR), Seismic hazard, Tectonic geomorphology
.....
Neotectonics of the eastern Matano fault in Sulawesi, Indonesia

Adi Patria1,2,*, Danny Hilman Natawidjaja1, Mudrik Rahmawan Daryono1, Muhammad Hanif1,
Anggraini Rizkita Puji1, Hiroyuki Tsutsumi3
1
2
Graduate School of Science and Engineering, Doshisha University, Kyotanabe, Japan
2
Department of Environmental Systems Science, Doshisha University, Kyotanabe, Japan
*corresponding author’s email : adip006@brin.go.id
ABSTRACT
Eastern Indonesia is a seismically active and rapidly deforming region due to the convergence of the Australian, Eurasian,
Philippine Sea, and Pacific plates. The Matano fault in Sulawesi is a part of a fast-slipping left-lateral strike-slip fault system,
extending from Papua to Sulawesi in eastern Indonesia. The Matano fault is the southeastern continuation of the Palu-Koro
fault which slipped during the 2018 Mw 7.5 earthquake. It is capable of generating large (Mw ≥7) earthquakes. However, the
Matano fault is less studied, unlike the Palu-Koro fault. This study aims to uncover the neotectonics of the eastern portion of
the Matano fault which may help in understanding the seismic hazard posed by the fault. Based on our investigation, the eastern
Matano fault slips slower than the Palu-Koro fault, corresponding to the fact that it has a longer recurrence interval of surface-
rupturing earthquakes. This slip-rate difference plausibly causes extensional deformation in the area north and east of the
Matano and Palu-Koro faults, respectively. The Matano fault may form two seismic gaps that potentially host the next large
earthquakes in the near future.
Keywords: Active fault, Earthquake geology, Matano fault, Neotectonics, Seismic hazard
21

Mataram Fault Section – New Active Fault Section Crosses East-West in The Centre of
Yogyakarta City
Mudrik R. Daryono1, Astyka Pamumpuni2, Dadan D. Wardhana1, Danny H. Natawidjaja1

1
Research Centre for Geological Disaster, National Research and Innovation Agency (BRIN)
2
Engineering Geology Department, Bandung Institute Technology (ITB)
*corresponding author’s email : mudr001@brin.go.id
ABSTRACT
Yogyakarta city has experienced two devastating earthquake disasters of the recent 2006 Mw6.4 (?) and the penultimate event
of the June 10, 1867 Mw7.7. The active fault that is thought to be responsible for the two earthquake events is the Opak Fault.
However, the Opak Fault has still not thoroughly well mapped yet. The lack of high-resolution image, dense vegetation cover,
high sediment flux from Merapi volcano, and human activities eroding the original landscape are challenges to study the on-
land fault in this area. Our recent study, however, indicate that the Opak Fault is not the only active fault that can cause major
disaster, but there is another fault strand in the area. We mapped the fault using the best available data of DEMNAS assisted
by more detailed DEM and Orthophoto, developed from drone survey, and we also conducted Earth Resistivity Tomography
(ERT) survey. We found the inferred new active fault is oriented East-West crosses the center of the high-populated city. The
fault runs parallel and close to the famous Mataram irrigation system (Selokan Mataram). The fault is likely to be the
continuation of previously recognized Dengkeng Fault, east of the Opak fault. The fault strand is indicated by morphological
lineaments and a few steam offsets. Our ERT 2D sections have been able to reveal the fault zone in several locations along
the inferred fault line. The fault line cut the Young Quaternary rocks, and hence it is an active fault. Further studies are needed
to get further details of this newly recognized active fault, such as conducting paleoseismological studies, detailed
seismological survey, geodetic GPS measurement, and acquiring high-resolution image acquisition survey (LiDAR – Light
Detection and Ranging Survey).
Keywords: Active Fault, Yogyakarta, Mataram Fault
.....
Geological slip rate study on Renun Segment of Sumatran Fault in Dolok Sanggul area,
North Sumatra using UAV photogrammetry and DEMNAS
Mudrik R. Daryono1, Muhammad Hanif1, Dadan Dany Wardhana1, Danny Hilman
Natawidjaja1
1
*corresponding author’s email : mudr001@brin.go.id
ABSTRACT
Renun Segment is located in the center of the Sumatran Fault, parallel to the Famous Toba Lake Volcano. The fault cuts the
thick and vast Young Toba Tuffs (YTT) that was deposited ~74,000 years ago during the super volcano eruptions. It shows a
series of a noticeable dextral strike-slip offsets of the river valleys that deeply incised the tuffs. Our study is to determine the
geological slip rate in this location by measuring river-valley offsets in the Dolok Sanggul area using UAV photogrammetry
and DEMNAS images. Detail morphology shows a fault lineament, monocline, and an anticline structure. Two rivers about a
hundred meters deep cross perpendicular to the fault line have 585±83 meters dextral offset. This amount is only a half of the
other dextral offsets on the northern and southern parts of the site. Hence, the slip rate is also a half slower, which is
approximately 7±1 mm/yr. This indicates that that the other half of the slip budget is accommodated by Samosir Fault in the
center of Toba Lake.
Keywords: Renun Segment, Sumatran Fault, River Offset, Slip rate
22
Correlation of Stress reconstruction, Correlation Dimension (Dc), and Seismic

Moment Rate Loading along the Sumatra subduction margin
Wahyu Triyoso*1, Dina A. Sarsito2 and Sigit Sukmono1
1
Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Jalan Ganesha 10, Bandung 40132, INDONESIA
2
Faculty of Earth Sciences and Technology ITB, Institut Teknologi Bandung, Jalan Ganesha 10, Bandung 40132, INDONESIA
*corresponding author’s email : wtriyoso@gmail.com
ABSTRACT
This study intends to find the possible correlation of the previous stress reconstruction result with the correlation dimension
(DC) and its correlation with the present-day strain loading along the Sumatra subduction margin. The purpose is to understand
the possible earthquake potential hazard that might be happened in the future. It is motivated based on the suggestion of the
finding based on our previous result. The areas with relatively high DC coincide with high seismic moment loading rates.
Therefore, it implies high tectonic stress loading that could pose the risk of producing significant earthquake hazards in the
future. The main reason causes the use of DC is that the Frequency Magnitude Distribution (FMD) parameter is affected by
several factors, and thus, to ﬁnd mechanisms that generate the stress perturbations are difﬁcult. By referring to our previous
study of the two seismotectonic parameters, the relationships of (DC = 2.851 - 1.272b) are used. The result shows that we can
characterize the possibility of earthquake potential hazards existence or seismic gaps existence. It is consistent with our
previous outcome based on the spatial correlation of maximum shear strain and correlation dimension. This study could better
understand the possible seismic hazard in the future since we can find the possible correlation between the high DC with the
possible present-day strain loading since we incorporate the pre-sent-day shallow crustal dynamic data.
Keywords: Stress Reconstruction, Correlation Dimension, Seismic Moment Rate Loading, Possible Future Earthquake Potential Hazard.
.....
23
Zoom 2 – Advance geophysics (1) Session 2

Laboratory Measurement of Velocity Dispersion in Sedimentary and Igneous Rocks
Fatkhan*1, Cysela, R. Y.1, Setiawan, T.1

1
Geophysical Engineering, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Bandung, 40132,
Indonesia
*corresponding author’s email: fatkhan@gf.itb.ac.id
ABSTRACT
Laboratory measurement was carried out to determine the elastic properties of igneous and sedimentary rocks. Sound velocity
measurements determined using contact transducer with frequency bandwidths 100 kHz to 1 MHz. At high frequencies,
the fluid in the compliant cracks will become isolated causing the rock to be stiffer. The behaviour of fluid in the cracks makes
the elastic moduli of the rock frequency dependent. Velocity dispersion is observed only for longitudinal/compressional waves.
The compressional velocities will be dispersive, and the accompanying attenuation will be frequency dependent. It was found
that the velocity has a linear relationship with the frequency.
Keywords:.
.....
A Freeware Tool for Masw Data Processing In Seira Software To Determine Shear
Wave Velocity
Putu Pradnya Andika,1,2* Tedi Yudistira,1* Ruhul Firdaus,2 Imam Catur Priambodo,3
1
Geophysical Engineering, Faculty of Mining and Petroleum Engineering, ITB, Indonesia
2
Geophysical Engineering, ITERA, Indonesia
3
Volcanology Survey Indonesia (PVMBG), Indonesia
*corresponding author’s email : putuandhika12@gmail.com
ABSTRACT
Multichannel Analysis of Surface Waves (MASW) data processing has been widely used to obtain the seismic velocity (Vs)
structure by utilizing surface wave recordings. The Vs profile is one of the fundamental parameters in the geotechnical field.
The Vs model can provide information about soil characteristics and can be used to design and predict dynamic interaction of
soil structures due to dynamic loading. However, the application of the method is very limited to intensively funded projects
due to the lacking of user-friendly freeware implementing the MASW technology. Here, we present our latest development of
such freeware. We have embedded a shear-wave dispersion analysis and Vs inversion module on the previous developed
freeware called SEIRA. We implement the damped least-squared inversion (DLS) method for solving the inverse problem. After
validating our program using synthetic models, we also challenge our development to solve a real case problem.
Keywords:. MASW method, dispersion curve, damped least-squares inversion
24

One Dimensional Magnetotelluric Data Inversion Modeling Using
Convolutional Neural Network
Muhammad Iqbal Khatami1*, Hendra Grandis1
1
Geophysical Engineering, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Indonesia
*corresponding author’s email : work.khatami@gmail.com
ABSTRACT
The magnetotelluric method is a geophysical method commonly used to map subsurface resistivity. The subsurface’s true
resistivity is generated by inversion of the magnetotelluric data. Inversions carried out using conventional methods such as linear
and global approaches have several limitations including long computation time, a lot of iterations, models trapped in local
minimums, and the need for an initial model guess. To overcome the drawbacks, this paper proposes to invert one-dimensional
magnetotelluric data using one of the deep learning methods, the convolutional neural network, which is heavily inspired by the
human nervous system. This method starts by training the network with large amounts of data. The trained network is then used
for inversion by receiving input in the form of apparent resistivity data and generating true resistivity and thickness values
instantly. This method has been tested on synthetic data with curves of type A, H, K, and Q. The inversion results show that the
convolutional neural network could approach the true resistivity and thickness values with a fairly small error and extremely
fast computation time without initial model guess and iteration.
Keywords:. convolutional neural network, inversion, magnetotelluric
.....

Numerical Example to Estimate Groundwater Level Changes using Time-lapse
Microgravity Data
Eko Januari Wahyudi,1* Setianingsih,1 Dasapta Erwin Irawan,1 and R. M. Rachmat Sule1
1
Institut Teknologi Bandung - Indonesia
*corresponding author’s email: ekojw@gf.itb.ac.id
ABSTRACT
The study using time-lapse microgravity (TLM) data can be helpfull to understand near surface variation caused by hidrological
variation. Shallow groundwater variation can be highly correlated with density variation as time-lapse microgravity source,
because the changes are very close to the surface. For several days (or weeks or months), we can observe density variation in
the near surface caused by the changes of water saturated zones. The time-lapse of density variation for several days can be
simplified as groundwater table changes. The purpose of this work is to study time-lapse microgravity application for
groundwater level estimation. Using numerical example, we utilize lateral variation of time-lapse microgravity signal and also
surface topography. In order to do the estimation groundwater level for second TLM survey, we need the groundwater level data
from first survey (and vice versa). After numerical example, we applied the algorithm in the study area of (5 x 5) km square in
the Blora regency. The shallow groundwater level (0 to 14 meters below the surface) are observed after the second (2016)
microgravity survey. From this work, we constrained the model using second (2016) groundwater level data. After iterative
calculation of using TLM anomaly (value in the range of -80 to +80 microGal) with surface topography variation (in the range
90 to 230 meters), we get the estimation the first (2014) groundwater level variation in the range of -6 to 6 meters from the
second (2016) groundwater level data.
Keywords:. time-lapse microgravity, groundwater level, hidrology, iterative calculation.
25
Preliminary Results and Interpretation of Gravity Survey Over Pre-Tertiary

Rock in Northern Kebumen
Faridz Nizar Ahmady1*, Djoko Santoso1, Susanti Alawiyah1
1
Applied Geophysics and Exploration, Institut Teknologi Bandung, Bandung, 40132, Indonesia
*corresponding author’s email: faridnyzer@itb.ac.id
ABSTRACT
The geology of Karangsambung has interesting tectonic setting to study due to complex geological structures. These structures
involve the presence of both pre-Tertiary and Tertiary rock in the area. The pre-Tertiary rock is related with the evidence of
paleo-subduction of Indo-Australian plate beneath Eurasian plate that occurred in Cretaceous era. Both geological and
geophysical approaches is necessity to conduct to help better understanding the subsurface structure in the area. Thus, this paper
will mainly focus on the data acquisition, processing and preliminary interpretation of gravity data in Karangsambung area. The
gravity survey had been conducted in 2017-2018 and had collected 818 data with averaged inter-station spacing of 300 m. The
interpretation will be conducted in respect with current geological interpretation. Based on the analysis, the result shows that
high anomaly in the center of research area indicates the tectonic activity involving Luk-Ulo Melange Complex (LMC) that
occurred in the Cretaceous era rather than the hypotheses of diabasic intrusion. Meanwhile, low anomaly in the northern area is
interpreted to be related with basement configuration of North Serayu Basin.
Keywords:. Karangsambung, tectonic, gravity, Luk-ulo, melange
.....
26
Zoom 3 – Passive seismic (2) Session 2
TUE-06 Zoom 3/ SEACG-001
Early Results of Modeling Subsurface Study in Lake Toba, North Sumatra Using
Receiver Function
Lio Julio,1* Nicolas Silaen,1 Wulan S Ratna,1 Indra R Siregar,1 Sandy K Suhardja,1
1
Universitas Pertamina, South Jakarta
*corresponding author’s email : lio.wae28@gmail.com
ABSTRACT
The island of Sumatra was produced by tectonic action between the Indo-Australian plate subducting beneath the Eurasian plate
to form Mount Toba, which had a large eruption several tens of thousands of years ago to form Lake Toba. The goal of this
study was to use the receiver function approach to calculate the depth variation of the Moho layer and identify the presence of a
low-velocity zone (LVZ) under Lake Toba, and then compare the results of the best receiver function method. Around Lake
Toba, there are 31 data stations from the Lake Toba (LT) project, with data sources from Geofon (Germany). Water level
(frequency domain deconvolution) and iterative time-domain deconvolution are the receiver function methods employed.
According to the findings, the depth of the Moho layer varies between 27.5 and 36.95 kilometers, and it thickens to the northeast
of Lake Toba. The average depth of Moho is roughly 31.05 km, and the Vp/Vs ratio obtained is around 1.86 on average. The
presence of a low-velocity zone was then detected at a depth of 15 km to 25 km in the A-A' section leading from northwest to
southeast on the east side of Lake Toba, and at a depth of 8 km to 22 km in the B-B' section leading from southwest to northeast
across Lake Toba, as indicated by low amplitude. The Investigator Fracture Zone (IFZ) beneath Lake Toba was identified as a
magma reservoir and source of volcanic activity. The iterative time-domain deconvolution approach was then used to illustrate
the result of the best receiver function in this investigation, which revealed that the arrival time of the Ps wave was more spiky
than the result of the water level.
Keywords:. Receiver Function, Low-Velocity Zone, Moho, Lake Toba, Vp/Vs

.....
TUE-07 Zoom 3/ SEACG-003
Early Forward RF Model In South Central Java
R. P. C. Tambunan1, S. K. Suharja2
1
-
1
Faculty of Exploration and Production Technology, University of Pertamina, Jl. Teuku Nyak Arief, Simprug, Kebayoran
2
*corresponding author’s email: putrichri08@gmail.com
ABSTRACT
Earth is a dynamic environment in which to exist. The crust, mantle, outer core, and inner core of the planet are the several
layers that make up the earth's internal structure. It's fascinating to investigate the presence of geological features in it, one of
which is the presence of the Moho layer. Between the lower crust and the upper mantle is the Moho layer, which is a geological
structure. The receiver function approach is one of the geophysical methods used to detect the presence of Moho. A considerable
variation in wave velocity while crossing the layer boundary characterizes this layer. The Moho layer was found in this study
utilizing the Receiver Function technique from Central Java. During the months of May to October 2004, earthquakes with a
magnitude of 5.6 to 9.5 were recorded. The arrival time of the P and Ps waves is calculated using this data. The results of picking
the data are then piled at numerous sites to establish the Ps wave's arrival time, and a Moho layer can then be created. Further
research using the Zhu-kanomori approach and modeling to characterize the Moho layer below the surface is needed to follow
up on the presence of the Moho layer below the surface.
Keywords:. Moho, Receiver Function, Central Java
27
Extracting Rayleigh wave phase velocities from ambient noise cross-correlation:

application to Jakarta seismic network
Rexha V. Ry1,2*, Phil R. Cummins1,3, Babak Hejrani1,3, Sri Widiyantoro2,4
1
Research School of Earth Sciences, The Australian National University, Australia
2
3
Geoscience Australia
4
Faculty of Engineering, Maranatha Christian University, Bandung, Indonesia
*corresponding author’s email : rexha.ry@anu.edu.au
ABSTRACT
For the past decade, cross-correlations of seismic ambient noise have been exploited widely to model the shallow- to deep-
structure of Earth’s interior through tomographic inversions. The stack of cross-correlations between a 2-station pair, i.e., the
noise correlation function (NCF), represents Green’s function and comprises the information of the subsurface structure between
those stations. In practice, NCF is then analyzed to reconstruct the group and phase velocity of surface waves. While extracting
group velocity is relatively straightforward by picking dispersion curves through frequency-time analysis, the procedure for
extracting phase velocity is more elusive, and somehow, it is usually avoided. This study presents a case for temporary seismic
networks deployed in the Jakarta basin; we applied a two-step approach to obtain a phase velocity dispersion curve based on the
relationship between NCF’s spectrum and the Bessel function. First, we calculated an average dispersion curve for our network
following Prieto et al. (2009) approach. Then, we used this average on each pair to select a reliable dispersion curve in the
frequency domain (Ekstrom et al., 2009). For all station pairs in the network, we successfully estimate ~2660 Rayleigh wave
phase velocity dispersion curves from vertical NCFs that, in further processing, they are usable for 2D tomographic inversions.
These results of phase velocities enrich the information on subsurface structure around the Jakarta basin, and we suggest that
the procedure can be conducted for other cases. It makes it possible to invert for shear-wave velocity (Vs) from Rayleigh wave
phase velocity or jointly invert for Vs using both phase and group velocity. Moreover, joint inversion is beneficial for reducing
the model’s ambiguity; hence, the result may represent a more accurate model.
Keywords:. ambient noise, cross-correlation, Rayleigh wave phase velocity, Bessel function.
.....
Preliminary Result of Shear-wave Velocity Estimation in Central Part of Sulawesi

using Microtremor Data
Sari Ratna Dewi1*, Zulfakriza2, Sri Widiyantoro2
1
Geophysical Engineer Undergraduate Program, Faculty of Mining and Petroleum Engineering, Bandung Institute of
Technology, Bandung, Indonesia
2
Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Bandung Institute of Technology,
Bandung, Indonesia
*corresponding author’s email : sratna518@gmail.com
ABSTRACT
Central part of Sulawesi lie on active tectonic region and higher seismicity events. Several significant earthquakes occurred in
this region, such as 1927 with a magnitude of 6.5 in Donggala, 1969 with a magnitude of 6 in the city of Tambu, and 2018 in
Palu with a magnitude of 7.5. Investigation of site effect, seismic amplification and Vs30 in Central Sulawesi is important to
identify. We used 33 seismic stations across central part of Sulawesi to estimate dominant frequency, seismic amplification and
shear wave velocity. Microtremor is a ground vibration due to the interference of body waves and surface waves that generated
by human activity and natural source. We apply the HVSR method to produce horizontal and vertical curve and allow us to
identify the structure of V_s. The dominant frequency and amplification factor in this region with a range of 0.55 – 18.42 Hz
and 1.19 – 6.0. Then, we apply the Neighborhood Algorithm to invers the HVSR curve to obtain the profile of V_s to depth .
The value of V_s in range of 74.44 – 3045.77 m/s for a depth range up to 818.56 meters. The value of Vs correspond to the
classification of the soil layer and the soil type is determined by averaging V_s to a depth of 30 meters. We interpreted the
central part of Sulawesi region is dominated by medium soil (SD) and very dense soil (SC).
Keywords:. microtremor, HVSR, amplification, predominant frequency, V_s
28
Modeling Crustal Structure in West Sumatera Using Teleseismic Receiver

Function
Wulan Sri Ratna,1* Lio Julio,1 Indra Rivaldi S.,1 M. Aldi,1 Mohamad Ramdhan,2 Sandy
Kurniawan Suhardja,3 Iktri Madrinovella,3 Ida Herawati,3 Bayu Pranata,4
1,3
Universitas Pertamina, Jl. Teuku Nyak Arief, RT.7/RW.8, Simprug, Kec. Kby. Lama, Kota Jakarta Selatan, DKI Jakarta
12220, Indonesia
2
National Research and Innovation Agency (BRIN), Kompleks BRIN, Jl. Sangkuriang, Kompleks BRIN, Kecamatan Coblong,
Kota Bandung, Jawa Barat 40135, Indonesia
4
Badan Meteorologi Klimatologi dan Geofisika (BMKG), Jl. Angkasa 1 No. 2, RW. 10,Gn.Sahari Sel., Kecamatan
Kemayoran, Kota Jakarta Pusat, DKI Jakarta 10610, Indonesia
*corresponding author’s email : wulansriratna572@gmail.com
ABSTRACT
Indonesia is located between three significant plate confluences that play a role in creating a subduction zone. The Semangko
Fault, the Mentawai Islands, a chain of non-volcanic forearc islands, and active faults that induce earthquakes and landslides
are all produced by subduction zones on numerous Indonesian islands, one of which is Sumatra Island. The receiver function
approach is developed to identify the discontinuity limits of the Earth's structures, such as the crust and mantle, from
teleseismic events. The information in this study comes from a temporal and permanent station in West Sumatra. The
earthquake data was taken from teleseismic earthquakes from the three components of seismic stations using criteria arc
distances between 30o and 90o and magnitudes ranging from 6.0 to 10.0. The Semangko Fault and active volcanic mountains
include the research region, affected by a subduction zone. This research aims to determine how crustal thickness is distributed
and what effect regional geological structure has on crustal thickness. This work used the receiver function approach and the
Zhu-Kanamori method to estimate crustal thickness. According to the early result study, the estimated range of crustal
thickness in the research region is 27.11–32.48 km, and the range of Vp/Vs ratio values is 1.70–2.19. The direction of crustal
thickness thickening in the area was likewise East-South.
Keywords:. Receiver Function, Crustal Thickness, Zhu Kanamori, Teleseismic, Vp/Vs ratio
.....
29
Application Of Seismic Multiattribute for Predicting The Distribution Of Reservoar

Property In “A” Field, Kutai Basin
P D Wardaya1, A Novianto2, A K Pratiknyo3, H Hamdalah2, Y Rizkianto4, A M Handini5

1
Pertamina Upstream Research and Technology Innovation
2
Geophysical Engineering Lecturer, UPN “Veteran” Yogyakarta
3
Petroleum Engineering Lecturer, UPN “Veteran” Yogyakarta
4
Geological Engineering Lecturer, UPN “Veteran” Yogyakarta
5
Geophysical Engineering Student, UPN “Veteran” Yogyakarta
*corresponding author’s email : -
ABSTRACT
Exploration in deep water area often has problems, this is due to the poor quality of seismic data and the limited data exploration
wells. This causes reservoir modeling to be difficult, so we need a method that can be used with limited data, that is Multi-
attribute method. This is a statistical method that uses more than one attribute to predict several physical properties of the earth
by looking for the relationship between logs and seismic data at the well location. This analysis is divided into 3rd step, that is
single attribute analysis (bivariate geostatistics), multi-attribute linear regression analysis and probabilistic neural network
(PNN) analysis. Bivariate geostatistics and multi-attribute linear regression utilize a linear relationship between physical
parameters of log data and seismic attributes, while probabilistic neural networks utilize non-linear relationships. This research
will predict the distribution of petrophysical parameter values including shale volume, porosity, and incompressibility in the
case study of the upper Miocene sandstone layer, deepwater area of Kutai Basin. The distribution of the results shows that there
is a potential zone in the southern of the study area in the slope channel fan with a low acoustic impedance value of 4000 to
7500 m/s*gr/cm3, the value of shale volume is less than 45%, the porosity value is more than 10%, and the incompressibility
value is less than 22 Gpa*g/cc.
Keywords:.
.....
Delineation Of Gas Sand Reservoir Prospect in The “A” Field, Deep Water Area
Of Kutai Basin
P D Wardaya1, A Novianto2, A K Pratiknyo3, H Hamdalah2, Y Rizkianto4, A M Handini5
1
Pertamina Upstream Research and Technology Innovation
2
Geophysical Engineering Lecturer, UPN “Veteran” Yogyakarta
3
Petroleum Engineering Lecturer, UPN “Veteran” Yogyakarta
4
Geological Engineering Lecturer, UPN “Veteran” Yogyakarta
5
Geophysical Engineering Student, UPN “Veteran” Yogyakarta
ABSTRACT
Potential hydrocarbons in Kutai basin is not only in the delta area, but extends to the deepwater area of East Kalimantan. Results
of the integration with petrophysical, rock physics, inversion and seismic multiattributes methods aim to delineate the gas
reservoir potential zone in the research area. The reservoir rocks consist of graywacke to sublitharenite of Upper Miocene
sandstone which has an average shale volume content of 18%, porosity 20%, and water saturation 60%. The fluid filling the
reservoir is gas which is reflected in the Vp/Vs ratio range of 1.58 to 1.9, the Poisson ratio value of 0.17 to 0.3 GPa, and the
incompressibility value of 8 to 22 Gpa*g/cc. The distribution potensial reservoir zone is indicated by the low acoustic impedance
anomaly value of 4000 to 8000 m/s*gr/cm3, high RMS attribute anomaly value of 2000 to 10000, and distribution of low shale
volume, high porosity, and low incompressibility values will be predicted using multi-attribute linear regression and probabilistic
neural network methods which is spread over the slope fan channel in the southern part of the study area.
Keywords:. -
30
Low Frequency Seismic Extrapolation in Full Waveform Inversion (FWI) using U-

Net
Asido Saputra Sigalingging,1,2*, Sonny Winardhi,1, Ekkal Dinanto,1
1
2
*corresponding author’s email : asido.saputra@gmail.com
ABSTRACT
Full Waveform Inversion (FWI) modeling is dependent on many factors, namely the initial model, source wavelet, and low
frequency of seismic data. The lack of initial model and low frequency data can affect the result of FWI modelling due to cycle
skipping problems. Low frequency data is one of the crucial problems that must be tackled. The loss of low-frequency data can
remove the trend of geological models. To deal with that, low-frequency data is recon-structed using deep learning methods.
We use a 1D U-Net algorithm to automatically extrapolate low frequency data from bandlimited Common Shot Gather (CSG)
seismic data in the time domain without pre-processing steps. The bandlimited seismic data is the input in deep learning, and
the algorithm predicts low-frequency seismic data as the output. The U-Net model was tested and validated with seismic
synthetic data. The results show that U-Net can extrapolate low frequencies with high accuracy.
Keywords:. Deep learning, U-Net, Full Wave-form Inversion (FWI), Low-Frequency, Seismic Imaging.
.....
Data Conditioning for Seismic Reservoir Characterization of Thin Bedded Fluvio-

Deltaic Reservoirs
Andri Hendriyana1), M. Noor Alamsyah2), Dona Sita Ambarsari1), Istifani Ferucha1), Diah Ayu
Widianingsih2), Sigit Sukmono1)
1)
Geophysical Engineering, Institute of Technology Bandung (ITB)
2)
PetroChina International Jabung Ltd.
ABSTRACT
Reservoir characterization using seismic data inversion is an important tool to make quantitative interpretation of reservoir
properties. An accurate estimate of reservoir properties could lead to better reservoir delineation, reduce uncertainty and
exploration risk, as well as better estimate of net pay. However, seismic inversion requires a good quality of seismic data. Then,
seismic data conditioning is necessary. The objective of this study is to define the best data conditioning workflow for seismic
reservoir characterization of thin-bedded fluvio-deltaic reservoirs in AHY field, Indonesia. In the legacy data, the stratigraphic
and structural features associated with the deposition of the studied reservoir are not clear due to the presence of noises and low
signal-noise ratio. The applied workflow for pre-stack gather consists of angle muting, bandpass filtering, super-gather stacking,
random noise attenuation and trim static. The results show that the workflow has successfully removed random noise and residual
moveout as results of inaccurate picked velocity. Furthermore, signal enhancement was employed to the partial angle stack
consisting of structure oriented mean filtering and bandwidth matching. After post-stack enhancement, the signal-to-noise ratio
was increased, resulting in clearer and crispier imaged of the stratigraphic and structural features. Overall, the conditioned
seismic data facilitate much better seismic attribute, inversion and AVO analysis of the thin-bedded fluvio-deltaic reservoirs.
Keywords:. -
31
3D Check shot Modeling for Time-Depth Conversion

Jefri Bambang Irawan1, Wahyu Triyoso2
1
2
Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Bandung Institute of Technology,
Indonesia
*corresponding author’s email : jefribambang@gmail.com
ABSTRACT
Time-depth conversion is one of the keys to successful drilling with a specific target. However, if the drilling target is not defined
correctly, it can cause significant material losses. For example, several cases found that the actual hydrocarbon depth was deeper
than the expected target. And often, the well will be assumed to be a dry well even though it still has hydrocarbon reserves. To
overcome this problem, in this study, we intend to develop 3D check shot modeling in which check shot and formation marker
data are referred to and taken at the well to obtain a precise time-depth relationship. However, not every well has check shot
data, so limited data can be used to do velocity modeling for the time-depth conversion process. In development fields with
dozens or even hundreds of wells, the time-to-depth conversion process becomes a challenge where it is not enough only to use
one single velocity to do it. Still, the velocity model must accommodate heterogeneity laterally so the depth structure map can
match well formation markers. On the other hand, seismic data can indirectly measure subsurface velocities vertically and
horizontally, so this research uses a method by combining well check shot data with seismic semblance velocity data to build a
velocity model that is calibrated by well data. In this study, the time-to-depth conversion process with the calibrated velocity
model has an average residual of less than 5-meters. In contrast, if only using check shot data, the average residual obtained is
25-meters.
Keywords:. time-depth conversion, 3D check shot modeling, well formation marker, velocity model
.....
32
Zoom 5 – Volcanology (1) Session 2
Geophysical Overview over the Paleo-Volcanic-Caldera-Rims Anomaly, from

Synthetically Modelling to the Field Cases
Accep Handyarso1*
1
Research Center for Geological Resources, National Research and Innovation Agency
Jalan Sangkuriang, Bandung 40135, Indonesia
*corresponding author’s email : accep.handyarso@brin.go.id
ABSTRACT
The concept of Superimposed-Volcanism, which was discovered in 2006 by Prof. Sutikno Bronto, implies the possibilities of
Superimposed-Paleo-Volcanic-Calderas or Calderas-Complex. Calderas are a crucial feature in any volcanic environment and
are commonly the prospect site of geothermal and mineralization. The Indonesian archipelagic area has several magmatic belts
with different periods. However, only the Tertiary-Magmatic-Arc is interpreted as the Tertiary-Volcanic-Arc, which refers to the
recent active volcanoes. We infer many older calderas exist in the region, but it hasn't been exposed yet. A comprehensive
understanding of regional geology will guide the geological resources explorations. Gravity and Magnetic are versatile tools
offering rapid and wide-coverage survey areas for delineating the causative bodies. Geophysical anomalies overview due to any
Volcanic-Calderas is needed, especially for Gravity and Magnetic methods. There are unique Gravity and Magnetic anomalies
found in Central Java, Indonesia. The interpretation processes invoke the inverse-modeling mechanism to examine the geological
situation. An inverse-modeling approximation was conducted to simulate the anomaly pattern due to a lateral offset
(displacement) and ring-shape causative body. We proposed the area to be the Paleo-Volcanic-Caldera. The existing fault zone
manifested as the weak zone where the intrusive rocks breakthrough and raise volcanoes as a Tecto-Volcano-System in the area.
Keywords: Paleo-Volcanic-Calderas, Gravity, Magnetic, Inverse-Modeling, Indonesia
.....
Analyses of Explosion Earthquakes at Sinabung Volcano: Characteristics of Waveform,

Spectra, and Energy
Aulia Dhita Pritaningtyas1*, Mohammad Hasib2*, Estu Kriswati2, Tati Zera1
1
Department of Physics, Faculty of Science and Technology, Syarif Hidayatullah State Islamic University Jakarta (UIN)
2
*corresponding author’s email : auliadhiza77@gmail.com
ABSTRACT
We investigate the characteristic of explosion earthquakes associated with avalanche at Sinabung Volcano located in Karo
Highlands, Karo District, North Sumatra, Indonesia during October 1 st to December 31 th , 2014. We use 3 short-period
seismometers which are Sukanalu (SKN), Laukawar (KWR), dan Kebayaken (KBY). All of them is located about 3-9 Km
away from the active crater. Our results show the characteristic of these waveforms occurs more than 60 seconds (c.a. 70 – 120
seconds), and most of them have dominant frequency is about 1-4 Hz. Our result may be different with the characteristic of
explosion earthquakes in general which shows short time occurrence (c.a. 10-40 seconds). We also found that the energy of
explosion earthquake changes with lapse time. Our results suggest the explosion earthquakes at Sinabung Volcano during
October 1 st to December 31 th , 2014 may have a different waveform generation and/or mechanism to the general explosion
earthquake at different volcano.
Keywords: Sinabung Volcano, explosion earthquake, avalanche, waveform, spectra, energy.
33
Geological Characteristics Identification of The Western Tangkuban Parahu

Volcanic Area Using H/V Analysis
Vigo Denis Radestian1*, Zulfakriza2, Andri Dian Nugraha2, Aditya Lesmana1

1
Geophysical Engineering Study Program, Faculty of Mining and Petroleum Engineering Bandung Institute of Technology, Bandung,
Indonesia
2
Global Geophysics Research Group, Faculty Mining and Petroleum Engineering Bandung Institute of Technology, Bandung, Indonesia
*corresponding author’s email : vigoradestian@gmail.com
ABSTRACT
Tangkuban Parahu Volcano is one of the active volcanoes in West Java that located about 15 km northern part of Bandung City.
Various methods have been carried out to identify geological characteristics in the Tangkuban Parahu Volcano Area such as
tephra-stratigraphy and surface geology studies. For further understanding of the geological characteristics at Tangkuban Parahu
Volcano, especially the western part, we used H/V method on microtremor data recorded by 17 temporary seismic stations across
the western part and around the volcano crater. We estimate the H/V curve for each station, then conduct the inverse modelling
to obtain the VS. We perform selection of data time range, demeaning, detrending, filtering, selection of STA/LTA, and
elimination of poor data. Then, base on H/V curve allow us to estimate the dominant frequency, amplification factor, and through
inverse modelling, VS30 values can be obtained for each station. We apply the Inverse Distance Weighted (IDW) interpolation
method to figure up the general profile of dominant frequency, amplification factor, and VS30 map. We obtain the dominant
frequency in range of 0.93-9.7 Hz, amplification factor in range of 2.4-11.2, and VS30 in range of 196,04-323,65 m/s. We
interpreted the results were correspond to volcanic tuff around the volcano crater and more compact geological unit that comes
from the older volcanoes in the western part of study area.
Keywords: Tangkuban Parahu, microtremor, geological characteristics, H/V analysis, Inverse Distance Weighted
.....
The Interpretation of Subsurface Structure of Banda Api Volcano Based on Gravity

Modelling
Faridz Nizar Ahmady1*, Djoko Santoso1, Fauzan Ardana2
1
2
Undergraduate Program, Institut Teknologi Bandung, Bandung, 40132, Indonesia
*corresponding author’s email : faridnyzer@itb.ac.id
ABSTRACT
Indonesia has 127 active volcanoes along subduction zone, one of them is Mount Banda Api. As active volcano which is
surrounded by the sea, volcanic eruption may trigger destructive tsunami in the adjacent coastline. Several geological studies
had been conducted to understand the history of volcanic eruption based on the surface mapping. However, in order to further
interpretation regarding subsurface activity of the volcano geophysical data is needed. Therefore, the purpose of this study was
to understand the subsurface structure of Mount Banda Api based on the interpretation of gravity data. The Complete Bouguer
anomaly are acquired by PSDMBP in 2016. The complete Bouguer anomaly range is 130-200 mgal with low values located
on the Banda Api volcano and high anomalies on Neira Island and Lonthoir Island. Then 2D spectral analysis is conducted to
estimate the depth of regional dan residual anomaly and determine the optimal Moving Average (MA) windows to conduct
regional-residual separation. Furthermore, the residual anomaly and background density of 2.63 g/cc was used to interpret the
subsurface structure beneath the volcano by conducting forward modelling 2,5 D. The resulting model shows that the low
anomaly in the study is estimated as a response to the presence of magma at the bottom of Mount Banda Api, while the high
anomaly is a response from rocks that have an old age. In the geological model, five faults were interpreted. Two faults indicate
the collapse of the Ancient Neira Volcano, resulting in the structure of the Neira Caldera.
Keywords: Banda Api, volcano, Banda, Neira, gravity, modelling
34
Digitization of Analog Seismograms using Machine Learning Methods
Martanto Martanto1, Andri Dian Nugraha2, David P. Sahara2, Zulfakriza Zulfakriza2, Devy Kamil
Syahbana1, Corentin Caudron3
1
Center for Volcanology and Geological Hazards Mitigation, Geological Agency, Ministry of Energy and Mineral Resources, Indonesia
2
Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Bandung, Indonesia
3
Université Libre de Bruxelles, Beligum
*corresponding author’s email : martanto@live.com
ABSTRACT
Indonesia is the country with the most volcanoes in the world, to provide warnings to communities in volcanic disaster-prone
areas, a fast decision-making process is needed and supported by the speed of monitoring data processing. To achieve this goal,
monitoring data needs to be acquired locally and processed meaningfully in real time. Volcano monitoring uses geophysical,
geochemical and geological methods, with efficient monitoring benefiting from multidisciplinary approaches. By providing an
unrivalled resolution in both time and space, volcano seismology constitutes the backbone of many volcano monitoring programs
worldwide. However, this method still has disadvantages, especially in Indonesia, whose monitoring still uses analog
seismograms. Volcanic signal detection and classification are still often carried out manually by operators. This manual process
is extremely time consuming and prone to errors and subjectivity. Identical volcano monitoring data, observed in very similar
settings, are described in different terminologies, and interpreted in such diverse ways that the resulting advice leads to opposing
conclusions. Recent advances in Machine Learning (ML) approaches provide schemes to automatically detect and classify
volcano-seismic events in an unbiased way. To approach this method, digitization of the seismogram is needed and carried out
by scanning the seismogram into a digital. Having merged digitised paper seismograms with existing digital data, we will have
built a new archive of continuous seismic data. The Convolutional neural network (CNN) method is used in digitised paper
seismograms to predict analog seismogram data (amplitude and time). This methodology has the potential to help speed up the
analysis process and understand volcanic processes more deeply.
Keywords: volcanology, volcano-tectonic, machine learning, seismogram
....
35
Zoom 1 – Tomography (1) Session 3
WED-01 Zoom 1 / SEACG-005

The early model of tomography in Eastern Indonesia using FMTOMO
Elta Purnama Wulan1, Sinta Nur Asyidah2, Yusrina Alfiani Haqi3
1
Geophysical Engineering Study Program
2
Pertamina University
ABSTRACT
The Indo-Australian Plate, Eurasian Plate, Pacific Plate, and Philippine Plate collide in eastern Indonesia, making this a
seismically active zone. Tectonics in Eastern Indonesia is extension and subduction. The Banda – Seram segment in Eastern
Indonesia results from the collision between three major plates, namely the Eurasian Plate, Australian Plate, and Pacific Plate,
causing this segment to have quite complex tectonics. There are also tectonic mechanisms in the Banda – Seram segment, such
as rollback, detachment, and slab tear. Tectonic tools and active tectonic motions generate partial melting, back-arc thrust, and
subduction. This study aims to find out how the tectonic structures in Eastern Indonesia, such as partial melting, subduction
patterns, and others, by using tomographic modeling and knowing the mitigation of natural disasters in the form of earthquakes
in this area based on the results of imaging the structure of seismic waves in the form of P waves. This study uses earthquake data
from the BMKG earthquake catalog. The data is processed using two software, namely hypoDD and fast marching. It wares
discovered that a subduction slab structure is associated with high seismic wave velocities up to a depth of 550 km based on the
results of this data processing. It is indeed visible in the high seismic wave velocity anomaly, which hares discovered to be a
back-arc thrust zone at a depth of 100 kilometers. This study demonstrated the relationship between subduction zones and their
influence on high seismic and magmatic activity in eastern Indonesia.
Keywords: Fast Marching Method, Banda-Seram Segment, Eastern Indonesia, Partial Melting, Back-arc Thrust
.....

2D and 3D Modelling of Molucca Sea Double Subduction Zone Based on Primary Waves
Travel Time Seismic Tomography Inversion
Gendis Ashella Nareswari Pratita1, Indriati Retno Palupi1, C. Prasetyadi2

1
Jurusan Teknik Geofisika, Fakultas Teknologi Mineral, Universitas Pembangunan Nasional “Veteran” Yogyakarta
2
Jurusan Teknik Geologi, Fakultas Teknologi Mineral, Universitas Pembangunan Nasional “Veteran” Yogyakarta
*corresponding author’s email : gendis.ashella@gmail.com
ABSTRACT
The Molucca Sea is one of the areas with a high level of seismicity in Indonesia due to activity between three plates, The
Eurasian Plate (Sangihe Microplate) and the Philippine Sea Plate (Halmahera Microplate), and Molucca Sea Microplate. The
activities of these three plates configure a unique subduction zone. The Eurasian Plate (Sangihe Microplate) which continues to
move towards the East and the Philippine Sea Plate (Halmahera Microplate) continues to move Westward, causing the Molucca
Sea Microplate to be pressed into the earth and subduct in two opposite directions at the same time forming an inverted “U.”
The research uses earthquake event data from as many as 3663 events and the coordinates of earthquake stations as many as 9
earthquake stations in the last 5 years period with a magnitude of 3 - 8 Mw and a depth of 3 - 630 km. Then the seismic
inversion of travel time tomography was be used to make 2D and 3D models of the double subduction zone in The Molucca
Sea.
The distribution of the high ΔVp anomaly with a value of 3.5% to 4% is associated with a solid medium, like the Sangihe and
Halmahera plates or the subduction slab of the Molucca Sea Plate. While the low ΔVp anomaly is 1.8% to 2% associated with
weak or crushed medium and the presence of thermal fluids such as magma or partial melting (the presence of volcanoes). From
the results of 2D and 3D modelling, it can be seen that the subduction under the Sangihe Microplate is deeper (627.2 km) with
an average subduction angle of 45.292o compared to the subduction under the Halmahera Microplate (280 km) with an average
subduction angle of 35,836o.
Keywords: Molucca Sea, Subduction, Travel Time Tomography, 2D and 3D Models, ΔVp
36

Early Results of Surface Tomography Using Automated Surface-Wave Phase-Velocity
Measuring System (ASWMS) in Hawaii
Kamaludin Nursal*, S.K. Suhardja
Faculty of Exploration and Production Technology, Pertamina University, Jl.Teuku Nyak Arief, Simprug, Kebayoran Lama, Jakarta
12220, Indonesia
ABSTRACT
The Hawaiian Islands are a well-known place to conduct seismicity research. Geologically, occurring earthquakes in this region
are unique since it is caused by the presence of a hot spot rather than the common fault activity. The research was conducted
using a tomography method called Surface Wave Tomography. The goals of this study are to produce tomographic images of
the Hawaiian Islands and to test the accuracy of the Automated Surface-Wave Phase-Velocity Measuring System (ASWMS).
ASWMS employs a cross-correlation process to calculate the phase delay between stations, then Eikonal equation to invert the
slowness vector to obtain apparent phase velocities, and lastly Helmholtz equation to correct the amplitude to obtain structural
phase velocities. The teleseismic data used in this study were collected from stations located throughout the Hawaiian Islands
between 2004 and 2009. This study discovered a low-velocity anomaly beneath the Hawaiian Islands with a value of 3.8 – 3.9
km/s, indicating the presence of a low-velocity body. The increasing velocity change is interpreted as the formation of a
Mohorovicic discontinuity as a result of magmatic underplating at a depth of 33 km. The decrease in velocity with increasing
depth indicates an increase in temperature associated with the path that magma takes from the Earth's mantle to the Earth's crust
beneath the Hawaiian Islands.
Keywords: ASWMS, Phase Velocity, Hawaii, Eikonal, Helmholtz
.....

Early Result of Imaging 3D Seismic Velocity Structure In Central Java Using Double-
Difference Tomography
Ryandika Saputra1*, Sandy Kurniawan Suhardja2, Iktri Madrinovella3,Mohammad Ramdhan4

1
Universitas Pertamina, Jakarta
2
BMKG
*corresponding author’s email : ryandikasaputra0@gmail.com
ABSTRACT
Central Java is a part of the Sunda Arc and has relatively high seismicity due to the subduction zone between the Indo-
Australian Plate and the Eurasian Plate. This study uses double-difference tomography to image the P and S wave 3D seismic
velocity structures associated with tectonic patterns due to subduction zones. The data used comes from the earthquake catalogs
of BMKG, BPPTKG, DOMERAPI, and MERAMEX, with a recording period from May 2004 to December 2020. The number
of earthquakes that successfully relocated was 1930 from 1937 earthquakes recorded by 285 recording stations. The results of
hypocentre relocation show that seismic activity in Central Java is relatively high and originates from geological structures,
including subduction zones, back-thrust, and the Opak Fault. The tomogram results of the P and S wave velocity models
support each other. The high seismic velocity anomaly is associated with the subduction zone of the Indo-Australian Plate to
the Eurasian Plate with a maximum resolution of up to 100 km depth. Low seismic velocity anomalies are found at the Merapi-
Lawu Anomaly (MLA), the Modern Volcanic Arc (MVA), the Sumbing-Sundoro-Dieng volcanic complex, the Banyumas-
Cilacap sedimentary basin, and the Kendeng Basin. A low-velocity anomaly is found at a depth of about 40-50 km and a depth
of 100 km, which is associated with the process of slab dehydration and partial melting. The inversion results show the impact
of the Indo-Australian Plate subduction and the Eurasian Plate on volcanic activity, seismicity, and the geological structures
developed in Central Java.
Keywords: Double-Difference Tomography, Hypocenter Relocation, P and S Wave Velocity Structures, Tectonic, Central Java
Preliminary Result of Local Earthquake Tomography Around Palu Koro Fault using BMKG
Seismic Network
37
H Jayadi1,2*, B J Santosa3, A D Nugraha4, D D Warnana5, R Jamroni6, P Supendi7, Z Zulfakriza4, F

Muttaqy8, G Rachman9, A S Putra8, I U Meidji10
1
Doctoral Program Department of Physics, Faculty of Science and Data Analytics, Institut Teknologi Sepuluh Nopember, Surabaya 60111,
Indonesia
2
Geological Engineering Study Program, Engineering Faculty, University of Tadulako, Palu 9411, Indonesia
3
Department of Physics, Faculty of Science and Data Analytics, Institut Teknologi Sepuluh Nopember, Surabaya 60111, Indonesia
4
Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Bandung 40132,
Indonesia
5
6
Agency for Meteorology, Climatology and Geophysics (BMKG), Makassar 90231, Indonesia
7
Agency for Meteorology, Climatology and Geophysics (BMKG), Jakarta 10720, Indonesia
8
Graduate Program of Geophysical Engineering, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Bandung
40132, Indonesia
9
Physics Education-FKIP, Universitas of Pattimura, Molucca 97233, Indonesia
10
Physics Study Program, Faculty of Natural Science and Mathematics, Universitas Negeri Gorontalo, Gorontalo 96128, Indonesia
*corresponding author’s email : harsanoj@gmail.com
ABSTRACT
The Central Sulawesi region has very complicated complex geological and tectonic conditions. One of the things that became
a topic in these tectonic conditions was the existence of the Palu Koro Fault, which crossed the Palu City area so that it often
resulted in solid and destructive earthquakes. Previous research has relocated the earthquake and mapped the coastline due to
the quake around the Palu Koro fault area. In this study, we used data from local earthquakes around the room. Data on the
arrival time of P from the local earthquake network belonging to the Meteorology, Climatology and Geophysics Agency
(BMKG) for the period 2010 to 2019. The seismic tomography method was applied to determine the structure of the seismic
velocity and, at the same time, locate the hypocenter. Using the inversion tomography method, we used a 1-D initial velocity
structure model, and the results obtained were evaluated using a checkerboard test, derivative weight sum, and ray hit count.
The initial results of the tomographic inversion show a relatively high seismic velocity in the west and east of the Palu Koro
fault area. These areas are Mendala West Sulawesi and Mendala East Sulawesi (Poso Fault) and Wekuli Fault. The low-speed
anomaly is around the Palu Koro fault area. The relocated hypocenter also clustered around the Palu Koro fault area due to the
high-speed thrust. For the next activity, we will determine the structure of the S wave velocity, which is expected to get better
geological information in its interpretation.
Keywords: -
38
Zoom 2 – Advance Geophysics (2) Session 3

Applying Map Filters and Inversion Model for Gravity Data on Matlab-based Application
M Tsaqif Ghifari,1*, Fawwaz Byru Fitrianto,2 Indriati Retno Palupi,3 Wiji Raharjo,3
1
Undergraduate Program of Geophysical Engineering, Faculty of Tecnology Mineral, Universitas Pembangunan Nasional “Veteran”
Yogyakarta, Indonesia
2
Postgraduate Program of Geophysical Engineering, Faculty of Tecnology Mineral, Universitas Pembangunan Nasional “Veteran”
Yogyakarta, Indonesia
3
Geophysical Engineering, Faculty of Tecnology Mineral, Universitas Pembangunan Nasional “Veteran” Yogyakarta, Indonesia
*corresponding author’s email : tsaqif.gm@gmail.com
ABSTRACT
In general, there are two types of processes using the gravity method, there is the basic process and the advanced process. This
basic process includes Gravity Meter reading conversion until the gravity corrections process. While the advanced process is a
process to put an edge on the geological appearance of the research area that is already in the form of a map. The advanced
process includes map filtering and inversion modeling. This advanced processing of gravity data is very useful, one of which
is to separate regional anomalies and local anomalies, while inversion modeling is useful for knowing the approximate location
of the anomaly depth. What will be applied to the software is advanced processing, because in its application this processing is
required in every geophysical research and requires data processing support tools. However, for now, gravity data processing
software, especially for the filtering and data modeling process, is available as a licensed software. Where not everyone can
access it because of limited costs and those who are still in the learning stage in lectures need practicing to process the
geophysical data. The purpose of implementing this filtering process and inversion modeling is to provide access to gravity
data processing for free and open source so that it can make it easier for students to process the data during practice in the field,
and in the future this software can be the basis for development and learning in the field of geophysical computing. This research
was made using a descriptive method with a complete step and explanation of the creating and applying of the geophysical data
processing into a software, so that the processes inside it can be known. A literature study was also conducted to determine the
basic formula that will be used as a reference in calculating gravity data in order to obtain the expected processing and modeling
results. After that, the final step is to apply all of that in a program to create software based on MATLAB. The final result of
this research is a software containing map filtering and inverse modeling for gravity data which can be used for data processing.
Keywords: Gravity, Data Processing, Software, Filtering, Inverse Modelling
.....

Tomographic Inversion Code Analysis using Regularization Iteratively Reweighted Norm (IRN)
Belista Ajeng Natasya*, Ryandika Syahputra, Elta W, Kamaludin Nursal, Aulia Rizki, Nicolas
Silaen, Sandy K Suhardja
Geophysical Engineering Department, Pertamina University
*corresponding author’s email : ntasya.lh@gmail.com
ABSTRACT
The island of Java was formed by the collision activity of the Indo-Australian plate and the Eurasian plate, which increased the
activity of several volcanoes in the Java Island area, particularly Central Java. One of the most active volcanoes in Central Java
and the world is Mount Merapi.). In the tomographic inversion process, ill-condition problems are often found in matrix
making. Therefore, a regularization technique is needed to obtain a consistent result. The regularization technique can be done
using IR Tools In this study, inversion techniques are used in the form of Least Square and Iteratively Reweighted Norm (IRN).
This study aims to analyze the tomographic inversion code with the regularization technique using the Iteratively Reweighed
Norm (IRN). In this study, two types of data were used: synthetic and real data. Using synthetic data, a Tomographic code was
applied and compared with the Least Square and Iteratively Reweighed Norm (IRN) methods. The results show that using the
Least Square method has better results where the resulting checkboard is more than adequate to reconstruct the initial velocity
and amplitude model, which is described more clearly and strongly. Based on the results of the real data inversion, there is a
low Vp anomaly (red) in the lower area of Mount Merapi and Mount Merbabu. Low seismic velocity anomalies are associated
with magmatic zones, partial melting, or magma reservoirs. However, a high seismic velocity anomaly (blue) is also associated
with the subduction zone of the Indo-Australian Plate to the Eurasian Plate.
Keywords: IRN, LSQR, Merapi, velocity
39
An Attempt to Invert Tomography Delay Time Problem With Hybrid Lsqr
Aulia Rizkiyana1*, Ryandika Syahputra1, Elta W 1, Kamaludin Nursal1, Nicolas Silaen1 , Belista
Ajeng N 1 ,Sandy K Suhardja1, 2
1
Faculty of Exploration and Production Technology, Pertamina University, Jl. Teuku Nyak Arief, Simprug, Kebayoran Lama, Jakarta
12220, Indonesia
2
*corresponding author’s email : auliarizki408@gmail.com
Indonesia is situated in Southeast Asia's active tectonic zone, influenced by the movement of four plates: the Eurasian Plate
in the south, the Eurasian Plate in the north, the Philippine Plate, and the Pacific Plate in the northeast. The Tomography
approach was used in this work to obtain a stable solution, which was achieved through regularization. Least Square and
Hybrid Lsqr are the inversion methods employed. This research aimed to use the IR Tools Package Irhybrid LSQR to
investigate inversion tomography with regulatory approaches in 3D. In this study, Matlab is used to perform inversion
tomography. Two types of data, real and synthetic, are employed at this stage. The study's genuine data came from earthquake
catalogues from BMKG, MARAMEX, DOMERAPI, and BPPTKG, with depths ranging from 0-658 km and recorded
between 20 December 2013 - 6 September 2020, with an area boundary of 110.1 – 110.8. LS -7.8–7.2 BT. The red colour
represents the magma reservoir area according to the inversion data. Hybrid LSQR is an iterative approach that is considered
stable enough to be used in tomography. When applied to tomographic coding and compared to the most widely used
inversion method, Least Square, the results obtained in the least Square have a higher value than Hybrid LSQR.
Keywords: Central Java, Hybrid, LSQR P and S waves, Tomography
.....

Ultrasonic Tomography of Core Plug Using Fresnel Volume Approach
Hafidz Naufal Aryan,1* Teuku Abdullah Sanny,1 Dona Sita Ambarsari,1

1
Geophysical Engineering, Faculty of Mining and Petroleum Engineering, Institute of Technology Bandung, Indonesia
*corresponding author’s email : hafidznaufalaryan@email.com
ABSTRACT
Tomography was a non-destructive method for investigating the internal structure of an object, usually used to find internal
anomalies caused by differences in the physical parameters. The fresnel volume tomography method was an alternative method
for reconstructing the image of an object using travel times, where this method did not use ray path in its calculations. In
calculating the fresnel zone, the finite difference method which was the solution to the equation of wave propagation was used.
The frequency of the waves was also considered, here we used ultrasonic waves with a frequency of 1 Mhz for inverse modelling
process. The reconstruction algorithm we used was Modified Simultaneous Iterative Reconstruction Technique for fresnel
volume. The application of the inversion of synthetic data resulted an estimation of velocity model that have good imaging
quality and similarity with the synthetic model.
Keywords: Tomography, Ultrasonic, Fresnel, Core
.....
40

Constrained Global Optimization of 1D Horizontal Electric Dipole Time Domain
Electromagnetic Method (HED TDEM) using Multi-Objective Particle Swarm Optimization
(MOPSO)
I. G. P. Fadjar Soerya Djaja*, Enjang J. Mustopa and Wahyu Srigutomo

Physics of Earth and Complex Systems, Physics Department, Faculty of Mathematics and Natural Sciences, Institut Teknologi
Bandung, Jl. Ganesa No.10, Bandung, Indonesia
*corresponding author’s email : fadjarsoerya@gmail.com
ABSTRACT
Similar to other geophysical methods, the inversion of TDEM data also suffers from inherent problem where there are many
models that can describe the observed data. This nonuniqueness is generally solved by adding a constraint function multiplied
by a constant called Lagrange multiplier to the minimized objective function. The consequence of this is the dependence of the
model obtained on the value of the Lagrange multiplier used, so finding a suitable Lagrange multiplier is a challenge in itself.
In this study, rather than minimizing one objective function consisting of a misfit function plus a constraint function, we utilized
MOPSO to minimize the two objective functions simultaneously. MOPSO utilizes the concept of Pareto efficiency to obtain the
approximation of Pareto front, which is a set of models that are Pareto optimal. We used the log natural value of the misfit
function and the commonly used roughness function as our objective functions. MOPSO succeeded in mapping the Pareto front
approximation of the two objective functions that we used. From the set of models of the approximated Pareto front, we then
calculated the median and both lower and upper quartile of the models to obtain the most optimum model and the boundaries of
uncertainties of the model.
Keywords: Inversion, Optimization, TDEM, MOPSO, PSO
.....

Identification of Andesite Distribution in Pasrepan Sub-District, Pasuruan District Using
Schlumberger Configuration Resistivity Method
Flavio Da Conceição Martins1, Almer Harahap1*, Mariyanto1, Nita Ariyanti1, Mochammad Fauzan
Dwiharto2, Ahmad Wafi Irsyad2
1
Department of Physics, Faculty Science and Data Analytics, Institut Teknologi Sepuluh Nopember, Surabaya 60111, Indonesia
2
CV. Safari Alam Teknik, Surabaya, Indonesia
*corresponding author’s email : almerharahap19@gmail.com
ABSTRACT
This research was conducted to determine the distribution of andesite layers and the depth and thickness of the layers in
Ampelsari Village, Pasrepan Sub-District, Pasuruan District, East Java Province. The benefits of this research can provide an
understanding of the subsurface section based on the resistivity value obtained from the results of the resistivity geoelectrical
measurement method, and as primary data for identification of andesite rock distribution. The research method uses a vertical
electrical sounding (VES) Schlumberger configuration to obtain a vertical distribution of values to the subsurface with 9
measurement points. The measured data is then input into the IP2win software and inverted until the error value is < 10. Based
on the results and interpretation of the data, it can be concluded that the andesite rocks in the research location, precisely in
Ampelsari Village, Pasrepan Sub-District, Pasuruan District, have a resistivity value distribution of 253-1048 Ω.m with a depth
of 14.8-59.1 m and a thickness of 10-45.6 m.
Keywords: Schlumberger Configuration, Resistivity, Andesit
41
Identification of Aquifer Layer Using Geoelectrical Resistivity Method With

Schlumberger Array in Bumi Harapan Village, New Capital City of Indonesia
Fahmi Saputra,1 Febrian Dedi Sastrawan, S.Si., M.Sc.,2 Meidi Arisalwadi, S.Si., M. Si.,3
1
Institut Teknologi Kalimantan, Indonesia
2
3
*corresponding author’s email : Febrian.dedi@lecturer.itk.ac.id
ABSTRACT
The Movement of National Capital Of Indonesia (IKN) from Jakarta, on Java Island to East Kalimantan induce additional
residents in the area, in particular Bumi Harapan Village, Sepaku sub-district. The increasing of additional residents will give
an impact to the water use on the environmental of resident. Water use at the Bumi Harapan village dominated by shallow water
well excavation, which is located near from the surface by the average range from 4-9 meters, while in the dry season the well
will experience drought. Based on the reasons, the research was conducted to perform any information about the distribution of
aquifer layer as preliminary survey on the area, based on depth and productivity of the aquifer layer. Geo-electrical Resistivity
Method with Schlumberger array was run at 5 measurement points. This method generally can be used to understand the
subsurface lithology based on the distribution of rock resistivity values. Data were obtained from the acquisition field, which is
consisting of potential differences, electricity current and apparent resistivity rock. Field data was processed by using forward
and invers modeling default from software progress that produce log resistivity models which is used to explain the subsurface
lithology. The analysis result shows that lithology structure of subsurface layer were identify topsoil on the surface, sand, sandy
clay, and clay with by the resistivity values range is 1.91 – 228.18 Ohm-m which is located in 69.69 – 112.32 m from the surface
on each measurement point. Based on the results are obtained at fifth measurement points, the resistivity value around 228.18
Ohm-m below the surface is interpreted as sand layer located in 112.32 m depth. The formation factor in this research using 18
well points with analysis results of the nearest well from the fifth measurement point is well point 17. The formation factor at
the well point 17 reveal the value is 4.7, indicated the type of aquifer soil was classified as good productivity.
Keywords: Aquifer, Formation Factor, Schlumberger Configuration, Lithology, Resistivity.
.....

Characterization of Agricultural Land Using Microtremor Approach
Edlyn Yoadan Nathania, 1,2* T. A. Sanny,1*

1
2
*corresponding author’s email : edlynathania02@gmail.com
ABSTRACT
Indonesia is a strategic area for agriculture because it is in the tropics and very suitable for developing agriculture compared
to subtropical regions and other regions in the world. Agricultural productivity is strongly correlated with the soil physical
characteristics. Generally, the characterization of agricultural land is carried out in a traditional way, but in recent decades
geophysical methods have been used in agriculture mainly to determine the characteristics of agricultural soils based on their
physical properties. In this study, geophysical measurements have been carried out using the microtremor method to obtain
information about the structure of the soil layer in the Cilembu sweet potato farmland which is integrated with geoelectric data
and other conventional soil characterization methods.
Keywords: agricultural geophysics, microtremor, soil characterization, subsoil structure
.....
42
WED-03 Zoom 3/ SEACG-048

Water Pipeline Leak Identification using Ground Penetrating Radar Method
Audi Dian Bakti,1, Mariyanto,2 Eko Minarto,1 Tigar Ade Atma Wijaya2,*
1
Department of Physics, Faculty Science and Data Analytics, Institut Teknologi Sepuluh Nopember, Surabaya 60111, Indonesia
2
Surabaya 60111, Indonesia.
*corresponding author’s email : wijayatigar@gmail.com
ABSTRACT
The Ground Penetrating Radar (GPR) method, commonly called Georadar, as a non-destructive geophysical investigation
method that uses electromagnetic waves. This method is widely used to determine conditions under the surface such as mapping
structures and subsurface lithology or buried objects. In this study, the main objective is GPR measurement using Geoscanners
GSSI 270 MHz antenna is conducted to identify and extract features from GPR images of the pipeline leakage. The GPR data
were initially processed using matGPR based on MATLAB. GPR data on each measuring line is processed first, a series of
filters and gaining have been conducted such as Adjust Signal Position, Remove DC, Dewow, Inverse Amplitude Decay,
Remove Global Background, Karhunen – Loeve Filter, then time to depth conversion with adjusted parameters. These results
show that the utilization filter can highlight the anomalies of the pipeline leakage, in an intuitive manner, which can be identified
by non-highly qualified personnel, using the 2D models. From the result of parameter analysis in processing data, found that
the GPR data can be interpreted which show the underground pipeline leak in the study area with the presence of diffraction
signatures(anomalies).
Keywords: Ground Penetrating Radar, matGPR, Pipeline
.....
WED-04 Zoom 3/ SEACG-051

Application of Crosshole Seismic Tomography to Inferred Cavities Condition
Fatkhan1*, Riskiray Ryannugroho2, Fernando Lawrens1, Andri Dian Nugraha3, Mohammad Rachmat
Sule1
1
Exploration and Engineering Seismology Research Grup, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung,
Ganesha 10, Bandung 40132, Indonesia
2
Geophysical Engineering, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Ganesha 10, Bandung 40132,
Indonesia
3
Global Geophysics Group, Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Ganesha 10, Bandung 40132,
Indonesia
*corresponding author’s email : dr.fatkhan@gmail.com
ABSTRACT
The construction of the Rajamandala hydroelectric power plant civil works has experienced critical problems that arise in the
implementation process, with the occurrence of cracks cavities conditions in the headrace tunnel locations. We applied the
crosshole seismic tomography to detect important heterogeneities and mechanical proprieties of the formations between two
boreholes in this case inferred cracks cavaties condition. The seismic source is inside the borehole "a sparker" and the receivers
"hydrophones" are in an adjacent one, the sparker is lowered into the borehole one step at the time. Seismic waves, propagating
in the tunnel wall rock, spread widely and are reflected and refracted when encountering interfaces between rocks with different
acoustic impedances. Reflected waves returning to the receivers are recorded. After the processing of the first arrivals of the P
waves we obtain a cross section of the seismic velocities in between the two wellbores. With the crosshole seismic tomography,
the basis could be provided to the tunnel construction and parameter adjustments to guarantee the construction safety. Our result
shown that there is a weak zone above the tunnel with low Vp zone (~1.2 km/s) may be related to weak zone (may be fracture,
unconsolidated rock, and fluid-filled rock or landslide caving). Checker-board resolution test is also conducted in order to
determine the tomogram areas that are reliable to be interpreted. Our challenges are we have poor resolution due to not enough
raypath. The Crosshole seismic tomography can effectively and safely guide the excavation of the tunnel section working surface
in combination with reconstructed images and excavation technology.
Keywords:. Tunnel; Caving; Seismic Image; Geological Conditions Reflected Waves, Crosshole Seismic Tomography
43
Reconstruction of Retrogressive Landslides Using Integrated Technique of Geoelectric Resistivity

and Borehole Logging in Sumampir, Purbalingga Regency, Central Java
Imam A. Sadisun,1,2,* Indra A. Dinata,1 Gunawan Wibisono,3

1
Applied Geology Research Group, Faculty of Earth Sciences and Technology, Institut Teknologi Bandung, Indonesia
2
Research Centre for Disaster Mitigation, Institut Teknologi Bandung, Indonesia
3
Public Works and Spatial Planning Department of Purbalingga Regency, Central Java, Indonesia
*corresponding author’s email : iasadisun@itb.ac.id
ABSTRACT
The landslide occurred in Sumampir of Purbalingga Regency has reached 120 m in length and 50 m in width. This landslide
is suspected to be a retrogressive landslide, which was manned by a landslide associated with the cliffs of the Brangkal River.
This landslide is interpreted as progressive development of a landslide associated with Brangkal river bank failures, which
continues to enlarge to the dimensions as measured currently. The main active landslide has a steep scarp of along 60 m with
a height of up to 8 m, with several transverse tensile cracks were found in between the scarp and above the crown. These
tensile cracks have caused subsidence on the side of the village road that crosses it. In order to determine the existence of slip
surface that have developed since the first landslide formed until the last one, an integrated technique of geoelectric resistivity
and borehole logging has been applied. This integrated technique was then used to reconstruct the landslide progression
retrogressively towards the top. The results of the geoelectric resistivity analyses indicated that there are four stages of slip
surface formation. The first slip surface is located at 80 m south of the village road, which is defined to be the initial of a
retrogressive landslide. The second slip surface located at 40 m from the village road is possible to be a former landslide
damaged the prior village road. The third slip surface is formed at a distance of 15 m from the village road, which is a landslide
considered to have caused the bridge to collapse. The fourth slip surface is located at the side of the village road, which is the
last landslide causing the concrete slab of the bridge moved 15 m to the south. The average depth of the slip surface is about
8 m. This depth has been confirmed using borehole logging data, by obtaining a relatively hard weathered claystone layers at
a depth of 6-8 m.
Keywords: retrogressive landslide, geoelectric resistivity, borehole logging, slaking, claystone
.....
Geophysical Prospecting of Contamination Pathways in Zeida Abandoned Mine (Morocco)
Mohamed Rouai,1* Abdelilah Dekayir,1 El Mehdi Benyassine,1 Jean C. Parisot,2 Ahmed Lachhab,3
1
Moulay Ismail University, Faculty of Sciences, BP. 11201, Zitoune, Meknes, Morocco
2
CEREGE, Europôle de l’Arbois - BP 80 – 13545, Aix En Provence, France
3
Susquehanna University, Earth and Environmental Sciences, Selinsgrove, USA
*corresponding author’s email : m.rouai@umi.ac.ma
ABSTRACT
The abandoned Zeida mine is located near the city of Midelt, in High Moulouya watershed between the Middle and the High
Atlas of Morocco. After its mining activity, huge tailings are left and stored either in large mining pit lakes with clay-
sandstone substratum or directly on granite bedrock.
The high contents of lead and arsenic in these tailings have transformed them into potential sources of pollution that disperse
and could reach the aquifers through faults and fractures.
In this work, the main goal is to identify the pathways of potential contaminated waters to the local aquifers and Moulouya
River.
Geophysical surveys, including mainly Electrical Resistivity Tomography (ERT), were carried out over the tailings, and
directly on the substratum outside the tailings.
The ERT profiles have successfully shown the location of geological fractures, some of which extend throughout the upper
formation to depths reaching the fractured granite.
The obtained results have shown that pollutants may percolate through fractures and reach the complex hydrological system
connecting groundwater, ponds and the river.
Keywords: Abandoned mine, Tailings, Geophysical prospecting, Zeida, Morocco.
44
Zoom 4 – Seismic Hazards (1) Session 3
Early Analysis of Vulnerability Study in Lake Toba
Dzikri N Rochim1*, Nicolas Silaen1, Riyan Yusami1 Sandy K Suhardja1, Soni Satiawan1, Dicky A
Zaky1, M. Ramdhan2
1
Faculty of Exploration and Production Technology, University of Pertamina, Jl. Teuku Nyak Arief, Simprug, Kebayoran Lama, Jakarta
12220, Indonesia
2
*corresponding author’s email : Dzdood11@gmail.com
ABSTRACT
Lake Toba is located on the Sumatra Island, Indonesia. The subduction of the Indo-Australian plate dominates the tectonic
setting in this region underneath the Eurasian plate. This area has active seismicity, especially from the subduction process
in the West. Most rocks in Lake Toba are tuff sediments. This study's objective is to map the vulnerability index that can be
used to identify high damage areas in the region if an intense event occurs. We use 26 stations 3-component seismometers
installed around Lake Toba. We gathered good quality recording data for one day on each station. Next, we divided our data
into 1 hour and performed the time to frequency domain using the fast Fourier transform. Lastly, we divided the horizontal
and vertical components and extracted some parameters such as dominant frequency value and its amplitude. The dominant
frequency value ranged from 0.18 to 27.6 Hz, H/V amplification was found to be between 2.84 and 11.9, and shear wave
velocity (Vs30) was found to be between 240 and 760 m/s, and the soil vulnerability index (Kg) was found to be between
0.42-43.4. The study area's highest vulnerability is in the southeast, whereas the study area's lowest vulnerability is located
in the West and Northeast.
Keywords: Amplification, Frequency, Toba, Vulnerability
.....
Development of a Land Subsidence Fragility Curve in Bandung Basin: a Preliminary Result
Irwan Gumilar1,*, Riantini Virtiana2, Teguh P. Sidiq1, Brian Bramanto1, Ramdan Sudarsono1, Titis
Idam Haq2
1
Geodesy Research Group, Institut Teknologi Bandung (ITB), Indonesia
2
Remote Sensing and Geographic Information Science Research Group, Institut Teknologi Bandung (ITB), Indonesia
*Corresponding author’s email: igumilar@itb.ac.id
ABSTRACT
Land subsidence has been historically detected and observed using the Global Positioning System (GPS) observations since
the year 2000. At present, the Interferometric Synthetic Aperture Radar (InSAR) method is used to monitor the ongoing
subsidence in the Bandung Basin. Due to its significant subsidence rate, reaching up to -20 cm/year, much evidence of its
impact has been found in this area, e.g., damaged buildings. This study aims to develop a land subsidence fragility curve for
damaged buildings, which is very beneficial for estimating economic losses and developing a risk map of land subsidence.
Field surveys were carried out in areas that experienced large land subsidence. More than 700 buildings were surveyed to
estimate the damage to buildings due to land subsidence. Damage classification is divided into four categories: heavy, medium,
low, and no subsidence. The classification is based on subsidence slope and cracks in the building. Survey results showed that
damage occurred in all these categories. The modeling of the fragility curve shows the tendency of the increasing subsidence
slope to cause the damage will be more severe. However, this study only shows the preliminary result of a land subsidence
fragility curve in the Bandung Basin and will be further investigated for optimal results.
Keywords: GPS, InSAR, fragility curve, damaged building, Bandung Basin
45
How does A Fast-Cascading Earthquake in a Multi-Scale Fracture Network Affect Ground

Motions?
Kadek Hendrawan Palgunadi1, Alice-Agnes Gabriel2,4, Dmitry Garagash3, Thomas Ulrich2, Paul
Martin Mai1
1
Physical Science and Engineering, King Abdullah University of Science and Technology, Thuwal, Saudi Arabia
2
Department of Earth and Environmental Sciences, Geophysics, Ludwig-Maximilians-Universität München, Munich, Germany
3
Dalhousie University, Department Civil Resource Engineering, Halifax, Canada
4
Institute of Geophysics and Planetary Physics, Scripps Institution of Oceanography, University of California, San Diego, CA, USA
ABSTRACT
A cascading earthquake has been shown to develop within a damage zone, involving a multiscale fracture network. A zigzagging
rupture progression among fractures produces complicated seismic wave propagation, generating unpredictable ground motion
distribution. Previous works are restricted to ground motion distribution made by a single or only a few fault segments without
accounting for the contribution of a fracture network. Therefore, the effect of a multiscale fracture network on the resulting
ground motion signatures remains unknown. These raise questions: (i) what is the contribution of each fault segment to the
overall ground motions in the near-field and far-field regions? (ii) How peak ground acceleration (PGA) and peak ground
velocity (PGV) distribution are affected by a cascading earthquake?
We conduct fully-coupled 3D dynamic rupture simulations with seismic wave propagation using open-source software SeisSol.
For far-field simulation, we use kinematic rupture simulation. The model comprises two conjugate fracture families striking
N20E and N120E with nearly vertical dips embedded with a listric fault.
Based on the simulations, each slipped fracture contributes to the overall ground motion at different times, producing higher
frequency content than the non-cascading rupture. Moreover, differences in the arrival time for each fracture create a coda-wave
signature for both near-field and far-field regions, although the model is homogeneous and elastic-isotropic. The PGA and PGV
distribution of a cascading earthquake significantly contrasts with the non-cascading earthquake. The cascading rupture
produces PGA and PGV distribution evenly to almost in all azimuths with higher values toward the directivity direction of the
fracture families. In contrast, the ground motion of the non-cascading rupture dominates the hanging-wall side of the listric fault.
Therefore, accounting for the contribution of a favorably oriented fracture network during seismic hazards assessment is
essential to estimate ground motion distribution near the fault damage zone or subsurface reservoirs.
Keywords: -
....

Calculation of Logarithmic Ratio of Energy-Seismic Moment, Rupture Duration, and Mwpd for
the Tsunami Early Warning System in Bali and Nusa Tenggara
Muh. Soekarno Saputra Rahman1,2*, Nanang T. Puspito3, Pepen Supendi2

1
Master Program in Earth Science, Faculty of Earth Sciences and Technology, Institut Teknologi Bandung, Indonesia
2
Agency for Meteorology, Climatology, and Geophysics (BMKG), Indonesia
3
*corresponding author’s email : soekarnotpi@gmail.com
ABSTRACT
Bali and Nusa Tenggara regions have high seismic activities due to subduction between the Indo-Australian Plate and the
Eurasian Plate, as well as several active faults around those regions. Subduction zone and active faults on the ocean floor could
trigger an earthquake thus generating a tsunami. This research aims to analyze the logarithmic ratio of the seismic moment (Θ),
rupture duration (Tdur), and Mwpd on strong to major earthquakes in the Bali and Nusa Tenggara regions. The data used in this
study are waveforms from the Indonesian Agency for Meteorology, Climatology, and Geophysics (BMKG) seismic stations
which recorded Mw > 6 earthquakes that occurred in Bali and Nusa Tenggara between 2010 and 2020. Our results indicate the
earthquakes that could generate tsunamis have Θ < -4.2, Tdur ≥ 111.65 s, and Mwpd ≥ 7.0, while earthquakes that could not
generate tsunamis have Θ ≥ -4.3, Tdur ≤ 101.47 s, and Mwpd ≤ 7.0. The difference between Mwpd from this study and Mw released
by BMKG is about 0.2 with a correlation value of 40%. Θ, Tdur, and Mwpd values could become supporting parameters for the
tsunami early warning system in Bali and Nusa Tenggara.
Keywords: -
46
....
Seismic Hazard Study of the 2018 Situbondo Earthquake for Situbondo and Its
Surrounding Cities
Sindi Hajah Patimah1, Wahyu Triyoso2, Endra Gunawan2
1
2
Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Bandung Institute of Technology, Indonesia
ABSTRACT
The 10 October 2018 Situbondo earthquake occurred around a city with a relatively high population density. These cities include
Surabaya, Probolinggo, Malang, Madiun, Kediri, and Sidoarjo. The Seismic Hazard Function (SHF) can provide the level of
seismic hazard information at each of these locations. Therefore, it can be used as a basis for hazard mitigation in these cities,
both structurally and non-structurally. This earthquake’s thrust fault was successfully modeled using the coseismic Global
Positioning System (GPS) data. This study aims to estimate the seismic hazard curve in Situbondo city and its surroundings
using the Probabilistic Seismic Hazard Analysis (PSHA) method. First, we calculated the fault slip rate using the empirical
function, and then the seismic hazard function was estimated. This result is then projected to the cities around Situbondo city
cities with high population density. The result suggests that Probolinggo and Surabaya city have higher seismic hazard potential
because the distance between these cities to the earthquake’s epicenter is closer. This result shows that in 500 years, with a
Probability Exceedance (PE) of 10%, at least an earthquake with Peak Ground Acceleration (PGA) of 5.2 gals in Probolinggo
and 1.2 gals in Surabaya will be felt.
Keywords: 2018 Situbondo earthquake, GPS, PSHA, SHF
47
Zoom 5 – Geomagnetic Session 3

Geomagnetic Data Analysis to Identify Subsurface Condition of Gondang Area, Bojonegoro
Yolanda Mustika Bohal1, Mariyanto Mariyanto1, Ihsan Fermantono1,*

1
*corresponding author’s email : fermantonoihsan@gmail.com
ABSTRACT
Manifestations of hot springs and mud pools in Gondang area, Bojonegor, East Java indicates certain geological condition.
Those manifestations are Banyu Kuning Hot Spring, Selo Gajah Hot Spring, and Jari Mud Burst. This study was observed near
Mount Pandan and conducted to find out determine distribution of magnetic fields and magnetic anomalies as well as
distribution patterns of potential geothermal energy in the study area based on its magnetism and to know the rock structure
based on susceptibility value in the manifestation area. The data were collected with 94 measuring points randomly dan
dominantly in manifestation area. The collected data have been corrected by IGRF correction and diurnal correction, then data
analysis was processed by using Surfer 11 to get contour map and Oasis Montaj to get the 3D models. The processed result
revealed that the value of the residual anomaly is -232.1 nT to 159.4 nT. With this residual anomaly, 2D and 3D models showed
low susceptibility distribution values in the range of -0.0298 to -0.0135 cgs around the manifestation area. While the high
susceptibility value has a value range of 0.0114 to 0.0466 cgs which is considered as an intrusion rock.
Keywords: Manifestation, Geomagnetic Anomaly, Gondang-Bojonegoro, Susceptibility.
.....

Characterization of Caldera Lakes’ Sediments Through Rock Magnetic and Geochemical
Analyses. Case Study: Lakes Buyan and Tamblingan in The Island of Bali
Ni Komang Tri Suandayani, Satria Bijaksana, Darharta Dahrin, Irwan Iskandar , Putu Billy
Suryanata, Silvia Jannatul Fajar
Faculty of Mining and Petroleum Engineering, Institut Teknologi Bandung, Bandung, Indonesia
ABSTRACT
Environmental changes as recorded in lake sediments might be due natural as well as anthropogenic causes. Combined rock
magnetic and geochemical analyses have been used intensively to identify anthropogenic as well as natural (lithogenic)
components in lake sediments. In this study this combined analyses were carried out on sediments from two adjacent caldera
lakes in Bali namely Lake Buyan and Lake Tamblingan. Geochemical analyses include measurements of heavy metals as
well as REE (Rare Earth Elements) content while magnetic analyses include magnetic susceptibility measurement. Surface
sediments were collected from 20 sampling points in Lake Buyan and from 16 sampling points in Lake Tamblingan. Rock
samples were also collected from 8 outcrops in the vicinity of these two lakes. Water samples were also collected from 3
different locations in each lake. Preliminary results would be presented.
Keywords: Lake Buyan, Lake Tamblingan, lithogenic component, anthropogenic component, Bali.
48

REEs Content and Magnetic Properties in Placer Deposits from Central Kalimantan
Gusfrimanuel, Satria Bijaksana, Putu Billy Suryanata, Silvia Jannatul Fajar
ABSTRACT
Central Kalimantan has been known to contain placer deposits that are rich in valuable heavy minerals such as zircon,
ilmenite, and rutile. These heavy minerals are probably associated with rare Earth elements or REEs which are currently
considered as strategic resources due to their uses in futuristic products such as magnets, batteries, electronic components
etc. In this study sieved samples taken from placer deposits from Central Kalimantan, that include gold mine tailings and
laterites from weathered granite, are subjected to geochemical as well as magnetic analyses to identify correlation between
heavy minerals, REEs content and magnetic properties such as magnetic susceptibility. Such correlation has been shown
earlier in riverine and lacustrine sediments elsewhere.
Keywords: REE, place deposits, magnetic susceptibility, heavy minerals, Central Kalimantan
.....

Inverse modeling using geomagnetic data of field camp geophysics in Karangsambung (2005-
2019)
Eko Januari Wahyudi1*

1
Institut Teknologi Bandung - Indonesia
*Corresponding author: ekojw@gf.itb.ac.id
ABSTRACT
Geomagnetic data from field camp surveys in Karangsambung 2005-2019 were collected and merged. Before 2020 (the
restriction of Covid-19 pandemic), there are 4713 and 5984 data of field and base observation in the study area. The purpose of
works in this paper are presenting the data collection, simple data processing, and simple calculation of inverse modeling.
Fifteen surveys from each year of field camp data acquisition are merged with 2017 survey as reference. The data observation
from proton magnetometer are corrected with geomagnetic regional field and diurnal variation (using available data from base
observation). Geomagnetic model for 2.5D is calculated using SW-NE slice sections with topographic variation in the study
area. The program for inverse model calculation was built to recover distribution of magnetization contrast from surface
geomagnetic anomaly. The interpretation of subsurface model should be able to be analyzed and correlated with rock
susceptibility and geological surface maps in the study area. From this work, the values of geomagnetic anomaly map are shown
in the range of -1600 to 1600 nanoTesla, and the inverse modeling conducted over the Volcanic Breccia in the part of Bukit
Brujul area, Basalts and Diabas between Geopark Karangsambung-Parangan-Sadang area, and Clay Breccia in the Sadang-
Totogan area.
Keywords: inverse modeling, geomagnetic, magnetization, Karangsambung
49

Magnetic Basement Modelling around Misool Island Raja Ampat Regency with MagB_Inv
Syaiful Apri Kurniawan, Satria Bijaksana, Silvia Jannatul Fajar, Putu Billy Suryanata
ABSTRACT
Misool Island is in the southwest of Papua, administratively included in the Raja Ampat Regency, West Papua Province. Misool
Island is bounded by the Sorong Fault and the Seram Trough. The limestone clusters form hills in the south and lowlands in the
north of Misool Island, resulting from uplifting bedrock characterized by regional anticlines in the southern part of Misool
Island. Judging from its stratigraphic arrangement, Misool Island has several hydrocarbons and mineral economic potentials.
Modelling of subsurface structures is significant for hydrocarbon and economic mineral exploration activities. The magnetic
method is a commonly used method to determine subsurface structures. By exploiting the natural magnetism of rocks, magnetic
bedrock can be mapped. In this study, inversion modelling was carried out to obtain a magnetic basement model from several
sedimentary basins located in the waters of Misool Island. The inversion modelling is carried out with the MagB_Inv, which
can model the crystalline magnetic basement and magnetic anomalies quickly and with high accuracy. This program processes
magnetic data, input parameters, frequency, and convergence criteria by utilizing the Fourier transform to produce a basement
relief interface. The data used are magnetic anomaly data and bathymetry data of Misool Island. From the results of the inversion
modelling, it can be concluded that the average depth of the magnetic basement in the waters of Misool Island is 5.96 km. The
basement is Malihan Ligu rock which belongs to the metamorphic rock and igneous rock. 3D modelling shows that the relief
of the basement in the waters of Misool Island has basins and ridges. These basins include the East Obi Basin, Salawati Basin,
and the Kofiau Basin.
Keywords: Misool Island, magnetic method, MagB_Inv, magnetic basement, West Papua
50
Zoom 1 – Seismology Session 4

Preliminary Study Of Hyposenter Relocation Distribution Using Double-Difference
Method In Bayah Dome Area, Lebak, Banten
Agung Sabtaji1, Hartanto1, Fitri Afiadi1
1
Indonesian Agency for Meteorological, Climatological, and Geophysic Regional II, Abdul Gani Street No.5, South Tangerang, 15412,
Indonesia.
ABSTRACT
Based on BMKG records, the Bayah Dome during the period 2009-2017 had a very low frequency of earthquake activity at very
shallow depths. However, in mid-2018 and late 2020, there was an increase in very shallow depth seismic activity in Bayah
Dome, namely around Bayah and Cilograng districts. The cause of the earthquake source in the area has not been mapped properly
so an approach is needed with an analysis of the distribution of earthquake events resulting from relocation from earthquake
events that have occurred. The relocation method used in this paper uses the double-difference method and uses the earthquake
arrival time catalog data from the BMKG Region II South Tangerang. The results of this relocation showed an improvement in
the depth of the earthquake, which was dominated by a depth of 10 km to be more varied in depth. In addition, it was also found
that there was a clear delineation or lineament trending northeast – southwest to the east of Bayah and also the alleged delineation
to the west of Bayah. So it can be assumed that this delineation is an active fault that has not been previously mapped around the
Bayah Dome.
Keywords:
.....

Preliminary Results: Determination of Microseismic Event Locations on Anisotropic
Medium Using Time Reverse Modeling
Syamsurizal Rizala1, Awali Priyono2, Andri Dian Nugraha2, Mochamad Apri3, Mochamad Agoes
Moelyadi4, Jamhir Safani5
1
Graduate Program of the Geophysical Engineering Department, Faculty of Mining and Petroleum Engineering, Institute of Technology
Bandung, Jalan Ganesha No.10, Bandung 40132, Indonesia
2
Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institute of Technology Bandung, Jalan Ganesha
No.10, Bandung 40132, Indonesia
3
Industrial and Financial Mathematics Research Group, Institute of Technology Bandung, Jalan Ganesha No.10, Bandung 40132,
Indonesia.
4
Department of Aerospace and Aeronautical Engineering, Institute of Technology Bandung, Jalan Ganesha No.10, Bandung 40132,
Indonesia.
5
Department of Geophysicsl Engineering, Faculty of Earth Sciences and Technology, Universitas Halu Olleo, Kampuss Bumi
Tridharma Anduonohu Kendari, Indonesia
ABSTRACT
Determining the microseismic event location is crucial in various fields of science such as hazard mitigation, exploration of new
fossil energy sources, and others. However, in determining the source location, several problems arise, namely the determination
of the source location that is not appropriate due to limited data. To determine the exact location of the event requires a lot of
microseismic recording data. We developed a time reverse modeling method for high-order elastic waves. The data used is
synthetic data that is generated from forward modeling which seems to originate a source that is located in subsurface with a depth
of 1,300 m. The seismic velocity model used is a layered seismic velocity model with the assumption that every layers is
anisotropic elastic medium. Data from the wavefield recording on the surface is propagated back to the source. We assume that
the subsurface layer is heterogeneous isotropic. From the results of the time reverse modeling, it was found that the microseismic
event was at a depth of 1,300 m.
Keywords: Time reverse modeling, microseismic source location, forward propagate, anisotropic medium.
51
Zoom 1 – Seismology Session 4
Hypocenter relocation of eastern Manggarai earthquake 2022 and its aftershocks based
on BMKG seismic network
Mohamad Ramdhan1*, Priyobudi2, Adityo Mursityanto3, Kadek Hendrawan Palgunadi4, Andersen
Lasbayakin Panjaitan2, Jajat Jatnika2
1
Research Center for Geological Disaster, National Research and Innovation Agency (BRIN), Kompleks BRIN, Jl. Sangkuriang,
Kecamatan Coblong, Kota Bandung, Jawa Barat 40135, Indonesia
2
Indonesian Agency for Meteorology, Climatology, and Geophysics (BMKG), Jl. Angkasa I No.2, Kemayoran, Jakarta 10720,
Indonesia
3
Geophysics Station of Bandung (BMKG Bandung), Jl. Cemara No.66, Pasteur, Kec. Sukajadi, Kota Bandung, Jawa Barat 40161,
Indonesia
4
Physical Science and Engineering, King Abdullah University of Science and Technology, Thuwal, Saudi Arabia
*corresponding author’s email: mramdhan123@gmail.com
ABSTRACT
The moderate earthquake (Mw 5.9) struck eastern Manggarai and its surroundings on February 21, 2022, at 12:36:01 UTC.
Even though it was classified as a moderate earthquake, it did produce quite many aftershocks. They could be observed as the
post-seismic activity, which was still well recorded until March 31, 2022. Those hypocenter parameters were updated by
applying the double-difference method. The results show the strike's west-east orientation and the dip's south direction. Those
two fault parameters are linear with the focal mechanism released by the Indonesian Agency for Meteorology, Climatology, and
Geophysics (BMKG). Most depth parameters of earthquakes after hypocenter relocation occurred at a depth of shallower than
25 km, in which the depth of mainshock is 21.5 km. The distribution of aftershocks in this study shows that the earthquake
rupture propagates from the down-dip to the up-dip direction. The relocation hypocenter results depict that most aftershocks
occurred in the complex splay faults, which is different from the main fault plane in which the co-seismic activity triggered
those aftershocks. In addition, this study shows that the system fault consists of some segment faults. Our study benefits
earthquake disaster mitigation, especially in mapping the segment faults of the Flores backarc thrust.
Keywords: hypocenter, relocation, aftershocks, eastern Manggarai
.....
52
Early Analisis on Teleseismic Shear Wave Splitting in Central Java

Eka W. A. R. P., Cicilia1*, Rita Saputri1, Sandy K. Suhardja1
1
Faculty of Technology Exploration and Production, Pertamina University, Indonesia
*corresponding author’s email: ciciliaeka10@gmail.com
ABSTRACT
Central Java is an area that has a high level of seismic activity from earthquakes on land or sea caused by the orthogonal
subduction of the Indo-Australian plate under the Eurasian plate. An understanding of the anisotropy medium is extremely
useful for predicting earthquakes and mapping fractures on Earth. Shear wave splitting is a method that is commonly used in
geophysics and in the interpretation of seismic anisotropy. Analysis was carried out to find and shift data particles and give
information related to delay time and direction polarization fast component from wave slide. The purpose of this study was to
determine the dynamics of plate tectonics in Central Java as well as to interpret and analyze the results of teleseismic data
processing using the shear wave splitting method. The Java region has a seismometer station with teleseismic events and range
time, as well as other parameters that have been determined, making this location study selected in the Central Java region. The
data used is IRIS recording data with 78 SKKS phase teleseismic earthquakes and 75 SKS phase teleseismic earthquakes in the
range of 01 January 2009 to 31 December 2019. From the results of data processing using the method of shear wave splitting
with bandpass filter, 0.0025– 0.25 Hz is obtained as the direction of plate movement in the North-East direction.
Keywords: Anisotropic Medium, Central Java, IRIS, Shear Wave Splitting, Telleseismic
.....
Early Model of Tomography in Timor Using FMTOMO

Nabiilah Rahmawati1*, Sandy K Suhardja1
1
Geophysical Engineering, Faculty of Exploration and Production Technology, University of Pertamina, Jl. Teuku Nyak Arief,
Simprug, Kebayoran Lama, Jakarta 12220, Indonesia
*corresponding author’s email: rahmawatinab@gmail.com
ABSTRACT
The Indo-Australian plate, the Eurasian plate, and the Pacific plate all meet in the eastern half of Indonesia, creating a
complicated tectonic system. Seismic and volcanic activity are caused by the interaction of these plates. This study aims to
analyze the tectonic system in the Timor Area based on the P wave velocity using the ray tracing fast-marching tomographic
method. This study uses two methods, i.e., the double difference method for the hypocenter relocation process and the fast-
marching ray tracing tomography method. The data used in this study came from the International Seismological Center (ISC)
from 1996 to 2018 and the Incorporated Research Institutions for Seismology (IRIS) from 2014 to 2016. The two datasets are
earthquake relocation data for 5970 events (ISC) and 375 events (IRIS). The relocation results show changes in the fixed depth
at shallow depths and the seismicity pattern moves northward to form a subduction pattern due to the subduction of the Indo-
Australian plate towards Eurasia. On tomography results, there are several structural patterns. The low velocity anomaly
imaged with a shallow depth of < 30 km, it is associated with a volcanic arc, while at a depth of 70 -120 km, it is believed to
be partial melting as a source of magma associated with a volcanic arc. High velocity anomalies imaged in a north-south
direction at a depth of ~30 km are associated with back arc thrust (Flores Thrust and Wetar Thrust). The presence of a
subduction slab is imaged by a continuous high velocity anomaly with a depth of up to ~600 km in a south-north direction. In
addition, there is a high velocity anomaly wedged between two low velocity anomalies with a depth of ~80 km, which is
believed to be a forearc sliver.
Keywords: Eastern Half Of Indonesia, Tectonic Sytem, Hypocenter Relocation, Ray Tracing Fast Marching Tomography
53
Early Modelling of Seismic Velocity Structure P, S, and Vp/Vs Ratio in Bengkulu

Region
Nicolas Silaen*, Yusrina A Haqi, Ryandika Syahputra, Amira Saqina, Dzikri N Rochim, M F
Firmansyah, Belista A Natasya, Riyan Yusami, Aulia Rizkiyana, Sandy K Suhardja
Faculty of Exploration and Production Technology, Pertamina University, Jl. Teuku Nyak Arief, Simprug, Kebayoran Lama, Jakarta
12220, Indonesia
*Corresponding author: nicolasilaen@gmail.com
ABSTRACT
Bengkulu is located in the southern region of Sumatra, where the Indo-Australian and Eurasian plates are subducted. As a
consequence of the subduction, seismicity in the area increases, as well as volcanic activity and the appearance of structures
such as the Sumatran fault and the Mentawai fault. Seismicity data from the International Seismological Center (ISC) and the
Meteorology, Climatology, and Geophysics Agency (BMKG) were used to reconstruct seismic velocity Vp, Vs, and Vp/Vs
using seismic tomography in Bengkulu and the surrounding areas. The earthquake data collection period was from January 1,
2010, to December 31, 2018. 19 earthquake stations recorded 1721 earthquakes. LOTOS-12 is used to do a tomography
inversion, starting with a 1D velocity model (AK-135) and Wadati's Vp/Vs ratio, and iteratively inverting for Vp, Vs, and
Vp/Vs ratios. The bending algorithm is used to calculate the arrival time, and the LSQR method is used to invert it. The
Checkerboard Resolution Test (CRT) yielded Tomography results. The results revealed the existence of zones with low
velocity anomalies and low Vp / Vs ratios, indicating the presence of weak zones caused by Sumatra faults, Mentawai faults,
and other local faults. Low-velocity anomalies with high Vp/Vs ratios are related with volcanic activity, whereas high-velocity
anomalies are associated with the existence of subduction plates and compact rocks. The existence of low-velocity anomalies
in the forearc at depths of 30-50 km coincides with fluid circulation caused by slab dehydration.
Keywords: Bengkulu, CRT, Tomography, Vp, Vs, Vp/Vs
.....
54
The Correlation of Land Survey and Satellite Gravity Data, Study Case: Timor Leste
Area
Ahmad Toifur1*, Ziyad Fakhri Kunadi1, Susanti Alawiyah1,2, and Indra Gunawan1,2
1
Geophysical Engineering Department, Institut Teknologi Bandung, Indonesia
2
Applied and Exploration Geophysics Research Group, Institut Teknologi Bandung, Indonesia
*corresponding author’s email : ahmad.toifur1@gmail.com
ABSTRACT
The gravity method has been widely used in tectonic studies and the determination of geological structure. Besides land surveys,
gravity data can also be measured by using satellites. However, the gravity data obtained from the satellite is usually not as good
as the land survey. For example, the satellite data provided by Topex has a resolution of 1-minute grids, while the land survey
could have a higher resolution. It means the gravity data obtained from the satellite cannot be as detailed as the land survey. To
understand the relationship between both data, we calculate the correlation on the CBA, regional, and residual map. We also
compared the maps qualitatively with geological data. During our investigation in the Timor Leste area, it was found that the
correlation value between the CBA, regional, and residual from both data (land survey and satellite) is 0.71, 0.76, and 0.06
sequentially. The CBA and regional of both gravity data have the best agreement, but the residual nearly doesn’t have any
correlation. This result is because the satellite data has 1-minute grids in resolution, which means the spacing between the gravity
data is roughly about 1800 m. With that spacing, the wavelength that can be observed is more than 3600 m, so anomalies in
lower wavelengths usually associated with shallow depths are not good enough to be appropriately measured. However, the
longer wavelengths will be sufficiently provided by satellite. The long-wavelength of gravity data can be used to investigate the
regional anomalies, such as basement and general geologic mapping. Our results also show that the regional gravity anomaly
has high compatibility with geological data.
Keywords: Gravity method, satellite, land survey, correlation, gravity anomaly
.....
Predictive Modeling of Land Subsidence Due to Groundwater Level Decline in

Gedebage District, Bandung, Indonesia
Dwi Sarah1*, Eko Soebowo1, Arifan Jaya Syahbana1, Nugroho Aji Satriyo1, Zulfahmi1, Wahyudin1
1
Research Center for Geological Disaster, National Research and Innovation Agency (BRIN), Indonesia
*corresponding author’s email : sarahpr28@gmail.com
ABSTRACT
The Gedebage district in Bandung city has experienced a vast physical development in recent years. Previous geodetic studies
showed that the area undergoes high-rate land subsidence. The subsurface soil of the Gedebage area has remarkably high
compressibility, void ratio, and water content. A thick sequence of soft organic clay of 20-27 m exists, making the area prone
to subside. As the district continues to develop, it is important to predict the future rate of land subsidence in the area. This
paper aims to analyze the contribution of groundwater exploitation to the land subsidence rate. A combination of 1-dimensional
Terzaghi consolidation analytical analysis and numerical modeling was employed in this study. Groundwater level data up to
2015 was used, and hydrostatic condition was assumed to occur in 1986. Results show that the Gedebage area has subsided
as much as 161 cm since 1986. Modeling results are consistent with geodetic survey results. It is predicted that the land
subsidence will slow down in the next 40 years, provided that the groundwater level remains stable.
Keywords: land subsidence, rate, Gedebage, groundwater, clay, soft
55
Application of Vertical Electrical Sounding Method Using Damped Least-Square

Inversion for Interpretation of Leachate in the Bakung Landfill, Lampung
Rizka1,2, Djoko Santoso2, Warsa2*
1
Geophysical Engineering, Institut Teknologi Sumatera, South Lampung 35365, Indonesia
2
*corresponding author’s email : warsa@gf.itb.ac.id
ABSTRACT
Landfill is a source of potentially hazardous pollution. Solid landfill produces leachate which contains various components
such as contaminants, concentration, and amount of leachate depending on the type of waste disposed of in the landfill. Leachate
is a substance that has seeped from the decay of waste, dissolved or suspended materials originating from the landfill. The
presence of leachate can be identified through a geophysical survey, namely the DC Resistivity geoelectrical method. In this
study, leachate identification was carried out in the Bakung Landfill, Lampung. Measurements were performed using the
vertical electrical sounding method with a Schlumberger configuration to attain preferable depth resolution. The application of
inverse damped least square was used for the inversion modeling. The damped Least Squares modeling was tested synthetic
model which have resistivity variations. Then Damped Least-Squares was applied in the Bakung Landfill. The results obtained
are in the form of a 1D resistivity model map which is used to determine the depth and thickness of the layer that is likely to
be leachate.
Keywords:. vertical electrical sounding, least-square inversion, Bakung landfill
....
56
Early Model on Teleseismic Study in Merapi Merbabu

Sirilus C. Sibagariang*, Sandy Kurniawan Suhardja, Sinta N Asyda, Cecilia E.P, M. Kholil
Universitas Pertamina; Jln. Teuku Nyak Arief, Simprug, Kebayoran Lama, Jakarta Selatan, DKI Jakarta, Indonesia, 12220
*corresponding author’s email: sirilusbagariang02@gmail.com
ABSTRACT
Merapi and Merbabu are two areas in Central Java with a high level of seismicity in Indonesia. Merapi-Merbabu was formed
due to activity between two continental plates, namely the Eurasian plate and the Indo-Australian plate. The movement of
these plates led to the formation of tectonic elements. This study observes seismic anisotropy in the subduction zone using
the shear wave splitting method. The shear wave splitting method is a study to analyze the anisotropy level of a medium. The
data used in this study is 3-component teleseismic data recorded by stations available through the GFZ (Geo Forschungs
Zentrum) network. The data for this study was collected between October 2013 and January 2015. The results of processing
the SWS method are the back azimuth parameter, which shows the direction of polarization and delay time. The results
obtained are different rotation patterns that describe the direction of rapid polarization of anisotropy. This is most likely
influenced by differences in age and the direction of absolute plate movement in Java. The direction of rapid anisotropy
polarization in the Merapi-Merbabu area is dominant in the NS direction, but there are several stations in the Merapi area
which almost point to the EW. This is most likely influenced by differences in age and the direction of absolute plate
movement on the island of Java. The direction of rapid anisotropy polarization in the Merapi-Merbabu area is dominant in
the NS direction, but there are several stations in the Merapi area which almost point to the EW. This is most likely influenced
by differences in age and the direction of absolute plate movement on the island of Java. The direction of rapid anisotropy
polarization in the Merapi-Merbabu area is dominant in the NS direction, but there are several stations in the Merapi area
which almost point to the EW. This is most likely influenced by differences in age and the direction of absolute plate
movement on the island of Java.
Keywords: Subduction, Anisotropy, Shear Wave Splitting, Back Azimuth, Delay Time
....
Preliminary Results on Receiver Function Study in Mt Merapi, Central Java, Indonesia

Ramadanti Fitria and Suhardja Sandy .K
Teknik Geofisika Universitas Pertamina
*corresponding author’s email : fitriamuryati@gmail.com
ABSTRACT
The Receiver Function method is currently widely used in subsurface study because it is simple to obtain and relatively
inexpensive. The Receiver Function method is used to determine the boundaries of discontinuities in earth structures such
as the crust and mantle by relying on passive seismic waves recorded at three component stations. This Receiver Function
study focuses on Mount Merapi, with the objective of imaging subsurface conditions, particularly the presence of a magma
reservoir, using synthetic velocity modeling, which is then processed into Synthetic Receiver Function curves. For the
Receiver Function, the synthetic curve is then correlated with the observation curve. The correlation between the synthetic
RF curve and the RF observation curve was tested in this study until the best correlation was found. The Receiver Function
response must first be analyzed to facilitate the fabrication of the Synthetic Receiver Function model, beginning with the
simplest model and progressing to include a sediment layer with a thickness of 1-3 km, crustal layer, low velocity zone or
low velocity zone interpreted as a magma reservoir, moho discontinuity layer, and the contrast of seismic wave velocity
values, particularly Vs. According to the results of synthetic receiver function curve modeling at the closest station to
Mount Merapi, as well as stations to the west and east of Mount Merapi, the moho depth is estimated to be between 26
and 29 km, with the depth decreasing to the south or towards Mount Merapi. For variations in the depth of the low velocity
zone to the west of Merapi, which is between 5 and 10 km, and the pattern is becoming shallower to the south or towards
Mount Merapi. The depth of the low velocity zone in the east is discovered to be 8-12 km.
Keywords: Central Java, Mount Merapi, Synthetic Receiver Function Curve, Low Velocity Zone, Magma Reservoir
57
Volcanostratigraphy and Identification of Geological Structure of the

Karaha-Sadakeling Colcano Complex Tasikmalaya Regency, West Java
Wilfridus FS Banggur1*, MN Kartadinata2, Sofyan Primulyana2, Estu Kriswati1, Azrie Ezziat Putera
Muhammad3, Gradita Trihadi2
1
National Research and Innovation Agency
2
Centre for Volcanology and Geological Hazard Mitigation
3
Geology Department, Faculty of Earth Science and Technology,Bandung Institute of Technology
*corresponding author’s email: horizondrone99@gmail.com
ABSTRACT
The purpose of this study was to determine the rock sediment and rock distribution both laterally and vertically as well as the
geological structure that developed in the Karaha-Sadakeling Volcanic Complex. Field research and rock sampling as well as
petrographic observation of rock thin sections and analysis of micro and major elements of the Karaha-sadakeling rock product
have been carried out. Volcanostratigraphy, the rocks of the Karaha-Sadakeling Volcanic Complex, respectively, from old to
young, are Khuluk Cakrabuana, Khuluk Sadakeling, Khuluk Karaha, and Khuluk Telagabodas. The petrographic and
geochemical analysis of rocks shows the evolution of magma that occurred at the Karaha-Sadakeling Crater Complex in the
form of fractional crystals, self-mixing due to convection currents, mixing of primitive magmas and decompression during the
eruption which was formed in the subduction zone with alkaline theolytic-calc magma series. Geological structures that
developed in the form of lineaments, fault structures and crater structures. The main fault line is a fault that passes through the
Malangbong area, in a northwest-southeast direction. Volcanic cones that grow in a relatively north-south direction are thought
to have formed connected to the main fault line.
Keywords: Volcanostratigraphy, Karaha-Sadakeling, Structure, geochemical analysis
....
58
Semeru Volcano Eruption Activity in 2021 Based on Seismic and Deformation
Monitoring Data
Kristianto1, Ahmad Basuki1, Heruningtyas Desi P.1, Devy Kamil Syahbana1
1
Center for Volcanology and Geologic Hazard Mitigation (CVGHM), Jl. Diponegoro 57, Bandung 40122, Indonesia

ABSTRACT
Semeru is a volcano with the highest peak on the island of Java and is one of the volcanoes in Indonesia that has a potential
threat of pyroclastic flow. The pyroclastic flow event on December 4, 2021, reached a distance of 16 km from the crater and
resulted in the loss of 51 people. The volcanic activity monitoring system consists of 4 seismic stations, deformation stations
(2 Tiltmeter and 3 GPS stations), and 2 CCTV.
An increased eruptive activity occurred in December 2021 with a series of pyroclastic flows reaching a maximum distance of
16 km on December 4, 2021. This event was triggered by an excess pressure at shallow depths a few hours before the pyroclastic
flows occurred. Deep Volcanic Earthquakes (VA) are relatively increased after the collapse of the lava dome or pyroclastic
flows which possibly caused by the sudden decrease in hydrostatic pressure of the rock mass around the magma pocket and
thus triggers gas expansion.
Deformation monitoring using Tiltmeter at Argosuko and Jawar stations indicates inflation of a deep source since 15 August
2021. Since 7 November 2021 Tiltmeter shows deep source deflation patterns which indicate magma movement from the deep
reservoir to a shallower reservoir. After November 21, inflation of the deep source was observed again indicating an increase
in magma supply from the deep pocket.
The results of the pressure source modeling from GPS vector data in the period before the eruption 1 November – 4 December
2021 showed an anomaly in the form of an increased volume at a depth of > 1.7 km of about 0.84 million/m3. Between
December 5 and December 31, 2021, after the December 4 pyroclastic flow, the deformation modeling indicates the transition
of the pressure source from a depthof 1.7 km to 8 km. In addition, we also observed a decrease in volume of 5.6 million/m3.
The deformation pattern at shallow depths showed a deflationary pattern indicating a decrease in magma overpressure.
Keywords: Semeru Volcano, Pyroclastic Flow, Seismic, and Deformation
.....

Observing Near-Real Time Volcanoes Deformation Using InSAR – Time Series of
LiCSBAS in Indonesia
Erlangga Ibrahim Fattah1,4*, Asep Saepuloh1, Andri Dian Nugraha2, Kristianto3, Ahmad Basuki3
1
Teknik Geologi, Fakultas Ilmu dan Teknologi Kebumian, Institut Teknologi Bandung, Indonesia.
2
Teknik Geofisika, Fakultas Teknologi Pertambangan dan Perminyakan, Institut Teknologi Bandung, Indonesia.
3
Pusat Vulkanologi Mitigasi Bencana Geologi.
4
Teknik Geofisika, Institut Teknologi Sumatera, Indonesia.
*corresponding author: erlanggaif@gmail.com
ABSTRACT
The Synthetic Aperture Radar (SAR) sensors onboard satellite are leading in advanced to monitor physical change of volcano
edifice such as the deformation. Currently, few day temporal resolution imageries provided by Sentinel-1 SAR makes possible
to a near-real time monitoring of volcanoes world widely. , The revisit time for certain region has been improved from twelve
to six days by constellation of Sentinel-1A and B. The Sentinel-1 constellation produces large amount of the SAR images
valuable for monitoring the deformation and hazard mitigation of volcanoes. However, processing the Interferometric SAR
(InSAR) for large amount of SAR data is very much time consuming and necessity of high performance computer. Overcoming
the problem, the COMET (Centre of the Observation and Modelling of Earthquakes, Volcanoes, and Tectonics) creates a
program termed as Looking into Continent from SAR (LiCSAR) to process and calculate the interferogram containing
deformation efficiently and integrate the process into the COMET server. Interferogram images produced by LiCSAR system
have been used to analyze the deformation InSAR time series using Looking into Continent using Small Baseline (LiCSBAS)
software for mapping deformation rate. In this study, we present the LiCSAR observation for volcanic activities at Mts.
Sinabung and Agung located at North Sumatra and Bali, respectively. We have observed the activities of both volcanoes from
2019 to 2022. We have analyzed 641 and 477 interferograms of Mt. Sinabung and Mt. Agung, respectively. The LicSBAS
results showed that the deformation rate of Mt. Sinabung varies for each sector. The east sector showed deflation with velocity
about 12 mm/yr. Meanwhile, the north and west sectors show inflation with velocity about 3.5 - 4 mm/yr. Contrary, the
deformation of Mt. Agung relatively constant with inflation at the summit about 1 - 2 mm/yr.
Keywords: InSAR time series, LiCSAR, LiCSBAS, Mount Sinabung, Mount Agung.
59
Ambient Seismic Noise Cross – correlation of the Tangkuban Parahu Volcano, West Java
Indonesia: Preliminary Result
Aditya Lesmana1*, Awali Priyono2, Andri Dian Nugraha2, Zulfakriza2, Shindy Rosalia2, Ahmad
Basuki3, Kristianto3
1
Geophysical Engineering Study Program, Faculty of Mining and Petroleum Engineering Bandung Institute of Technology, Bandung,
Indonesia
2
Global Geophysics Research Group, Faculty Mining and Petroleum Engineering Bandung Institute of Technology, Bandung, Indonesia
3
Center for Volcanology and Geological Hazard Mitigation, Bandung, Indonesia
*corresponding author: lesmana.aditya47@gmail.com
ABSTRACT
The Tangkuban Parahu Volcanic Complex (TPVC) is one of an active stratovolcano located in West Java, Indonesia. To deeper
understand the recent internal seismic structure of the TPVC, we used ambient seismic noise recorded by 63 temporary seismic
stations deployed around the TPVC during October 2021 until February 2022. The resulting data is a three-component seismic
recording (vertical, north-south, and east-west). In this study, we computed the cross-correlations of the vertical component of
continuous records of ambient noise data to extract empirical Green’s functions. Empirical Green’s function were obtained from
the daily data series and stacked the day-by-day cross-correlation data into one inter-station cross-correlation. The Empirical
Green’s Function is seen at the band period 0.5-10 s. These functions provide detailed images of Rayleigh wave group velocity
for different periods. The group velocity of Rayleigh waves obtained varies from 1.09 km.s – 3.07 km/s. Low group velocity
might be indicated the presence of a weak zone and high group velocity might be indicated volcanic rocks produced by
volcanoes.
Keywords: Tangkuban Parahu, cross-corelation, temporary station, ambient noise seismic
.....
60

Untitled

Uploaded by

Copyright:

Available Formats

Untitled

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Untitled

Uploaded by

Copyright:

Available Formats

SEACG 2022 Program Book

Welcome to SEACG 2022 Head of Geophysical Engineering

Dr. Ir. Fatkhan, M.T.

<< 11.50 - 12.45 SESSION BREAK >>

TIME 12.45 – 13.45 (UTC+7)

<< 13.45 - 14.00 COFFEE BREAK >>

SCIENTIFIC PARALLEL SESSION II

<< 10.15 - 10.30 COFFEE BREAK >>

Wednesday, 10 August INVITED SESSION IV

TIME 12.45 – 14.15 (UTC+7)

The 4th Southeast Asian Conference on Geophysics

SEACG Room Layout

WED, 10 AUG 2022 | Session 3

WED, 10 AUG 2022 | Session 4

TUE-01 Zoom 1 / SEACG-004

Investigation of Postseismic Deformation of The 2018 Lombok Earthquake Sequence

*corresponding author’s email :

*corresponding author’s email : ajengmarina7@gmail.com

Moment Tensor Inversion Implementation In Determining Focal Mechanism Solution

*corresponding author’s email : diann.kusumawati03@gmail.com

TUE-04 Zoom 1 / SEACG-058

*corresponding author’s email : indraandradinata@gmail.com

TUE-01 Zoom 2 / SEACG-030

*corresponding author’s email : amirasaqina1407@gmail.com

Keywords: Cross-Correlation, Phase-velocity, Surface-wave tomography

TUE-02 Zoom 2 / SEACG-16

M.Aldi1*, Nicolas Silaen1, Dayyan Dhaifullah1, Sandy K Suhardja1,2

*corresponding author’s email : aldi.mhd07@gmail.com

TUE-03 Zoom 2 / SEACG-039

*corresponding author’s email : riyanyusami03@gmail.com

Sediment thickness in the eastern part of Kalimantan from high-frequency receiver

*corresponding author’s email : hadda9@gmail.com

Keywords: Receiver function, Kalimantan, Sediment, H-κ stacking

TUE-01 Zoom 3 / SEACG-054

*corresponding author’s email : rayhanfr1@gmail.com

Keywords: Georadar, Old Tram Tracks, Dielectric Permittivity

TUE-02 Zoom 3 / SEACG-032

Steam-Vent Identification Using Frequency Charasteristic Method In Ground

*corresponding author’s email : rarifpambudi@gmail.com

Keywords: Ground Penetrating Radar, frequency characteristic, steam vent, FDTD

TUE-03 Zoom 3 / SEACG-078

Detecting Subsurface Utilities with Multi-frequency and Multi-antenna Configuration

*corresponding author’s email: wahyudi@gf.itb.ac.id

Keywords: GPR, multi-frequency, multi-configuration antenna, PPSDM Geominerba’s Test Site

TUE-04 Zoom 3 / SEACG-079

*corresponding author’s email : wahyudi@gf.itb.ac.id

Keywords: Impedance inversion, Ground-Penetrating Radar, GPR, groundwater level

TUE-01 Zoom 4 / SEACG-017

Saturated Fluid Variation Effect on Seismic Anisotropy: A Case Study of Deep-Water

*corresponding author’s email : madaniyaoktariena@outlook.com

Fault Delineation and Geometrical Modeling of the Geological Carbon

Keywords: Climate-Change, Decarbonization, Geological-Carbon-Storage, Fault, Capacity, Gravity

TUE-03 Zoom 4 / SEACG-035

*corresponding author’s email : rahmanurhazlita00@gmail.com

Keywords: Anisotropy, Core plug, Porosity, Reservoir, Sandstone

TUE-01 Zoom 5 / SEACG-052

*corresponding author’s email : rahmiazizah05@gmail.com

Keywords: VLF-EM, NA-MEMD, Selogajah

TUE-02 Zoom 5 / SEACG-068

*corresponding author’s email : farrasfadhil5@gmail.com