首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 29 毫秒
1.
Summary. Using nine IDA records for the Indonesian earthquake of 1977 August 19, we have formed an optimal linear combination of the records and have measured the frequency and Q of 0 S 0 and 1 S 0. The frequency was measured using the moment ratio method. The attenuation was measured by the minimum width method and by the time-lapse method. The frequency and attenuation were measured simultaneously by varying them to obtain a best fit to the data. A 2000-hr stack, the sum of nine individual records, for 0 S 0 gave a frequency of 0.814664 mHz±4 ppm. The values for the Q of 0 S 0 for the three different methods of measurement were 5600,5833 and 5700, respectively. The error in the estimates of Q -1 is about 5 per cent for the minimum power method. For 1 S 0 a 300-hr stack yielded a frequency of 1.63151 mHz±30 ppm. The values of Q for this mode were 1960, 1800 and 1850, respectively, with an error in Q -1 of about 12 per cent for the minimum power method.  相似文献   

2.
The deformation at the core–mantle boundary produced by the 2004 Sumatra earthquake is investigated by means of a semi-analytic theoretical model of global coseismic and postseismic deformation, predicting a millimetric coseismic perturbation over a large portion of the core–mantle boundary. Spectral features of such deformations are analysed and discussed. The time-dependent postseismic evolution of the elliptical part of the gravity field ( J 2) is also computed for different asthenosphere viscosity models. Our results show that, for asthenospheric viscosities smaller than 1018 Pa s, the postseismic J 2 variation in the next years is expected to leave a detectable signal in geodetic observations.  相似文献   

3.
Summary. The statistical capability of the m b: M s discriminant for the discrimination of earthquake and explosion populations is examined by application of discriminant functions to a group of 83 explosions and 72 earthquakes in Eurasia. Equations are derived for the probability that an event is an earthquake or an explosion. The positive sign of DIS in the decision index equation, DIS i = 34.3383 – 11.9569 mb t + 7.1161 M si , indicates that the event i is an earthquake. Its negative sign indicates that event i is an explosion. The probability of correct classification for an event, P i , is related to its DIS i value, by P i = [1-exp (DIS i )]−1, where a large, positive DIS indicates a high probability that an event is an earthquake and a large, negative DIS indicates a high probability that an event is an explosion. The discrimination line M s = 1.680 m b– 4.825, or m b= 0.595 M s+ 2.872 very successfully separates the explosion population from the earthquake population. The points on this line have an equal chance of being an earthquake or an explosion; moreover, for any event, the distance parallel to the M s-axis from the point representing that event in the m b: M s plane to this line is a measure of the probability for the correct classification of that event.  相似文献   

4.
We explore a practical approach to earthquake early warning in southern California by determining a ground-motion period parameter  τ c   and a high-pass filtered displacement amplitude parameter Pd from the initial 3 s of the P waveforms recorded at the Southern California Seismic Network stations for earthquakes with M > 4.0. At a given site, we estimate the magnitude of an event from  τ c   and the peak ground-motion velocity ( PGV ) from Pd . The incoming three-component signals are recursively converted to ground acceleration, velocity and displacement. The displacements are recursively filtered with a one-way Butterworth high-pass filter with a cut-off frequency of 0.075 Hz, and a P -wave trigger is constantly monitored. When a trigger occurs,  τ c   and Pd are computed. We found the relationship between  τ c   and magnitude ( M ) for southern California, and between Pd and PGV for both southern California and Taiwan. These two relationships can be used to detect the occurrence of a major earthquake and provide onsite warning in the area around the station where onset of strong ground motion is expected within seconds after the arrival of the P wave. When the station density is high, the methods can be applied to multistation data to increase the robustness of onsite early warning and to add the regional warning approach. In an ideal situation, such warnings would be available within 10 s of the origin time of a large earthquake whose subsequent ground motion may last for tens of seconds.  相似文献   

