Imaging mantle transition zone thickness with SdS-SS finite-frequency sensitivity kernels

We invert differential SdS-SS traveltime residuals measured from stacked waveforms and finite-frequency sensitivity kernels for topography on the 410- and 660-km discontinuities. This approach yields higher resolution images of transition zone thickness than previous stacking methods, which simply average/smooth over topographic features. Apparent structure measured using simple stacking is highly dependent upon the bin size of each stack. By inverting for discontinuity topography with a variety of bin sizes, we can more accurately calculate the true structure. The inverted transition zone model is similar to simple stack models with an average thickness of 242 km, but the lateral variations in thickness are larger in amplitude and smaller in scale. Fast seismic velocities in 3-D mantle models such as SB4L18 correlate with areas of thicker transition zone. The elongated curvilinear regions of thickened transition zone that occur near subduction zones are narrow and high amplitude, which suggests relatively little lateral spreading and warming of subducted lithosphere within the transition zone. The anomalously thin transition zone regions are laterally narrow, and not broadly continuous. If these variations in transition zone thickness are interpreted as thermal in nature, then this model suggests significant temperature variations on small lateral scales.     

Inversion for laterally heterogeneous upper mantle S-wave velocity structure using iterative waveform inversion

In contrast to previous work, which treats the Earth's lateral heterogeneity as an infinitesimal perturbation to a spherically symmetrical starting model, we conduct iterative linearized waveform inversion for the Earth's laterally heterogeneous structure. We use the Direct Solution Method (DSM) (Geller et al. 1990a) to calculate synthetic seismograms and their partial derivatives for a laterally heterogeneous earth model. We invert surface-wave data from the IDA and GEOSCOPE networks. We expand the lateral heterogeneity of rigidity in spherical harmonics up to angular order number 8 and use three parameters to specify the depth dependence of each harmonic, giving us a total of 240 unknowns. The short-wavelength lateral heterogeneity (s = 4, 6 and 8) in the deeper part of the upper mantle obtained by this study differs significantly from M84A. The relative improvement in the variance reduction as compared with model M84A is about 20 per cent for the IDA data and more than 100 per cent for the GEOSCOPE data.     

Slab low-velocity layer in the eastern Aleutian subduction zone

Local earthquakes in the vicinity of the Alaskan Peninsula's Shumagin Islands often produce arrivals between the main P and S arrivals not predicted by standard traveltime tables. Based on traveltime and polarization, these anomalous arrivals appear to be from P -to- S conversions at the surface of the subducted Pacific Plate beneath the recording stations. The P -to- S conversion occurs at the top of a low-velocity layer which extends to at least 150 km depth and is 8 ˜ 2 per cent slower than the overlying mantle. The slab is ˜ 7 per cent faster than the mantle. The low-velocity layer contains the foci of the earthquakes in the upper plane of the double seismic zone and confines PS ray paths to lie within it. These observations indicate that layered structures persist to positions well past the surface location of the volcanic front. Reactions forming high-pressure minerals do not yield slab-like velocities until beyond the point that subduction zone magma genesis occurs. If the subducted oceanic crust forms the layer, it is subducted essentially intact.     

The three-dimensional shear wave structure in the mantle by overtone waveform inversion - I. Radial seismogram inversion

Summary. The three-dimensional (3-D) shear wave structure of the mantle, down to the depth of about 900 km, is obtained by inverting waveforms of radial component seismograms. Radial component seismograms contain large amplitude overtone signals which circle the Earth as wave packets and are sometimes called X1, X2, X3, … We use data which contain R1, X1 and X2 and filtered between 2 and 10mHz. It is shown that, unless each seismogram is weighted, all seismograms are not fitted uniformly. Only data from large earthquakes are fitted and the final velocity anomalies are biased by the small number of large earthquake data. Resolution is good at shallow depths, becomes worse in the intermediate depth range between about 400 and 500 km and then becomes better at greater depth ranges (600–900km). Even though we use only spheroidal mode data, velocity anomalies in the shallow structure show excellent correlation with the age of the surface rocks of the Earth. In the deeper regions, between about 600 and 900km, South America shows a fast velocity anomaly which may indicate the slab penetration beyond 700 km there. Another region which shows a fast velocity anomaly is the Mariana trench, but other subduction regions do not show such features.     

Whole mantle SH velocity model constrained by waveform inversion based on three-dimensional Born kernels

A whole mantle SH velocity model is obtained by using a unique data set and techniques. Body and surface waveforms including major and multi-orbit phases are used as a data set and are inverted by using 3-D Born kernels. The resultant model, SH18CE, reveals the different natures of the two major upwelling systems: the strong low velocity anomalies beneath Africa extend for more than 1000 km from the core–mantle boundary (CMB), whereas those beneath the Pacific are restricted to 300–400 km from the CMB. The results also show the variable natures of stagnant slabs on the 670 discontinuity around Japan: the depths of the strongest high velocity anomalies within the stagnant slabs are different region by region, which is consistent with the detailed delay time tomography model in this area.     

ScS wave splitting of deep earthquakes around Japan

ScS wave splitting of five deep earthquakes in subduction zones near Japan is investigated using horizontal seismograms recorded al JMA stations. For each earthquake, we clearly observe uniform ScS wave splitting in all stations over Japan, especially for the events located south of Honshu in 1982, 1984 and 1993. However, the directions of fast-polarized waves of these events differed by a maximum of about 50° from one another. The orientation of fast-polarized waves in the 1982 event was NNW-SSE; those in the two later events WNW-ESE. We also recognize this discrepancy in the results of the analysis of the 1971 Sea of Okhotsk event reported by Fukao (1984). The Sakhalin Islands event in 1990 reveals a linear particle motion without such a change in direction of the second arrivals, implying no anisotropy. These observations are interpreted as indicating an anisotropic region within the slab near the earthquake sources but not beneath the receivers, since the orientations of fast-polarized waves recorded at each station are not common to all the earthquakes. Furthermore, we consider that anisotropy exists non-uniformly within the slab. The event in 1982, which occurred in almost the same area as those in 1984 and 1993, showed a fast direction different from the events in 1984 and 1993. The 1982 event was 179 km deep, but the two later events were at 398 km and 360 km, respectively. The fast direction observed from the 1982 event is parallel to the fossil plate motion, whereas those from the events in 1984 and 1993 are parallel to the compression axis within the subducting slab. The depth of 400 km is a phase boundary, where olivine changes to β spinel. We consider that the most likely cause of the change in anisotropy direction is the re-orientation of crystals associated with the phase change of olivine to β spinel due to subduction of the slab.