首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
We analyzed receiver function of teleseismic events recorded at twelve Indonesian-GEOFON (IA-GE) broadband stations using nonlinear Neighbourhood Algorithm (NA) inversion and H-k stacking methods to estimate crustal thickness, V p /V s ratios and S-wave velocity structure along Sunda-Banda arc transition zone. We observed crustal thickness of 34–37 km in Timor Island, which is consistent with the previous works. The thick crust (> 30 km) is also found beneath Sumba and Flores Islands, which might be related to the arc-continent collision causing the thickened crust. In Timor and Sumba Islands, we observed high V p /V s ratio (> 1.84) with low velocity zone that might be associated with the presence of mafic and ultramafic materials and fluid filled fracture zone. The high V p /V s ratio observed at Sumbawa and Flores volcanic Islands might be an indication of partial melt related to the upwelling of hot asthenosphere material through the subducted slab.  相似文献   

2.
Accurate determination of seismic velocity of the crust is important for understanding regional tectonics and crustal evolution of the Earth. We propose a stepwise joint linearized inversion method using surface wave dispersion, Rayleigh wave ZH ratio (i.e., ellipticity), and receiver function data to better resolve 1D crustal shear wave velocity (v S) structure. Surface wave dispersion and Rayleigh wave ZH ratio data are more sensitive to absolute variations of shear wave speed at depths, but their sensitivity kernels to shear wave speeds are different and complimentary. However, receiver function data are more sensitive to sharp velocity contrast (e.g., due to the existence of crustal interfaces) and v P/v S ratios. The stepwise inversion method takes advantages of the complementary sensitivities of each dataset to better constrain the v S model in the crust. We firstly invert surface wave dispersion and ZH ratio data to obtain a 1D smooth absolute v S model and then incorporate receiver function data in the joint inversion to obtain a finer v S model with better constraints on interface structures. Through synthetic tests, Monte Carlo error analyses, and application to real data, we demonstrate that the proposed joint inversion method can resolve robust crustal v S structures and with little initial model dependency.  相似文献   

3.
We made an attempt to assess the shear wave velocity values V S and, to a lesser extent, the V P values from ambient noise recordings in an array configuration. Five array sites were situated in the close proximity to borehole sites. Shear wave velocity profiles were modeled at these five array sites with the aid of two computational techniques, viz. spatial autocorrelation (SPAC) and H/V ellipticity. Out of these five array sites, velocity estimates could be reliably inferred at three locations. The shear wave velocities estimated by these methods are found to be quite consistent with each other. The computed V S values up to 30 m depth are in the range from 275 to 375 m/s in most of the sites, which implies prevalence of a low velocity zone at some pocket areas. The results were corroborated by evidence of site geology as well as geotechnical information.  相似文献   

4.
It is noticed that few geophysical studies have been carried out to decipher the crustal structure of southwestern part of the Northeast India comprising of Tripura fold belt and Bengal basin as compared to the Shillong plateau and the Brahmaputra basin. This region has a long history of seismicity that is still continuing. We have determined first-order crustal features in terms of Moho depths (H) and average VP/VS ratios (κ) using H-κ stacking technique. The inversion of receiver functions data yields near surface thick sedimentary layer in the Bengal basin, which is nearly absent in the Shillong plateau and Tripura fold belt. Our result suggests that the crust is thicker (38–45 km) in the Tripura fold belt region with higher shear-wave velocity in the lower crust than the Shillong plateau. The distribution of VP/VS ratio indicates heterogeneity throughout the whole region. While low to medium value of Poisson’s ratio (1.69–1.75) indicates the presence of felsic crust in the Shillong plateau of the extended Indian Archean crust. The medium to high values of VP/VS ratio (> 1.780) in the Bengal basin and the Tripura fold belt region represent mafic crust during the formation of the Bengal delta and the Tripura fold belt creation in the Precambrian to the Permian age. The depth of the sediments in the Bengal basin is up to 8 km on its eastern margin, which get shallower toward its northeastern and southeastern margins.  相似文献   

