Interseismic deformation along subduction zones: studying the 2010 Maule earthquake by means of 2D finite element model
Abstract
We implemented a two-dimensional finite element model that simulates the accumulation of crustal deformation due to the tectonic loading on a locked subduction fault and applied this model to study the seismic cycle of the Mw8.8 2010 Maule (Central Chile) earthquake. Our goal is to gain insigth into the fundamental factors controling elastic strain build-up and release in subduction zones and to evaluate different approaches proposed for modeling surface deformation as observed by GPS-based crustal velocities. By applying the finite element technique we developed a linear elasticity solver that allows us to assess a realistic plate geometry, rheology and relative velocity of subducting plate in a coupled seismic zone. Constraining parameters such as convergence velocity as well as the geometry of the subduction zone are supported by independent geophysical data so we concentrate on the influence of mechanical slab thickness, variations in the updip and downdip limit, rheology of upper ans lower plates and velocity reduction at the base of the subducted plate. We have introduced idealized geometric models, noting that our numerical solution reproduce the analytical solution for an elastic half-space and that the surface displacement field obtained for a curved fault and non-zero slab thickness model mimics the predictions of a simple backslip model when the slab thickness tends to zero. We compared model predictions with GPS observations in a EW profile crossing the Maule earthquake rupture area in an attempt for determining the parameters of the seismogenic zone most suitable for this region. Our results, that consider a realistic geometry and uniform convergence velocity, suggest little influence of the subducting plate thickness for the same downdip limit and the fit to observations is achieved only after reducing the velocity at the base of the subducted slab by 10% with respect to convergence in a localized region below the trench axis. We will present these novel results that expand our comprehension of seismo-tectonic processes at subduction zones and discuss possible causes for this mechanical behaviour.
- Publication:
-
EGU General Assembly Conference Abstracts
- Pub Date:
- April 2013
- Bibcode:
- 2013EGUGA..1512613C