首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
By inversion analysis of the baseline changes and horizontal displacements observed with GPS (Global Positioning System) during 1990–1994, a high-angle reverse fault was detected in the Shikoku-Kinki region, southwest Japan. The active blind fault is characterized by reverse dip-slip (0.7±0.2  m yr−1 within a layer 17–26  km deep) with a length of 208±5  km, a (down-dip) width of 9±2  km, a dip-angle of 51°±2° and a strike direction of 40°±2° (NE). Evidence from the geological investigation of subfaults close to the southwestern portion of the fault, two historical earthquakes ( M L=7.0, 1789 and 6.4, 1955) near the centre of the fault, and an additional inversion analysis of the baseline changes recorded by the nationwide permanent GPS array from 18 January to 31 December 1995 partially demonstrates the existence of the fault, and suggests that it might be a reactivation of a pre-existing fault in this region. The fact that hardly any earthquakes ( M L>2.0) occurred at depth on the inferred fault plane suggests that the fault activity was largely aseismic. Based on the parameters of the blind fault estimated in this study, we evaluated stress changes in this region. It is found that shear stress concentrated and increased by up to 2.1 bar yr−1 at a depth of about 20  km around the epicentral area of the 1995 January 17  Kobe earthquake ( M L=7.2, Japan), and that the earthquake hypocentre received a Coulomb failure stress of about 5.6 bar yr−1 during 1990–1994. The results suggest that the 1995  Kobe earthquake could have been induced or triggered by aseismic fault movement.  相似文献   

2.
Summary. Attenuation of earthquake intensities with epicentral distance was studied by analysing the intensity data for 39 earthquakes in the United States. Attenuation of MM intensity ( I ) with distance (Δ km) follows a simple relation of the type log I = log I 0 - m Δ, where I 0 is the intensity at the epicentre and m is a constant. Slope m is found to be inversely proportional to the square of the focal depth. Intensity attenuation pattern in the United States in general can be represented by a unified relation I/I 0= exp [-(0.8999/ h 2+ 0.0014)Δ] where 16km ≤ h ≤ 60km. Intensities were calculated with the help of this equation and a good agreement with the observed intensities were found. A comparative study has also been made between the attenuation relations applicable to India and the United States.  相似文献   

3.
Summary . Vertical component Rayleigh-wave amplitudes from 1461 shallow earthquakes recorded in the distance range 0–150° are analysed to separate the effects of earthquake size, epicentral distance (Δ) and recording station.
The estimated decay of amplitude with distance has the form of a theoretical curve for the decay of Rayleigh waves with distance if the assumption is made that the decay due to dispersion for the data analysed is that of an Airy phase. Writing the decay due to anelastic attenuation as exp (- k Δ), k is estimated to be 0.676/rad over the whole range of distance. If the distance effects are represented by a straight line of the form h log Δ+ constant, h is estimated to be 1.15. The calibration function for computing M s derived from the estimated distance effects is very similar to that of Marshall & Basham.
Station effects on Rayleigh-wave amplitudes though statistically significant are small, and can probably be ignored in the computation of M s.
Comparing the estimated surface-wave magnitudes (earthquake size) obtained in this study with the long and short period body-wave magnitudes ( m LPb and m SPb respectively) obtained by Booth, Marshall & Young for the same earthquake shows that m LPb is about equal to M s over the magnitude range of interest (˜4.0–7.0). The m LPb and Ms relationship shows that the greater the long-period energy radiated by an earthquake the smaller proportionately is the short-period energy.  相似文献   

4.
Summary. A low-velocity zone (LVZ) shadow on a refraction profile may contain coherent energy. If the LVZ contains reflectors, the corresponding pre-critical reflection branches will lie within the geometrical shadow. Such energy is observed in data and can, via inversion or synthetic seismogram modelling, potentially provide information on the location and character of the main structural features (those with significant reflectivity) within the LVZ.  相似文献   

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

6.
Summary. The asymptotic properties of spheroidal mode dispersion at high frequency for fixed phase velocity are related to the intercept times τβ( p ) for P and S waves. If the mode eigenfrequency and the ratio of horizontal to vertical displacement at the surface for the mode are known τα( p ) and τβ( p ) may be separately estimated. If discontinuities exist in the velocity model then 'solotone' effects occur, in frequency at fixed slowness, and in τα( p ), τβ( p ) estimated from the mode dispersion as a function of slowness. The coupling of P and S waves in the spheroidal modes means that the interaction of P waves with upper-mantle discontinuities affects also the estimates of the S wave τβ( p ) values for which the corresponding turning points lie in the lower mantle. The asymptotic formalism also shows that sharp pulses formed by superposition of spheroidal modes correspond to multiple PS reflections.
A study of τα( p ), τβ( p ) estimates derived from spheroidal modes with periods from 45–50s, calculated for model 1066B, shows that even in the presence of strong upper-mantle discontinuities the errors in intercept time are only about one-tenth of a period. The asymptotic properties may there-for provide a useful means of estimating intercept times from modes with a few seconds period as a supplement to travel-time methods.  相似文献   

