Seismic velocity structure of the Juan de Fuca and Gorda plates revealed by a joint inversion of ambient noise and regional earthquakes

SUMMARY Seismic body wave traveltime tomography and surface wave dispersion tomography have been used widely to characterize earthquakes and to study the subsurface structure of the Earth. Since these types of problem are often significantly non-linear and have non-unique solutions, Markov chain Monte Carlo methods have been used to find probabilistic solutions. Body and surface wave data are usually inverted separately to produce independent velocity models. However, body wave tomography is generally sensitive to structure around the subvolume in which earthquakes occur and produces limited resolution in the shallower Earth, whereas surface wave tomography is often sensitive to shallower structure. To better estimate subsurface properties, we therefore jointly invert for the seismic velocity structure and earthquake locations using body and surface wave data simultaneously. We apply the new joint inversion method to a mining site in the United Kingdom at which induced seismicity occurred and was recorded on a small local network of stations, and where ambient noise recordings are available from the same stations. The ambient noise is processed to obtain inter-receiver surface wave dispersion measurements which are inverted jointly with body wave arrival times from local earthquakes. The results show that by using both types of data, the earthquake source parameters and the velocity structure can be better constrained than in independent inversions. To further understand and interpret the results, we conduct synthetic tests to compare the results from body wave inversion and joint inversion. The results show that trade-offs between source parameters and velocities appear to bias results if only body wave data are used, but this issue is largely resolved by using the joint inversion method. Thus the use of ambient seismic noise and our fully non-linear inversion provides a valuable, improved method to image the subsurface velocity and seismicity.

Download Full-text

Seismic velocity structure at the southern termination of the 2016 Kumamoto Earthquake rupture, Japan

Earth Planets and Space ◽

10.1186/s40623-020-01276-1 ◽

2020 ◽

Vol 72 (1) ◽

Author(s):

Yasuhira Aoyagi ◽

Haruo Kimura ◽

Kazuo Mizoguchi

Keyword(s):

Fault Zone ◽

Velocity Structure ◽

Seismic Velocity ◽

Earthquake Rupture ◽

2016 Kumamoto Earthquake ◽

Seismic Velocity Structure ◽

Kumamoto Earthquake ◽

Seismogenic Layer ◽

The 2016 Kumamoto Earthquake ◽

Tectonic Line

Abstract The earthquake rupture termination mechanism and size of the ruptured area are crucial parameters for earthquake magnitude estimations and seismic hazard assessments. The 2016 Mw 7.0 Kumamoto Earthquake, central Kyushu, Japan, ruptured a 34-km-long area along previously recognized active faults, eastern part of the Futagawa fault zone and northernmost part of the Hinagu fault zone. Many researchers have suggested that a magma chamber under Aso Volcano terminated the eastward rupture. However, the termination mechanism of the southward rupture has remained unclear. Here, we conduct a local seismic tomographic inversion using a dense temporary seismic network to detail the seismic velocity structure around the southern termination of the rupture. The compressional-wave velocity (Vp) results and compressional- to shear-wave velocity (Vp/Vs) structure indicate several E–W- and ENE–WSW-trending zonal anomalies in the upper to middle crust. These zonal anomalies may reflect regional geological structures that follow the same trends as the Oita–Kumamoto Tectonic Line and Usuki–Yatsushiro Tectonic Line. While the 2016 Kumamoto Earthquake rupture mainly propagated through a low-Vp/Vs area (1.62–1.74) along the Hinagu fault zone, the southern termination of the earthquake at the focal depth of the mainshock is adjacent to a 3-km-diameter high-Vp/Vs body. There is a rapid 5-km step in the depth of the seismogenic layer across the E–W-trending velocity boundary between the low- and high-Vp/Vs areas that corresponds well with the Rokkoku Tectonic Line; this geological boundary is the likely cause of the dislocation of the seismogenic layer because it is intruded by serpentinite veins. A possible factor in the southern rupture termination of the 2016 Kumamoto Earthquake is the existence of a high-Vp/Vs body in the direction of southern rupture propagation. The provided details of this inhomogeneous barrier, which are inferred from the seismic velocity structures, may improve future seismic hazard assessments for a complex fault system composed of multiple segments.

Download Full-text