Displacement, strain and stress fields due to shear and tensile dislocations in a viscoelastic half-space

Okada (1992) provided expressions for the displacement and strain fields due to a finite rectangular source in an elastic, homogeneous and isotropic half-space. Starting with these results, we applied the correspondence principle of linear viscoelasticity to derive the quasi-static displacement, strain and stress fields in a viscoelastic, homogeneous and isotropic half-space. We assume that the medium deforms viscoelastically with respect to both the shear and the normal stresses but keeps a constant bulk modulus; in particular, the shear modulus relaxes as Maxwell fluid. We presented the viscoelastic effect on displacement, displacement gradient and stress fields, for a choice of parameter values. The viscoelastic effect due to the sudden dislocation reaches a limit value after about 10 times the Maxwell time. The expressions obtained here provide tools for the study of viscoelastic relaxation of lithosphere associated with seismic and volcanic phenomena.     

Scaling relations of earthquakes and aseismic deformation in a damage rheology model

We perform analytical and numerical studies of scaling relations of earthquakes and partition of elastic strain energy between seismic and aseismic components using a thermodynamically based continuum damage model. Brittle instabilities occur in the model at critical damage level associated with loss of convexity of the strain energy function. A new procedure is developed for calculating stress drop and plastic strain in regions sustaining brittle instabilities. The formulation connects the damage rheology parameters with dynamic friction of simpler frameworks, and the plastic strain accumulation is governed by a procedure that is equivalent to Drucker–Prager plasticity. The numerical simulations use variable boundary forces proportional to the slip-deficit between the assumed far field plate motion and displacement of the boundary nodes. These boundary conditions account for the evolution of elastic properties and plastic strain in the model region. 3-D simulations of earthquakes in a model with a large strike-slip fault produce scaling relations between the scalar seismic potency, rupture area, and stress drop values that are in good agreement with observations and other theoretical studies. The area and potency of the simulated earthquakes generally follow a linear log–log relation with a slope of 2/3, and are associated with stress drop values between 1 and 10 MPa. A parameter-space study shows that the area-potency scaling is shifted to higher stress drops in simulations with parameters corresponding to lower dynamic friction, more efficient healing, and higher degree of seismic coupling.