Influence of the Thickness of the Overriding Plate on Convergence Zone Dynamics

The important role played by the upper plate in convergence zones dynamics has long been underestimated but is now more and more emphasized. However, the influence of its thickness and/or strength on orogenic systems evolution remains largely unknown. Here we present results from 3D thermo-mechanical numerical simulations of convergence zones (including oceanic subduction followed by continental subduction/collision), in which we vary the rheological profile of the overriding plate (OP). For this, we systematically modify the crustal thickness of the overriding lithosphere and the temperature at the Moho to obtain a thermal thickness of the overriding lithosphere ranging from 80 to 180 km. While all models share a common global evolution (i.e., slab sinking, interaction between slab and the 660 km discontinuity, continental subduction/collision, and slab breakoff), they also highlight first-order differences arising from the variations in the OP strength (thermal thickness). With a thin/weak OP, slab rollback is favored, the slab dip is low, the mantle flow above the slab is vigorous, and the trench migrates at a high rate compared to a thick/strong OP. In addition, slab breakoff and back-arc basin formation events occur significantly earlier than in models involving a thick OP. Our models therefore highlight the major role played by the thickness/strength of the OP on convergence zone dynamics and illustrate its influence in a quantitative way.

Horizontal mantle flow controls subduction dynamics

It is generally accepted that subduction is driven by downgoing-plate negative buoyancy. Yet plate age –the main control on buoyancy– exhibits little correlation with most of the present-day subduction velocities and slab dips. “West”-directed subduction zones are on average steeper (~65°) than “East”-directed (~27°). Also, a “westerly”-directed net rotation of the lithosphere relative to the mantle has been detected in the hotspot reference frame. Thus, the existence of an “easterly”-directed horizontal mantle wind could explain this subduction asymmetry, favouring steepening or lifting of slab dip angles. Here we test this hypothesis using high-resolution two-dimensional numerical thermomechanical models of oceanic plate subduction interacting with a mantle flow. Results show that when subduction polarity is opposite to that of the mantle flow, the descending slab dips subvertically and the hinge retreats, thus leading to the development of a back-arc basin. In contrast, concordance between mantle flow and subduction polarity results in shallow dipping subduction, hinge advance and pronounced topography of the overriding plate, regardless of their age-dependent negative buoyancy. Our results are consistent with seismicity data and tomographic images of subduction zones. Thus, our models may explain why subduction asymmetry is a common feature of convergent margins on Earth.

Compositional mantle layering revealed by slab stagnation at ~1000-km depth

Improved constraints on lower-mantle composition are fundamental to understand the accretion, differentiation, and thermochemical evolution of our planet. Cosmochemical arguments indicate that lower-mantle rocks may be enriched in Si relative to upper-mantle pyrolite, whereas seismic tomography images suggest whole-mantle convection and hence appear to imply efficient mantle mixing. This study reconciles cosmochemical and geophysical constraints using the stagnation of some slab segments at ~1000-km depth as the key observation. Through numerical modeling of subduction, we show that lower-mantle enrichment in intrinsically dense basaltic lithologies can render slabs neutrally buoyant in the uppermost lower mantle. Slab stagnation (at depths of ~660 and ~1000 km) and unimpeded slab sinking to great depths can coexist if the basalt fraction is ~8% higher in the lower mantle than in the upper mantle, equivalent to a lower-mantle Mg/Si of ~1.18. Global-scale geodynamic models demonstrate that such a moderate compositional gradient across the mantle can persist can in the presence of whole-mantle convection.

Structural control on the Tohoku earthquake rupture process investigated by 3D FEM, tsunami and geodetic data

The 2011 Tohoku earthquake (M_w = 9.1) highlighted previously unobserved features for megathrust events, such as the large slip in a relatively limited area and the shallow rupture propagation. We use a Finite Element Model (FEM), taking into account the 3D geometrical and structural complexities up to the trench zone, and perform a joint inversion of tsunami and geodetic data to retrieve the earthquake slip distribution. We obtain a close spatial correlation between the main deep slip patch and the local seismic velocity anomalies, and large shallow slip extending also to the North coherently with a seismically observed low-frequency radiation. These observations suggest that the friction controlled the rupture, initially confining the deeper rupture and then driving its propagation up to the trench, where it spreads laterally. These findings are relevant to earthquake and tsunami hazard assessment because they may help to detect regions likely prone to rupture along the megathrust, and to constrain the probability of high slip near the trench. Our estimate of ~40 m slip value around the JFAST (Japan Trench Fast Drilling Project) drilling zone contributes to constrain the dynamic shear stress and friction coefficient of the fault obtained by temperature measurements to ~0.68 MPa and ~0.10, respectively.