*Ryosuke AZUMA1, Fumiaki TOMITA2, Takeshi IINUMA1, Ryota HINO1, Motoyuki KIDO1, Hiromi FUJIMOTO1, Yukihito OSADA1, Yusaku OHTA3, Ikuko WADA1
(1.IRIDeS, Tohoku Univ., 2.Graduate School of Science, Tohoku Univ., 3.RCPEVE, Tohoku Univ.)
Keywords:ocean bottom geodetic observation
The development of a technique for GPS-acoustic (GPS/A) geodetic observations has enabled us to understand the slip distribution of the 2011 Tohoku-oki earthquake. However, there remains an issue with the precision of GPS/A measurement is still lower by two orders of magnitude than that of on-land GPS measurement due to problems in observations and data processing methods. In this study, we focus on the problem for determination of traveltimes of acoustic signals obtained from GPS/A measurement.The conventional approach for determining the two-way traveltime of observed acoustic signals is to determine the maximum peak of the cross-correlation waveforms between the transmitted and returned signals. However, the maximum peak often differs from the true peak due to the distortion in the correlation waveform which depends on the relative spatial geometry of the ship and station. These misread traveltimes have been re-read manually so far. Such procedure is no longer applicable for processing vast array of data obtained at newly installed over 20 GPS/A stations after the 2011 Tohoku-oki earthquake. The aim of this study is to develop fully automated algorithms for analyzing GPS/A data with high precision. We introduce here two algorithms.1) We read the maximum peak in the observed correlogram and then deconvolve it by the synthetic correlogram. Then, we apply the same operation to the deconvolved waveform. This procedure is iterated until the correlation coefficient decreases lower than a pre-defined threshold. A true traveltime is defined as the fastest traveltime during the iterations.2) We classify the observed correlograms into several groups based on their similarity through cluster analyses and choose a master waveform in each group. Then evaluate the traveltime residual between the maximum peak and the true peak in the observed correlogram. Thus obtained residual is applied as the correction value of each clustered group.We also use a seismic data analysis tool to visually inspect whether above algorithms work properly. We confirmed that the both new methods properly correct for misreadings in the current method, which sometimes amount to several hundred micro seconds. This corresponds roughly to a 0.3-m difference in the slant range. Therefore, with the new algorithms, significant improvement in the estimation of the station location is expected. However, both methods have to be assigned an arbitrary value as a threshold. Further analyses are needed to determine arbitrary threshold values and to construct fully automated algorithms.