# Relocation of hypocenters from DOMERAPI and BMKG networks: a preliminary result from DOMERAPI project

Mohamad Ramdhan, Sri Widiyantoro, Andri Dian Nugraha, Jean-Philippe Métaxian, Asep Saepuloh, Said Kristyawan, Andry Syaly Sembiring, Agus Budi Santoso, Antoine Laurin, Ahmad Ali Fahmi

<figure><img src="image_1.png" /></figure>

► To cite this version:

Mohamad Ramdhan, Sri Widiyantoro, Andri Dian Nugraha, Jean-Philippe Métaxian, Asep Saepuloh, et al.. Relocation of hypocenters from DOMERAPI and BMKG networks: a preliminary result from DOMERAPI project. Earthquake Science, 2017, 30 (2), pp.67-79. 10.1007/s11589-017-0178-3. hal-04085169

HAL Id: hal-04085169

https://hal.science/hal-04085169

Submitted on 29 Apr 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Relocation of hypocenters from DOMERAPI and BMKG networks: a preliminary result from DOMERAPI project

Mohamad Ramdhan · Sri Widiyantoro · Andri Dian Nugraha · Jean-Philippe Métaxian · Asep Saepuloh · Said Kristyawan · Andry Syaly Sembiring · Agus Budi Santoso · Antoine Laurin · Ahmad Ali Fahmi

Received: 21 November 2016/Accepted: 7 March 2017/Published online: 4 May 2017
© The Author(s) 2017. This article is an open access publication

**Abstract** Merapi volcano located in central Java, Indonesia, is one of the most active stratovolcanoes in the world. Many Earth scientists have conducted studies on this volcano using various methods. The geological features around Merapi are very attractive to be investigated because they have been formed by a complex tectonic process and volcanic activities since tens of millions of years ago. The southern mountain range, Kendeng basin and Opak active fault located around the study area resulted from these processes. DOMERAPI project was conducted to understand deep magma sources of the Merapi volcano comprehensively. The DOMERAPI network was running from October 2013 to mid-April 2015 by deploying 46 broad-band seismometers around the volcano. Several steps, i.e., earthquake event identification, arrival time picking of P and S waves, hypocenter determination and hypocenter relocation, were carried out in this study. We used Geiger's method (Geiger 1912) for hypocenter determination and double-difference method for

hypocenter relocation. The relocation result will be used to carry out seismic tomographic imaging of structures beneath the Merapi volcano and its surroundings. For the hypocenter determination, the DOMERAPI data were processed simultaneously with those from the Agency for Meteorology, Climatology and Geophysics (BMKG) seismic network in order to minimize the azimuthal gap. We found that the majority of earthquakes occurred outside the DOMERAPI network. There are 464 and 399 earthquakes obtained before and after hypocenter relocation, respectively. The hypocenter relocation result successfully detects some tectonic features, such as a nearly vertical cluster of events indicating a subduction-related backthrust to the south of central Java and a cluster of events to the east of Opak fault suggesting that the fault has an eastward dip.

**Keywords** Earthquake · Relocation · DOMERAPI · BMKG · Tectonic feature

<footnote>M. Ramdhan · S. Kristyawan · A. S. Sembiring
Study Program of Earth Sciences, Faculty of Earth Science and Technology, Institute of Technology Bandung, Jalan Ganesa No. 10, Bandung 40132, Indonesia
e-mail: mramdhan123@gmail.com</footnote>

<footnote>M. Ramdhan · S. Kristyawan · A. S. Sembiring
Agency for Meteorology, Climatology and Geophysics (BMKG), Jakarta, Indonesia</footnote>

<footnote>S. Widiyantoro · A. D. Nugraha (✉)
Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, Institute of Technology Bandung, Jalan Ganesa No. 10, Bandung 40132, Indonesia
e-mail: nugraha@gf.itb.ac.id</footnote>

<footnote>S. Widiyantoro
Research Center for Disaster Mitigation, Institute of Technology Bandung, Jalan Ganesa No. 10, Bandung 40132, Indonesia</footnote>

