New insight into the velocity and anisotropy structures of the subduction zone in northern Sumatra

In this study, we conducted seismic tomographic inversions to investigate the velocity and anisotropy structures of northern Sumatra, using 9774 P-wave and 8405 S-wave arrivals from regional earthquakes. Isotropic P-wave velocity, isotropic S-wave velocity, P-wave azimuthal anisotropy, and P-wave radial anisotropy models were generated using eikonal equation-based traveltime tomography methods. The study identified low-velocity zones beneath the Toba and Sinabung volcanoes, potentially indicating the presence of magma reservoirs. Furthermore, low-velocity anomalies above the subduction slab were detected, which were likely caused by the dehydration of the slab and interpreted as channels of upwelling flow. The tomographic results revealed a trench-parallel high-velocity belt in the uppermost mantle, representing the subducting slab of the India-Australian plate. The trench-parallel fast velocity directions in the slab suggested that the subducted oceanic slab retains its frozen-in anisotropy formed at the mid-ocean ridge, or that the anisotropy is induced by the lattice-preferred orientation of the B-type olivine. Negative radial anisotropy in the mantle wedge was observed, reflecting hot upwelling flows and transitions of olivine fabrics in the presence of water due to slab dehydration. The results also indicated a multilevel magma plumbing system beneath the Toba Caldera. In summary, the results of this study provided new insights into the structure and dynamic processes of the northern Sumatra subduction zone.