Electrical conductivity in the mantle transition zone beneath eastern Central Asian Orogenic Belt revealed by geomagnetic signals

In the eastern segment of the Central Asian Orogenic Belt (CAOB), there is widespread volcanic magma activity. However, there is still considerable controversy over the formation mechanisms and material sources of these volcanoes. The mantle transition zone (MTZ), as a necessary channel for the upward and downward movement of mantle material and energy exchange may provide crucial constraints on the dynamic mechanisms of volcanic activity. This paper intends to obtain the deep structure beneath the eastern CAOB based on the geomagnetic depth sounding (GDS) method. First, the data of geomagnetic observatories in the study region are collected and processed, and the C-response curves are obtained by the bounded influence remote reference processing method (BIRRP). Then, the staggered grid finite difference method is used for forward modeling, and the finite memory quasi-Newton method based on L1-norm is used for three-dimensional (3-D) inversion. After that, 3-D inversion is carried out in spherical coordinates. Finally, the electrical conductivity model is obtained. The inversion model shows that there are two high conductivity anomalies in the MTZ beneath the Mongol-Okhotsk suture. Combined with the geological background of the structural domain, and constrained by the spatiotemporal variations in magmatism, we speculate that the high conductivity anomaly bodies are the stagnant oceanic crust material of the Okhotsk Ocean or the delaminated island arc accretionary wedge. The sinking slab or the detached lithosphere residual descending into the lower MTZ causes the upwelling of hot mantle material, forming widely distributed volcanic rocks on both sides of the Mongol-Okhotsk suture.

Slab Geometry and Upper Mantle Flow Patterns in the Central Mediterranean From 3D Anisotropic P-Wave Tomography

We present the first three-dimensional (3D) anisotropic teleseismic P-wave tomography model of the upper mantle covering the entire Central Mediterranean. Compared to isotropic tomography, it is found that including the magnitude, azimuth, and, importantly, dip of seismic anisotropy in our inversions simplifies isotropic heterogeneity by reducing the magnitude of slow anomalies while yielding anisotropy patterns that are consistent with regional tectonics. The isotropic component of our preferred tomography model is dominated by numerous fast anomalies associated with retreating, stagnant, and detached slab segments. In contrast, relatively slower mantle structure is related to slab windows and the opening of back-arc basins. To better understand the complexities in slab geometry and their relationship to surface geological phenomenon, we present a 3D reconstruction of the main Central Mediterranean slabs down to 700 km based on our anisotropic model. P-wave seismic anisotropy is widespread in the Central Mediterranean upper mantle and is strongest at 200-300 km depth. The anisotropy patterns are interpreted as the result of asthenospheric material flowing primarily horizontally around the main slabs in response to pressure exerted by their mid-to-late Cenezoic horizontal motion, while sub-vertical anisotropy possibly reflects asthenospheric entrainment by descending lithosphere. Our results highlight the importance of anisotropic P-wave imaging for better constraining regional upper mantle geodynamics.

The closure of the Vardar Ocean (the western domain of the northern Neotethys) from the early Middle Jurassic to the Paleocene time, based on the surface geology of eastern Pelagonia and the Vardar zone, biostratigraphy, and seismic-tomographic images of the mantle below the Central Hellenides

Seismic tomographic images of the mantle below the Hellenides indicate that the Vardar Ocean probably had a composite width of over 3000 km. From surface geology we know that this ocean was initially located between two passive margins: Pelagonian Adria in the west and Serbo-Macedonian-Eurasia in the east. Pelagonia was covered by a carbonate platform that accumulated, during Late Triassic to Early Cretaceous time, where highly diversified carbonate sedimentary environments evolved and reacted to the adjacent, converging Vardar Ocean plate. We conceive that on the east side of the Vardar Ocean, a Cretaceous carbonate platform evolved from the Aptian to the Maastrichtian time in the forearc basin of the Vardar supra-subduction volcanic arc complex. The closure of the Vardar Ocean occurred in one episode of ophiolite obduction and in two episodes of intra-oceanic subduction. 1. During the Middle Jurassic time a 1200-km slab of west Vardar lithosphere subducted beneath the supra-subduction, 'Eohellenic', arc, while a 200-km-wide slab obducted onto Pelagonia between the Callovian and Valanginian times. 2. During the Late Jurassic through to the Cretaceous time a 1700-km-wide slab subducted beneath the evolving east Vardar-zone arc-complex. Pelagonia, the trailing edge of the subducting east-Vardar Ocean slab, crashed and underthrust the Vardar arc complex during the Paleocene time and ultimately crashed with Serbo-Macedonia. Since the late Early Jurassic time, the Hellenides have moved about 3000 km toward the northeast while the Atlantic Ocean spread.