5.
We investigated ground response for Baku (Azerbaijan) from two earthquakes of magnitude M6.3 occurred in Caspian Sea (characterized as a near event) and M7.5 in Shamakhi (characterized as a remote extreme event). S-wave velocity with the average shear wave velocity over the topmost 30 m of soil is obtained by experimental method from the V P values measured for the soils. The downtown part of Baku city is characterized by low VS30 values (< 250 m/s), related to sand, water-saturated sand, gravel-pebble, and limestone with clay. High surface PGA of 240 gal for the M7.5 event and of about 190 gal for the M6.3 event, and hence a high ground motion amplification, is observed in the shoreline area, through downtown, in the north-west, and in the east parts of Baku city with soft clays, loamy sands, gravel, sediments.  相似文献   

6.
The evaluation and design of stone column improvement ground for liquefaction mitigation is a challenging issue for the state of practice. In this paper, a shear wave velocity-based approach is proposed based on the well-defined correlations of liquefaction resistance (CRR)-shear wave velocity (V s)-void ratio (e) of sandy soils, and the values of parameters in this approach are recommended for preliminary design purpose when site specific values are not available. The detailed procedures of pre- and post-improvement liquefaction evaluations and stone column design are given. According to this approach, the required level of ground improvement will be met once the target V s of soil is raised high enough (i.e., no less than the critical velocity) to resist the given earthquake loading according to the CRR-V s relationship, and then this requirement is transferred to the control of target void ratio (i.e., the critical e) according to the V s-e relationship. As this approach relies on the densification of the surrounding soil instead of the whole improved ground and is conservative by nature, specific considerations of the densification mechanism and effect are given, and the effects of drainage and reinforcement of stone columns are also discussed. A case study of a thermal power plant in Indonesia is introduced, where the effectiveness of stone column improved ground was evaluated by the proposed V s-based method and compared with the SPT-based evaluation. This improved ground performed well and experienced no liquefaction during subsequent strong earthquakes.  相似文献   

7.
Pressure effect on the electrical conductivity of San Carlos olivine was investigated by the newly installed electrical conductivity measurement system at China University of Geosciences. Electrical conductivity of San Carlos olivine aggregates was measured up to 12 GPa and 1475 K using the Walker-type multi-anvil apparatus equipped with eight WC cubes as the second-stage anvils. The pressure generation against applied load for the experimental assemblage was examined by phase transition of Bi,quartz, forsterite under different P-T conditions. To check the data validity of this new system, electrical conductivities of the serpentinites and talc samples were measured. The results are consistent with the published data of the same samples. Electrical conductivity(σ) of the San Carlos olivine aggregates and temperature(T) satisfy the Arrhenian formula: σ=σ0exp[.(ΔE+PΔV)/kT].The pre-exponential factor(σ0), activation energy(ΔE) and activation volume(ΔV) yield value of 7.74 S/m, 0.85 eV and 0.94cm3/mol, respectively. Electrical conductivities of the San Carlos olivine aggregates decline with increasing pressure at same temperatures. The negative pressure effect can be interpreted by strain energy model of defect energy together with the lattice deformation. In addition, the electrical conductivity-depth 1-D profile of the upper mantle was constructed based on our results and some assumptions. The calculated profile is concordant with the geophysical observation at the depth of 180–350 km beneath Europe, which indicates that the upper mantle beneath Europe might be dry.  相似文献   

