Computational modelling of evolving fault systems based on a multi-mechanism damage rheology

In seismically active regions, faults nucleate, propagate, and form networks that evolve over time. Progressive strain localisation and periodic fault system re-configuration alter the geometrical complexity and internal structure of fault systems. To better simulate long-term crustal deformation including fault nucleation and propagation (i.e. creation of new frictional surfaces) a numerical model should incorporate concepts of fracture mechanics in a numerical framework that enables representation of evolving interfaces. For this we incorporate the damage rheology model developed by Lyakhovsky and Ben-Zion (Lyakhovsky et al. 1997, 2001) in the numerical framework of Escript (Gross et al., 2007). The progress of our work and initial application of our numerical models will be presented and compared with previous results. Incorporation of damage rheology in Escript is also used as an example of the applicability of this numerical framework for a wide range of Earth science models and disciplines