Tectono-stratigraphic correlations between Northern Evvoia, Skopelos and Alonnisos, and the postulated collision of the Pelagonian carbonate platform with the Paikon forearc basin (Pelagonian-Vardar zones, Internal Hellenides, Greece)

The Pelagonian stratigraphy of the Internal Hellenides consists of a Permo-Triassic basement and an Upper Triassic and Jurassic carbonate platform formation that has been overthrust by the Eohellenic ophiolite sheet during the Early Cretaceous. Intensive erosion, during the Cretaceous, removed most of the ophiolite and parts of the Jurassic formation. It is hypothesised that uplift and erosion of eastern Pelagonia was triggered by the break-off of the subducted oceanic leading edge of the Pelagonian plate. An investigation of the rocks that succeed the erosional unconformity shows that they constitute a shear-zone that is tectonically overlain by Cretaceous platform carbonates. Geochemical analyses of the shear-zone rocks substantiate that they are of mid-oceanic ridge and island arc provenience. Eastern Pelagonia collided with a Cretaceous carbonate platform, probably the Paikon forearc basin, as the Almopias ocean crust subducted beneath that island-arc complex. The Cretaceous platform, together with a substrate of sheared-off ocean floor mélange, overthrust eastern Pelagonia as subduction continued, and the substrate was dynamically metamorphosed into cataclastic rocks, mylonite, phyllonite and interpreted pseudotachylite. This complex of Cretaceous platform rocks and a brittle-ductile shear-zone-substrate constitute the here named Paikon-Palouki nappe, which was emplaced during Early Palaeocene. The Paikon-Palouki nappe did not reach Evvoia. Seismic tomographic models of the Aegean region apparently depict images of two broken-off ocean-plate-slabs, interpreted as Almopias-lithosphere-slabs. It is concluded that the western Almopias slab began to sink during the Early Cretaceous, while the eastern Almopias slab broke off and sank after the Paikon-Palouki nappe was emplaced in the Early Palaeocene.

Southwest Pacific Absolute Plate Kinematic Reconstruction Reveals Major Cenozoic Tonga-Kermadec Slab Dragging

Tectonic plates subducting at trenches having strikes oblique to the absolute subducting plate motion undergo trench-parallel slab motion through the mantle, recently defined as a form of "slab dragging." We investigate here long-term slab-dragging components of the Tonga-Kermadec subduction system driven by absolute Pacific plate motion. To this end we develop a kinematic restoration of Tonga-Kermadec Trench motion placed in a mantle reference frame and compare it to tomographically imaged slabs in the mantle. Estimating Tonga-Kermadec subduction initiation is challenging because another (New Caledonia) subduction zone existed during the Paleogene between the Australia and Pacific plates. We test partitioning of plate convergence across the Paleogene New Caledonia and Tonga-Kermadec subduction zones against resulting mantle structure and show that most, if not all, Tonga-Kermadec subduction occurred after ca. 30 Ma. Since then, Tonga-Kermadec subduction has accommodated 1,700 to 3,500 km of subduction along the southern and northern ends of the trench, respectively. When placed in a mantle reference frame, the predominantly westward directed subduction evolved while the Tonga-Kermadec Trench underwent ~1,200 km of northward absolute motion. We infer that the entire Tonga-Kermadec slab was laterally transported through the mantle over 1,200 km. Such slab dragging by the Pacific plate may explain observed deep-slab deformation and may also have significant effects on surface tectonics, both resulting from the resistance to slab dragging by the viscous mantle.

The Collaborative Seismic Earth Model: Generation 1

We present a general concept for evolutionary, collaborative, multiscale inversion of geophysical data, specifically applied to the construction of a first-generation Collaborative Seismic Earth Model. This is intended to address the limited resources of individual researchers and the often limited use of previously accumulated knowledge. Model evolution rests on a Bayesian updating scheme, simplified into a deterministic method that honors today's computational restrictions. The scheme is able to harness distributed human and computing power. It furthermore handles conflicting updates, as well as variable parameterizations of different model refinements or different inversion techniques. The first-generation Collaborative Seismic Earth Model comprises 12 refinements from full seismic waveform inversion, ranging from regional crustal- to continental-scale models. A global full-waveform inversion ensures that regional refinements translate into whole-Earth structure.

