Application of repeated passive source travel time tomography to reveal weak velocity changes related to the 2011 Tohoku-Oki M sub(w) 9.0 earthquake

• Published in: Journal of geophysical research. Solid earth Vol. 121; no. 6; pp. 4408 - 4426
• Main Authors: Koulakov, Ivan, Gladkov, Valery, El Khrepy, Sami, Al-Arifi, Nassir, Fathi, Ismail Husain
• Format: Journal Article
• Language: English
• Published: 01-06-2016
• Subjects: Algorithms; Amplitudes; Anomalies; Earth; Earthquakes; Rupture; Tomography; INW, Japan, Honshu, Miyagi Prefect., Tohoku
• ISSN: 2169-9313; 2169-9356
• DOI: 10.1002/2016JB013002

Temporal changes of seismic velocities may provide important information on the processes that occur inside the Earth. However, using body wave data with passive sources faces the problem of an uneven distribution of rays, which may cause artifacts with stronger amplitudes than the actual velocity changes in the Earth. We propose an algorithm for the selection of similar data sets in different time periods that minimize the artifacts related to variable data distributions. In this study, we used the data of the Japan Meteorological Agency for several years before and after the M sub(w) 9.0 Tohoku-Oki event that occurred on 11 March 2011. We performed careful testing using different synthetic models, showing that the selected data subsets allow detecting weak velocity changes with amplitudes above 0.2%. The analysis of the experimental data revealed important features associated with the stress and deformation distributions after the megathrust event. In the upper crust, we found a large zone along the coast with significant P velocity increase likely caused by compression of crustal rocks. This zone was cut by several elongated anomalies with local velocity decrease coinciding with the limits of the maximum slip area. These anomalies possibly mark the areas of major ruptures and deformations after the Tohoku-Oki earthquake. In the coupling zone at a depth of 40km, we observe a velocity decrease in the area of the M sub(w) 7.7 aftershock representing strong fracturing in the focal zone. Beneath the volcanic arc, we observe significant (up to 0.5%) decrease of P velocity but less prominent S velocity changes. Key Points * We present a repeated tomography algorithm based on selecting similar data subsets * Several tests proved the capacity of the algorithm to reveal weak temporal changes * In the Tohoku earthquake area, seismic velocity changes are associated with rupture processes