Geology topics

Andri Dian Nugraha

Publications and source records attributed to Andri Dian Nugraha.

4 recordsLinked to original sources

USGS2020

Relocated aftershocks and background seismicity in eastern Indonesia shed light on the 2018 Lombok and Palu earthquake sequences

High seismicity rates in eastern Indonesia occur due to the complex interaction of several tectonic plates which resulted in two deadly, destructive earthquake sequences that occurred in Lombok Island and the city of Palu, Sulawesi in 2018. The first sequence began in July with an M w 6.4 event near Lombok, culminating in an M w 7.0 event 8 d later. This was then followed by a nearby M w 6.9 event 12 d later. Approximately 1000 km to the northeast, a separate sequence began several weeks later near Palu where an M w 7.5 event occurred that triggered a tsunami. In this study, we present hypocentre relocations for both earthquake sequences as well as all other regional earthquakes in eastern Indonesia. The relocations were performed using a teleseismic double-difference relocation method and arrival times for P and S waves from stations at local, regional, and teleseismic distances. The catalogue and phase data were taken from the Agency for Meteorology, Climatology and Geophysics (BMKG) of Indonesia and the International Seismological Centre (ISC) for the period of April 2009 through November 2018. The relocated catalogue provides an improved view of seismicity in eastern Indonesia over the study period, sharpening locations and interpretations of seismogenic features throughout the region. In the Lombok area, the relocated earthquakes clearly show a backarc thrust to the north of the Sunda-Banda Arc transition zone. The relocated aftershocks show that the destructive M w 7.0 and M w 6.9 earthquakes of the Lombok sequence ruptured two different regions: The M w 7.0 earthquake propagated westward, whereas the M w 6.9 earthquake propagated eastward. The entire sequence of Lombok earthquakes was most likely started by the M w 6.4 event as the initial event or foreshock, which then triggered backarc thrusts on both sides. Several weeks later and far to the northeast, the M w 7.5 Palu earthquake occurred along the Palu-Koro Fault, filling a seismic gap that had not ruptured in an M w 6.0 event or larger since at least 1900. The distribution of aftershocks indicates that the northern part of the Palu-Koro Fault has lower relative seismicity rates than the southern part at shallow depths, and that off fault aftershocks are mostly located to the east of the Palu-Koro Fault.

Lombok, Palu↗

USGS2019

Joint 3-D tomographic imaging of Vp, Vs and Vp/Vs and hypocenter relocation at Sinabung volcano, Indonesia from November to December 2013

We conducted travel time tomography using P- and S-wave arrival times of volcanic-tectonic (VT) events that occurred between November and December 2013 to determine the three-dimensional (3D) seismic velocity structure (Vp, Vs, and Vp/Vs) beneath Sinabung volcano, Indonesia in order to delineate geological subsurface structure and to enhance our understanding of the volcanism itself. This was a time period when phreatic explosions became phreatomagmatic and then magma migrated to the surface forming a summit lava dome. We used 4846 VT events with 16,138 P- and 16,138 S-wave arrival time phases recorded by 6 stations for the tomographic inversion. The relocated VTs collapse into three clusters at depths from the surface to sea level, from 2 to 4 km below sea level, and from 5 to 8.5 km below sea level. The tomographic inversion results show three prominent regions of high Vp/Vs (~ 1.8) beneath Sinabung volcano at depths consistent with the relocated earthquake clusters. We interpret these anomalies as intrusives associated with previous eruptions and possibly surrounding the magma conduit, which we cannot resolve with this study. One anomalous region might contain partial melt, at sea level and below the eventual eruption site at the summit. Our results are important for the interpretation of a conceptual model of the “plumbing system” of this hazardous volcano.

Sinabung Volcano↗

USGS2019

3-D seismic tomographic study of Sinabung Volcano, Northern Sumatra, Indonesia, during the inter-eruptive period October 2010-July 2013

We estimated 3D P-wave velocity (Vp), S-wave velocity (Vs), and Vp/Vs tomographic structures in and around the Sinabung Volcano area, Northern Sumatra, Indonesia during the intereruptive period between October 2010 and July 2013 using high-quality P- and S-wave phases from > 600 local volcano tectonic (VT) earthquakes and more than 5000 phase picks. We used the iteritive damped-least-squares method of SIMULPS12 to simultaneously invert for the velocity structure and hypocenter relocations. We find that the majority of the VT earthquakes locate beneath Sinabung volcano to depths of about 8 km, with distal hypocenters dominantly to the north and northwest of the summit. We find six anomalous regions with the most notable including an area of low seismicity, high Vp/Vs and low Vp at depths of around 5 km slightly south of the summit. We believe this anomaly most likely represents hot material with some partial melt that may be the source for the 2010 phreatic eruptions and continued unrest throughout the study period. We also find a region below the summit with high seismicity, high Vp, high Vs and low Vp/Vs that we believe likely represents old intrusive material that has been fractured by the processes of the 2010 phreatic eruptions. Other shallow anomalous values of Vp, Vs and Vp/Vs are consistent with the known geology of Sinabung: anomalies of high Vp and Vs likely represent old, well-consolidated rock, and are observed to the SW and NE of Sinabung; while low Vp, low Vs and high Vp/Vs likely represent increased fracture densities, temperatures, or gas contents in geothermal features to the west of the summit. Finally we compare our results to those of the tomographic study during the October to November 2013 eruptive period in Nugraha et al. (this issue) and find evidence that suggests magma rose from the lower region imaged in our study to within a few kilometers of the surface.

Journal of Volcanology and Geothermal Research↗

USGS2018

Hypocenter relocation along the Sunda arc in Indonesia, using a 3D seismic velocity model

The tectonics of the Sunda arc region is characterized by the junction of the Eurasian and Indo‐Australian tectonic plates, causing complex dynamics to take place. High‐seismicity rates in the Indonesian region occur due to the interaction between these tectonic plates. The availability of a denser network of seismometers after the earthquakes of M w "> M w 9.1 in 2004 and M w "> M w 8.6 in 2005 supports various seismic studies, one of which regards the precise relocation of the hypocenters. In this study, hypocenter relocation was performed using a teleseismic double‐difference (DD) relocation method (teletomoDD) combining arrival times of P and S waves from stations at local, regional, and teleseismic distances. The catalog data were taken from the Agency of Meteorology, Climatology, and Geophysics (BMKG) of Indonesia, and the International Seismological Centre (ISC) for the time period of April 2009 to May 2015. The 3D seismic‐wave velocity model with a grid size 1 ° × 1 ° "> 1°×1° was used in the travel‐time calculations. Relocation results show a reduction in travel‐time residuals compared with the initial locations. The relocation results better illuminate subducted slabs and active faults in the region such as the Mentawai back thrust and the outer rise in the subduction zone south of Java. Focal mechanisms from the Global Centroid Moment Tensor catalog are analyzed in conjunction with the relocation results, and our synthesis of the results provides further insight into seismogenesis in the region.

Sunda arc↗