Shallow earthquakes in a viscoelastic shear zone with depth-dependent friction and rheology

Summary. Most crustal earthquakes of the world are observed to occur within a seismogenic layer which extends from the Earth's surface to a depth of a few tens of kilometres at most. A model is proposed in which the shear zone along a transcurrent plate margin is represented as a viscoelastic medium with depth-dependent power-law rheology. A frictional resistance linearly increasing with depth is assumed on a vertical transcurrent fault within the shear zone. Such a model is able to reproduce a continuous transition from the brittle behaviour of the upper crust to the ductile behaviour at depth. Assuming that the shear zone is subjected to a constant strain rate from the opposite motions of the two adjacent plates, it is found that there exists a maximum depth H below which tectonic stress can never reach the frictional threshold: this may be identified as the maximum depth of earthquake nucleation. The value of H is consistent with observations for plausible values of the model parameters. The stress evolution in the shear zone is calculated in the linear approximation of the constitutive equation. A change in rigidity with depth, which is also introduced in the model, may reproduce the high vertical gradient of shear stress, which has been measured across the San Andreas fault, and the fact that most earthquakes are nucleated at some depth in the seismogenic layer. A crack which drops the ambient stress to the dynamic frictional level is then introduced in the model. To this aim, a crack solution is employed without a stress singularity at its edges, which is compatible with a frictional stress threshold criterion for fracture. A constraint on the vertical friction gradient is obtained if such cracks are assumed to be entirely confined within the seismogenic layer.     

The influence of earthquakes on the Chandler wobble during 1977–1983

Summary. A direct calculation is made of the effect on the Chandler wobble of 1287 earthquakes that occurred during 1977–1983. The hypocentral parameters (location and origin time) and the moment tensor representation of the best point source for each earthquake as determined by the 'centroidmoment tensor' technique were used to calculate the change in the Chandler wobble's excitation function by assuming this change is due solely to the static deformation field generated by that earthquake. The resulting theoretical earthquake excitation function is compared with the 'observed' excitation function that is obtained by deconvolving a Chandler wobble time series derived from LAGEOS polar motion data. Since only 7 years of data are available for analysis it is not possible to resolve the Chandler band and determine whether or not the theoretical earthquake excitation function derived here is coherent and in phase with the 'observed' excitation function in that band. However, since the power spectrum of the earthquake excitation function is about 56 dB less than that of the 'observed' excitation function at frequencies near the Chandler frequency, it is concluded that earthquakes, via their static deformation field, have had a negligible influence on the Chandler wobble during 1977–1983. However, fault creep or any type of aseismic slip that occurs on a time-scale much less than the period of the Chandler wobble could have an important (and still unmodelled) effect on the Chandler wobble.