Seismic Structure of Upper-Mantle Discontinuities Based on High-Frequency Triplication Data

Constraining seismic properties of the 410- and 660-km discontinuities which delineate the mantle transition zone, is crucial in understanding the mantle composition and convection dynamics. One approach to studying the transition zone is to use “triplicated” arrivals of seismic data. When a wavefield encounters an abrupt increase in wave speed with depth, three different phases, direct, reflected, and transmitted phases, result, and they can arrive at the same distance. In the traveltime curve, the three arrivals appear to make a bowtie shape (Figure 1 Top Left), hence the name triplication. Triplication pattern and its location, i.e., travel times of three types of waves, provide constraints on the depth and velocity jump at the discontinuity.

Figure 1: Predicted and observed data in T-D domain and τ-p domain. (Top) The bowtie shaped travel time curve becomes unwrapped as an “ω” shaped τ-p curve after Radon transform. Each phase is separated from one another with different signs of curvature: (+) for prograde and (-) for retrograde. (Bottom) Radon transform of noisy data shows clear τ-p curve and brings out the information from the data.

One of the challenging components in using triplication data, however, is to identify the three individual phases, since they arrive close in time and overlap with each other (Figure 1 Botton Left). In order to separate the three phases, we apply Radon transform to short-period seismograms recorded by a dense array of stations. This approach unwraps the triplication pattern, and brings out the high-frequency information that is not easily accessible in the original form of data (Figure 1 Right). Subsequent modeling of the unwrapped data allows velocity jump, depth, and width of the discontinuities to be obtained (Figure 2).

Figure 2: Investigated parameters. Velocities at the top (V1) and bottom (V2), depth (D), and width (W) of the discontinuity

This method is applied to study the transition zone around the Kuril subduction zone, a region northeast of Japan. We take advantage of the High-Sensitivity Seismograph Network in Japan that consists of more than 700 stations. These stations provide dense sampling in distance that allows us to capture the triplication pattern. The wave speeds immediately above and below the 410- and 660-km discontinuities as well as their depths and widths are constrained. In general, both discontinuities are depressed compared to the global average, and exhibit finite widths. The width estimates have implications on effects such as the existence of water or melt, and garnet transformations occurring at similar depths as the post-spinel transition.

Publications:

Park, S., Ishii, M., 2019, Double 660-km Discontinuity as a Selective Barrier to Mantle Flow, under review in Nature.

Park, S., Ishii, M., 2019, Patchy Melt Distribution Inferred from Missing 410-km Discontinuity, in prep.

Park, S., Ishii, M., 2019, A Novel Approach Utilizing High-Frequency Triplication Data to Study Local Variations in the Upper-Mantle Transition Zone, in prep.