Constraining Earth's nonlinear mantle viscosity using plate-boundary resolving global inversions.

Variable viscosity in Earth's mantle exerts a fundamental control on mantle convection and plate tectonics, yet rigorously constraining the underlying parameters has remained a challenge. Inverse methods have not been sufficiently robust to handle the severe viscosity gradients and nonlinearities (arising from dislocation creep and plastic failure) while simultaneously resolving the megathrust and bending slabs globally. Using global plate motions as constraints, we overcome these challenges by combining a scalable nonlinear Stokes solver that resolves the key tectonic features with an adjoint-based Bayesian approach. Assuming plate cooling, variations in the thickness of continental lithosphere, slabs, and broad scale lower mantle structure as well as a constant grain size through the bulk of the upper mantle, a good fit to global plate motions is found with a nonlinear upper mantle stress exponent of 2.43 [Formula: see text] 0.25 (mean [Formula: see text] SD). A relatively low yield stress of 151 [Formula: see text] 19 MPa is required for slabs to bend during subduction and transmit a slab pull that generates asymmetrical subduction. The recovered long-term strength of megathrusts (plate interfaces) varies between different subduction zones, with South America having a larger strength and Vanuatu and Central America having lower values with important implications for the stresses driving megathrust earthquakes.

Southward expanding plate coupling due to variation in sediment subduction as a cause of Andean growth.

Growth of the Andes has been attributed to Cenozoic subduction. Although climatic and tectonic processes have been proposed to be first-order mechanisms, their interaction and respective contributions remain largely unclear. Here, we apply three-dimensional, fully-dynamic subduction models to investigate the effect of trench-axial sediment transport and subduction on Andean growth, a mechanism that involves both climatic and tectonic processes. We find that the thickness of trench-fill sediments, a proxy of plate coupling (with less sediments causing stronger coupling), exerts an important influence on the pattern of crustal shortening along the Andes. The southward migrating Juan Fernandez Ridge acts as a barrier to the northward flowing trench sediments, thus expanding the zone of plate coupling southward through time. Consequently, the predicted history of Andean shortening is consistent with observations. Southward expanding crustal shortening matches the kinematic history of inferred compression. These results demonstrate the importance of climate-tectonic interaction on mountain building.

A rapid burst in hotspot motion through the interaction of tectonics and deep mantle flow.

Volcanic hotspot tracks featuring linear progressions in the age of volcanism are typical surface expressions of plate tectonic movement on top of narrow plumes of hot material within Earth's mantle. Seismic imaging reveals that these plumes can be of deep origin--probably rooted on thermochemical structures in the lower mantle. Although palaeomagnetic and radiometric age data suggest that mantle flow can advect plume conduits laterally, the flow dynamics underlying the formation of the sharp bend occurring only in the Hawaiian-Emperor hotspot track in the Pacific Ocean remains enigmatic. Here we present palaeogeographically constrained numerical models of thermochemical convection and demonstrate that flow in the deep lower mantle under the north Pacific was anomalously vigorous between 100 million years ago and 50 million years ago as a consequence of long-lasting subduction systems, unlike those in the south Pacific. These models show a sharp bend in the Hawaiian-Emperor hotspot track arising from the interplay of plume tilt and the lateral advection of plume sources. The different trajectories of the Hawaiian and Louisville hotspot tracks arise from asymmetric deformation of thermochemical structures under the Pacific between 100 million years ago and 50 million years ago. This asymmetric deformation waned just before the Hawaiian-Emperor bend developed, owing to flow in the deepest lower mantle associated with slab descent in the north and south Pacific.

The dynamics of plate tectonics and mantle flow: from local to global scales.

Plate tectonics is regulated by driving and resisting forces concentrated at plate boundaries, but observationally constrained high-resolution models of global mantle flow remain a computational challenge. We capitalized on advances in adaptive mesh refinement algorithms on parallel computers to simulate global mantle flow by incorporating plate motions, with individual plate margins resolved down to a scale of 1 kilometer. Back-arc extension and slab rollback are emergent consequences of slab descent in the upper mantle. Cold thermal anomalies within the lower mantle couple into oceanic plates through narrow high-viscosity slabs, altering the velocity of oceanic plates. Viscous dissipation within the bending lithosphere at trenches amounts to approximately 5 to 20% of the total dissipation through the entire lithosphere and mantle.

Reconstructing Farallon plate subduction beneath North America back to the Late Cretaceous.

Using an inverse mantle convection model that assimilates seismic structure and plate motions, we reconstruct Farallon plate subduction back to 100 million years ago. Models consistent with stratigraphy constrain the depth dependence of mantle viscosity and buoyancy, requiring that the Farallon slab was flat lying in the Late Cretaceous, consistent with geological reconstructions. The simulation predicts that an extensive zone of shallow-dipping subduction extended beyond the flat-lying slab farther east and north by up to 1000 kilometers. The limited region of flat subduction is consistent with the notion that subduction of an oceanic plateau caused the slab to flatten. The results imply that seismic images of the current mantle provide more constraints on past tectonic events than previously recognized.

Seismological support for the metastable superplume model, sharp features, and phase changes within the lower mantle.

Recently, a metastable thermal-chemical convection model was proposed to explain the African Superplume. Its bulk tabular shape remains relatively stable while its interior undergoes significant stirring with low-velocity conduits along its edges and down-welling near the middle. Here, we perform a mapping of chemistry and temperature into P and S velocity variations and replace a seismically derived structure with this hybrid model. Synthetic seismogram sections generated for this 2D model are then compared directly with corresponding seismic observations of P (P, P(C)P, and PKP) and S (S, S(C)S, and SKS) phases. These results explain the anticorrelation between the bulk velocity and shear velocity and the sharpness and level of SKS travel time delays. In addition, we present evidence for the existence of a D" triplication (a putative phase change) beneath the down-welling structure.

Sharp sides to the African superplume.

Beneath southern Africa is a large structure about 1200 kilometers across and extending obliquely 1500 kilometers upward from the core-mantle boundary with a shear velocity reduction of about 3%. Using a fortuitous set of SKS phases that travel along its eastern side, we show that the boundary of the anomaly appears to be sharp, with a width less than 50 kilometers, and is tilted outward from its center. Dynamic models that fit the seismic constraints have a dense chemical layer within an upwardly flowing thermal structure. The tilt suggests that the layer is dynamically unstable on geological time scales.