The mechanisms of shallow earthquakes and the monitoring of a comprehensive test ban

The ability of seismological criteria to identify earthquakes from underground explosions depends partly on the orientation of the earthquake source. Well-determined double-couple moment tensor solutions for a large number of earthquakes have been published in the Harvard centroid moment tensor (CMT) and United Slates Geological Survey (USGS) catalogues. Statistical analyses of these catalogues indicate that the distribution of the orientation of earthquake mechanisms is not random. The distribution of the T axes shows significant clustering around the downward vertical, indicating that a larger number of earthquake mechanisms radiate compressional P -wave energy to teleseismic distances from near the maximum of the radiation pattern than is predicted if earthquake sources are randomly oriented double couples. The clustered T axes correspond to compressional dip-slip mechanisms, and it is this type of mechanism which is believed to cause both the m _b: M _s (the ratio of body-wave to surface-wave magnitude) and first-motion criteria to misidentify an earthquake as an explosion.     

Centroid moment tensor inversion of low-frequency seismic spectra using Green's functions for aspherical earth models

We present a new method for centroid moment tensor (CMT) inversion, in which we employ the Green's function computed for aspherical earth models using the Direct Solution Method. We apply this method to CMT inversion of low-frequency seismic spectra for the 1994 Bolivia and 1996 Flores Sea deep earthquakes. The estimated centroid locations agree well with those obtained by multiple-shock analyses using body-wave data. This shows that it is possible to obtain reliable CMT solutions by analyses of low-frequency seismic spectra using accurate Green's functions computed for present 3-D earth models.     

Determination of the polynomial moments of the seismic moment rate density distribution with positivity constraints

We present a new formulation of the inverse problem of determining the temporal and spatial power moments of the seismic moment rate density distribution, in which its positivity is enforced through a set of linear conditions. To test and demonstrate the method, we apply it to artificial data for the great 1994 deep Bolivian earthquake. We use two different kinds of faulting models to generate the artificial data. One is the Haskell-type of faulting model. The other consists of a collection of a few isolated points releasing moment on a fault, as was proposed in recent studies of this earthquake. The positions of 13 teleseismic stations for which P - and SH -wave data are actually available for this earthquake are used. The numerical experiments illustrate the importance of the positivity constraints without which incorrect solutions are obtained. We also show that the Green functions associated with the problem must be approximated with a low approximation error to obtain reliable solutions. This is achieved by using a more uniform approximation than Taylor's series. We also find that it is necessary to use relatively long-period data first to obtain the low- (0th and 1st) degree moments. Using the insight obtained into the size and duration of the process from the first-degree moments, we can decrease the integration region, substitute these low-degree moments into the problem and use higher-frequency data to find the higher-power moments, so as to obtain more reliable estimates of the spatial and temporal source dimensions. At the higher frequencies, it is necessary to divide the region in which we approximate the Green functions into small pieces and approximate the Green functions separately in each piece to achieve a low approximation error. A derivation showing that the mixed spatio-temporal moments of second degree represent the average speeds of the centroids in the corresponding direction is given.     

Rayleigh-wave dispersion function for a transversely isotropic layered spherical earth using the Thomson-Haskell matrix method

We consider the two coupled differential equations of the two radial functions appearing in the displacement components of spheroidal oscillations for a transversely isotropic (TI) medium in spherical coordinates. Elements of the layer matrix have been explicitly written—perhaps for the first time—to extend the use of the Thomson-Haskell matrix method to the derivation of the dispersion function of Rayleigh waves in a transversely isotropic spherical layered earth. Furthermore, an earth-flattening transformation (EFT) is found and effectively used for spheroidal oscillations. The exponential function solutions obtained for each layer give the dispersion function for TI spherical media the same form as that on a flat earth. This has been achieved by assuming that the five elastic parameters involved vary as r ^p and that the density varies as r ^p-2, where p is an arbitrary constant and r is the radial distance. A numerical illustration with p = - 2 shows that, in spite of the inhomogeneity assumed within layers, the results for spherical harmonic degree n , versus time period T , obtained here for the Primary Reference Earth Model (PREM), agree well with those obtained earlier by other authors using numerical integration or variational methods. The results for isotropic media derived here are also in agreement with previous results. The effect of transverse isotropy on phase velocity for the first two modes of Rayleigh waves in the period range 20 to 240 s is calculated and discussed for continental and oceanic models.