The effect of obliquity on temperature in subduction zones: insights from 3D numerical modeling
Abstract. In subduction zones the geotherm is thought to vary as a function of the subduction rate and the age of the subducting lithosphere. Along a single subduction zone the rate of subduction can strongly vary due to changes in the angle between the trench and the plate convergence vector, namely the subduction obliquity. We currently observe such a configuration all around the Pacific (e.g. Marianna, Chile, Aleutians). Subduction obliquity is also supposed in the geological record of Western and Central Turkey. In order to investigate this effect, we designed and computed simple thermo-kinematic finite element 3D numerical models. We prescribe the trench geometry by means of a simple mathematical function and compute the mantle flow in the mantle wedge only by solving the equation of mass and momentum conservation. We then solve the energy conservation equation until steady-state is reached. We analyse the results (i) in terms of mantle wedge flow with emphasis on the trench-parallel component and (ii) in terms of temperature along the plate interface by means of maps and depths-temperature path at the interface. We show that the effect of the trench curvature on the geotherm is substantial. A small obliquity yields a small but not negligible trench parallel mantle flow leading to differences of 50 °C along strike of the model. With increasing obliquity, the trench parallel component of the velocity consequently increases and the temperature variation can be as important as 200 °C along strike. This can even be larger with varying plate velocity. Finally, we discuss the implication of our simulations for the ubiquitous oblique systems that are observed on Earth, the limitation of our modeling approach and the significance for the geological record with an emphasis on the case study of Western and Central Turkey.