Subduction dynamics and plate coupling along seismically active margin: the role of slab steps

Subduction dynamics and plate coupling along seismically active margin: the role of slab steps

The negative buoyancy of the slab drives subducting plate and trench motions, influencing tectonic stresses at convergent margins. Variations in slab depth along the trench affect plate coupling and seismicity. We investigated these effects using 3-D subduction models, focusing on slab depth variations and their impact on force balance. Results show convergence velocity is primarily controlled by slab mass and subducting plate rheology, while trench movement is sensitive to slab depth heterogeneity. The highest stress occurs around slab steps, aligning with deformation patterns observed along the Sunda margin. This study highlights the crucial role of trench-parallel forces in plate dynamics.

Mantle flow and upper plate deformation due to slab step

Mantle flow and upper plate deformation due to slab step

Research Tags

None specified.

Associated Publication

Publication information not available.

Abstract

The negative buoyancy of the slab primarily controls the subducting plate and trench motions, and tectonic stresses around the convergent margins. Lateral variations in negative buoyancy associated with varying slab depth along strike must affect plate and margin motions, and, most importantly, have an impact on the stress acting across the margin, thereby setting the context for plate coupling, tectonics and present-day seismicity. Here, we investigated these interactions in 3-D subduction numerical models, focusing on along-trench variations in the subduction depth and the resulting perturbations to the force balance. While we focus on the steps in the slab depth, we additionally test the role of subducting plate, i.e., cohesion and viscosity, and upper plate properties, i.e., thickness, viscosity, and cohesion. The results show that the magnitude of convergence velocity only depends on the integrated slab mass and rheology of the subducting plate when the upper plate is thin. Instead, the trench retreat/advance is sensitive to the heterogeneity in the slab depth, and a complex pattern arises atop the slab step, with a characteristic length of ~500 km from the slab depth perturbation. The remaining parameters of the subducting and upper plate mainly affect the magnitude of the trench velocities with minor influence on the pattern. The highest deformation/stress in the upper plate is observed around the slab step due to the rigidity of the plate, causing mutual perturbation between the deep and shallow slab portions, and the lateral flow around the step. These results are compared with the observations along the Sunda margin, where similar slab depth variations are found. The upper plate deformation in the model shows remarkable compatibility with the observed distributions of compression and extension of the Andaman- Sumatra-Java segments. Our study indicates that trench-parallel forces, arising from the natural variations of slab depth, exert a first-order control on the plate coupling and deformation along convergent margins and should not be neglected.

Graphic abstract

Compute Tags

None specified.

Software

Software information not available.

Model Setup

Numerical model setup. The parameter d denotes the depth of the long slab (i.e., 660 km, green plane), and Δ𝑑 represents step length (i.e., the difference between the long and short slab). The long and short slabs extend from 0-2000 and 2000-4000 km in the Y-direction. The trench is located at X = 2000 km.

Numerical model setup. The parameter d denotes the depth of the long slab (i.e., 660 km, green plane), and Δ𝑑 represents step length (i.e., the difference between the long and short slab). The long and short slabs extend from 0-2000 and 2000-4000 km in the Y-direction. The trench is located at X = 2000 km.

Dataset (NCI catalogue):
https://thredds.nci.org.au/thredds/catalog/nm08/MATE/gollapalli-2022-sunda-subduction/catalog.html

Dataset existing identifier:
10.25914/vw5k-q269

Dataset notes: Output data mainly contains .h5 and .xdmf files.

Model files (NCI catalogue):
https://thredds.nci.org.au/thredds/catalog/nm08/MATE/gollapalli-2022-sunda-subduction/catalog.html

Model files existing identifier:
10.5281/zenodo.7022845

Model files notes: Code and inputs for computational model

Source repository:
https://github.com/ModelAtlasofTheEarth/gollapalli-2022-sunda-subduction

Citation

Gollapalli, T. (2024). Subduction dynamics and plate coupling along seismically active margin: the role of slab steps [Data set]. AuScope, National Computational Infrastructure. https://doi.org/vw5k-q269

Licence

CC-BY-4.0

Funders

  • No response