首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
Seismic sources with observable glut moments of spatial degree two   总被引：1，自引：0，他引：1  
Let ζΛ and r Λ. be the hypocentral position and time of an extended indigenous seismic source. Backus showed that the force moment tensors of the source, Γ( m +1, n )(ζΛ, r Λ), determine and are determined by the motion which the source produces. For small m + n , only the long-period motion is relevant. The glut moment tensor Λ( m,n ) (ζΛ, r Λ.) can be calculated uniquely from γ( m +1, n )(ζΛ r Λ) only if m = 0 or m = 1. The tensor G =Λ(2,0)(ζΛ) gives the spatial variance tensor WΛ of the source, and WΛ. roughly describes the size, shape and orientation of the source region. Therefore the failure of the observed F =Γ(3,0)(ζΛ) to determine G uniquely is of seismological interest. In the present paper we show that F determines G uniquely if we assume the source to be a simple straight line source (SSLS) or an ideal fault in an isotropic medium with isotropic prestress (IFIMIP). We give tests on F which determine whether it can come from a SSLS, from an IFIMIP or from a simple plane surface source (SPSS). If we assume the source to be a SPSS then knowing F and the fault plane determines G to within an unknown scalar multiple of a certain tensor tangent to the fault plane. Moreover F determines the fault plane uniquely unless F can come from a SSLS. If it can, then F determines this virtual source line uniquely, and F permits the fault plane to be any plane containing the virtual source line.  相似文献   

2.
Summary. We construct a catalogue of all possible elementary point sources for static deformation in an elastic solid. the familiar double-couples, CLVD's, centres of compression and dilatation, etc., are all members of the complete catalogue. the sources are classified according to the rank of the seismic moment tensor, and according to the weight (or order) of the irreducible tensor representation of the 3-D rotation group. These sources can be classified as belonging to one of three general classes. the static excitation functions are calculated for an infinite, homogeneous, isotropic medium for all these sources. We show that, except for sources belonging to these three general classes, all other sources — which are numerous for the tensors of high rank — are null static sources. That is, sources that do not produce any static displacement outside of the source region. Due to the presence of null sources, an inversion of the static deformation data is non-unique. the expansion of the equivalent-force tensors and the stress glut tensors (or seismic moment tensors) into a set of the symmetric trace-free source tensors is proposed. the sources corresponding to seismic moment tensors of the second, third and fourth ranks are considered in more detail. We identify the third-rank sources with rotational dislocations or disclinations.  相似文献   

3.
  相似文献   

4.
Summary. Parameters pertaining to the kinematics of a finite source are usually estimated by fitting specific fault models to the data. On the other hand, these parameters, including source location, are also contained in the moment tensors of higher degree. In this paper, the seismic response is represented in terms of 20 source parameters which are related to components of the moment tensors; they are also related to the parameters of fault models, as will be demonstrated for a number of 'classical' models. A linearized inversion for the moment tensor shows that with real data, or with realistic synthetic data, the results are not necessarily physically meaningful, unless constraints are imposed. The constraints are precisely those appearing as a priori assumptions in the conventional methods of source analysis; it is thus possible to investigate the impact of these assumptions. We will discuss in particular the assumption of a general deviatoric point source (not necessarily a double couple) versus that of a plane fault in finite sources. Although at this stage experience with practical performance of the new method is limited, it is suggested that in the appropriate circumstances constrained inversion for the seismic moment tensors offers a viable alternative to estimate kinematic source parameters.  相似文献   

5.
The first-order statistical moment of the seismic moment tensor   总被引：2，自引：0，他引：2  
Summary. If a complex earthquake is assumed to be a set of individual, randomly oriented elementary pure double couple sources, the solution for the seismic moment of the complex event projected on the mean trend of the fault will perforce be comprised of sources of both double couple and compensated linear vector dipole (CLVD) types. We investigate the statistical properties of these two components of seismic sources in terms of the invariants of the seismic moment tensor of a realistic set of synthetic earthquakes. It is very likely that the size of the CLVD component is two to three orders of magnitude smaller than that of the double couple component.  相似文献   

