Significant earthquakes on the island of Sumatra, Indonesia, have predominantly been earthquakes with a thrust mechanism that occurred due to the subduction process and seismotectonics near coastal cities of West and South Sumatra, which could be affected by earthquakes triggered by these seismic sources. We compared the Seismic Hazard Function (SHF) of two coastal cities of Sumatra: Bengkulu and Padang. The results showed that the SHF of Bengkulu is higher than that of Padang. Estimated earthquake hazards are presented in the form of seismic hazard maps expressed as the PGA of 10% rate of exceedance probability in 50 years. In estimating the seismic potential in Sumatra, the seismic moment rate was jointly estimated from the smoothed mean seismicity rate and the pre-seismic subduction surface strain rate model. In this study, the island of Sumatra was chosen as a master model for Seismic Hazard Analysis (SHA). The motivation for choosing Sumatra for the SHA was because of the large body of complete historical earthquake data of the North Western Sunda Arc. The SHF is calculated based on a magnitude range of 6.0 to 9.0 during 50 years with the radius distance from the source less than or equal to 100 km.

No abstract

The island of Lombok in Indonesia is located between the Indo-Australian and Eurasian subduction trenches and the Flores back-arc thrust, making it vulnerable to earthquakes. On 29 July 2018, a significant earthquake Mw 6.4 shook this region and was followed by series of major earthquakes (Mw>5.8) on 5, 9, and 19 August, which led to severe damage in the northern Lombok area. In this study, we attempt to reveal the possible cause of the sequences of the 2018 Lombok earthquakes based on aftershock monitoring data. Twenty stations were deployed to record earthquake waveform data from 4 August to 9 September 2018. In total, 3259 events were identified using 28,728 P- and 20,713 S-wave arrival times during the monitoring. The aftershock hypocenters were determined using a nonlinear approach and relocated using double-difference method. The moment magnitude (Mw) of each event was determined by fitting the displacement spectrum amplitude using a Brune-type model. The magnitudes of the aftershocks range from Mw 1.7 to 6.7. The seismicity pattern reveals three clusters located in the Flores oceanic crust, which fit well with the occurrences of the four events with Mw>6. We interpret these events as the main rupture area of the 2018 Lombok earthquake sequence. Furthermore, an aseismic zone in the vicinity of Rinjani extending toward the northwestern part of Lombok was observed. We propose that the crust in this area has elevated temperatures and is highly fractured thus inhibiting the generation of large earthquakes. The aseismic nature is therefore an artifact of the detection threshold of our network (Mw 4.6).

Magma is fed to a volcano through a complex "plumbing" system that involves not only shallow structures beneath the volcano edifice, but also deep structures and processes within the underlying crust and upper mantle. This paper summarizes seismic experiments carried out over many years at Gunung Merapi in Central Java. These have resolved the 3D seismic velocity structure of the Merapi edifice, and provided a 3D structural image of the lithosphere and subduction zone beneath Central Java. Earthquake locations reveal that with distance from the trench, the dip of the subducting slab steepens from nearly horizontal (0-150 km), through 45 degrees (150-250 km), to 70 degrees (>250 km). The slab appears as a 30 km thick double layer of seismicity in a depth range of 80 km to 150 km, and it can be identified seismically to a depth of more than 600 km. The active volcanoes of Merapi, Sumbing, and Lawu are located at the edge of a large low velocity body that extends from the upper crust to the upper mantle beneath Central Java. Shear wave signals recorded above this anomaly are strongly attenuated compared to neighboring areas. The anomalous body 2 has a detected volume of >50,000 km 3 and a decrease in P and S velocities relative to adjacent regions of up to 30%. The resulting Vp/Vs ratio of up to 1.9 is unusually high for lower crust. Additionally, the anomaly extends along a 45 degree-slope downward from beneath the volcanic arc and meets the slab at 100 km depth. We interpret this sloping anomaly as a pathway for fluids and partial melts. Increased seismicity is observed at depths of ~100 km, possibly as a result of dehydration of the subducting slab with related fluid releases causing partial melting of overlying mantle material. The large velocity reduction and high Vp/Vs ratio in the region are consistent with an increase in temperature, a reduction of shear strength, and the presence of fluids or melts of 13 to 25 vol. %. The detected strong anomaly beneath Central Java is unique in size and amplitude compared to other subduction zones. The geophysical evidence suggests that this segment of the arc has a high magma flux and is thus capable of developing even larger shallow crustal reservoirs and more voluminous explosive eruptions in the future.

The Agung-Batur Volcanic Complex (ABVC), part of the Sunda volcanic arc, is the source of some of the most hazardous volcanic activity in Indonesia. The ABVC has undergone many small (VEI 1-2) eruptions since historical records began in the early 19th century, but Mt. Agung has experienced much larger (VEI 5) eruptions, both in the modern (1963) and historical (1843) eras, as well as several times during the past 2000-3000 years. The 1963 eruption caused more than 1000 deaths, and a more recent eruption in 2017 caused the evacuation of 140,000 people. Delineating the magma structure beneath ABVC is an important first step in understanding the physics of these eruptions. This paper presents the first local-scale study of Rayleigh wave group velocity structure and the seismic velocity structure beneath the ABVC using ambient seismic noise tomography. Seismic data were collected using 25 seismometers deployed across the ABVC during early January to March 2019. The local seismic network provides good resolution beneath both Mt. Agung and Mt. Batur. We obtained 158 Rayleigh Green's Functions, extracted by cross correlating noise simultaneously recorded at available station pairs. We used sub-space inversion to calculate group velocity at different periods and to estimate the lateral variations in group velocity for given periods. 2-D tomographic maps obtained from the inversion of the group velocity of Rayleigh waves clearly showed some pronounced velocity anomalies beneath the ABVC. We applied the Neighbourhood Algorithm (NA) technique to invert the Rayleigh wave dispersion curves to obtain shear wave velocity (Vs) vs. depth profiles. These profiles indicate a low Vs of about 1 km/s underlying the volcanic complex between Mt. Agung and Mt. Batur at depths up to 2 km, which we suggest is due to a combination of low-Vs volcanic deposits as well as a shallow hydrothermal fluids system associated with magma fluids and/or gases produced by magma intrusion at depths >7 km.