Double-difference tomography of P- and S-wave velocity structure beneath the western part of Java, Indonesia

West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults. In this study, double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java. To infer the geometry of the structure beneath the study area, precise earthquake hypocenter determination was first performed before tomographic imaging. For this, earthquake waveform data were extracted from the regional Meteorological, Climatological, Geophysical Agency (BMKG) network of Indonesia from South Sumatra to Central Java. The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected, the key features being events of magnitude > 3, azimuthal gap < 210° and number of phases > 8. A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters. The hypocenter locations were then relocated using double-difference tomography (tomoDD). A significant reduction of travel-time (root mean square basis) and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java. Double-difference tomography was found to give a clear velocity structure, especially beneath the volcanic arc area, i.e., under Mt Anak Krakatau, Mt Salak and the mountains complex in the southern part of West Java. Low velocity anomalies for the P and S waves as well as the v_P/v_S ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.