6.
Summary. Linear-programming methods are powerful and efficient tools for objectively analysing seismic focal mechanisms and are applicable to a wide range of problems, including tsunami warning and nuclear explosion identification. The source mechanism is represented as a point in the six-dimensional space of moment-tensor components. Each observed polarity provides an inequality constraint, linear with respect to the moment tensor components, that restricts the solution to a half-space bounded by a hyperplane passing through the origin. The intersection of these half-spaces is the convex set of all acceptable solutions. Using linear programming, a solution consistent with the polarity constraints can be obtained that maximizes or minimizes any desired linear function of the moment tensor components; the dilatation, the thrust-like nature, and the strike-slip-like nature of an event are examples of such functions. The present method can easily be extended to fit observed seismic-wave amplitudes (either signed or absolute) subject to polarity constraints, and to assess the range of mechanisms consistent with a set of measured amplitudes.  相似文献   

7.
We investigate the effect of laterally varying earth structure on centroid moment tensor inversions using fundamental mode mantle waves. Theoretical seismograms are calculated using a full formulation of surface wave ray theory. Calculations are made using a variety of global tomographic earth models. Results are compared with those obtained using the so-called great-circle approximation, which assumes that phase corrections are given in terms of mean phase slowness along the great circle, and which neglects amplitude effects of heterogeneity. Synthetic tests suggest that even source parameters which fit the data very well may have large errors due to incomplete knowledge of lateral heterogeneity. The method is applied to 31 shallow, large earthquakes. For a given earthquake, the focal mechanisms calculated using different earth models and different forward modelling techniques can significantly vary. We provide a range of selected solutions based on the fit to the data, rather than one single solution. Difficulties in constraining the dip-slip components of the seismic moment tensor often produce overestimates of seismic moment, leading to near vertical dip-slip mechanisms. This happens more commonly for earth models not fitting the data well, confirming that more accurate modelling of lateral heterogeneity can help to constrain the dip-slip components of the seismic moment tensor.  相似文献   

8.
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.  相似文献   

9.
Summary. The long wavelength radiation patterns of P - and S -waves are determined for an elastic prestressed earth model. The prestress is treated as a perturbation of an isotropic medium. Non-zero first-order corrections are found to the eigenvalues and eigenvectors of the deviatoric seismic moment tensor. The radiation pattern given by a purely tangential dislocation is equivalent to a double couple contaminated by a dipolar term linear in the perturbation.  相似文献   

10.
Cluster analysis of seismic moment tensor orientations   总被引：1，自引：0，他引：1  
This paper demonstrates that well-known methods of cluster analysis and multivariate data analysis are useful for geodynamic interpretation of seismic moment tensors. To use these methods, moment tensors are expressed as vectors in a 6-D space. These are vectors in a rigorous sense, rather than an arbitrary set of ordered numbers, because a dot product can be defined that is independent of the coordinate system. In this vector space, non-isotropic moment tensors are a 5-D linear subspace and normalized moment tensors are unit vectors, or points on a unit sphere. Distance along a great circle of the unit sphere satisfies reasonable requirements for any measure of the difference between normalized moment tensors. In regions with a few isolated sets of orientations, cluster analysis based on the great circle distance identifies the same groups of earthquakes that a seismologist would. Figures based on principal component analysis and discriminant analysis illustrate orientation clustering better than equal area projections of moment tensor principal axes. In one case where clusters have been claimed to exist, orientations appear to be continuously distributed and no evidence is found for separate populations of moment tensors.  相似文献   

11.
The two-point correlation function of the seismic moment tensor   总被引：2，自引：0，他引：2  
Summary. We use the invariants of the two-point correlation function of the seismic moment to investigate the degree of irregularity of an earthquake fault, i.e. to study the rapidity with which a complex fault changes its direction of orientation. The two-point correlation function is a fourth-order tensor which has three scalar invariants in the isotropic case. Although the accuracy of present-day catalogues of fault plane solutions is rather low for our purpose, nevertheless the invariants of these correlation tensors confirm the generally  相似文献   

