Abstract

The surface deformation associated with the 2010 M-w 8.8 Maule earthquake in Chile was recorded in great detail before, during and after the event. The high data quality of the continuous GPS (cGPS) observations has facilitated a number of studies that model the postseismic deformation signal with a combination of relocking, afterslip and viscoelastic relaxation using linear rheology for the upper mantle. Here, we investigate the impact of using linear Maxwell or power-law rheology with a 2D geomechanical-numerical model to better understand the relative importance of the different processes that control the postseismic deformation signal. Our model results reveal that, in particular, the modeled cumulative vertical postseismic deformation pattern in the near field (< 300 km from the trench) is very sensitive to the location of maximum afterslip and choice of rheology. In the model with power-law rheology, the afterslip maximum is located at 20-35 km rather than > 50 km depth as suggested in previous studies. The explanation for this difference is that in the model with power-law rheology the relaxation of coseismically imposed differential stresses occurs mainly in the lower crust. However, even though the model with power-law rheology probably has more potential to explain the vertical postseismic signal in the near field, the uncertainty of the applied temperature field is substantial, and this needs further investigations and improvements.