Models of postseismic deformation after megaearthquakes: the role of various rheological and geometrical parameters of the subduction zone
Abstract. The postseismic deformations following subduction megaearthquakes are characterized by a horizontal velocity which, once non-dimensionalized by the coseismic displacement, increases with distance to the trench then presents an almost constant value for distances between 500 and 1500 km. The vertical velocity features a strong narrow peak on the trenchward side of the volcanic arc. Subsidence is observed in the far-field.
In order to understand better the implications of these observations, the influence of the geometry of low viscosity regions in subduction zones on the postseismic deformations is analyzed using a 2-D finite element model with viscoelastic rheologies. The slab dip in the top 80 km Θtop, and deeper Θbottom and the locking depth all have a limited impact on the ratio of horizontal postseismic velocity over coseismic displacement. The smaller Θbottom, the smaller the amplitude of the predicted vertical velocity in the middle-field (200–500 km from the trench). The presence of the slab at asthenospheric depths affects very significantly both the horizontal and vertical velocities. Models with an 80 km thick lithosphere, where the relaxation occurs only in the asthenosphere, are characterized by a trenchward horizontal velocity decreasing very moderately in the middle-field and an uplift maximum on the continental side of the volcanic arc, at odds with the observations. A low viscosity channel (LVCh) over the deep parts of the subduction interface or a low viscosity wedge (LVW) have a considerable impact on the middle-field horizontal and vertical velocities: the trenchward horizontal velocities are very significantly increased while the vertical velocities are characterized by strong uplift over the deep parts of the subduction interface. In the case of a low viscosity wedge, a marked subsidence further away from the trench, on the continent side of the volcanic arc is predicted. While the low-viscosity wedge affects little the far-field horizontal velocities, the LVCh increases them significantly. The thicknesses of the lithosphere and the asthenosphere also have a strong impact on both the middle-field and the far-field velocities. The larger they are, the further from the trench are the maxima of the ratio of the postseismic over coseismic horizontal displacement and of the far-field subsidence. 3-D modeling with a geometry as precise as possible of the various zones with postseismic creep associated with each megaearthquake is necessary to derive more precise conclusions. However, the 2-D modeling results obtained here, compared with postseismic data, point towards lithospheres and asthenospheres surprisingly similar in various areas of the world, with thicknesses around 70 and 200 km respectively and towards the presence of a LVW and/or a LVCh. The systematic description of the role of each parameter presented here will facilitate the choice of the parameters to vary in 3-D models.
O. Trubienko et al.
O. Trubienko et al.
O. Trubienko et al.
9 citations as recorded by crossref.
- Deep Transient Slow Slip Detected by Survey GPS in the Region of Atacama, Chile E. Klein et al. 10.1029/2018GL080613
- Global coseismic deformations, GNSS time series analysis, and earthquake scaling laws L. Métivier et al. 10.1002/2014JB011280
- Constraints on Transient Viscoelastic Rheology of the Asthenosphere From Seasonal Deformation K. Chanard et al. 10.1002/2017GL076451
- Afterslip and viscoelastic relaxation model inferred from the large-scale post-seismic deformation following the 2010Mw8.8 Maule earthquake (Chile) E. Klein et al. 10.1093/gji/ggw086
- Stress‐driven relaxation of heterogeneous upper mantle and time‐dependent afterslip following the 2011 Tohoku earthquake Y. Hu et al. 10.1002/2015JB012508
- Deep postseismic viscoelastic relaxation excited by an intraslab normal fault earthquake in the Chile subduction zone L. Bie et al. 10.1016/j.tecto.2017.07.012
- Transient Deformation in California From Two Decades of GPS Displacements: Implications for a Three‐Dimensional Kinematic Reference Frame E. Klein et al. 10.1029/2018JB017201
- Afterslip following the 2007 Mw 8.4 Bengkulu earthquake in Sumatra loaded the 2010 Mw 7.8 Mentawai tsunami earthquake rupture zone L. Tsang et al. 10.1002/2016JB013432
- The Geodetic Signature of the Earthquake Cycle at Subduction Zones: Model Constraints on the Deep Processes R. Govers et al. 10.1002/2017RG000586