5.
Absolute S-velocity estimation from receiver functions   总被引：2，自引：0，他引：2  
We present a novel method to recover absolute S velocities from receiver functions.
For a homogeneous half-space the S velocity can be calculated from the horizontal slowness and the angle of surface particle motion for an incident P wave. Generally, the calculated S velocity is an apparent half-space value which depends on model inhomogeneity and P -waveform. We therefore, suggest to calculate such apparent half-space S velocities from low-pass filtered (smoothed) receiver functions using a suite of filter-parameters, T . The use of receiver functions neutralize the influence of the P -waveform, and the successive low-pass filterings emphasize the variation of S velocity with depth.
We apply this   V S ,app.( T )  technique to teleseismic data from three stations: FUR, BFO and SUM, situated on thick sediments, bedrock and the Greenland ice cap, respectively. The observed   V S ,app.( T )  curves indicate the absolute S velocities from the near surface to the uppermost mantle beneath each station, clearly revealing the different geological environments. Application of linearized, iterative inversion quantify these observations into   V S ( z )  models, practically independent of the S -velocity starting model. The obtained models show high consistency with independent geoscientific results. These cases provide also a general validation of the   V S ,app.( T )  method.
We propose the computation of   V S ,app.( T )  curves for individual three-component broad-band stations, both for direct indication of the S velocities and for inverse modelling.  相似文献   

6.
We evaluate the stress field in and around the southern Korean Peninsula with focal mechanism solutions, using the data collected from 71 earthquakes ( ML = 1.9–5.2) between 1999 and 2004. For this, the hypocentres were relocated and well-constrained fault plane solutions were obtained from the data set of 1270 clear P -wave polarities and 46 SH / P amplitude ratios. The focal mechanism solutions indicate that the prevailing faulting types in South Korea are strike-slip-dominant-oblique-slip faultings with minor reverse-slip component. The maximum principal stresses (σ1) estimated from fault-slip inversion analysis of the focal mechanism solutions show a similar orientation with E–W trend (269°–275°) and low-angle plunge (10°–25°) for all tectonic provinces in South Korea, consistent with the E–W trending maximum horizontal stress (σHmax) of the Amurian microplate reported from in situ stress measurements and earthquake focal mechanisms. The directions of the intermediate (σ2) and minimum (σ3) principal stresses of the Gyeongsang Basin are, however, about 90 deg off from those of the other tectonic provinces on a common σ2–σ3 plane, suggesting a permutation of σ2 and σ3. Our results incorporated with those from the kinematic studies of the Quaternary faults imply that NNW- to NE-striking faults (dextral strike-slip or oblique-slip with a reverse-slip component) are highly likely to generate earthquakes in South Korea.  相似文献   

7.
Summary. One method to determine the depths of sources of anomalies in the Earth's gravity field is to plot log [(2 n + 1) σ n ] versus n where σ n is the n th term in the amplitude spectrum of the Earth's gravitational potential. This procedure assumes that the amplitude spectrum of the anomalous density variations does not vary with n. Such an assumption may not apply to the Earth.  相似文献   

8.
The blockage of the L g wave by crustal barriers such as continental margins and graben structures has long been recognized as providing a very useful tool for mapping large-scale lateral crustal variations along the propagation path. Numerical simulation of L g -wave propagation in complex anelastic media using the pseudospectral method provides insight into the nature of the propagation process using both snapshots of the wavefield and synthetic seismograms. A variety of 2-D structures have been investigated, including the influence of sediments, crustal thickness and attenuation.
Thick sedimentary basins covering a graben structure can have a major influence, since they remove L g energy by generating P conversion and scattering–the principal mechanisms for strong L g attenuation across a graben. The reduction of the L g energy is reinforced by anelastic attenuation in the sediments as well as the influence of the gradually thinning crustal waveguide associated with an elevated Moho.
The extinction of L g in a sequence of explosions fired across the central graben of the North Sea can be simulated by numerical calculations for the structure derived from refraction experiments.  相似文献   

9.
Summary. Tidal gravity measurements have been made at six sites in Britain with two nulled LaCoste and Romberg Earth tide gravitymeters. The M 2 observations from these and two further sites are compared with calculations of the tidal loading from the seas around the British Isles and the major oceans. Models of the M 2 marine tides are convolved with Green's functions for appropriate radially stratified Earth models. The differences between the M 2 observations and the theoretical calculations are less than 0.6 μ gals and it is shown that these differences contain further information concerning the errors in the marine tide models. The M 2 marine tides on the north-west European continental shelf are reasonably well known and this allows a useful test of the feasibility of using tidal gravity measurements for the inverse ocean tide problem in areas where the ocean tides are less well known. The differential gravity loading signal between pairs of gravity stations is shown to be important for considerations of the uniqueness and accuracy of the inverse problem. M 2 tidal gravity loading maps for the British Isles and Europe have been produced which are of use in making corrections to various geodetic measurements.  相似文献   