8.
Attenuation characteristics in the New Madrid Seismic Zone (NMSZ) are estimated from 157 local seismograph recordings out of 46 earthquakes of 2.6?≤?M?≤?4.1 with hypocentral distances up to 60 km and focal depths down to 25 km. Digital waveform seismograms were obtained from local earthquakes in the NMSZ recorded by the Center for Earthquake Research and Information (CERI) at the University of Memphis. Using the coda normalization method, we tried to determine Q values and geometrical spreading exponents at 13 center frequencies. The scatter of the data and trade-off between the geometrical spreading and the quality factor did not allow us to simultaneously derive both these parameters from inversion. Assuming 1/R 1.0 as the geometrical spreading function in the NMSZ, the Q P and Q S estimates increase with increasing frequency from 354 and 426 at 4 Hz to 729 and 1091 at 24 Hz, respectively. Fitting a power law equation to the Q estimates, we found the attenuation models for the P waves and S waves in the frequency range of 4 to 24 Hz as Q P?=?(115.80?±?1.36) f (0.495?±?0.129) and Q S?=?(161.34?±?1.73) f (0.613?±?0.067), respectively. We did not consider Q estimates from the coda normalization method for frequencies less than 4 Hz in the regression analysis since the decay of coda amplitude was not observed at most bandpass filtered seismograms for these frequencies. Q S/Q P?>?1, for 4?≤?f?≤?24 Hz as well as strong intrinsic attenuation, suggest that the crust beneath the NMSZ is partially fluid-saturated. Further, high scattering attenuation indicates the presence of a high level of small-scale heterogeneities inside the crust in this region.  相似文献   

9.
Solutions of P-SV equations of motion in a homogeneous transversely isotropic elastic layer contain a factor exp(±ν j z), where z is the vertical coordinate and j?=?1, 2. For computing Rayleigh wave dispersion in a multi-layered half space, ν j is computed at each layer. For a given phase velocity (c), ν j becomes complex depending on the transversely isotropic parameters. When ν j is complex, classical Rayleigh waves do not exist and generalised Rayleigh waves propagate along a path inclined to the interface. We use transversely isotropic parameters as α H , β V , ξ, ? and η and find their limits beyond which ν j becomes complex. It is seen that ν j depends on ? and η, but does not depend on ξ. The complex ν j occurs when ? is small and η is large. For a given c/β V , the region of complex ν j in a ? -η plane increases with the increase of α H /β V . Further, for a given α H /β V , the complex region of ν j increases significantly with the decrease of c/β V . This study is useful to compute dispersion parameters of Rayleigh waves in a layered medium.  相似文献   

10.
For studying the structure of the lithosphere in southern Ukraine, wide-angle seismic studies that recorded the reflected and refracted waves were carried out under the DOBRE-4 project. The field works were conducted in October 2009. Thirteen chemical shot points spaced 35–50 km apart from each other were implemented with a charge weight varying from 600 to 1000 kg. Overall 230 recording stations with an interval of 2.5 km between them were used. The high quality of the obtained data allowed us to model the velocity section along the profile for P- and S-waves. Seismic modeling was carried out by two methods. Initially, trial-and-error ray tracing using the arrival times of the main reflected and refracted P- and S-phases was conducted. Next, the amplitudes of the recorded phases were analyzed by the finite-difference full waveform method. The resulting velocity model demonstrates a fairly homogeneous structure from the middle to lower crust both in the vertical and horizontal directions. A drastically different situation is observed in the upper crust, where the V p velocities decrease upwards along the section from 6.35 km/s at a depth of 15–20 km to 5.9–5.8 km/s on the surface of the crystalline basement; in the Neoproterozoic and Paleozoic deposits, it diminishes from 5.15 to 3.80 km/s, and in the Mesozoic layers, it decreases from 2.70 to 2.30 km/s. The subcrustal V p gradually increases downwards from 6.50 to 6.7–6.8 km/s at the crustal base, which complicates the problem of separating the middle and lower crust. The V p velocities above 6.80 km/s have not been revealed even in the lowermost part of the crust, in contrast to the similar profiles in the East European Platform. The Moho is clearly delineated by the velocity contrast of 1.3–1.7 km/s. The alternating pattern of the changes in the Moho depths corresponding to Moho undulations with a wavelength of about 150 km and the amplitude reaching 8 to 17 km is a peculiarity of the velocity model.  相似文献   

