Early Results of Time Domain Receiver Function Data Processing in Mt Merapi and Mt Merbabu

Abstract

The receiver function method is a method to image the earth subsurface by utilizing Ps conversion waves that are formed due to the velocity boundary. In general, the receiver function method estimates depth of structures such as the mantle-crust boundary by deconvoluting the vertical component from the horizontal component. Typical receiver function data processing is done in the frequency domain where the deconvolution process can be seen as a division between two components. In this study, we tried to reprocess the data using a deconvolution technique in time domain, popularly known as iterative time-domain deconvolution. The principle of iterative time domain deconvolution consists of iterative cross-correlation between the horizontal and vertical component. We use data from the DOMERAPI seismic station network located in the vicinity of Mt Merapi and Mt Merbabu. Mt Merapi is one of the most active volcanoes in the world with frequent eruptions and located at the ring of fire chain volcano in Indonesia. Note that the previous receiver function study in this region showed complex signals at some stations that may be related to sediment at shallow sediment and possible layers of low velocity zone that interfering main signal for a crust-mantle boundary. Our current results show iterative time domain RFs have clearer and smoother signal than the frequency domain that help interpreting the waveform signals. We estimate a range of crust thickness between 26-31 km near Mt Merapi. However, we noticed that iterative time domain calculation requires longer computation time and input signal.