10.
Some comments on the descriptions of the polarization states of waves   总被引：4，自引：0，他引：4  
Summary. Procedures are developed for specifying the polarization characteristics of n -dimensional waves, and in particular three-dimensional waves of geophysical interest. We show that when a wave is in a pure state or is totally polarized, all the polarization information can be represented by a single vector u in an n -dimensional unitary space. Simple measures of the degree of polarization of the wave are constructed from the characteristic equation of the spectral matrix S . These measures are functions only of the scalar invariants of S and consequently S need not be diagonalized. If S represents a purely polarized wave, the unitary vector which contains the polarization information about the wave can be obtained directly from S using any 2 n – 1 equations of n 2 possible equations. By multiplying by a phase-factor this unitary vector can be written in the form u = r 1+ i r 2 where r 1 and r 2 are orthogonal vectors in a real space. For an elliptically polarized wave, r 1 and r 2 locate the major and minor axes of the ellipse, and the ellipticity is given by the ratio of their magnitudes. The polarization parameters of ULF magnetic waves at the Earth's surface are computed from one set of five equations ( n = 3) and compared with parameters calculated using established techniques.  相似文献   

11.
The tidal dynamics of the Irish and Celtic Seas   总被引：1，自引：0，他引：1  
Summary. Current meter data collected over periods of more than 14 day from the Irish and Celtic Seas are harmonically analysed and presented in maps of tidal stream information. Making use of the analysed current data, and by constructing time series of frictional and inertial stresses which are also harmonically analysed, harmonic constituents of the surface tidal slopes at current meter stations are obtained. Using these with data collected from offshore tide gauges, and in conjunction with coastal tide data, cotidal maps are drawn with some confidence for M 2, S 2, O 1 and K 1, the M 2 chart resolving the discrepancy which exists between the different charts of the Celtic Sea already produced. Cotidal maps for M 3 and M 4 are also presented.
The mean over a tidal cycle of the energy flux for M 2, S 2 and O 1 is also presented in the form of the total energy flux in these constituents which crosses different sectional lines. A flux of 44 × 106 kW is observed to enter the Celtic Sea from the Atlantic and this is compared with previous estimates. An energy budget is also performed for M 2, including all the effects of astronomical forcing and Earth tides to enable comparison to be made between the true energy inflow and the estimated frictional dissipation. Finally, comparison is made between the mean of the instantaneous energy flux and the sum of the energy fluxes associated with the major harmonics.  相似文献   

12.
Summary. Based on accurately located 23 very shallow earthquakes ( h = 1–14 km) in northern and central Greece by portable networks of seismic stations and by the joint epicentre method, the travel times of the Pn -waves from the foci of these earthquakes to the sites of 54 permanent stations in the Balkan region have been determined. The travel times of Pn -waves in the central and eastern part of the area (eastern Greece, south-eastern Yugoslavia, the Aegean Sea, Bulgaria, southern Romania, western Turkey) fit a straight line very well with the Pn velocity equal to 7.9 ± 0.1 km s-1. On the contrary, the travel times of Pn -waves to stations in the western part of the area (Albania, western Greece) do not fit this curve because the Pn -waves travelling to these stations are delayed by more than 1 s due to the thicker crust under the Dinarides–Hellenides mountain range. Time delays for Pn -waves have been calculated for each permanent station in the Balkan area with respect to the mean travel-time curve of these waves in the central and eastern part of the area. Corrections of the travel times for these delays contribute very much to the improvement of the accuracy in the location of the shallow earthquakes in the Aegean and surrounding area.  相似文献   

