Subduction initiation is a critical link in the plate tectonic cycle. Intraoceanic subduction zones can form along transform faults and fracture zones, but how subduction nucleates parallel to mid-ocean ridges, as in e.g., the Neotethys Ocean during the Jurassic, remains a matter of debate. In recent years, extensional detachment faults have been widely documented adjacent to slow-spreading and ultraslowspreading ridges where they cut across the oceanic lithosphere. These structures are extremely weak due to widespread occurrence of serpentine and talc resulting from hydrothermal alteration, and can therefore effectively localize deformation. Here, we show geochemical, tectonic, and paleomagnetic evidence from the Jurassic ophiolites of Albania and Greece for a subduction zone formed in the western Neotethys parallel to a spreading ridge along an oceanic detachment fault. With 2-D numerical modeling exploring the evolution of a detachment-ridge system experiencing compression, we show that serpentinized detachments are always weaker than spreading ridges. We conclude that, owing to their extreme weakness, oceanic detachments can effectively localize deformation under perpendicular far-field forcing, providing ideal conditions to nucleate new subduction zones parallel and close to (or at) spreading ridges. Direct implication of this, is that resumed magmatic activity in the forearc during subduction initiation can yield widespread accretion of suprasubduction zone ophiolites at or close to the paleoridge. Our new model casts the enigmatic origin of regionally extensive ophiolite belts in a novel geodynamic context, and calls for future research on three-dimensional modeling of subduction initiation and how upper plate extension is associated with that.

[1] The way individual faults and rift segments link up is a fundamental aspect of lithosphere extension and continental break-up. Little is known however about the factors that control the selection of the different modes of rift interaction observed in nature. Here we use state-of-the-art large deformation 3D numerical models to examine the controls on the style and geometry of rift linkage between rift segments during extension of crustal brittle-ductile coupled systems. We focus on the effect of viscosity of the lower layer, the offset between the rift basins and the amount of strain weakening on the efficiency of rift linkage and rift propagation and the style of extension. The models predict three main modes of rift interaction: 1) oblique to transform linking graben systems for small to moderate rift offset and low lower layer viscosity, 2) propagating but non linking and overlapping primary grabens for larger offset and intermediate lower layer viscosity, and 3) formation of multiple graben systems with inefficient rift propagation for high lower layer viscosity. The transition between the linking (Mode 1) and non-linking mode (Mode 2) is controlled by the trade-off between the rift offset, the strength of brittle-ductile coupling, and the amount of strain weakening. The mode transition from overlapping non-connecting rift segments (Mode 2) to distributed deformation (Mode 3) is mainly controlled by the viscosity of the lower layer and can be understood from minimum energy dissipation analysis arguments.

Abstract. ASPECT (Advanced Solver for Problems in Earth's ConvecTion) is a massively parallel finite element code originally designed for modeling thermal convection in the mantle with a Newtonian rheology. The code is characterized by modern numerical methods, high-performance parallelism and extensibility. This last characteristic is illustrated in this work: we have extended the use of ASPECT from global thermal convection modeling to upper-mantle-scale applications of subduction.Subduction modeling generally requires the tracking of multiple materials with different properties and with nonlinear viscous and viscoplastic rheologies. To this end, we implemented a frictional plasticity criterion that is combined with a viscous diffusion and dislocation creep rheology. Because ASPECT uses compositional fields to represent different materials, all material parameters are made dependent on a user-specified number of fields.The goal of this paper is primarily to describe and verify our implementations of complex, multi-material rheology by reproducing the results of four well-known two-dimensional benchmarks: the indentor benchmark, the brick experiment, the sandbox experiment and the slab detachment benchmark. Furthermore, we aim to provide hands-on examples for prospective users by demonstrating the use of multi-material viscoplasticity with three-dimensional, thermomechanical models of oceanic subduction, putting ASPECT on the map as a community code for high-resolution, nonlinear rheology subduction modeling.

No abstract

Surface processes and inherited structures are widely regarded as factors that strongly influence the evolution of mountain belts. The first-order effects of these parameters have been studied extensively throughout the last decades, but their relative importance remains notoriously difficult to assess and document. We use lithospheric scale plane-strain thermomechanical model experiments to study the effects of surface processes and extensional inheritance on the internal structure of contractional orogens and their foreland basins. Extensional inheritance is modeled explicitly by forward modeling the formation of a rift basin before reversing the velocity boundary conditions to model its inversion. Surface processes are modeled through the combination of a simple sedimentation algorithm, where all negative topography is filled up to a prescribed reference level, and an elevation-dependent erosion model. Our results show that (1) extensional inheritance facilitates the propagation of basement deformation in the retro-wedge and (2) increases the width of the orogen; (3) sedimentation increases the length scale of both thin-skinned and thick-skinned thrust sheets and (4) results in a wider orogen; (5) erosion helps to localize deformation resulting in a narrower orogen and a less well-developed retro-wedge. A comparison of the modeled behaviors to the High Atlas, the Pyrenees, and the Central Alps, three extensively studied natural examples characterized by different degrees of inversion, is presented and confirms the predicted controls of surface processes and extensional inheritance on orogenic structure.