On the consistency of seismically imaged lower mantle slabs

The geoscience community is increasingly utilizing seismic tomography to interpret mantle heterogeneity and its links to past tectonic and geodynamic processes. To assess the robustness and distribution of positive seismic anomalies, inferred as subducted slabs, we create a set of vote maps for the lower mantle with 14 global P-wave or S-wave tomography models. Based on a depth-dependent threshold metric, an average of 20% of any given tomography model depth is identified as a potential slab. However, upon combining the 14 models, the most consistent positive wavespeed features are identified by an increasing vote count. An overall peak in the most robust anomalies is found between 1000-1400 km depth, followed by a decline to a minimum around 2000 km. While this trend could reflect reduced tomographic resolution in the middle mantle, we show that it may alternatively relate to real changes in the time-dependent subduction flux and/or a mid-lower mantle viscosity increase. An apparent secondary peak in agreement below 2500 km depth may reflect the degree-two lower mantle slow seismic structures. Vote maps illustrate the potential shortcomings of using a limited number or type of tomography models and slab threshold criteria.

Slab detachment under the Eastern Alps seen by seismic anisotropy

We analyze seismic anisotropy for the Eastern Alpine region by inspecting shear-wave splitting from SKS and SKKS phases. The Eastern Alpine region is characterized by a breakdown of the clear mountain-chain-parallel fast orientation pattern that has been previously documented for the Western Alps and for the western part of the Eastern Alps. The main interest of this paper is a more detailed analysis of the anisotropic character of the Eastern Alps, and the transition to the Carpathian-Pannonian region. SK(K)S splitting measurements reveal a rather remarkable lateral change in the anisotropy pattern from the west to the east of the Eastern Alps with a transition area at about 12°E. We also model the backazimuthal variation of the measurements by a vertical change of anisotropy. We find that the eastern part of the study area is characterized by the presence of two layers of anisotropy, where the deeper layer has characteristics similar to those of the Central Alps, in particular SW-NE fast orientations of anisotropic axes. We attribute the deeper layer to a detached slab from the European plate. Comparison with tomographic studies of the area indicates that the detached slab might possibly connect with the lithosphere that is still in place to the west of our study area, and may also connect with the slab graveyard to the East, at the depth of the upper mantle transition zone. On the other hand, the upper layer has NW-SE fast orientations coinciding with a low-velocity layer which is found above a more-or-less eastward dipping high-velocity body. The anisotropy of the upper layer shows large-scale NW-SE fast orientation, which is consistent with the presence of asthenospheric flow above the detached slab foundering into the deeper mantle.

Subsidence and uplift by slab-related mantle dynamics: a driving mechanism for the Late Cretaceous and Cenozoic evolution of continental SE Asia?

Continental SE Asia is the site of an extensive Cretaceous–Paleocene regional unconformity that extends from Indochina to Java, covering an area of c. 5 600 000 km2. The unconformity has previously been related to microcontinental collision at the Java margin that halted subduction of Tethyan oceanic lithosphere in the Late Cretaceous. However, given the disparity in size between the accreted continental fragments and area of the unconformity, together with lack of evidence for requisite crustal shortening and thickening, the unconformity is unlikely to have resulted from collisional tectonics alone. Instead, mapping of the spatial extent of the mid–Late Cretaceous subduction zone and the Cretaceous–Paleocene unconformity suggests that the unconformity could be a consequence of subduction-driven mantle processes. Cessation of subduction, descent of a northward dipping slab into the mantle, and consequent uplift and denudation of a sediment-filled Late Jurassic and Early Cretaceous dynamic topographic low help explain the extent and timing of the unconformity. Sediments started to accumulate above the unconformity from the Middle Eocene when subduction recommenced beneath Sundaland.

Australia–SE Asia collision: plate tectonics and crustal flow

