The composition of the Earth

We relocate nearly 100,000 events that occurred during the period 1964 to 1995 and are well-constrained teleseismically by arrival-time data reported to the International Seismological Centre (ISC) and to the U.S. Geological Survey9s National Earthquake Information Center (NEIC). Hypocenter determination is significantly improved by using, in addition to regional and teleseismic P and S phases, the arrival times of PKiKP, PKPdf , and the teleseismic depth phases pP, pwP , and sP in the relocation procedure. A global probability model developed for later-arriving phases is used to independently identify the depth phases. The relocations are compared to hypocenters reported in the ISC and NEIC catalogs and by other sources. Differences in our epicenters with respect to ISC and NEIC estimates are generally small and regionally systematic due to the combined effects of the observing station network and plate geometry regionally, differences in upper mantle travel times between the reference earth models used, and the use of later-arriving phases. Focal depths are improved substantially over most other independent estimates, demonstrating (for example) how regional structures such as downgoing slabs can severely bias depth estimation when only regional and teleseismic P arrivals are used to determine the hypocenter. The new data base, which is complete to about Mw 5.2 and includes all events for which moment-tensor solutions are available, has immediate application to high-resolution definition of Wadati-Benioff Zones (WBZs) worldwide, regional and global tomographic imaging, and other studies of earth structure.

Global teleseismic earthquake relocation with improved travel times and procedures for depth determination

Seismic tomography models of the three-dimensional upper mantle velocity structure of the Mediterranean-Carpathian region provide a better understanding of the lithospheric processes governing its geodynamical evolution. Slab detachment, in particular lateral migration of this process along the plate boundary, is a key element in the lithospheric dynamics of the region during the last 20 to 30 million years. It strongly affects arc and trench migration, and causes along-strike variations in vertical motions, stress fields, and magmatism. In a terminal-stage subduction zone, involving collision and suturing, slab detachment is the natural last stage in the gravitational settling of subducted lithosphere.

http://www.mantleplumes.org/WebDocuments/Wortel2000.pdf

Subduction and Slab Detachment in the Mediterranean-Carpathian Region

Seismic tomography based on P-wave travel times and improved earthquake locations provides further evidence for mantle-wide convective flow. The use of body waves makes it possible to resolve long, narrow structures in the lower mantle some of which can be followed to sites of present-day plate convergence at the Earth's surface. The transition from subduction-related linear structures in the mid-mantle to long-wavelength heterogeneity near the core-mantle boundary remains enigmatic, but at least some slab segments seem to sink to the bottom of the mantle.

/pdf/evidence-for-deep-mantle-circulation-from-global-tomography-3unqcjsa05.pdf

Evidence for deep mantle circulation from global tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al. [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [Spakman and Bijwaard, 1998]. Second, we exploit the global data set of Engdahl et al. [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used (P, pP, and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/98JB02467

Closing the gap between regional and global travel time tomography

The seismic tomography problem does not have a unique solution, and published tomographic images have been equivocal with regard to the deep structure of subducting slabs. An improved tomographic method, using a more realistic background Earth model and surface-reflected as well as direct seismic phases, shows that slabs beneath the Japan and Izu Bonin island arcs are deflected at the boundary between upper and lower mantle, whereas those beneath the northern Kuril and Mariana arcs sink into the lower mantle.

/pdf/tomographic-imaging-of-subducted-lithosphere-below-northwest-1g3v35uqe6.pdf

Robert Engdahl

Papers

Tomographic imaging of subducted lithosphere below northwest Pacific island arcs