11.
To determine the shear wave velocity structure and predominant period features of T?naztepe in ?zmir, Turkey, where new building sites have been planned, active–passive surface wave methods and single-station microtremor measurements are used, as well as surface acquisition techniques, including the multichannel analysis of surface waves (MASW), refraction microtremor (ReMi), and the spatial autocorrelation method (SPAC), to pinpoint shallow and deep shear wave velocity. For engineering bedrock (V s > 760 m/s) conditions at a depth of 30 m, an average seismic shear wave velocity in the upper 30 m of soil (AVs30) is not only accepted as an important parameter for defining ground behavior during earthquakes, but a primary parameter in the geotechnical analysis for areas to be classified by V s30 according to the National Earthquake Hazards Reduction Program (NEHRP). It is also determined that Z1.0, which represents a depth to V s = 1000 m/s, is used for ground motion prediction and changed from 0 to 54 m. The sediment–engineering bedrock structure for T?naztepe that was obtained shows engineering bedrock no deeper than 30 m. When compared, the depth of engineering bedrock and dominant period map and geology are generally compatible.  相似文献   

12.
  相似文献   

13.
Over the past decade, the H-κ stacking technique of Zhu and Kanamori (J Geophys Res 105:2969–2980, 2000) has become a standard tool to determine the crustal thickness H and the bulk crustal vP/vS ratio κ from teleseismic receiver functions. It is obvious that unfavorable noise conditions as well as a complex 3D velocity structure can severely hamper the interpretation of receiver-function data. However, we observe that ambiguities can even arise from a simple 1D layered velocity structure which raises a high potential for misinterpretations. To analyze the feasibility and basic limitations of the H-κ stacking method, we conduct a series of tests based on synthetic data. The impact of different given elementary parameters, related either to the velocity structure or to the data processing, is evaluated in a series of eight individual tests. We deliberately exclude complications such as 3D structural variations and/or noise to show that even a simple 1D velocity structure, involving, e.g., an additional inter-crustal discontinuity, can have significant consequences for the interpretation of the results. However, our modeling suggests that more complex crustal structures may lead to even less reliable results. Additionally, our tests illustrate that time shifts of the maxima in the H-κ domain due to the superposition and merging of individual phases can lead to significantly overestimated vP/vS ratios. In general, the depth to the Moho (or other discontinuities of interest) is less significantly affected. Our tests indicate the necessity to accurately check delay times derived from the maxima of the H-κ stacks against corresponding phases in the receiver functions. Repeating the stacking with varied weighting factors and filter ranges can help to reduce the ambiguities and to avoid possible misinterpretation.  相似文献   

14.
We consider a transversely isotropic medium with vertical axis of symmetry (VTI). Rayleigh wave displacement components in a homogeneous VTI medium contain exp(±krjz), where z is the vertical coordinate, k is the wave number, and j?=?1, 2; rj may be considered as depth-decay factor. In a VTI medium, rj can be a real or imaginary as in an isotropic medium, or it can be a complex depending on the elastic parameters of the VTI medium; if complex, r1 and r2 are complex conjugates. In a homogeneous VTI half space, Rayleigh wave displacement is significantly modified with a phase shift when rj changes from real to complex with variation of VTI parameters; at the transition, the displacement becomes zero when r1?=?r2. In a liquid layer over a VTI half space, Rayleigh waves are dispersive. Here, also Rayleigh wave displacement significantly modified with a phase shift when rj changes from real to complex with a decrease of period. At very low period, phase velocity of Rayleigh waves becomes less than P-wave velocity in the liquid layer giving rise to Scholte waves (interface waves). The amplitudes of Scholte waves with a VTI half space are found to be significantly larger than those with an isotropic half space.  相似文献   

