Imaging of a magma system beneath the Merapi Volcano complex, Indonesia, using ambient seismic noise tomography

• Published in: Geophysical journal international Vol. 226; no. 1; pp. 511 - 523
• Main Authors: Yudistira, T, Metaxian, J-P, Putriastuti, M, Widiyantoro, S, Rawlinson, N, Beauducel, F, Zulfakriza, Z, Nugraha, A D, Laurin, A, Fahmi, A A, Budi-Santoso, A
• Format: Journal Article
• Language: English
• Published: Oxford University Press 01.07.2021; Oxford University Press (OUP)
• Subjects: Sciences of the Universe; Tomography; Seismic noise; Asia; Volcanic structure
• ISSN: 0956-540X; 1365-246X
• DOI: 10.1093/gji/ggab104

SUMMARY Mt Merapi, which lies just north of the city of Yogyakarta in Java, Indonesia, is one of the most active and dangerous volcanoes in the world. Thanks to its subduction zone setting, Mt Merapi is a stratovolcano, and rises to an elevation of 2968 m above sea level. It stands at the intersection of two volcanic lineaments, Ungaran–Telomoyo–Merbabu–Merapi (UTMM) and Lawu–Merapi–Sumbing–Sindoro–Slamet, which are oriented north–south and west–east, respectively. Although it has been the subject of many geophysical studies, Mt Merapi's underlying magmatic plumbing system is still not well understood. Here, we present the results of an ambient seismic noise tomography study, which comprise of a series of Rayleigh wave group velocity maps and a 3-D shear wave velocity model of the Merapi–Merbabu complex. A total of 10 months of continuous data (October 2013–July 2014) recorded by a network of 46 broad-band seismometers were used. We computed and stacked daily cross-correlations from every pair of simultaneously recording stations to obtain the corresponding inter-station empirical Green's functions. Surface wave dispersion information was extracted from the cross-correlations using the multiple filtering technique, which provided us with an estimate of Rayleigh wave group velocity as a function of period. The group velocity maps for periods 3–12 s were then inverted to obtain shear wave velocity structure using the neighbourhood algorithm. From these results, we observe a dominant high velocity anomaly underlying Mt Merapi and Mt Merbabu with a strike of 152°N, which we suggest is evidence of old lava dating from the UTMM double-chain volcanic arc which formed Merbabu and Old Merapi. We also identify a low velocity anomaly on the southwest flank of Merapi which we interpret to be an active magmatic intrusion.