7.
The earthquake swarm that struck Shadwan Island at the entrance of the Gulf of Suez in 2001 August included 408 events. Almost all of these events (94 per cent) were microearthquakes and only 6 per cent had small measurable magnitudes  (5.0 > M L≥ 3.0)  . Most of the earthquakes were weak and followed each other so closely in time that they could not be identified at more distant stations. The fault plane solutions of the strongest events of the swarm show almost identical focal mechanisms, predominately normal faulting with a significant sinistral strike-slip component for nodal planes trending NW–SE. A comparison with the mechanisms of the 1969 and 1972 events which took place 20 km north of the swarm epicentral region shows similarities in faulting type and orientation of nodal planes. The azimuths of T -axes determined from focal mechanisms in this study are oriented in the NNE–SSW direction. This direction is consistent with the present-day stress field derived from borehole breakouts in the southern Gulf of Suez and the last phase of stress field changes in the Late Pleistocene, as well as with recent GPS results.
The source parameters of the largest  ( M L≥ 3.0)  events of the 2001 August Shadwan swarm have been estimated from the P -wave spectra of the Egyptian National Seismograph Network (ENSN). Averaging of the values obtained at different stations shows relatively similar source parameters, including a fault length of  0.65 ≤ L ≤ 2 km  , a seismic moment of  7.1 × 1012≤ Mo ≤ 3.0 × 1014 N m  and a stress drop of  0.4 ≤Δσ≤ 10  bar.  相似文献   

8.
It has long been known that S waves on seismograms of local explosions are often accompanied by strong low-frequency, low-velocity, Rg surface wave trains, often significantly diminished for earthquakes. We utilize this fact to construct a new formal discriminator between earthquakes and explosions by measuring the S -surface-wave group velocity. The method is based on analysing the velogram; that is, the display of the envelope of ground motion versus group velocity V = R/T , where R is the epicentral distance and T  the traveltime. We examine the distribution of seismic energy in time and space using envelopes of records from the Israel Seismic Network (ISN), from which we compute the velograms and observe differences in the velograms of quarry blasts and earthquakes. The data include 143 seismic events occurring in three areas (Galilee, Dead Sea, and Gilad) monitored by the ISN; the magnitude range is M L = 1.0–2.8 at distances of 15–310 km. From the velograms we measure the group velocity, V m s , within the 1–4 km s1 range at which the velogram reaches its maximum for each available station. The resulting V m s (R) function is closely fitted by the empirical relationship a + b  ln  R , with a and b coefficients varying from event to event. A simple linear function c = b + 0.33a at a threshold C = 0.69 completely separates ( a,b ) pairs for the 67 Galilee events, and, for the 76 remaining events, one earthquake and four explosions are wrongly classified. After data validation and application of the Fisher linear discriminator, adapted to the events from Galilee, only two misidentified events remain for the whole data set.  相似文献   

9.
The 2003 August 21 Fiordland earthquake ( M L7.0, M W7.2) was the largest earthquake to occur in New Zealand for 35 yr and the fifth of M6+ associated with shallow subduction in Fiordland in the last 15 yr. The aftershocks are diffuse and do not distinguish between the two possible main shock fault planes implied by the Harvard CMT solution, one corresponding to subduction interface thrusting and the other corresponding to steeply seaward dipping thrusting. The distinction is important for calculating the induced stress changes on the overlying Alpine Fault which has a history of very large earthquakes, the last possibly in 1717. We have relocated the aftershocks, using data from temporary seismographs in the epicentral region and the double difference technique. We then use the correlation between aftershock hypocentres and regions of positive changes in Coulomb Failure Stress (CFS) due to various candidate main shock fault planes to argue for concentrated slip on the shallow landward dipping subduction interface. Average changes in CFS on the offshore segments of the Alpine Fault are then negative, retarding any future large events. In our models the change in CFS is evaluated on faults of optimal orientation in the regional stress field as determined by inversion of P -wave polarities.  相似文献   