13.
Summary. The inverse gravity potential problem consists in the determination of the form and the density of the body by its exterior gravity potential. We describe two similar classes of bodies for which this problem has a unique constructive solution.
(1) The first class contains the cylindrical bodies with finite length, arbitrary form of section and ρ( R , ø, z) =ρ1( z )ρ2( R , ø) density distribution, where z is the cylindrical coordinate; R , ø are the polar coordinates in a section plane. This class is important for prospecting geophysics in that it allows us to determine in a unique and constructive way, the function ρ1( R , ø), the length, form and orientation of the cylinder if we know the function ρ1( z ) and the exterior potential. The classical moment problem of functions is the basis for the solution of this problem.
(2) The analogous problem for the class of the spherical cylinders, or bodies bounded by arbitrary similar sections of two different concentric spheres and the radial lateral surface, appears when bodies of planetary size are studied. (An example of these bodies would be the Moon mascons.) The density distribution of these cylinders is ρ(τ, θ, ø) =ρ1(τ)ρ2(θ, ø) where τ, θ, ø are the spherical coordinates. The function ρ1(θ, ø), length and form of spherical sections can be uniquely determined by exterior potential if we know the function ρ1(τ). We propose a new constructive method for harmonic continuation of the gravity potential into the region containing the perturbing masses for the solution of the problem.  相似文献   

14.
Focal mechanisms of small earthquakes with magnitudes of about 3 in the SE Brazilian shield are calculated using S / P amplitude ratios. Low attenuation ( Q p from 400 to 800) in the shield upper-crustal layers allowed sharp S arrivals to be recorded up to distances of 100 km. Besides P -wave polarities, SH -wave first motions were also used to constrain the nodal-plane orientations. Normal and reverse faulting mechanisms with strike-slip components were found. The inversion of four mechanisms to estimate the stress tensor indicated a strike-slip stress regime with roughly E–W-orientated σ 1 and N–S σ 3. Both the orientations and the shape factor ( φ =0.7) of the inverted stress are in excellent agreement with theoretical predictions for that part of Brazil from the driving-force model of Coblentz & Richardson (1996) . Good agreement with the nature of the stress, as well as its orientation, was also found for the model of Meijer (1995) . Both of these theoretical models include spreading stresses along the continent/ocean lithospheric transition. Because the earthquakes are more than 300 km from the continental shelf they should not be affected by the local flexural forces caused by sediment load in the marginal basins. The agreement between observed and theoretical stresses then confirms the importance of continental spreading forces in modelling intraplate stresses.  相似文献   

15.
The influence of barometric-pressure variations on gravity   总被引：2，自引：0，他引：2  
Summary. The superconducting gravimeter has been used to measure the influence of barometric pressure on gravity in the frequency range 0.1–10 cycles/day. These measurements show that the incoherent barometric fluctuations are the major cause of random fluctuations in local gravity and account for much of the 'noise' on our gravimeter records. A simple model has been constructed which adequately explains the response of gravity to the local pressure fluctuations. These measurements also show a response to the global atmospheric tides at S 1 and S 2 which is much larger than the response to local fluctuations. Although this behaviour is most likely due to the response of the world-wide oceans to the atmospheric tides, no theoretical model has yet been constructed.  相似文献   

