Effects of slab geometry and obliquity on the interplate thermal regime associated with the subduction of three-dimensionally curved oceanic plates

We investigated the relationships among slab geometry, obliquity, and the thermal regime associated with the subduction of oceanic plates using a three-dimensional (3D) parallelepiped thermal convection model. Various models with convex and concave slab shapes were constructed in the numerical simulation, and the temperature and mantle flow distributions were calculated. The results revealed that when the slab dip angle increases, or the obliquity of subduction becomes steeper, the interplate temperature decreases remarkably. Cooler (warmer) zones on the plate interface were identified from the modeling where there was a larger (smaller) subduction angle. Consequently, the interplate temperature distribution is partly controlled by the true subduction angle (TSA), which is a function of the slab dip angle and the obliquity of subduction. The rate of change of the interface temperature for the TSA was 10–50 °C (10°< TSA < 20°) at depths ranging from (TSA – 10) × 5 km to 60 + (TSA – 10) × 5 km for a flat slab after a subduction history of 7 Myrs. The along-arc slab curvature affects the variation in TSA. The slab radius also appeared to influence the radius of induced mantle flow.