Postseismic deformation due to subcrustal viscoelastic relaxation following dip-slip earthquakes

Abstract

The deformation of the Earth following a dip-slip earthquake is calculated using a three layer rheological model and finite element techniques. The three layers are an elastic upper lithosphere, a standard linear solid lower lithosphere, and a Maxwell viscoelastic asthenosphere-a model previously analyzed in the strike-clip case (Cohen, 1981, 1982). Attention is focused on the magnitude of the postseismic subsidence and the width of the subsidence zone that can develop due to the viscoelastic response to coseismic reverse slip. Detailed analysis for a fault extending from the surface to 15 km with a 45 deg dip reveals that postseismic subsidence is sensitive to the depth to the asthenosphere but is only weakly dependent on lower lithosphere depth. The greatest subsidence occurs when the elastic lithosphere is about 30 km thick and the asthenosphere lies just below this layer (asthenosphere depth = 2 times the fault depth). The extremum in the subsidence pattern occurs at about 5 km from the surface trace of the fault and lies over the slip plane. In a typical case after a time t = 30 tau (tau = Maxwell time) following the earthquake the subsidence at this point is 60% of the coseismic uplift. Unlike the horizontal deformation following a strike slip earthquake, significant vertical deformation due to asthenosphere flow persists for many times tau and the magnitude of the vertical deformation is not necessarily enhanced by having a partially relaxing lower lithosphere