Open AccessStudy of plate deformation and stress in subduction processes using two‐dimensional numerical models
Author(s) -
Hassani Riad,
Jongmans Denis,
Chéry Jean
Publication year - 1997
Publication title -
journal of geophysical research: solid earth
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.67
H-Index - 298
eISSN - 2156-2202
pISSN - 0148-0227
DOI - 10.1029/97jb01354
Subject(s) - geology , asthenosphere , subduction , lithosphere , seismology , geophysics , trench , compression (physics) , lithospheric flexure , tectonics , materials science , layer (electronics) , composite material
Two‐dimensional finite element modeling is used to model subduction of an oceanic lithospheric plate beneath continental lithosphere. The subduction process is initiated along a preexisting inclined fault and continues until reaching 400 km of total convergence. The lithosphere is assumed to be underlain by an in viscid asthenosphere. Different rheological laws have been considered for the lithosphere, including elasticity and elastoplasticity. The modeling shows that both the stress system in the plates and the surface topography are strongly dependent on two main parameters: the density contrast between lithosphere and asthenosphere (Δρ = ρ L ‐ ρ A ) and the coefficient of friction along the subduction plane. Varying these two parameters allows explanation of the main characteristics of real subduction zones and results in two major regimes manifested by extension or compression in the arc‐back arc system. Extension and back arc rifting corresponds to a positive density contrast and a low coefficient of friction, while negative Δρ values and/or high friction leads to a compressional regime. The coexistence of trench arc compression and back arc tension is only possible for a coefficient of friction lower than 0.1. The results of the numerical experiments agree with those of experimental modeling conducted under similar physical assumptions.