<footnote>J.-P. Métaxian · A. Laurin · A. A. Fahmi
ISTerre, IRD R219, CNRS, Université de Savoie Mont Blanc, Le Bourget-du-Lac, France</footnote>

<footnote>A. Saepuloh
Faculty of Earth Sciences and Technology, Institute of Technology Bandung, Jalan Ganesa No. 10, Bandung 40132, Indonesia</footnote>

<footnote>A. B. Santoso
Center for Volcanology and Geological Hazard Mitigation, Geological Agency, Jalan Diponegoro No. 57, Bandung 40122, Indonesia</footnote>
# 1 Introduction

Merapi volcano is located in central Java and surrounded by densely populated areas so that its eruption effect will be hazardous for the people living around the volcano. Merapi has a relatively very short return period of eruptions, i.e., from 4 to 6 years with volcanic explosivity index (VEI) of 1–3 (Andreastuti et al. 2000; Voight et al. 2000). The eruption type of Merapi is dominated by pyroclastic flow caused by dome collapse. High VEIs (more than 4) with an explosive type occur in a longer return period, for instance the 2010 eruption event. Such an explosive type of eruption occurs every ±100 years (Surono et al. 2012). This eruption type has been controlled by deep magma sources with different mechanism (Costa et al. 2013). Based on the result of a petrological study (Costa et al. 2013), there are two magma bodies beneath the Merapi volcano, i.e., at ~30 and ~13 km depths. The DOMERAPI project was conducted partly to understand the internal structure and characteristics of the Merapi volcano.

Many earth scientists have conducted research on Merapi using various methods to understand the process related to physical, chemical and geological parameters beneath the volcano and its summit. Seismological approach can also be used to understand physical parameters beneath a volcano. For example, studies using MERAMEX (Merapi Amphibious Experiment) data have resulted in some publications explaining the relation between subduction zone and volcanic arc in central Java (Wagner et al. 2007; Koulakov et al. 2007; Rohadi et al. 2013; Haberland et al. 2014). These studies depict clearly a low velocity ascending from the subducted slab to the volcanic arc, but they do not show a detailed velocity structure of Merapi due to the lack of data coverage beneath the volcano.

Detailed studies of Merapi have been performed using volcano-ectonic (VT) earthquake data to analyze seismicity and focal mechanism of events beneath the volcano (Ratdomopurbo and Poupinet 2000; Hidayati et al. 2008; Budi-Santoso et al. 2013). They show that the maximum depth of VT earthquakes is about 5 km depth from the summit, so structures below this depth cannot be delineated.

In this study, the data from the BMKG network in the same period of the DOMERAPI project were incorporated in the hypocenter determination and relocation to maximize azimuthal coverage. Epicenters from the DOMERAPI network have an azimuthal gap average of more than 300°. The joint networks have successfully minimized the azimuthal gap to <250° (Ramdhan et al. 2015).

# 2 Data and method

The DOMERAPI network covers Merapi and Merbabu volcanoes and geological features, such as the southern

mountain range of central Java, part of the Kendeng basin and Opak fault (Fig. 1). The network consists of 53 sites of broad-band seismometers operated from October 2013 to mid-April 2015. The network consists of 46 seismometers at the beginning, and then seven seismometers are relocated to different sites to obtain better S/N ratio. The seismometer spacing average is about 4 km. This network design is potential to enhance seismic image resolution in the study area in unprecedented detail. A preliminary work on seismic ambient noise tomography was applied using the network to delineate subsurface structure beneath Merapi and Merbabu volcanoes (Trichandi et al. 2015).

Broad-band seismometers detect wide frequency signals, so they detect not only local events but also regional and teleseismic events as well. For this study, we selected regional events in the radius of 500 km from the Merapi volcano. Figure 2a, b shows examples of 3-component waveforms from local and regional events recorded by the DOMERAPI network. We identified events and picked P- and S-wave arrival times carefully. Events were determined by at least six DOMERAPI stations for P and S phases or P phase only. The P phase has much better onset quality than S leading to a less number of S phase than P.