16.
The traveltime perturbation equations for the quasi-compressional and the two quasi-shear waves propagating in a factorized anisotropic inhomogeneous (FAI) media are derived. The concept of FAI media simplifies considerably these equations. In the FAI medium, the density normalized elastic parameters a ijkl ( X i ) can be described by the relation a ijkl ( X i) = f 2( x i ) A ijkl, where A ijkl are constants, independent of coordinates x i and f 2( x i) is a continuous smooth function of x i . The types of anisotropy ( A ijkl ) and inhomogeneity [ f ( x i)] are not restricted. The traveltime perturbations of individual seismic body waves ( q P , qS 1 and qS 2) propagating in the FAI medium depend, of course, both on the structural pertubations [δ f 2( x i)] and on the anisotropy perturbations (δ A ijkl ), but both these effects are fully separated. The perturbation equations for the time delay between the two qS -waves propagating in the FAI medium are simplified even more. If the unperturbed (background) medium is isotropic, the perturbation of the time delay does not depend on the structural perturbations (δ f 2( x i) at all. This striking result, valid of course only in the framework of first-order perturbation theory, will simplify considerably the interpretation of the time delay between the two split qS -waves in inhomogeneous anisotropic media. Numerical examples are presented.  相似文献   

17.
We propose the use of ℓ1 regularization in a wavelet basis for the solution of linearized seismic tomography problems   A m = d   , allowing for the possibility of sharp discontinuities superimposed on a smoothly varying background. An iterative method is used to find a sparse solution m that contains no more fine-scale structure than is necessary to fit the data d to within its assigned errors.  相似文献   

18.
We invert surface-wave and geodetic data for the spatio-temporal complexity of slip during the M w =8.1 Chile 1995 event by simulated annealing. This quasi-global inversion method allows for a wide exploration of model space, and retains the non-linearity of the source tomography problem. Complex source spectra are obtained from 5 to 45 mHz from first- and second-orbit fundamental-mode Rayleigh waves using an empirical Green's function cross-correlation technique. Coseismic displacement vectors were measured at 10 GPS sites near Antofagasta. They are part of a French-Chilean experiment which monitors the Northern Chile seismic gap. The spectra, together with the geodetic data, are inverted for the moment distribution on a 2-D dipping fault, under the physical constraints of slip positivity and causality. Marginal a posteriori distributions of the model parameters are obtained from several independently inverted solutions. In general, features of the slip model are well resolved. Data are well fitted by a purely unilateral southward rupture with a nearly uniform velocity around 2.5–3.0 km s−1, and a total duration of 65 s. Several regions of moment release were imaged, one near the hypocentre, a major one 80 km south of it and a minor one 160 km south of it. The major patch of moment release seemed to have propagated to relatively shallow depths near the trench, 100 km SSW of the epicentre. The region of major slip is located updip of the 1987, M w =7.5 earthquake, suggesting a causal relationship. Most of the slip occurred updip of the hypocentre (36 km), but the entire coupled plate interface (20–40 km) ruptured during the Chile 1995 event.  相似文献   

19.
We study the tectonic setting and lithospheric structure of the greater Barents Sea region by investigating its isostatic state and its gravity field. 3-D forward density modelling utilizing available information from seismic data and boreholes shows an apparent shift between the level of observed and modelled gravity anomalies. This difference cannot be solely explained by changes in crustal density. Furthermore, isostatic calculations show that the present crustal thickness of 35–37 km in the Eastern Barents Sea is greater than required to isostatically balance the deep basins of the area (>19 km). To isostatically compensate the missing masses from the thick crust and deep basins and to adequately explain the gravity field, high-density material (3300–3350 kg m−3) in the lithospheric mantle below the Eastern Barents Sea is needed. The distribution of mantle densities shows a regional division between the Western and Eastern Barents and Kara Seas. In addition, a band of high-densities is observed in the lower crust along the transition zone from the Eastern to Western Barents Sea. The distribution of high-density material in the crust and mantle suggests a connection to the Neoproterozoic Timanide orogen and argues against the presence of a Caledonian suture in the Eastern Barents Sea. Furthermore, the results indicate that the basins of the Western Barents Sea are mainly affected by rifting, while the Eastern Barents Sea basins are located on a stable continental platform.  相似文献   

20.
This paper presents a method to invert underside-reflection ( P d P or S d S arrivals) data for lateral depth variations of upper-mantle discontinuities, combining traveltime and amplitude data. The method greatly improves the resolution of small-scale undulations obtained by existing imaging methods and does not suffer from the long-wavelength biases that are likely to be present in currently available models. Existing inversion methods account for the large size of the Fresnel zone of underside reflections, but not for its complexity, arising from the mini-max traveltime nature of PP- and SS -related waves. This neglect results in long-wavelength artefacts from small-scale undulations of the discontinuities, obscuring true long-wavelength depth variations. The inversion method presented in this paper uses a complex-valued sensitivity kernel, derived from the representation of underside reflections through a Kirchhoff integral formulation. The sensitivity kernel accounts for the varying sensitivity of the waveforms to discontinuity structure over the Fresnel zone. The method is applied to a large, synthetic data set. The data set consists of P d P amplitudes and traveltimes. The results show that the new inversion method resolves depth variations on a lateral scale that is smaller than the size of the Fresnel zone of individual underside reflections (but larger than the dominant wavelength), retaining the resolution of large-scale variations. The results presented here suggest that the discontinuity depth variations induced by slab penetration of the 670 discontinuity could be resolved by current broad-band P 670 P data sets.  相似文献   

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

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