15.
It is a common opinion that only crustal earthquakes can occur in the Crimea–Black Sea region. Since the existence of deep earthquakes in the Crimea–Black Sea region is extremely important for the construction of a geodynamic model for this region, an attempt is made to verify the validity of this widespread view. To do this, the coordinates of all earthquakes recorded by the stations of the Crimean seismological network are reinterpreted with an algorithm developed by one of the authors. The data published in the seismological catalogs and bulletins of the Crimea–Black Sea region for 1970–2012 are used for the analysis. To refine the coordinates of hypocenters of earthquakes in the Crimea–Black Sea region, in addition to the data from stations of the Crimean seismological network, information from seismic stations located around the Black Sea coast are used. In total, the data from 61 seismic stations were used to determine the hypocenter coordinates. The used earthquake catalogs for 1970–2012 contain information on ~2140 events with magnitudes from–1.5 to 5.5. The bulletins provide information on the arrival times of P- and S-waves at seismic stations for 1988 events recorded by three or more stations. The principal innovation of this study is the use of the original author’s hypocenter determination algorithm, which minimizes the functional of distances between the points (X, Y, H) and (x, y, h) corresponding to the theoretical and observed seismic wave travel times from the earthquake source to the recording stations. The determination of the coordinates of earthquake hypocenters is much more stable in this case than the usual minimization of the residual functional for the arrival time of an earthquake wave at a station (the difference between the theoretical and observed values). Since determination of the hypocenter coordinates can be influenced by the chosen velocity column beneath each station, special attention is focused on collecting information on velocity profiles. To evaluate the influence of the upper mantle on the results of calculating the velocity model, two different low-velocity and high-velocity models are used; the results are compared with each other. Both velocity models are set to a depth of 640 km, which is fundamentally important in determining hypocenters for deep earthquakes. Studies of the Crimea–Black Sea region have revealed more than 70 earthquakes with a source depth of more than 60 km. The adequacy of the obtained depth values is confirmed by the results of comparing the initial experimental data from the bulletins with the theoretical travel-time curves for earthquake sources with depths of 50 and 200 km. The sources of deep earthquakes found in the Crimea–Black Sea region significantly change our understanding of the structure and geotectonics of this region.  相似文献   

16.
We have analyzed the behavior of the F2 layer parameters during nighttime periods of enhanced electron concentration by the results of vertical sounding of the ionosphere carried out with five-minute periodicity in Almaty (76°55′ E, 43°15′ N) in 2001–2012. The results are obtained within the frameworks of the unified concept of different types of ionospheric plasma disturbances manifested as variations in the height and half-thickness of the layer accompanied by an increase and decrease of N m F2 at the moments of maximum compression and expansion of the layer. A good correlation is found between height h Am , which corresponds to the maximum increase, and layer peak height h m F, while h Am is always less than h m F. The difference between h Am and h m F linearly increases with increasing h m F. Whereas the difference is ~38 km for h m F = 280 km, it is ~54 km for h m F = 380 km. Additionally, the correlation is good between the increase in the electron concentration in the layer maximum ΔN m and the maximum enhancement at the fixed height ΔN; the electron concentration enhancement in the layer maximum is about two to three times lower than its maximum enhancement at the fixed height.  相似文献   

17.
In this study, continuous wavelet transform is applied to estimate the frequency-dependent quality factor of shear waves, Q S , in northwestern Iran. The dataset used in this study includes velocigrams of more than 50 events with magnitudes between 4.0 and 6.5, which have occurred in the study area. The CWT-based method shows a high-resolution technique for the estimation of S-wave frequency-dependent attenuation. The quality factor values are determined in the form of a power law as Q S (f)?=?(147?±?16)f 0.71?±?0.02 and (126?±?12)f 0.73?±?0.02 for vertical and horizontal components, respectively, where f is between 0.9 and 12 Hz. Furthermore, in order to verify the reliability of the suggested Q S estimator method, an additional test is performed by using accelerograms of Ahar-Varzaghan dual earthquakes on August 11, 2012, of moment magnitudes 6.4 and 6.3 and their aftershocks. Results indicate that the estimated Q S values from CWT-based method are not very sensitive to the numbers and types of waveforms used (velocity or acceleration).  相似文献   

