Strains and stresses within the earth constitute important precur- sors of earthquake. Therefore, the determination of the static defor- mation of the earth around surface faults is important for any scheme for the prediction of earthquake. Ben-Menahem and Singh (1976a) derived explicit expressions for the deformation of a homogeneous sphere due to an internal dislocation of arbitrary orientation and depth. Singh (1970) formulated the problem of the static deformation of a multilayered model of the earth in terms of layer matrices. This for- mulation is now being widely used to compute the residual strain, stress and tilt fields in half-space earth models. Singh and Garg (1985) applied the method of layer matrices to solve the two-dimensional problem of a long displacement dislocation in a multilayered half space. It is reasonable to consider the postseismic crustal movements to be controlled by a quasi-static process of the relaxation of stress changes produced by the sudden appearance of a fault in a viscoelastic medium. The correspondence principle of linear viscoelasticity has been used to obtain the quasi-static solutions from the solutions of the corresponding elastic problems (Singh and Rosenman 1974, Singh and Singh 1989a).