The Sundaland core of SE Asia is a heterogeneous assemblage of Tethyan sutures and Gondwana fragments. Its complex basement structure was one major influence on Cenozoic tectonics; the rifting history of the north Australian margin was another. Fragments that rifted from Australia in the Jurassic collided with Sundaland in the Cretaceous and terminated subduction. From 90 to 45 Ma Sundaland was largely surrounded by inactive margins with localized strike-slip deformation, extension and subduction. At 45 Ma Australia began to move north, and subduction resumed beneath Sundaland. At 23 Ma the Sula Spur promontory collided with the Sundaland margin. From 15 Ma there was subduction hinge rollback into the Banda oceanic embayment, major extension, and later collision of the Banda volcanic arc with the southern margin of the embayment. However, this plate tectonic framework cannot be reduced to a microplate scale to explain Cenozoic deformation. Sundaland has a weak thin lithosphere, highly responsive to plate boundary forces and a hot weak deep crust has flowed in response to tectonic and topographic forces, and sedimentary loading. Gravity-driven movements of the upper crust, unusually rapid vertical motions, exceptionally high rates of erosion, and massive movements of sediment have characterized this region.

Evidence for a chemical-thermal structure at base of mantle from sharp lateral P-wave variations beneath Central America

Compressional waves that sample the lowermost mantle west of Central America show a rapid change in travel times of up to 4 s over a sampling distance of 300 km and a change in waveforms. The differential travel times of the PKP waves (which traverse Earth's core) correlate remarkably well with predictions for S-wave tomography. Our modeling suggests a sharp transition in the lowermost mantle from a broad slow region to a broad fast region with a narrow zone of slowest anomaly next to the boundary beneath the Cocos Plate and the Caribbean Plate. The structure may be the result of ponding of ancient subducted Farallon slabs situated near the edge of a thermal and chemical upwelling.

Finite-frequency tomography reveals a variety of plumes in the mantle

We present tomographic evidence for the existence of deep-mantle thermal convection plumes. P-wave velocity images show at least six well-resolved plumes that extend into the lowermost mantle: Ascension, Azores, Canary, Easter, Samoa, and Tahiti. Other less well-resolved plumes, including Hawaii, may also reach the lowermost mantle. We also see several plumes that are mostly confined to the upper mantle, suggesting that convection may be partially separated into two depth regimes. All of the observed plumes have diameters of several hundred kilometers, indicating that plumes convey a substantial fraction of the internal heat escaping from Earth.

Seismic evidence for water deep in Earth's upper mantle

Water in the deep upper mantle can influence the properties of seismic discontinuities in the mantle transition zone. Observations of converted seismic waves provide evidence of a 20- to 35-kilometer-thick discontinuity near a depth of 410 kilometers, most likely explained by as much as 700 parts per million of water by weight.

Chemical and seismological constraints on mantle heterogeneity

Recent seismological studies that use scattered waves to detect heterogeneities in the mantle reveal the presence of a small, distributed elastic heterogeneity in the lower mantle which does not appear to be thermal in nature. The characteristic size of these heterogeneities appears to be ca. 8 km, suggesting that they represent subducted recycled oceanic crust. With this stimulus, old ideas that the mantle is heterogeneous in structure, rather than stratified, are reinterpreted and a simple, end-member model for the heterogeneity structure is proposed. The volumetrically largest components in the model are recycled oceanic crust, which contains the heat-producing elements, and mantle depleted of these and other incompatible trace elements. About 10% of the mantle's mass is made up of recycled oceanic crust, which is associated with the observed small-scale seismic heterogeneity. The way this heterogeneity is distributed is in convectively stretched and thinned bodies ranging downwards in size from 8 km. With the present techniques to detect small bodies through scattering, only ca. 55% of the mantle's small-scale heterogeneities are detectable seismically.

Mantle shear-wave tomography and the fate of subducted slabs

A new seismic model of the three-dimensional variation in shear velocity throughout the Earth's mantle is presented. The model is derived entirely from shear bodywave travel times. Multibounce shear waves, core-reflected waves and SKS and SKKS waves that travel through the core are used in the analysis. A unique aspect of the dataset used in this study is the use of bodywaves that turn at shallow depths in the mantle, some of which are triplicated. The new model is compared with other global shear models. Although competing models show significant variations, several large-scale structures are common to most of the models. The high-velocity anomalies are mostly associated with subduction zones. In some regions the anomalies only extend into the shallow lower mantle, whereas in other regions tabular high-velocity structures seem to extend to the deepest mantle. The base of the mantle shows long-wavelength high-velocity zones also associated with subduction zones. The heterogeneity seen in global tomography models is difficult to interpret in terms of mantle flow due to variations in structure from one subduction zone to another. The simplest interpretation of the seismic images is that slabs in general penetrate to the deepest mantle, although the flow is likely to be sporadic. The interruption in slab sinking is likely to be associated with the 660 km discontinuity.