Delineation of Upper Crustal Structure Beneath the Island of Lombok, Indonesia, Using Ambient Seismic Noise Tomography

We have successfully conducted the first ambient noise tomography on the island of Lombok, Indonesia using local waveform data observed at 20 temporary stations. Ambient noise tomography was used to delineate the seismic velocity structure in the upper crust. The waveform data were recorded from August 3rd to September 9th, 2018, using short-period and broadband sensors. There are 185 Rayleigh waves retrieved from cross-correlating the vertical components of the seismograms. We used frequency-time analysis (FTAN) to acquire the interstation group velocity from the dispersion curves. Group velocity was obtained for the period range of 1 s to 6 s. The group velocity maps were generated using the subspace inversion method and Fast Marching Method (FMM) to trace ray-paths of the surface waves through a heterogeneous medium. To extract the shear wave velocity (Vs) from the Rayleigh wave group velocity maps, we utilize the Neighborhood Algorithm (NA) method. The 2-D tomographic maps provide good resolution in the center and eastern parts of Lombok. The tomograms show prominent features with a low shear velocity that appears up to 4 km depth beneath Rinjani Volcano, Northern Lombok, and Eastern Lombok. We suggest these low velocity anomalies are associated with Quaternary volcanic products, including the Holocene pyroclastic deposits of Samalas Volcano (the ancient Rinjani Volcano) which erupted in 1257. The northeast of Rinjani Volcano is characterized by higher Vs, and we suggest this may be due to the presence of igneous intrusive rock at depth

Seismic Structure Beneath the Molucca Sea Collision Zone from Travel Time Tomography Based on Local and Regional BMKG Networks

The Molucca Sea Plate, and Sangihe and Halmahera plates have a complex tectonic setting and interact to create the Molucca Sea Collision Zone. We re-picked 1647 events recorded from 2010 to 2017 from 32 of The Agency for Meteorology, Climatology, and Geophysics (BMKG) stations and obtained P- and S-arrivals of ~17,628 phases. Hypocenter locations were determined using the software NonLinLoc. Then, we performed a travel time tomography in order to image the subsurface and approximate the Molucca Sea Plate subduction angle beneath Sulawesi’s north arm, the relationship subduction zone and volcanic activity in Halmahera, and depth comparison of the Molucca Sea Plate. Our results show that (i) high Vp, high Vs, and low Vp/Vs are associated with the Molucca Sea Plate beneath Sulawesi’s north arm, and the approximate subduction angle is ~40°. (ii) Low Vp, low Vs, and high Vp/Vs beneath the northern and southern Halmahera Volcanic Arc are associated with a possible magma source. Volcanoes in the north have experienced eruptions and are dormant in the south. This group of volcanoes is connected by partial melting above the Molucca Sea Plate subducts to the east. (iii) High Vp, high Vs, and low Vp/Vs are interpreted as double subduction of the Molucca Sea Plate. It is submerged deeper in the north compared with the south, which is nearer to the surface

Seismic Structure Beneath the Molucca Sea Collision Zone from Travel Time Tomography Based on Local and Regional BMKG Networks

The Molucca Sea Plate, and Sangihe and Halmahera plates have a complex tectonic setting and interact to create the Molucca Sea Collision Zone. We re-picked 1647 events recorded from 2010 to 2017 from 32 of The Agency for Meteorology, Climatology, and Geophysics (BMKG) stations and obtained P- and S-arrivals of ~17,628 phases. Hypocenter locations were determined using the software NonLinLoc. Then, we performed a travel time tomography in order to image the subsurface and approximate the Molucca Sea Plate subduction angle beneath Sulawesi’s north arm, the relationship subduction zone and volcanic activity in Halmahera, and depth comparison of the Molucca Sea Plate. Our results show that (i) high Vp, high Vs, and low Vp/Vs are associated with the Molucca Sea Plate beneath Sulawesi’s north arm, and the approximate subduction angle is ~40°. (ii) Low Vp, low Vs, and high Vp/Vs beneath the northern and southern Halmahera Volcanic Arc are associated with a possible magma source. Volcanoes in the north have experienced eruptions and are dormant in the south. This group of volcanoes is connected by partial melting above the Molucca Sea Plate subducts to the east. (iii) High Vp, high Vs, and low Vp/Vs are interpreted as double subduction of the Molucca Sea Plate. It is submerged deeper in the north compared with the south, which is nearer to the surface

Shear wave splitting of the 2018 Lombok earthquake aftershock area, Indonesia

Abstract Lombok is one of the islands in the transitional zone from the Sunda Arc to the Banda Arc, Indonesia. In the mid-2018, the island of Lombok was shaken by a series of strong earthquakes, started with a moment magnitude (Mw) 6.4 earthquake on July 29, 2018 followed by earthquakes on August 5 (Mw 7.0), August 9 (Mw 5.9), and August 19 (Mw 6.3 and 6.9). Some researchers suggested that this phenomenon occurred due to a segmentation rupture in the northern part of Lombok Island. This study aims to obtain information on the distribution of the Lombok earthquake fault zone 2018 and also to understand the character of seismic anisotropy around the Lombok earthquake fault zone 2018 through Shear Wave Splitting (SWS) study. Splitting, or S-wave separation, occurs when the S wave passes through an anisotropic medium. The S wave is split into fast and slow S waves with almost orthogonal polarizations and has parameters such as delay time and polarization direction of the fast S wave. To determine the SWS parameters, we used a Lombok earthquake aftershock data set recorded from 4 August to 9 September 2018, using 16 seismographic stations. The steps taken to obtain the SWS parameters are event selection, windowing using short time Fourier transform, and rotation-correlation process. The results of the SWS analysis indicate that the fast polarization directions probably have a linkage with the local fault system and the fault related to the Lombok earthquake fault zone

Seismic attenuation tomography from 2018 Lombok earthquakes, Indonesia

Local earthquake data was used to determine a three-dimensional (3D) seismic attenuation structure around the aftershock source region of the 2018 Lombok earthquake in Indonesia. The aftershocks were recorded by 13 seismic stations from August 4 to September 9, 2018. The selected data consist of 6,281 P-wave t∗ values from 914 events, which had good t∗ quality in at least four stations. Our results show that the two aftershock clusters northwest and northeast of Lombok Island have different attenuation characteristics. A low P-wave quality factor (low-Qp), low P-wave velocity (Vp), and high ratio of P-wave velocity and S-wave velocity (Vp/Vs), which coincide with a shallower earthquake (20 km) northeast of Lombok Island, might be associated with an area that lacks fluid content. The difference in fluid content between the northwest and northeast regions might be the cause of the early generation of aftershocks in the northwest area. The significant earthquake that happened on August 5, 2018, took place in a region with moderate Qp, close to the contrast of high and low-Qp and high Vp, which suggests that the earthquake started in a strong material before triggering the shallower aftershocks occurring in an area affected by fluid content. We also identified an old intrusive body on the northeast flank of the Rinjani volcano, which was characterized by a high-Qp, high-velocity, and a high Bouguer anomaly.Published versionThis study was supported by “Riset P3MI ITB 2020” awarded to AP. This study was partially supported by “Riset Penelitian Dasar Unggulan Perguruan Tinggi (PDUPT) BP-PTNBH Kemristek/BRIN 2019-2020,” World Class Research Kemristek/BRIN 2021-2022 awarded to AN, and was also partially supported by Center for Earthquake Science and Technology, Research Center for Disaster Mitigation, Institut Teknologi Bandung (CEST, PPMB, ITB)