Gondwanan flood basalts linked seismically to plume-induced lithosphere instability.

Delamination of the continental lithospheric mantle is well recorded beneath several continents. However, the fate of the removed continental lithosphere has been rarely noted, unlike subducted slabs reasonably well imaged in the upper and mid mantle. Beneath former Gondwana, recent seismic tomographic models indicate the presence of at least 5 horizontal fast-wavespeed anomalies at ~600 km depths that do not appear to be related to slab subduction, including fast structures in locations consistent with delamination associated with the Paraná Flood Basalt event at ~134 Ma and the Deccan Traps event at ~66 Ma. These fast-wavespeed anomalies often lie above broad slow seismic wavespeed trunks at 500 to 700 km depths beneath former Gondwana, with slow wavespeed anomalies branching around them. Numerical experiments indicate that delaminated lithosphere tends to stagnate in the transition zone and mid-mantle above a mantle plume where it shapes subsequent plume upwelling. For hot plumes, the melt volume generated during plume-influenced delamination can easily reach ~2 to 4 × 106 km3, consistent with the basalt eruption volume at the Deccan Traps. This seismic and numerical evidence suggests that observed high-wavespeed mid-mantle anomalies beneath the locations of former flood basalts are delaminated fragments of former continental lithosphere, and that lithospheric delamination events in the presence of subcontinental plumes induced several of the continental flood basalts associated with the multiple breakup stages of Gondwanaland. Continued upwelling in these plumes can also have entrained subcontinental lithosphere in the mid-mantle to bring its distinctive geochemical signal to the modern mid-ocean spreading centers that surround southern and western Africa.

New constraints on Cenozoic subduction between India and Tibet.

The type of lithosphere subducted between India and Tibet since the Paleocene remains controversial; it has been suggested to be either entirely continental, oceanic, or a mixture of the two. As the subduction history of this lost lithosphere strongly shaped Tibetan intraplate tectonism, we attempt to further constrain its nature and density structure with numerical models that aim to reproduce the observed history of magmatism and crustal thickening in addition to present-day plateau properties between 83°E and 88°E. By matching time-evolving geological patterns, here we show that Tibetan tectonism away from the Himalayan syntaxis is consistent with the initial indentation of a craton-like terrane at 55 ± 5 Ma, followed by a buoyant tectonic plate with a thin crust, e.g., a broad continental margin (Himalandia). This new geodynamic scenario can explain the seemingly contradictory observations that had led to competing hypotheses like the subduction of Greater India versus largely oceanic subduction prior to Indian indentation.

Seismic constraints and geodynamic implications of differential lithosphere-asthenosphere flow revealed in East Asia.

During the last 50 Ma, the East Asian continent has been a zone of massive continental collision and lithospheric deformation. While the consequences of this for Asian surface and lithospheric deformation have been intensively studied over the past 4 decades, the relationships between lithospheric deformation and underlying asthenospheric flow have been more difficult to constrain. Here we present a high resolution 3-D azimuthal anisotropy model for the northeastern Tibetan Plateau and its eastward continuation based on surface-wave tomography and shear-wave splitting measurements. This model shows that eastward lateral flow of asthenosphere beneath the northeastern Tibetan Plateau is being blocked by thick Ordos and Sichuan cratonic keels. The damming effect of these keels induces flow to first rotate around the Ordos keel and then transition into strong east-west flow beneath the thinner lithosphere that forms the lithospheric suture between the two cratonic keels. We further find that asthenosphere flow directions can differ from those of overlying lithosphere, with the asthenosphere neither being passively dragged by overlying lithosphere, nor being able to drag the overlying plate to mimic its subsurface flow. Finally, the region of eastward-channeled asthenospheric flow from Tibet underlies a belt of stronger intracontinental deformation in eastern China.

A strength inversion origin for non-volcanic tremor.

