期刊
EARTH AND PLANETARY SCIENCE LETTERS
卷 272, 期 3-4, 页码 698-708出版社
ELSEVIER
DOI: 10.1016/j.epsl.2008.06.009
关键词
2-D numerical modelling; subduction zone; solid-fluid coupling; stress regime; slab dynamics
We developed a new 2-D numerical approach to study solid-fluid coupling applied to subduction zones. The lithosphere is characterised by an elastic or elastoplastic behaviour and the asthenosphere by a homogeneous isoviscous fluid. The temperature effects are ignored and viscosity and density are constant in time. The solid and the fluid problem are discretised by the finite elements method (FEM). The same solid code used in Hassani et al. [Hassani, R., Jongmans, D., Chery, J., Study of plate deformation and stress in subduction processes using two-dimensional numerical models,J. Geophys. Res. (1997) 102 17951-17965.] has been used to compute the solution of the solid problem. The Stokes problem is solved by a direct solver with a stabilisation procedure. We used a very simple staggered coupling method where the fluid domain is regularly re-meshing. We observed numerical instabilities when the time step is not sufficiently small, especially when strong Coupling between the solid and the fluid occurs. We have tested different configurations where the lithosphere is elastic or elastoplastic and show how the slab geometry, the topography and the stress regime in the plates are affected by the Viscous resistance of the mantle. We observed that the asthenosphere viscosity is a fundamental parameter in the subduction process. For subduction with an extensional regime in the upper plate, we observe a linear decrease of the extensional stress as a function of the asthenospheric viscosity. (c) 2008 Elsevier B.V. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据