TY - ABST
T1 - A numerical reference model for themomechanical subduction
AU - Quinquis, Matthieu
AU - Chemia, Zurab
AU - Tosi, Nicola
AU - Buiter, Susanne
AU - Dolejš, David
PY - 2010/5/1
Y1 - 2010/5/1
N2 - Building an advanced numerical model of subduction requires choosing
values for various geometrical parameters and material properties, among
others, the initial lithosphere thicknesses, representative lithological
types and their mechanical and thermal properties, rheologies, initial
temperature profiles as well as model boundary conditions. Some of these
can be constrained by observations of present-day subduction zones, such
as lithosphere age and thickness, whereas others have a larger
uncertainty and require critical evaluation of geophysical and
experimental results. To test the model response to systematic
variations in input parameters, numerical studies often start from a
'reference' subduction model. However, the reference model often varies
between different numerical studies, making it difficult to compare
results directly. We aim therefore to define a numerical reference model
for thermomechanical subduction. This reference setup will facilitate
comparisons of a series of numerical models that focus on different
aspects of subduction, such as the effects of elasticity on the stress
distribution, the energetic impact of phase transformations or the
influence of devolatilization reactions. Our reference model represents
ocean-ocean convergence and describes initial geometries and
lithological stratification for a three-layered subducting slab and
overriding plate along with their respective flow laws and chemical
composition. It also includes kinematic and thermal boundary conditions,
and initial temperature distribution. We will show results of the
evolution and dynamics of the subduction reference model using different
numerical codes: a finite element code, SULEC, and two finite difference
codes, YACC and FDcon.
AB - Building an advanced numerical model of subduction requires choosing
values for various geometrical parameters and material properties, among
others, the initial lithosphere thicknesses, representative lithological
types and their mechanical and thermal properties, rheologies, initial
temperature profiles as well as model boundary conditions. Some of these
can be constrained by observations of present-day subduction zones, such
as lithosphere age and thickness, whereas others have a larger
uncertainty and require critical evaluation of geophysical and
experimental results. To test the model response to systematic
variations in input parameters, numerical studies often start from a
'reference' subduction model. However, the reference model often varies
between different numerical studies, making it difficult to compare
results directly. We aim therefore to define a numerical reference model
for thermomechanical subduction. This reference setup will facilitate
comparisons of a series of numerical models that focus on different
aspects of subduction, such as the effects of elasticity on the stress
distribution, the energetic impact of phase transformations or the
influence of devolatilization reactions. Our reference model represents
ocean-ocean convergence and describes initial geometries and
lithological stratification for a three-layered subducting slab and
overriding plate along with their respective flow laws and chemical
composition. It also includes kinematic and thermal boundary conditions,
and initial temperature distribution. We will show results of the
evolution and dynamics of the subduction reference model using different
numerical codes: a finite element code, SULEC, and two finite difference
codes, YACC and FDcon.
M3 - Conference abstract in journal
VL - 12
SP - 9439
JO - EGU General Assembly 2010, held 2-7 May, 2010 in Vienna, Austria, p.9439
JF - EGU General Assembly 2010, held 2-7 May, 2010 in Vienna, Austria, p.9439
ER -