Mantle Flow Induced by the Interplay of Downgoing Slabs Revealed by Seismic Anisotropy Beneath the Sula Block in Eastern Indonesia

Abstract

The North Sulawesi subduction zone is characterized by southward subduction of the Celebes Sea slab to a depth of ∼250 km, mainly overlying the Sangihe slab that subducts west from the Molucca Sea and penetrates the mantle transition zone. The palaeo‐subducted Sula slab dips northward and partially underlies both the Sangihe and Celebes Sea slabs. Adjacent subduction zones with horizontal overlapping subducting slabs in the upper mantle have unclear dynamic interactions. An extensive strike‐slip fault forms the western boundary of the active North Sulawesi subduction zone, providing an ideal setting to study mantle flow between overlapping slabs. We use local S‐wave and teleseismic S and SK(K)S waveform splitting analysis to measure seismic anisotropy in the northern Sulawesi region. Our observations reveal typical mantle wedge corner flow within the Sangihe subduction system. In the Gulf of Tomini, the observed trench‐oblique fast‐axis orientations above the Celebes Sea slab are likely a consequence of the interaction between two subducting slabs. The southernmost measurement with an E–W‐trending fast direction in the mantle wedge might be related to the subduction of the Sula slab. Furthermore, fault‐parallel fast‐axis orientations of anisotropy near the southern segment of the Palu‐Koro fault are attributed to large‐scale shearing across this lithospheric‐scale strike‐slip fault system. Overall, our observations suggest that the strain caused by lithospheric and asthenospheric deformation is mainly confined within the microplate, displaying a restricted flow pattern and localized effects due to the size of the plate boundaries, such as the Palu‐Koro fault.