We reveal the existence of a previously unknown fault that generated the Mw 7.3 Flores Sea earthquake, which occurred on 14 December 2021, approximately 100 km to the north of Flores Island, in one of the most complex tectonic settings in Indonesia. We use a double-difference method to relocate the hypocenters of the mainshock and aftershocks, determine focal mechanisms using waveform inversion, and then analyze stress changes to estimate the fault type and stress transfer. Our relocated hypocenters show that this earthquake sequence ruptured on at least three segments: the source mechanism of the mainshock exhibits dextral strike-slip motion (strike N72°W and dip 78° NE) on a west–east-trending fault that we call the Kalaotoa fault, whereas rupture of the other two segments located to the west and east of the mainshock (striking west-northwest and southeast, respectively) may have been triggered by this earthquake. The Coulomb stress change imparted by the rupture of these segments on nearby faults is investigated, with a focus on regions that experience a stress increase with few associated aftershocks. Of particular interest are stress increases on the central back-arc thrust just north of Flores and the north–south-striking Selayar fault in the northwest of our study region, both of which may be at increased risk of failure as a result of this unusual earthquake sequence.

The Central part of Sumatra is a region that has a high potential for earthquakes. This research intended to determine the crustal thickness of the earth, P and S wave velocity profiles, and vp/vs value in the Central part of Sumatra using stacking H-k and inversion techniques based on the analysis of receiver function. This study utilized teleseismic earthquake data with a distance of 30° to 90° from the station and magnitude more than 6 (M>6). The stations used in this study were 3 BMKG broadband stations located in 3 zones, the fore arc ridge zone (SISI), the volcanic zone (PLSI) and the back arc zone (TPRI). The crustal thickness varies in the fore arc ridge zone (SISI) estimated 17.8 km, volcanic zone (PLSI) reaches 29.7 km and the back arc zone (TPRI) reaches 34 km. The crustal thickness is quite thick under the PLSI and thicker beneath TPRI station. These possibly due to the effect of topography and isostatic compensation in the station. However, whether there is a correlation between crustal thickness and topography needs further research using more stations. The highest vp/vs value was found in the volcanic zone of 1.9, that might be associated with the presence partial melting beneath the station. Meanwhile, the vp/vs value in the back arc zone is 1.72, indicating a relatively more homogeneous structure.

Lampung region is seismically and volcanic active because located in subduction zone of Indo-Australian and Eurasian plate. We applied receiver function and stacking H-k analysis to estimate the crustal structure in Lampung region. We used teleseismic earthquake data (epicenter distance 30°-90°) and M>6 recorded at 3 seismic broadband stations owned by Agency for Meteorology Climatology and Geophysics (BMKG). Those stations are PSLI (located on Sebesi Island approximately 20 km from Anak Krakatau) represented volcanic arc zone, KASI (located on Kota Agung, Lampung) represented Sumatran Fault Zone and KLI (located on Kotabumi, Lampung) represented back-arc basin. Crustal thickness estimated at PSLI station 32-36 km, KASI station 36-40 km, and KLI station 30-36 km. Furthermore, in 3 stations P wave velocity estimated 4.1-11 km/s, S wave velocity 2.2-6.2 km/s, while vp/vs value estimated 1.7-2.05. We estimated Anak Krakatau volcano’s magma chamber beneath PSLI station in depth 16-30 km, Great Sumatran Fault structure in depth about 8-14 km beneath KASI station, and thick sediment layer about 4 km near surface beneath KLI station. This study result is expected to explain more detail crustal of Lampung region and can be useful for developing of BMKG’s seismic monitoring systems and other geophysical fields in future.

Tsunami modeling with the Manila Trench earthquake source, the Philippines is very necessary as one of the mitigation measures for the prospective nuclear power plant (PLTN) facility in West Kalimantan. In this preliminary study, tsunami modeling using the COMCOT of non linear and the nesting grid method was carried out to obtain an estimate of the height and arrival time of the tsunami at the closest tidal initial points of a nuclear power plant. The Manila Trench is considered as the closed potential sources of earthquake threats to the prospective nuclear power plant site in Gosong Beach, West Kalimantan. The source is divided into 33 multiple segments and accumulated to magnitude 9.1 The results of this study can be used as material for risk assessment for nuclear power plant site development. The results obtained are that the estimated tsunami height at Gosong Beach is around 0.3 – 0.5 meters on the beach with an arrival time of about 550 minutes after the earthquake at its center.

Aceh region has a very complex crustal structure from the forearc ridge to the backarc basin. This study aims to determine the velocity model of P and S waves and the depth of Moho discontinuity. This research was conducted using teleseismic earthquake data (30°-90° from the station) with M>6 from four seismic stations belonging to the BMKG in Aceh region. The stations are qualified based on the volcanic arc system zone. Furthermore, the velocity model determined by result of forward modelling, while the depth of the Moho layer estimated by migrated receiver function from time domain to the depth domain. At station SNSI that represented the forearc ridge zone, the depth of Moho is ±28 km, at station TPTI represent the forearc basin is ±16 km, while at zone with higher topography, namely volcanic arc zone represented by station KCSI, the Moho depth was identified at ±38 km, and the backarc basin represented by station LASI with ±40 km depth of Moho. This variation occurs because the composition of the earth’s layers below the station is diverse also different topography for each station.