<figure><img src="image_4.png" /><figcaption>**Fig. 1** Distribution of seismographic stations of the DOMERAPI network. The *inset picture* shows the boundary of DOMERAPI network area. The geological map was modified from Smyth et al. (2008). The 2006 earthquake location was taken from global CMT catalogue (Ekström et al. 2016)</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 2** a Examples of 3-component waveforms from local events recorded by the DOMERAPI network. b The same as a, but from regional event. c Shows the epicentral distribution of earthquakes and stations recording these two events. The left and right figures are close up the study region to show the station distribution recording the local and regional events, respectively</figcaption></figure>
b

<figure><img src="image_6.png" /><figcaption>Fig. 2 continued</figcaption></figure>

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /><figcaption>Fig. 2 continued</figcaption></figure>

The earthquake hypocenters were first determined using the Geiger’s method (Geiger 1912) implemented with the hypoellipse software (Lahr 1999). The method was also applied to determine earthquakes along Cimandiri and Baribis faults in West Java based on BMKG data (Supendi and Nugraha 2016). We relocated the hypocenters using a double-difference earthquake location algorithm (Waldhauser and Ellsworth 2000; Waldhauser 2001). This method was successfully applied to relocate hypocenters in the Indonesia region based on data from BMKG network (Ramdhan and Nugraha 2013; Cahyaningrum et al. 2015; Sabtaji and Nugraha 2015; Utama et al. 2015). The hypocenter relocation was conducted carefully to obtain precise hypocenters before conducting seismic travel time tomography for our future work.

## 3 Results and discussion

In this study, we incorporated events recorded by the BMKG seismometers outside the DOMERAPI network as shown in Fig. 3 to maximize azimuthal coverage. There are 33 broad-band seismometers selected from the BMKG network. XMIS is one of the Australian seismographic stations that was used in the processing in order to better constrain event locations, especially to the south of Java. The hypocenter determination and relocation hypocenter used a 1-D P velocity model for central Java (Koulakov et al. 2007) and $v_P/v_S$ ratio of 1.73 obtained from the Wadati diagram constructed using the combined DOMERAPI and BMKG data sets (see Wadati 1933).

There are 464 earthquakes determined from the data recorded by the DOMERAPI and BMKG networks. After hypocenter relocation, we end up with 399 events. There are some earthquakes that could not be relocated which may be due to poorly clustered events. A comparison between epicenter distributions from the DOMERAPI and BMKG networks before and after relocation is shown in Fig. 3. Most of the resulting hypocenters are located outside the DOMERAPI network, so they can be useful for imaging the deep structure of Merapi in our future work on seismic tomography. Our new data set indicates that there are still many events to the east of the Opak fault events after nearly a decade of the Yogyakarta earthquake in 2006. In addition, a thrust fault to the south of central Java, which may be related to a backthrust, is clearly delineated by the relocated events.

Epicenter shifts after relocation and a rose diagram depicting the number of relocations for 5° azimuth bins are given in Fig. 4. It is shown clearly that most of the shifts are perpendicular to the trench. The average epicenter shift is ~5 km, and the maximum shift is ~27 km. The majority of the epicenter shifts are northward which may be caused by the seismographic station distribution, in which almost all of the stations are located in Java Island. The shifting perpendicular to the trench has been interpreted as most likely due to an effect of a 3-D velocity model incorporating the slab, because teleseismic rays spend a lot of time in the slab (see, e.g., Pesicek et al. 2010). In this study, however, we used a lot of local data and 1-D velocity model, but we still obtained a similar result. So the shifting perpendicular to the trench seems to be caused by the station distribution, which is generally
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 3** a Epicentral distribution from the DOMERAPI and BMKG networks before hypocenter relocation. b The same as a, but after relocation. Symbol *yellow triangles* depict the location of seismographic stations of the BMKG network used to determine earthquake hypocenters along with the DOMERAPI stations</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Fig. 4** a Epicenter shifts owing to event relocations. Rectangles depict the location of vertical cross sections A–A' and B–B' shown in Figs. 6 and 8, respectively. b Rose diagrams showing number of relocations for 5° azimuth bins. The shifting values are normalized first to be plotted in the diagram. The maximum value in the rose diagram is 1.0, i.e., a normalized value. The increment scale of 0.2 in the radial direction is shown</figcaption></figure>

