Numerical modeling of induced subduction initiation: Insights from the oceanic plateau accretion

Abstract

Subduction initiation is a crucial process for plate tectonics. Two scenarios of induced subduction initiation–subduction polarity reversal and subduction transference–are primarily attributed to the accretion of oceanic plateaus. However, the factors governing these two distinct types of subduction initiation remain elusive. Here, we employ 2D thermomechanical numerical models to investigate the dynamics of induced subduction initiation during oceanic plateau accretion. Specifically, we explore the impact of the thickness, density, and rheological strength of the oceanic plateau on subduction initiation. Our results reveal three geodynamic regimes: continuous subduction, subduction transference, and subduction polarity reversal. Thin oceanic plateaus tend to subduct along with slabs into the deep mantle, while thick oceanic plateaus often accrete to the margin, triggering subduction initiation. The rheological strength and density of an oceanic plateau play a significant role in determining whether subduction transference or polarity reversal occurs. Increasing the viscosity and decreasing the density of the oceanic plateau promotes strain concentration in front of the oceanic plateau, resulting in subduction polarity reversal. Our model results provide insights into the dynamics of subduction transference in the Caroline Sea and polarity reversal in Solomon in the western Pacific Ocean.