18.
This paper deals with the interpretation of Bouguer gravity anomalies measured along a 250 km long Suhaitu-Etuokeqi gravity profile located at the transitional zone of the Alxa and Ordos blocks where geophysical characteristics are very complex. The analysis is carried out in terms of the ratio of elevation and Bouguer gravity anomaly, the normalized full gradient of a section of the Bouguer gravity anomaly (G h ) and the crustal density structure reveal that (1) the ratio of highs and lows of elevation and Bouguer gravity anomaly is large between Zhengyiguan fault (F4) and Helandonglu fault (F6), which can be explained due to crustal inhomogeneities related to the uplift of the Qinghai-Tibet block in the northeast; (2) the main active faults correspond to the G h contour strip or cut the local region, and generally show strong deformation characteristics, for example the Bayanwulashan mountain front fault (F1) or the southeast boundary of Alxa block is in accord with the western change belt of G h , a belt about 10 km wide that extends to about 30 km; (3) Yinchuan-Pingluo fault (F8) is the seismogenic structure of the Pingluo M earthquake, and its focal depth is about 15 km; (4) the Moho depth trend and Bouguer gravity anomaly variation indicates that the regional gravity field is strongly correlated with the Moho discontinuity.  相似文献   

19.
A method for determining medium quality factor is developed on the basis of analyzing the attenuation dispersion of the arrived first period P wave. In order to enhance signal to noise ratio, improve the resolution in measurement and reduce systematic error we applied the data resampling technique. The group velocity delay of P wave was derived by using an improved multi-filtering method. Based on a linear viscoelastic relaxation model we deduced the medium quality factor Q m, and associated error with 95% confidence level. Applying the method to the seismic record of the Xiuyan M=5.4 earthquake sequences we obtained the following result: (1) High Q m started to appear from Nov. 9, 1999. The events giving the deduced high Q m value clustered in a region with their epicenter distances being between 32 and 46 km to the Yingkou station. This Q m versus distance observation obviously deviates from the normal trend of Q m linearly increasing with distance. (2) The average Q m before the 29 Dec. 1999 M=5.4 earthquake is 460, while the average Q m between the M=5.4 event and the 12 Jan. 2000 M=5.1 earthquake is 391, and the average Q m after the M=5.1 event is 204.  相似文献   

20.
On the 27 June 2015, at 15:34:03 UTC, a moderate-sized earthquake of M w 5.0 occurred in the Gulf of Aqaba. Using teleseismic P waves, the focal mechanism of the mainshock was investigated by two techniques. The first technique used the polarities of the first P wave onsets, and the second technique was based on the normalized waveform modeling technique. The results showed that the extension stress has a NE orientation with a shallow southward plunge while the compression stress has a NW trend with a nearly shallow westward plunge, obtaining a strike-slip mechanism. This result agrees well with the typical consequence of crustal deformation resulting from the ongoing extensional to shear stress regime in the Gulf of Aqaba (NE-SW extension and NW-SE compression). The grid search method over a range of focal depths indicates an optimum solution at 15 ± 1 km. To identify the causative fault plane, the aftershock hypocenters were relocated using the local waveform data and the double-difference technique. Considering the fault trends, the spatial distribution of relocated aftershocks demarcated a NS-oriented causative fault, in consistence with one of the nodal planes of the focal mechanism solution, emphasizing the dominant stress regime in the region. Following the Brune model, the estimates of source parameters exhibited fault lengths of 0.29 ≤ L ≤ 2.48 km, moment magnitudes of 3.0 ≤ M w ≤ 5.0, and stress drops of 0.14 ≤ Δσ < 1.14 MPa, indicating a source scaling similar to the tectonic earthquakes related to plate boundaries.  相似文献   

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

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