10.
Seismic amplitude tomography for crustal attenuation beneath China   总被引:1,自引:0,他引:1  
Amplitude tomography reconstructs seismic attenuation directly from recorded wave amplitudes. We have applied the tomography to amplitude data reported in the 'Annual Bulletin of Chinese Earthquakes' and interpreted the regionally varying crustal attenuation in terms of tectonics. The seismic amplitudes were originally recorded for determining the M L and M S magnitudes. They generally correspond to the maximum amplitudes of the horizontal components of the short-period S waves and intermediate-period Rayleigh waves. Both sets of measurements are sensitive to crustal structure. The peak amplitudes from M L amplitudes spread spherically with significant dispersion and scattering. M S amplitudes show cylindrical spreading with little dispersion. Average crustal Q values for attenuation at 1 Hz are 737 and 505 for M L and M S, respectively, with substantial regional variations. Frequency dependence in the attenuation is also indicated. Regions with the lowest attenuation (high Q values) are beneath the south China Block, Sichuan Basin, Ordos Platform, the Daxinganling and the Korea Craton. These tend to be tectonically inactive regions, which are generally dominated by intrusive and cratonic rocks in the upper crust. Regions with the highest attenuation (low Q values) are beneath Bohai Basin, Yunnan, eastern Songpan-Ganzi Terrain, margins of the Ordos platform and the Qilian Shan. These are predominantly active basins, grabens and fold belts. The continental margin also highly attenuates both S and surface waves.  相似文献   

11.
The ratio of shear to compressional travel times from a local earthquake gives a good estimate of the average velocity ratio for the wave path, providing the origin time has been moderately well determined. Thus changes of velocity ratio can be identified and mapped from standard observatory data. The method is applied to the principal New Zealand earthquakes since 1964: Gisborne, 1966, M L= 6. 2; Seddon, 1966, M L= 6.0; Inangahua, 1968, M L= 7.1. The results suggest that monitoring could be achieved with a station spacing of 100 km. It appears that velocity change is essentially a rapid process and that the return to normality may not begin until the earthquake is imminent.  相似文献   

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

13.
Summary. The computational effectiveness of travel-time inversion methods depends on the parameterization of a 3-D velocity structure. We divide a region of interest into a few layers and represent the perturbation of wave slowness in each layer by a series of Chebyshev polynomials. Then a relatively complex velocity structure can be dcscribed by a small set of parameters that can be accurately evaluated by a linearized inversion of travel-time residuals. This method has been applied to artificial and real data at small epicentral distances and in the teleseismic distance range. The corresponding matrix equations were solved using singular value decomposition. The results suggest that the method combines resolution with computational convenience.  相似文献   

14.
Earthquake populations have recently been shown to have many similarities with critical-point phenomena, with fractal scaling of source sizes (energy or seismic moment) corresponding to the observed Gutenberg–Richter (G–R) frequency–magnitude law holding at low magnitudes. At high magnitudes, the form of the distribution depends on the seismic moment release rate M˙ and the maximum magnitude m max . The G–R law requires a sharp truncation at an absolute maximum magnitude for finite M˙ . In contrast, the gamma distribution has an exponential tail which allows a soft or 'credible' maximum to be determined by negligible contribution to the total seismic moment release. Here we apply both distributions to seismic hazard in the mainland UK and its immediate continental shelf, constrained by a mixture of instrumental, historical and neotectonic data. Tectonic moment release rates for the seismogenic part of the lithosphere are calculated from a flexural-plate model for glacio-isostatic recovery, constrained by vertical deformation rates from tide-gauge and geomorphological data. Earthquake focal mechanisms in the UK show near-vertical strike-slip faulting, with implied directions of maximum compressive stress approximately in the NNW–SSE direction, consistent with the tectonic model. Maximum magnitudes are found to be in the range 6.3–7.5 for the G–R law, or 7.0–8.2 m L for the gamma distribution, which compare with a maximum observed in the time period of interest of 6.1 m L . The upper bounds are conservative estimates, based on 100 per cent seismic release of the observed vertical neotectonic deformation. Glacio-isostatic recovery is predominantly an elastic rather than a seismic process, so the true value of m max is likely to be nearer the lower end of the quoted range.  相似文献   

