Geodynamics of divergent double subduction: 3‐D numerical modeling of a Cenozoic example in the Molucca Sea region, Indonesia

Abstract

AbstractGeological observations reveal existence of a unique form of plate subduction featuring subduction on both sides of one single oceanic plate, which is termed divergent double subduction (DDS). DDS may play an important role in facilitating tectonic processes like closure of oceanic basins, accretion and amalgamation of magmatic arcs, and growth of continents. However, this type of subduction has been largely a conceptual model and the geodynamics behind DDS are still poorly constrained. The Molucca Sea subduction zone in SE Asia has been considered as a Cenozoic example of DDS based on geophysical and geological data and provides an opportunity for detailed assessment of how DDS occurs. Here we present 3‐D numerical modeling with aims to reproduce the geodynamic processes of DDS. Several factors that may have important influences on the evolution of DDS are evaluated, including the geometry of the subducting plate, the order of subduction initiation on both sides, the far‐field boundary conditions and thickness of the overriding plates, and the negative buoyancy of the subducting plate. Our results reproduce the observed asymmetrical shape of the subducting Molucca Sea plate and the bending of Halmahera and Sangihe arcs and suggest that DDS is possible if effective escape of the slab‐trapped upper mantle overcomes the space problem, otherwise the slab‐trapped mantle may hinder the sustainability of subduction. We therefore conclude that DDS is associated with closure of narrow and short oceanic plate, and large‐scale double subduction is rare in nature probably owing to space problem.