12.
Summary. The method proposed by Mendiguren to determine the source parameters from free oscillation data is applied to the 1970 July 31 deep Colombian earthquake. The results indicate a source propagating horizontally for about 150 km along the lithosphere and cutting across its width. The slab behaves as a guide for source propagation. The horizontal propagation velocity is determined as 3.8 km/s. The intensity of the source grew proportionately to the second power of the propagation distance. This rate of source intensity growth may be interpreted either by a fan-shaped fault model or by a cone-shaped volume source. The average slip and stress drop are estimated as 360 cm and 300 bar for the fault model. For the volume source model the transformational shear strain and stress are estimated as 11 × 10−5 and 160 bar. There is no evidence of a double couple radiation preceding the P origin time. It is shown that the isotropic and deviatoric components of the moment tensor cannot be uniquely resolved when only observations of a single mode are available. It is observed that, from a statistical basis, the available 0 Sn data for Colombian shock can be equally well explained by a pure deviatoric source model or by a source model including an isotropic component. Numerical experiments indicate that the inclusion of higher mode data does not change this situation. But, on the other hand, numerical experiments show that the available data and the scheme used for the inversion would not result in a solution including an artificial implosive component if the actual source were pure deviatoric. If the departure from a pure deviatoric source is produced by noise, it has to be non-random, as it could be produced by lateral heterogeneities not included in the inversion scheme.  相似文献   

13.
FIR filter effects and nucleation phases   总被引：1，自引：0，他引：1  
The symmetric impulse response of linear phase Finite Impulse Response (FIR) filters most commonly used in modern seismic recording systems produces precursory signals to impulsive arrivals. These acausal filter-generated artefacts may result in misinterpretations of various onset properties. Prior to any onset interpretation, these effects have to be removed from the seismic record. This can be achieved without loss of bandwidth by post-filtration of the digital seismograms if the filter coefficients and the decimation ratios are known. We have analysed numerous signals from different instruments and sampling rates for precursory phases and found that—in contrast to commonly held beliefs—FIR-filter-related precursory signals are not always easy to recognize visually from their waveform signature. Furthermore, they can exhibit surprisingly similar properties to those reported for nucleation phases, although the majority of nucleation phases reported in the past have been obtained on instruments with a causal response. We demonstrate examples of filter-related precursory signals for events scanning nine orders of moment, from 1010 N m to 1019 N m. Surprisingly, the lower bound of the artefact durations as a function of seismic moment scales close to the cube root of the seismic moment. We interpret this as being caused by the fact that above a certain seismic moment, the attenuated source signal acts as a causal lowpass filter of a smaller bandwidth than the FIR filter. Assuming an ω-2 source model, constant stress drop and an empirical relationship between the maximum artefact duration and the cut-off frequency of the FIR filter, the artefact durations are expected to scale proportional to the 1/2.5 power of the seismic moment, in comparison to 1/3 as proposed for nucleation phases.  相似文献   

14.
Finite difference (FD) simulation of elastic wave propagation is an important tool in geophysical research. As large-scale 3-D simulations are only feasible on supercomputers or clusters, and even then the simulations are limited to long periods compared to the model size, 2-D FD simulations are widespread. Whereas in generally 3-D heterogeneous structures it is not possible to infer the correct amplitude and waveform from 2-D simulations, in 2.5-D heterogeneous structures some inferences are possible. In particular, Vidale & Helmberger developed an approach that simulates 3-D waveforms using 2-D FD experiments only. However, their method requires a special FD source implementation technique that is based on a source definition which is not any longer used in nowadays FD codes. In this paper, we derive a conversion between 2-D and 3-D Green tensors that allows us to simulate 3-D displacement seismograms using 2-D FD simulations and the actual ray path determined in the geometrical optic limit. We give the conversion for a source of a certain seismic moment that is implemented by incrementing the components of the stress tensor.
Therefore, we present a hybrid modelling procedure involving 2-D FD and kinematic ray-tracing techniques. The applicability is demonstrated by numerical experiments of elastic wave propagation for models of different complexity.  相似文献   