15.
Summary. In order to separate the scattering effect from the intrinsic attenuation, we need a multiple scattering model for seismic wave propagation in random heterogeneous media. In this paper, we apply radiative transfer theory to seismic wave propagation and formulate in the frequency domain the energy density distribution in space for a point source. We consider the cases of isotropic scattering and strong forward scattering. Some numerical examples are shown. It is seen that the energy density–distance curves have quite different shapes depending on the values of medium seismic albedo B 0s/(ηsa) where ηs is the scattering coefficient and ηa is the absorption coefficient of the medium. For a high albedo ( B > 0.5) medium, the energy–distance curve is of arch shape and the position of the peak is a function of the extinction coefficient of the medium ηesa. Therefore it is possible to separate the scattering effect and the absorption based on the measured energy density distribution curves.  相似文献   

16.
Analyses of relative P - and S -wave amplitudes of 15 selected earthquakes ( M L <2.3) from a seismic swarm, which occurred in May and June 1994 at the Eyjafjallajökull volcano in South Iceland, reveal similar radiation patterns, a thrust-type double-couple with an additional source component. All focal solutions have nearly vertical T -axes and horizontally oriented P -axes, with E-W-oriented nodal planes. The volume increase corresponding to an isotropic source component is estimated to be in the range of 24 m3. The temporal and spatial seismic pattern, small magnitude range, focal mechanisms and depth range of the Eyjafjallaökull earthquakes indicate vertical intrusion of magma into a confined region at the northern flank of the volcano.  相似文献   

17.
Summary. A method is outlined to determine the dynamic behaviour of a phase boundary in the Earth when non-uniform time-varying pressure and temperature conditions are assumed at the Earth's surface. An integral equation describing the phase boundary motion is derived and it is solved under a linearizing assumption. The solution is obtained in the form of a double integral transform. Short and long time-expansions of the solution can be obtained from series expansion and integration of the Laplace transform along a branch cut. The method is illustrated by considering a stepwise change in surface pressure conditions.
For short times, the solution exhibits the same type of time dependence (i.e. the first-order term is in t 1/2) as the one obtained in the one-dimensional case (i.e. uniform pressure perturbation at the Earth's surface).
For long times, it is shown that the time dependence of the phase boundary motion is almost identical to the one derived for the one- dimensional case if the wavenumber k L of the surface excitation is such that κ k 2Lτ≤ 1 (where τ is the relaxation time associated with the one-dimensional phase boundary motion and κ is the thermal diffusivity). If κ k 2Lτ > 1, then the relaxation time for the phase boundary motion in two dimensions is of the order of κ−1 k −2L.
When considering parameters that would be appropriate for a basalt to eclogite phase transition at Moho depth, the latter situation is met only when the load wavelength is smaller than 35 km.  相似文献   

18.
A cause of bias in estimates of distance and station terms used to compute seismic magnitudes ( m b) is demonstrated in a comparative study on the world amplitude–distance curve. This bias results from the presence of station thresholds in the measurement and reporting of seismic amplitudes and is analogous to that known to be present in routinely computed magnitudes. Distance terms corresponding to low amplitudes such as in the core shadow will tend to be positively biased as are station terms for stations with high thresholds.  相似文献   

19.
The relation between p-Δcurves for surface and deep focus sources is investigated in order to construct synthetic body wave seismograms for non-zero focal depths by the quantized ray theory algorithm. The transformation of a surface focus p-Δ curve into a deep focus p-Δ curve is denned in terms of that curve which corresponds to surface focus rays reflected from the depth at which the deep focus is located. By analogy with the geometry of the surface focus formulation, paths of integration to obtain absolute travel-time and velocity-depth curves can be denned in the p-Δ plane. Explicit inversion from deep focus data is possible only when the velocity-depth structure above the depth of focus is known. Through a comparison of short period quantized ray theory synthetic seismograms with similar Cagniard-de Hoop computations, it is shown that quantized ray theory can be used for accurate predictions of body wave amplitude behaviour corresponding to a wide range of focal depths.  相似文献   

20.
Summary. Group velocities for first and second higher mode Rayleigh waves, in the frequency range 0.8–4.8 Hz, generated from a local earthquake of magnitude 3.7 M L in western Scotland, are measured at stations along the 1974 LISPB line. These provide detailed information about the crustal structure west of the line. The data divide the region into seven apparently homogeneous provinces. Averaged higher mode velocity dispersion curves for each province are analysed simultaneously using a linearized inversion technique, yielding regionalized shear velocity profiles down to a depth of 17 km into the upper crust. Shear wave velocity is between 3.0 and 3.4 km s−1 in the upper 2 km, with a slow increase to around 3.8 km s−1. P -wave models computed using these results agree with profiles from the LISPB and LUST refraction experiments.  相似文献   

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

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