Non-volcanic tremor is a particularly enigmatic form of seismic activity. In its most studied subduction zone setting, tremor typically occurs within the plate interface at or near the shallow and deep edges of the interseismically locked zone. Detailed seismic observations have shown that tremor is composed of repeating small low-frequency earthquakes, often accompanied by very-low-frequency earthquakes, all involving shear failure and slip. However, low-frequency earthquakes and very-low-frequency earthquakes within each cluster show nearly constant source durations for all observed magnitudes, which implies characteristic tremor sub-event sources of near-constant size. Here we integrate geological observations and geomechanical lab measurements on heterogeneous rock assemblages representative of the shallow tremor region offshore the Middle America Trench with numerical simulations to demonstrate that these tremor events are consistent with the seismic failure of relatively weaker blocks within a stronger matrix. In these subducting rocks, hydrothermalism has led to a strength-inversion from a weak matrix with relatively stronger blocks to a stronger matrix with embedded relatively weaker blocks. Tremor naturally occurs as the now-weaker blocks fail seismically while their surrounding matrix has not yet reached a state of general seismic failure.

Transmogrification of ocean into continent: implications for continental evolution.

When continents collide, the typical embayments and protrusions along their rifted margins make it likely that fragments of seafloor will be trapped within the growing orogenic belt. These trapped seafloor fragments become preferential depocenters for marine and terrestrial sedimentation. After ∼0.5 Gy, the high radioactivity of their thick terrigenous sediment pile converts former seafloor into a unique form of continental crust and underlying lithosphere. We call this process transmogrification. Initially strong and low-lying basins that act as mechanically stronger blocks in the collisional orogeny will eventually warm, weaken, and thermoisostatically rise and will eventually transform into preferred sites for future continental rifting. In modern Asia, transmogrifying basins have induced the characteristic paired-mountain belt geomorphology associated with the assembly of this supercontinent, for example, the Himalaya/Tibet + Tian Shan surrounding the Tarim Basin that has greatly strengthened the East Asian Monsoon. The time-dependent temperature, uplift, and strength changes associated with transmogrification are relevant for improving our understanding of continental evolution, basin modeling, paleoclimate studies, and natural resources prospection.

Mechanism of progressive broad deformation from oceanic transform valley to off-transform faulting and rifting.

Oceanic transform faults (TFs) are commonly viewed as single, narrow strike-slip seismic faults that offset two mid-ocean ridge segments. However, broad zones of complex deformation are ubiquitous at TFs. Here, we propose a new conceptual model for the progressive deformation within broad zones at oceanic TFs through detailed morphological, seismic, and stress analyses. We argue that, under across-transform extension due to a change in plate motion, plate deformation occurs first along high-angle transtensional faults (TTFs) within the transform valleys. Off-transform normal faults (ONFs) form when across-transform deviatoric extensional stresses exceed the yield strength of the adjacent oceanic lithosphere. With further extension, these normal faults can develop into off-transform rift zones (ORZs), some of which can further develop into transform plate boundaries. We illustrate that such progressive complex deformation is an inherent feature of oceanic TFs. The new conceptual model provides a unifying theory to explain the observed broad deformation at global transform systems.

Extensional tectonics and two-stage crustal accretion at oceanic transform faults.

Oceanic transform faults are seismically and tectonically active plate boundaries1 that leave scars-known as fracture zones-on oceanic plates that can cross entire ocean basins2. Current descriptions of plate tectonics assume transform faults to be conservative two-dimensional strike-slip boundaries1,3, at which lithosphere is neither created nor destroyed and along which the lithosphere cools and deepens as a function of the age of the plate4. However, a recent compilation of high-resolution multibeam bathymetric data from 41 oceanic transform faults and their associated fracture zones that covers all possible spreading rates shows that this assumption is incorrect. Here we show that the seafloor along transform faults is systemically deeper (by up to 1.6 kilometres) than their associated fracture zones, in contrast to expectations based on plate-cooling arguments. Accretion at intersections between oceanic ridges and transform faults seems to be strongly asymmetric: the outside corners of the intersections show shallower relief and more extensive magmatism, whereas the inside corners have deep nodal basins and seem to be magmatically starved. Three-dimensional viscoplastic numerical models show that plastic-shear failure within the deformation zone around the transform fault results in the plate boundary experiencing increasingly oblique shear at increasing depths below the seafloor. This results in extension around the inside corner, which thins the crust and lithosphere at the transform fault and is linked to deepening of the seafloor along the transform fault. Bathymetric data suggest that the thinned transform-fault crust is augmented by a second stage of magmatism as the transform fault intersects the opposing ridge axis. This makes accretion at transform-fault systems a two-stage process, fundamentally different from accretion elsewhere along mid-ocean ridges.