15.
Stable inversions for complete moment tensors   总被引：5，自引：0，他引：5  
The seismic moment tensors for certain types of sources, such as volcanic earthquakes and nuclear explosions are expected to contain an isotropic component. Some earlier efforts to calculate the isotropic component of these sources are flawed due to an error in the method of Jost & Herrmann. We corrected the method after Herrmann & Hutchensen and found great improvement in the recovery of non-double-couple moment tensors that include an isotropic component. Tests with synthetic data demonstrate the stability of the corrected linear inversion method, and we recalculate the moment tensor solutions reported in Dreger et al. for Long Valley caldera events and Dreger & Woods for Nevada Test Site nuclear explosions. We confirm the findings of Dreger et al. that the Long Valley volcanic sources contain large statistically significant isotropic components. The nuclear explosions have strikingly anomalous source mechanisms, which contain very large isotropic components, making it evident that these events are not tectonic in origin. This indicates that moment tensor inversions could be an important tool for nuclear monitoring.  相似文献   

16.
Summary. We give a derivation of the displacement potentials and the wave equations which they satisfy. The derivation is similar to one given by Richards but is more general and yields explicit formulas for the source terms. This generality is retained when the moment tensor representation of the source is used. Formulas for the source terms are given in both spherical and cylindrical coordinate systems and are evaluated for the particular case of a point source with second order moment tensor.  相似文献   

17.
Summary New formulae are obtained for the displacement potentials and displacements due to a point source with moment tensor, and for a fault with prescribed slip. These formulae, unlike previous formulae, are invariant, i.e. they are valid in any coordinate system, not just Cartesian coordinates or orthogonal curvilinear coordinates.
Apart from their invariance, these formulae have other advantages: they are exceedingly simple; the expressions for P -motion are nearly the same as those for S -motion; and the separation of far field motion from final static displacement is automatic.  相似文献   

18.
Summary. The response of a stratified elastic medium can be conveniently characterized by the Green's tensor for the medium. For coupled seismic wave propagation ( P—SV or fully anisotropic), the Green's tensor may be constructed directly from two matrices of linearly independent displacement solutions. Rather simple forms for the Green's tensor can be found if each displacement matrix satisfies one of the boundary conditions on the seismic field. This approach relates directly to 'reflection matrix' representations of the seismic field.
For a stratified elastic half space the Green's tensor is used to give a spectral representation for coupled seismic waves. By means of a contour integration a general completeness relation is obtained for the 'body wave' and 'surface wave' parts of the seismic field. This relation is appropriate for SH and P–SV waves in an isotropic medium and also for full anisotropy.  相似文献   

19.
Summary. The flow of the mantle near subducted lithosphere applies a force (dynamic load) to the slab which has a normal component as well as a shear component. In addition, the greater density of the slab compared to the surrounding mantle results in a downward force (negative buoyancy) which also has a component normal to the slab. Because the slab is long compared to its thickness, the stresses in the slab due to the normal components of the loads are accurately estimated from the bending moment. We have used previous estimates of the loads to compute the bending moment in the slab as a function of distance down the slab. Our results indicate that the loads applied to the slab are large compared to its strength; bending moments due to normal loads will not be transmitted over the length of the slab. Instead, the slab will deflect so that the normal loads are balanced locally.  相似文献   

20.
Summary. The use of a moment tensor formulation in source mechanism inversions is becoming more widespread. However, this method does not yield results that correspond to a double-couple source, which is still the preferred source model. This is because the double couple is characterized by extra constraints on the general symmetric moment tensor. A method is presented of imposing these non-linear constraints while retaining the linear nature of the inversion. An application of this result shows (1) that the double-couple model fits the body-wave data as well as more general models if one explicitly considers the fewer number of free parameters in the double couple and (2) that the principal axes of the moment tensor of the unconstrained model do not necessarily correspond to the P and T axes of the double couple.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号