located on land. The horizontal errors (X and Y error) and
vertical errors (Z error) decrease significantly after
hypocenter relocation. The majority of errors for all
direction are <0.3 km. Figure 5 shows the error

distribution before and after hypocenter relocation. The
result of the hypocenter relocation is used to analyze the
seismotectonics in the study region as depicted by the
following figures.
<figure><img src="image_12.png" /><figcaption>**Fig. 5** a Plot of errors in X (east-west) direction before hypocenter relocation (red square) and after hypocenter relocation (blue square). b The same as a, but in Y (north-south) direction. c The same as a, but in Z direction</figcaption></figure>
There are two tectonic features detected clearly by the relocated hypocenters. They are: (1) a nearly vertical cluster of events that may be related to a backthrust to the south of the study region and (2) the active Opak fault to the south of Merapi, which is indicated by small events even after a decade of the Yogyakarta earthquake of May 26, 2006, at 22:54 UTC (Ekström et al. 2016).

A vertical section (A-A') crossing the Opak fault is shown in Fig. 6. The cluster of events beneath the top of slab, from a depth of ~40 km down to ~115 km, may be related to a double seismic zone as discussed by Koulakov et al. (2007), although we could not delineate it clearly. Notice the existence of fixed depths before hypocenter relocation, which was successfully overcome after the hypocenter relocation.

Figure 7 depicts the seismic activity occurring to the east of the Opak fault during the DOMERAPI project period. This fault has been interpreted to be related to the Yogyakarta earthquake occurring on May 26, 2006. The earthquake mechanism is sinistral strike slip with $M_W$ 6.3 (Nakano et al. 2006; Ekström et al. 2016). The earthquake occurred to the east of the Opak fault line. The earthquake source may come from a hidden fault because no geomorphic feature was well noticed before the earthquake (Tsuji et al. 2009). Geological condition shows that surrounding area beneath the hidden fault is dominated by sediment covers. Usually an earthquake related to an active fault occurs at the basement part of continental crust. For the case of a thick sediment cover, the slip that is created by an earthquake is not much stronger to develop surface rupture and make a geomorphological feature like surface faults. It is consistent with the geological condition, in which the hypocenter of the earthquake was located beneath an alluvial deposit (Smyth et al. 2008). Based on the relocation result, the depth of the majority of events to the east of the Opak fault is <25 km. This is in agreement with the crustal thickness in the area, i.e., around 35 km (Wölbern and Rümpker 2016). From the hypocenter depths, it is clear that these events indeed occurred in the basement part of continental crust. At this stage, we speculate that the cluster of events parallel to Opak fault is alternatively related to the fault. Our preferred interpretation is that the Opak fault is dipping eastward. However, a further investigation is still needed.

Figure 8 shows a vertical cross section B-B' depicting hypocenters before and after hypocenter relocation. It is

<figure><img src="image_13.png" /><figcaption>Fig. 6 a Vertical cross sections A-A' showing the distribution of hypocenters a before and b after relocation. Blue line depicts the slab 1.0 model from USGS (Hayes et al. 2012). The location of cross sections is shown in Fig. 4</figcaption></figure>

shown that the dip of the slab is quite steep, similar with the one, e.g., beneath Andaman (Pesicek et al. 2008) associated with the subduction of an old oceanic lithosphere. A striking feature observed after hypocenter relocation is the nearly vertical of events down to a depth of ~30 km. We interpret that this cluster of events represents a backthrust existing to the south of central Java. A backthrust, in which its displacement is generally in an opposite direction to that of the main thrust direction, is
<figure><img src="image_14.png" /><figcaption>**Fig. 7** Seismicity to the east of Opak fault based on the hypocenter relocation result using jointly the data from the DOMERAPI and BMKG networks. The 2006 earthquake location was taken from global CMT catalogue (Ekström et al. 2016)</figcaption></figure>
assumed to form as a result of layer-parallel shortening in a
late stage of thrust sequences (Xu et al. 2015). In the case
of Java, it is similar to that existing along Sumatra, which is

<figure><img src="image_15.png" /><figcaption>Fig. 8 The same as Fig. 6, but for vertical cross section B-B'</figcaption></figure>

more obvious (Nugraha et al. 2015). This observation
agrees well with the result of teletomoDD hypocenter
relocation using the BMKG catalogue from 2009 to 2014
(Nugraha et al. 2015). Our interpretation is also supported
by the focal mechanism of some selected events (Ekström
et al. 2016) shown in Fig. 9. In general, they have a
mechanism of a thrust fault that is related to the main thrust
and a backthrust.

4 Conclusions and further work

The hypocenter relocation result reveals some important
tectonic features. One of them is the observed nearly ver-
tical cluster of events to the south of central Java, which is
related to a backthrust fault. In addition, the relocated
events also depict a cluster of earthquakes parallel to the
Opak fault characterized by a contrast in wavespeed
(Zulfakriza et al. 2014), and a possible double seismic zone
in the study region.

For further work, we will use our hypocenter relocation
result to perform seismic travel time tomography in order
to construct 3-D vP, vS and vP/vS models. These models will
be used to understand the deep structure of the Merapi
volcano, which has not been detected clearly by previous
works due to lack of data coverage. As previous volcanic
studies (e.g., Lei et al. 2013), a multiscale seismic tomo-
graphic imaging will be carried out utilizing the
DOMERAPI and BMKG data and by incorporating the
data from MERAMEX as used in previous works (Wagner
et al. 2007; Koulakov et al. 2007; Rohadi et al. 2013;
Haberland et al. 2014) to reconstruct magma bodies of
Merapi and ascending low velocity features from the sub-
ducting slab to the volcanic arc and some other tectonic
features. We will also combine data from the Center for
Investigation and Technology Development of Geological
Disaster, Indonesia (BPPTKG), to obtain images,
<figure><img src="image_16.png" /><figcaption>**Fig. 9** Focal mechanism solutions obtained by events selected from the global CMT catalogue (Ekström et al. 2016) showing that most events have thrust mechanisms</figcaption></figure>

especially below the Merapi summit with unprecedented
detailed resolution.

**Acknowledgements** We are grateful to Institut de Recherche pour le Développement (IRD), France, for funding the DOMERAPI project; Center for Volcanology and Geohazard Mitigation as the main counterpart of the DOMERAPI project in Indonesia; the Agency for Meteorology, Climatology and Geophysics (BMKG) for providing the waveform data used in this study; and the Indonesia Endowment Fund for Education (LPDP) for granting a doctoral scholarship to MR. This study has been supported in part by the Indonesian Directorate General of Higher Education (DIKTI) research funding 2015–2016 and the Institut Teknologi Bandung (ITB) through a WCU research Grant 2016 awarded to SW.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

References

Andreastuti SD, Alloway BV, Smith IEM (2000) A detailed tephrostratigraphic framework at Merapi Volcano, Central Java, Indonesia: implications for eruption predictions and hazard assessment. J Volcanol Geotherm Res 100(1-4):51-67

Budi-Santoso A, Lesage P, Dwiyono S, Sumarti S, Subandriyo, Surono, Metaxian J-P (2013) Analysis of the seismic activity associated with the 2010 eruption of Merapi Volcano, Java. J Volcanol Geotherm Res 261:153-170

Cahyaningrum AP, Nugraha AD, Puspito NT (2015) Earthquake hypocenter relocation using double difference method in East Java and surrounding areas. In: AIP conference proceedings, New York, p 30021. doi:10.1063/1.4915029

Costa F, Andreastuti S, Bouvet de Maisonneuve C, Pallister JS (2013) Petrological insights into the storage conditions, and magmatic processes that yielded the centennial 2010 Merapi explosive eruption. J Volcanol Geotherm Res 261:209-235

Ekström G, Nettles M, Dziewonski AM (2016) Global Centroid Moment Tensor Project. http://www.globalcmt.org/. Accessed 7 Nov 2016

Geiger L (1912) Probability method for the determination of earthquake epicenters from the arrival time only. Bull St Louis Univ 8:56-71
Haberland C, Bohm M, Asch G (2014) Accretionary nature of the crust of Central and East Java (Indonesia) revealed by local earthquake travel-time tomography. *J Asian Earth Sci* 96:287–295

Hayes GP, Wald DJ, Johnson RL (2012) Slab1.0: a three-dimensional model of global subduction zone geometries: slab1.0 3D subduction geometry. *J Geophys Res Solid Earth*. doi:10.1029/2011JB008524

Hidayati S, Ishihara K, Iguchi M, Ratdomopurbo A (2008) Focal mechanism of volcano-tectonic earthquakes at Merapi volcano, Indonesia. *Indones J Phys* 19(3):75–82

Koulakov I, Bohm M, Asch G, Luhr B-G, Manzanares A, Brotopuspito KS, Fauzi, Purbawinata MA, Puspito NT, Ratdomopurbo A, Kopp H, Rabbel W, Shevkunova E (2007) P and S velocity structure of the crust and the upper mantle beneath central Java from local tomography inversion. *J Geophys Res* 112(B8):B08310. doi:10.1029/2006JB004712

Lahr JC (1999) Revised 2012, HYPOELLIPSE: a computer program for determining local earthquake hypocentral parameters, magnitude, and first-motion pattern. *US Geol Surv Openfile report* 99

Lei J, Xie F, Fan Q, Santosh M (2013) Seismic imaging of the deep structure under the Chinese volcanoes: an overview. *Phys Earth Planet Int* 224:104–123. doi:10.1016/j.pepi.2013.08.008

Nakano M, Kumagai H, Miyakawa K, Yamashita T, Inoue H, Ishida M, Aoi S, Morikawa N, Harjadi PJ (2006) Source estimates of the May 2006 Java earthquake. *EOS Trans Am Geophys Union* 87(45):493–494

Nugraha AD, Shiddiqi HA, Widiyantoro S, Ramdhan M, Wandono, Sutiyono, Handayani T, Nugroho H (2015) Preliminary results of teleseismic double-difference relocation of earthquakes around Indonesian archipelago region. In: 4th International symposium on earthquake and disaster mitigation 2014 (ISEDM 2014). AIP Publishing, New York, p 30002. doi:10.1063/1.4915010

Pesicek JD, Thurber CH, Widiyantoro S, Engdahl ER, DeShon HR (2008) Complex slab subduction beneath northern Sumatra. *Geophys Res Lett* 35(20):L20303. doi:10.1029/2008GL035262

Pesicek JD, Thurber CH, Zhang H, DeShon HR, Engdahl ER, Widiyantoro S (2010) Teleseismic double-difference relocation of earthquakes along the Sumatra–Andaman subduction zone using a 3-D model. *J Geophys Res* 115:B10303. doi:10.1029/2010JB007443

Ramdhan M, Nugraha AD (2013) Study of seismicity around Toba area based on relocation hypocenter result from BMKG catalogue. In: Padjadjaran international physics symposium 2013 (PIPS2013): contribution of physics on environmental and energy conservations. AIP Publishing, New York, pp 242–244. doi:10.1063/1.4820330

Ramdhan M, Nugraha AD, Widiyantoro S, Mètaxian J-P, Valencia AA (2015) Earthquake location determination using data from DOMERAPI and BMKG seismic networks: a preliminary result of DOMERAPI project. In: 4th International symposium on earthquake and disaster mitigation 2014 (ISEDM 2014). AIP Publishing, New York, p 30007. doi:10.1063/1.4915015

Ratdomopurbo A, Poupinet G (2000) An overview of the seismicity of Merapi volcano (Java, Indonesia), 1983–1994. *J Volcanol Geotherm Res* 100(1–4):193–214

Rohadi S, Widiyantoro S, Nugraha AD, Masturyono (2013) Tomographic imaging of P- and S-wave velocity structure beneath central Java, Indonesia: joint inversion of the MERAMEX and MCGA earthquake data. *Int J Tomogr Stat* 24(3):1–16

Sabtaji A, Nugraha AD (2015) 1-D seismic velocity model and hypocenter relocation using double difference method around West Papua region. In: AIP conference proceedings, New York. doi:10.1063/1.4915013

Smyth HR, Hall R, Nichols GJ (2008) Cenozoic volcanic arc history of East Java, Indonesia: the stratigraphic record of eruptions on an active continental margin. In: Special paper 436: formation and applications of the sedimentary record in arc collision zones. *Geol Soc Am* 436:199–222

Supendi P, Nugraha AD (2016) Preliminary result of earthquake hypocenter determination using hypoellipte around western Java region. In: The 5th international symposium on earthhazard and disaster mitigation: the annual symposium on earthquake and related geohazard research for disaster risk reduction. AIP Publishing, New York. doi:10.1063/1.4947370

Surono JP, Pallister J, Boichu M, Buongiorno MF, Budisantoso A, Costa F, Andreastuti A, Prata F, Schneider D, Clarisse L, Humaida H, Sumarti S, Bignami C, Griswold J, Carn S, Oppenheimer C, Lavigne F (2012) The 2010 explosive eruption of Java's Merapi volcano-a "100-year" event. *J Volcanol Geotherm Res* 241–242:121–135

Trichandi R, Yudistira T, Nugraha AD, Zulfakriza Z, Saygin E (2015) The preliminary results: seismic ambient noise Rayleigh wave tomography around Merapi volcano, central Java, Indonesia. In: AIP conference proceedings, New York. doi:10.1063/1.4915049

Tsuji T, Yamamoto K, Matsuoka T, Yamada Y, Onishi K, Bahar A, Meilano I, Abidin HZ (2009) Earthquake fault of the 26 May 2006 Yogyakarta earthquake observed by SAR interferometry. *Earth Planets Space* 61(7):e29. doi:10.1186/BF03353189

Utama MRJ, Nugraha AD, Puspito NT (2015) Seismicity studies at Moluccas area based on the result of hypocenter relocation using HypoDD. In: AIP conference proceedings, New York. doi:10.1063/1.4915030

Voight B, Constantine EK, Siswowidjoyo S, Torley R (2000) Historical eruptions of Merapi Volcano, Central Java, Indonesia, 1768–1998. *J Volcanol Geotherm Res* 100(1–4):69–138

Wadati K (1933) On the travel time of earthquake waves II. *Geophys Mag* 7:101–111

Wagner D, Koulakov I, Rabbel W, Lüehr B-G, Wittwer A, Kopp H, Bohm M, Asch G, MERAMEX Scientists (2007) Joint inversion of active and passive seismic data in Central Java. *Geophys J Int* 170(2):923–932

Waldhauser F (2001) HypoDD: a computer program to compute double- difference earthquake location. *US Geol Surv Openfile report*, pp 1–113

Waldhauser F, Ellsworth WL (2000) A double-difference earthquake location algorithm: method and application to the northern Hayward fault, California. *Bull Seismol Soc Am* 90(6):1353–1368

Wölbern I, Rümpker G (2016) Crustal thickness beneath Central and East Java (Indonesia) inferred from P receiver functions. *J Asian Earth Sci* 115:69–79

Xu S, Fukuyama E, Ben-Zion Y, Ampuero J-P (2015) Dynamic rupture activation of backthrust fault branching. *Tectonophysics* 644–645:161–183. doi:10.1016/j.tecto.2015.01.011

Zulfakriza Z, Saygin E, Cummins PR, Widiyantoro S, Nugraha AD, Luhr B-G, Bodin T (2014) Upper crustal structure of central Java, Indonesia, from transdimensional seismic ambient noise tomography. *Geophys J Int* 197(1):630–635. doi:10.1093/gji/ggu016