Causes and consequences of asymmetric lateral plume flow during South Atlantic rifting.

Volcanic rifted margins are typically associated with a thick magmatic layer of seaward dipping reflectors and anomalous regional uplift. This is conventionally interpreted as due to melting of an arriving mantle plume head at the onset of rifting. However, seaward dipping reflectors and uplift are sometimes asymmetrically distributed with respect to the subsequent plume track. Here we investigate if these asymmetries are induced by preexisting lateral variations in the thickness of continental lithosphere and/or lithospheric stretching rates, variations that promote lateral sublithospheric flow of plume material below only one arm of the extending rift. Using three-dimensional numerical experiments, we find that South Atlantic rifting is predicted to develop a strong southward asymmetry in its distribution of seaward dipping reflectors and associated anomalous relief with respect to the Tristan Plume that "drove" this volcanic rifted margin, and that the region where plume material drains into the rift should experience long-lived uplift during rifting-both as observed. We conclude that a mantle plume is still needed to source the anomalously hot sublithospheric material that generates a volcanic rifted margin, but lateral along-rift flow from this plume, not a broad starting plume head, is what controls when and where a volcanic rifted margin will form.

How and when plume zonation appeared during the 132 Myr evolution of the Tristan Hotspot.

Increasingly, spatial geochemical zonation, present as geographically distinct, subparallel trends, is observed along hotspot tracks, such as Hawaii and the Galapagos. The origin of this zonation is currently unclear. Recently zonation was found along the last ∼70 Myr of the Tristan-Gough hotspot track. Here we present new Sr-Nd-Pb-Hf isotope data from the older parts of this hotspot track (Walvis Ridge and Rio Grande Rise) and re-evaluate published data from the Etendeka and Parana flood basalts erupted at the initiation of the hotspot track. We show that only the enriched Gough, but not the less-enriched Tristan, component is present in the earlier (70-132 Ma) history of the hotspot. Here we present a model that can explain the temporal evolution and origin of plume zonation for both the Tristan-Gough and Hawaiian hotspots, two end member types of zoned plumes, through processes taking place in the plume sources at the base of the lower mantle.

Hybrid shallow on-axis and deep off-axis hydrothermal circulation at fast-spreading ridges.

Hydrothermal flow at oceanic spreading centres accounts for about ten per cent of all heat flux in the oceans and controls the thermal structure of young oceanic plates. It also influences ocean and crustal chemistry, provides a basis for chemosynthetic ecosystems, and has formed massive sulphide ore deposits throughout Earth's history. Despite this, how and under what conditions heat is extracted, in particular from the lower crust, remains largely unclear. Here we present high-resolution, whole-crust, two- and three-dimensional simulations of hydrothermal flow beneath fast-spreading ridges that predict the existence of two interacting flow components, controlled by different physical mechanisms, that merge above the melt lens to feed ridge-centred vent sites. Shallow on-axis flow structures develop owing to the thermodynamic properties of water, whereas deeper off-axis flow is strongly shaped by crustal permeability, particularly the brittle-ductile transition. About 60 per cent of the discharging fluid mass is replenished on-axis by warm (up to 300 degrees Celsius) recharge flow surrounding the hot thermal plumes, and the remaining 40 per cent or so occurs as colder and broader recharge up to several kilometres away from the axis that feeds hot (500-700 degrees Celsius) deep-rooted off-axis flow towards the ridge. Despite its lower contribution to the total mass flux, this deep off-axis flow carries about 70 per cent of the thermal energy released at the ridge axis. This combination of two flow components explains the seismically determined thermal structure of the crust and reconciles previously incompatible models favouring either shallower on-axis or deeper off-axis hydrothermal circulation.