首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
《International Geology Review》2012,54(10):1213-1225
P- and S-wave tomography of the upper mantle beneath the Cape Verde hotspot is determined using arrival-time data measured precisely from three-component seismograms of 106 distant earthquakes recorded by a local seismic network. Our results show a prominent low-velocity anomaly imaged as a continuous column <100 km wide from the uppermost mantle down to about 500 km beneath Cape Verde, especially below the Fogo active volcano, which erupted in 1995. The low-velocity anomaly may reflect a hot mantle plume feeding the Cape Verde hotspot.  相似文献   

2.
Long-period teleseismic P waves recorded at AAE (Addis Ababa, Ethiopia) and NAI (Nairobi, Kenya) show comparable Ps conversions on the radial component of ground motion. The timing and amplitude of the Ps conversions are modeled with synthetic seismograms to get an estimate of crustal thickness under both stations. Instrument response and effective source time function are removed from the data using a P-wave equalization procedure. The timing of the Ps conversions, relative to direct P, suggests that both stations have similar crustal thickness. Using constraints on crustal velocities determined by previous surface wave dispersion and travel-time studies, the Ps-P timing suggests a crust of 41 km thickness. This agrees reasonably well with previous crustal estimates. Tangential wave forms exhibit large amplitudes and are consistent at most backazimuths. However, these wave forms could not be explained with models containing simple planar dipping interfaces. The crustal thickness of 41 km taken in conjunction with the close proximity of the stations to the rift zone suggests that crustal thinning is localized to the rift itself.  相似文献   

3.
The Sinai Peninsula has been recognized as a subplate of the African Plate located at the triple junction of the Gulf of Suez rift, the Dead Sea Transform fault, and the Red Sea rift. The upper and lower crustal structures of this tectonically active, rapidly developing region are yet poorly understood because of many limitations. For this reason, a set of P- and S-wave travel times recorded at 14 seismic stations belonging to the Egyptian National Seismographic Network (ENSN) from 111 local and regional events are analyzed to investigate the crustal structures and the locations of the seismogenic zones beneath central and southern Sinai. Because the velocity model used for routine earthquake location by ENSN is one-dimensional, the travel-time residuals will show lateral heterogeneity of the velocity structures and unmodeled vertical structures. Seismic activity is strong along the eastern and southern borders of the study area but low to moderate along the northern boundary and the Gulf of Suez to the west. The crustal Vp/Vs ratio is 1.74 from shallow (depth ≤ 10 km) earthquakes and 1.76 from deeper (depth > 10 km) crustal events. The majority of the regional and local travel-time residuals are positive relative to the Preliminary Reference Earth Model (PREM), implying that the seismic stations are located above widely distributed, tectonically-induced low-velocity zones. These low-velocity zones are mostly related to the local crustal faults affecting the sedimentary section and the basement complex as well as the rifting processes prevailing in the northern Red Sea region and the ascending of hot mantle materials along crustal fractures. The delineation of these low-velocity zones and the locations of big crustal earthquakes enable the identification of areas prone to intense seismotectonic activities, which should be excluded from major future development projects and large constructions in central and southern Sinai.  相似文献   

4.
This study presents the crustal shear wave velocity structure and radial anisotropy along two linear seismic arrays across the North China Craton (NCC) from ambient noise tomography. About a half to one year long ambient noise data from 87 stations were used for obtaining the inter-station surface wave empirical Green's functions (EGFs) from cross-correlation. Rayleigh and Love dispersion curves within the period band 5–30 s were measured from the EGFs of the vertical and transverse components, respectively. These dispersion data were then used to determine the crustal shear wave velocity structure (VSV and VSH) and radial anisotropy (2(VSH ? VSV) / (VSH + VSV)) from point-wise linear inversion with constraints from receiver function analysis. Our results reveal substantial structural variations among different parts of the NCC. The Bohai Bay Basin in the eastern NCC is underlain by a thin crust (~ 30 km) with relatively low velocities (particularly VSV) and large positive radial anisotropy in the middle to lower crust. Such a crustal structure is no longer of a cratonic type and may have resulted from the widespread tectonic extension and intensive magmatism in this region since late Mesozoic. Beneath the Ordos Basin in the western NCC, the crust is relatively thicker (≥ 40 km) and well stratified, and presents a large-scale low velocity zone in the middle to lower crust and overall weak radial anisotropy except for a localized lower crust anomaly. The overall structural features of this region resemble those of typical Precambrian shields, in agreement with the long-term stability of the region. The crustal structure under the Trans North China Orogen (TNCO, central NCC) is more complicated and characterized by smaller scale velocity variations, strong positive radial anisotropy in the middle crust and rapid change to weak-to-negative anisotropy in the lower crust. These features may reflect complex deformations and crust–mantle interactions, probably associated with tectonic extension and magmatic underplating during the Mesozoic to Cenozoic evolution of the region. Our structural images in combination with previous seismic, geological and geochemical observations suggest that the Phanerozoic lithospheric reactivation and destruction processes may have affected the crust (especially the middle and lower crust) of the eastern NCC, and the effect probably extended to the TNCO, but may have minor influence on the crust of the western part of the craton.  相似文献   

5.
E.A. Hetland  F.T. Wu  J.L Song   《Tectonophysics》2004,386(3-4):157-175
During 1998–1999, we installed a temporary broadband seismic network in the Changbaishan volcanic region, NE China. We estimated crustal structure using teleseismic seismograms collected at the network. We detected a near surface region of strong anisotropy directly under the main volcanic edifice of the volcanic area. We modeled 109 receiver functions from 19 broadband stations using three techniques. First we used a “slant-stacking” method to model the principal crustal P reverberation phases to estimate crustal thickness and the average crustal P to S speed ratio (vp/vs), assuming an average P-wave velocity in the crust. We then estimated crustal S-wave velocity (vs) and vp/vs profiles by modeling stacked receiver functions using a direct search. Finally, we inverted several receiver functions recorded at stations closest to the main volcanic edifice using least squares to estimate vs velocity profiles, assuming a vp/vs value. The results from the three estimation techniques were consistent, and generally we found that the receiver functions constrained estimates of changes in wave speeds better than absolute values. We resolved that the crust is 30–39 km thick under the volcanic region and 28–32 km thick away from the volcanic region, with a midcrust velocity transition at about 10–15 km depth. We estimated that the average crust P-wave velocity is about 6.0–6.2 km/s surrounding the main volcanic region, while it is slightly lower in the vicinity of the main volcanic edifice. The estimates of vp/vs were more ambiguous, but we inferred that the bulk crustal Poisson's ratio (which is related to vp/vs) ranges between 0.20 and 0.30, with a suggestion that the Poisson's ratio is lower under the central volcanic region compared to the surrounding areas. We resolved low S-wave velocities (down to about 3 km/s) in the middle crust in the region of the main volcanic edifice. The low velocity anomaly extends from about 5–10 to 15–25 km below the surface, probably indicating a region of elevated temperatures. We were unable to determine if partial melt is present with the data we considered in this paper.  相似文献   

6.
The elastic properties of cordierite, a common volatile-bearing metamorphic mineral, were measured using Brillouin spectroscopy under ambient conditions. We obtain a bulk modulus of KS =129(1) GPa, and a shear modulus of G=54.0(4) GPa. The bulk modulus of cordierite is much larger than those of other crustal framework silicates (e.g., quartz and feldspars), but is similar to KS for denser upper mantle phases such as olivine. This is likely a result of the cordierite crystal structure, as suggested by a similarly high value of KS for minerals with closely related structures. Cordierite has an unusually high K/G ratio of about 2.4, and a Poisson’s ratio of 0.31,which may be a diagnostic seismic properties of areas in which cordierite-rich metamorphic rocks occur. The overall velocity anisotropy of cordierite is relatively low (<14%) in comparison with many other metamorphic minerals. Calculated velocities for a representative lower crustal rock suggest that cordierite is not likely to explain the high seismic anisotropy observed in some lower crustal sections. Cordierite would have a strong influence on the bulk seismic anisotropy only in rocks where it is present in large concentrations and has a strong preferred orientation. Although such rocks are known to occur, they are uncommon. Received: 23 Deceber 1997/ Revised, accepted: 12 October 1998  相似文献   

7.
This study is based on the seismic data collected as a result of explosions carried out during the 1976 and 1978 Deep Seismic Sounding (DSS) field operations in the Koyna region. These shots were exploded from twelve shot points by the National Geophysical Research Institute along the Guhagar-Chorochi and Kelsi-Loni profiles.Refraction studies of the records reveal a two-layered crust. The top layer consists 17 km of granite and the second layer 19 km of basalt, giving the average depth of the Moho as 36 km in the region. The velocities of the phases Pg, P* and Pn have been computed as 5.82 ± 0.01, 6.61 ± 0.05 and 8.23 ± 0.05 km/sec respectively and those of Sg, S* and Sn as 3.41 ± 0.00, 4.09 ± 0.07 and 4.60 ± 0.08 km/sec respectively. The shear wave velocity in the basement rock has been found to be lower in comparison with other regions of the peninsular India.In some cases reflections were recorded both from the Moho as well as from the intermediate layer. These reveal a crustal thickness of 39 km with 19 km of granitic and 20 km of basaltic layers.Coda signal durations from DSS explosions recorded by microearthquake seismographs indicate a lateral heterogeneity in the crust on either side of Karad in an east-west direction.  相似文献   

8.
Group velocity dispersion data of fundamental-mode Rayleigh and Love waves for 12 wave paths within southeastern China have been measured by applying the multiple-filter technique to the properly rotated three-component digital seismograms from two Seismic Research Observatory stations, TATO and CHTO. The generalized surface wave inversion technique was applied to these group velocity dispersion data to determine the S-wave velocity structures of the crust and upper mantle for various regions of southeastern China. The results clearly demonstrate that the crust and upper mantle under southeastern China are laterally heterogeneous. The southern China region south of 25°N and the eastern China region both have a crustal thickness of 30 km. The eastern Tibet plateau along the 100°E meridian has a crustal thickness of 60 km. Central China, consisting mainly of the Yangtze and Sino-Korean platforms, has a crustal thickness of 40 km. A distinct S-wave low-velocity layer at 10–20 km depth in the middle crust was found under wave paths in southeastern China. On the other hand, no such crustal low-velocity layer is evident under the eastern Tibet plateau. This low-velocity layer in the middle crust appears to reflect the presence of a sialic low-velocity layer perhaps consisting of intruded granitic laccoliths, or possibly the remnant of the source zone of widespread magmatic activities known to have taken place in these regions since the late Carboniferous.  相似文献   

9.
One in-line wide-angle seismic profile was conducted in 1990 in the course of the Southeastern China Continental Dynamics project aimed at the study of the contact between the Cathaysia block and the Yangtze block. This 380-km-long profile extended in NW–SE direction from Tunxi, Anhui Province, to Wenzhou, Zhejiang Province. Five in-line shots were fired and recorded at seismic stations with spacing of about 3 km along the recording line. We have used two-dimensional ray tracing to model P- and S-wave arrivals and provide constraints on the velocity structure of the upper crust, middle crust, lower crust, Moho discontinuity, and the top part of the lithospheric mantle. P-wave velocity, S-wave velocity and VP/VS ratio are mapped. The crust is 36-km thick on average, albeit it gradually thins from the northwest end to the southeast end (offshore) of the profile. The average crustal velocity is 6.26 km/s for P-waves but 3.6 km/s for S-waves. A relatively narrow low-velocity layer of about 4 km of thickness, with P- and S-wave velocities of 6.2 km/s and 3.5 km/s, respectively, marks the bottom of the middle crust at a depth of 23-km northwest and 17-km southeast. At the crust–mantle transition, the P- and S-wave velocity change quickly from 7.4 to 7.8 km/s (northwest) and 8.0 to 8.2 km/s (southeast) and from 3.9 to 4.2 km/s (northwest) and 3.9 to 4.5 km/s (southeast), respectively. This result implies a lateral contrast in the upper mantle velocity along the 140 km sampled by the profile approximately. The average VP/VS ratio ranges from 1.68–1.8 for the upper crust to 1.75 for the middle and 1.75–1.85 for lower crust. With the interpretation of the wide-angle seismic data, Jiangshan–Shaoxin fault is considered as the boundary between the Yangtze and the Cathaysia block.  相似文献   

10.
Based upon the deep seismic sounding profiles carried out in the Tengchong Volcano-Geothermal Area (TVGA), western Yunnan Province of China, a 2-D crustal P velocity structure is obtained by use of finite-difference inversion and forward travel-time fitting method. The crustal model shows that a low-velocity anomaly zone exists in the upper crust, which is related to geothermal activity. Two faults, the Longling–Ruili Fault and Tengchong Fault, on the profile extend from surface to the lower crust and the Tengchong Fault likely penetrates the Moho. Moreover, based on teleseismic receiver functions on a temporary seismic network, S-wave velocity structures beneath the geothermal field show low S-wave velocity in the upper crust. From results of geophysical survey, the crust of TVGA is characterized by low P-wave and S-wave velocities, low resistivity, high heat-flow value and low Q. The upper mantle P-wave velocity is also low. This suggests presence of magma in the crust derived from the upper mantle. The low-velocity anomaly in upper crust may be related to the magma differentiation. The Tengchong volcanic area is located on the northeast edge of the Indian–Eurasian plate collision zone, away from the eastern boundary of the Indian plate by about 450 km. Based on the results of this paper and related studies, the Tengchong volcanoes can be classified as plate boundary volcanoes.  相似文献   

11.
Seismic anisotropy of the crystalline crust: what does it tell us?   总被引:2,自引:0,他引:2  
The study of the directional dependence of seismic velocities (seismic anisotropy) promises more refined insight into mineral composition and physical properties of the crystalline crust than conventional deep seismic refraction or reflection profiles providing average values of P-and S-wave velocities. The alignment of specific minerals by ductile rock deformation, for instance, causes specific types of seismic anisotropy which can be identified by appropriate field measurements.
Vice versa , the determination of anisotropy can help to discriminate between different rock candidates in the deep crust. Seismic field measurements at the Continental Deep Drilling Site (KTB, S Germany) are shown as an example that anisotropy has to be considered in crustal studies. At the KTB, the dependence of seismic velocity on the direction of wave propagation in situ was found to be compatible with the texture, composition and fracture density of drilled crustal rocks.  相似文献   

12.
The main aim of the TOR project is to study the lithospheric–asthenospheric boundary structure under the Sorgenfrei–Tornquist Zone, across northern Germany, Denmark and southern Sweden. Relative arrival-time residuals of teleseismic P and S phases from 51 earthquakes, recorded by 150 seismic stations along the TOR array, were used to delineate the transition zone in the studied area. The effects of crustal structures were investigated by correcting the teleseismic residuals for travel-time variations in the crust based on a 3D crustal model derived from other data. The inversion was carried out for S phases. The results were then compared with the corresponding P-wave models. As expected, the derived models show that the relatively old and cold Baltic Shield has higher velocity at depth than the younger lithosphere farther South. The models show two sharp and distinct increases in depth to velocities which are low compared to our reference model, as we move from South to North. The location and sharpness of these boundaries suggests that the features resolved are, at least partially, compositional in origin, presumably related to mantle depletion. A sharp and steep subcrustal boundary is found roughly coincident with the southern edge of Sweden. This is below where the edge of the Baltic Shield is usually placed, based on surface geological evidence (the Sorgenfrei–Tornquist Zone). Another less significant transition is recognised more or less beneath the Elbe-lineament. Relatively high d(Vp / Vs) ratios under the central part of the profile (Denmark) indicate relatively low S-velocity in an area where a gravity high supports the hypothesis of extensive mafic intrusions.  相似文献   

13.
14.
A significant drop in seismic travel times of up to 1.0% occurred in the Charlevoix region between 1979 and 1980, possibly related to the M = 5.0 (Aug., 1979) earthquake in the vicinity. A travel-time drop of this magnitude could have been produced either by the closing of dry or saturated cracks in the upper crustal material or by the saturation of dry or partly saturated cracks. However, the anisotropy of travel-time changes in this area supports the view that this travel-time drop was caused by the closing of water-saturated aligned vertical cracks in the crustal material. Three different crack directions with respect to north were resolved: 0 ° or 90 ° in the Precambrian rocks underneath the St. Lawrence River, −18° or 72° in the shallow rocks (< 5 km) of the Charlevoix crater, region, and +35° or 125° in the Paleozoic cover rocks. Crack closure would require a decrease in the pore volume of the rocks which would be expected to produce subsidence in the Charlevoix area. Since repeated levellings restrict the vertical crustal motion during this time interval to less than 2 cm, we conclude that either the effective aspect-ratio of cracks is less than 0.0001 or the process of crack closure occurred in a number of unconnected regions. More specifically the crack deformation would have to occur in isolated inclusions less than 1 km in diameter and no deeper than about 6 km. The process of crack closure may have been triggered by the passage of seismic waves from the M = 5.0 earthquake.  相似文献   

15.
ABSTRACT

Strong seismic anisotropy is observed in many subduction zones. This effect is attributed partly to subducting oceanic crust that is transformed into blueschist facies rocks. Because blueschist facies constituents such as glaucophane, epidote, and phengite show strong anisotropic elasticity, seismic anisotropy in subducting oceanic crust can be attributed to the lattice preferred orientation (LPO) of these minerals. We studied the deformation fabrics and seismic properties of phengite-rich epidote–glaucophane schists from the Franciscan Complex of Ring Mountain, California. The samples are composed mainly of glaucophane, epidote, and phengite. Some samples contain abundant phengite, the maximum being 40%. The LPOs of glaucophane showed that the [001] axes are aligned subparallel to lineation, and both (110) poles and [100] axes are aligned subnormal to foliation. The epidote [001] axes are aligned subnormal to foliation, with both (110) and (010) poles aligned subparallel to lineation. The LPOs of phengite are characterized by the maxima of [001] axes subnormal to foliation, and both (110) and (010) poles and [100] axes are aligned in a girdle subparallel to foliation. The phengite showed substantially strong seismic anisotropy (AVP = 42%, max.AVS = 37%). The glaucophane schist with abundant phengite showed significantly stronger seismic anisotropy (AVP = 30%, max.AVS = 23%) than the epidote–glaucophane schist (AVP = 13%, max.AVS = 9%). When the subduction angle of phengite-rich glaucophane schist is considered, the polarization direction of the fast S-waves for vertically propagating S-waves changed to a nearly trench-parallel direction for the subduction angle of 45?60°, and the S-wave anisotropy became stronger for vertically propagating S-waves with increasing subduction angles. Our data showed that phengite-rich blueschist facies rock can therefore contribute to the strong trench-parallel seismic anisotropy occurring at the subducting oceanic crust and at the slab–mantle interface in many subduction zones.  相似文献   

16.
This study presents a detailed analysis of the seismic records of a strong explosion that occurred on 21 September 2001 at a chemical complex located south of Toulouse, France, and provoked important damages. The explosion, which is equivalent to a 3.4 magnitude earthquake, has been recorded at most of the stations of the National Seismological Network, as well as at a station under test at the ‘Observatoire Midi-Pyrénées’, 4.2 km away from the epicentre. The main seismic phases are interpreted using the known crustal structures, and a modelling with synthetic seismograms is performed. To cite this article: A. Souriau et al., C. R. Geoscience 334 (2002) 155–161.  相似文献   

17.
Earthquake source parameters and crustal Q are being estimated simultaneously through the inversion of S-wave displacement spectra from three-component recordings of ten local cratonic intraplate earthquakes from 3-6 broadband stations in the eastern Indian shield, wherein, an iterative Levenberg-Marquardt inversion technique is used. The estimated seismic moment (Mo) and source radii (r) vary from 7.4 x 1012 to 7.1 x 1014 N-m and 144.2 to 211.3 m, respectively, while estimated stress drops (Δσ) and multiplicative factor (Emo) values range from 0.11 to 4.13 MPa and 1.33 to 2.16, respectively. The corner frequencies range from 6.23 to 8.62 Hz while moment magnitudes vary from 2.44 to 3.57. The radiated seismic energy and apparent stresses range from 8.3 x 106 to 2.0 x 1010 Joules and 0.06 to 0.94 MPa, respectively, wherein the estimated corner frequencies and seismic moment satisfy the relation Mo ∞ f c –(3+ε) for ε = 12.7. Thus, the source scaling of these events clearly deviates from the self-similarity i.e. f–3. Estimated Zuniga parameters reveal that all selected events satisfy the partial stress drop model, which is in good agreement with the global observations. Our estimated crustal S-wave quality factors vary from 1091 to 4926 with an average of 3006, suggesting a less heterogeneous crustal structure underlying the study region.We also perform moment tensor inversion of five selected local events using ISOLA software, which reveals that the dominant deformation mode for the eastern Indian shield is left-lateral strike slip motion with minor normal dip-slip component on an almost vertical plane. This observation suggests that neotectonic vertical movements might have played a key role in generating these earthquakes. Our modeling also depicts that the seismically mildly active Singhbhum shear zone and Eastern Ghats mobile belt are characterized by the left-lateral strike motion while two events in the Chotanagpur half graben belt suggest a normal dip-slip motion along a south dipping plane. A north-south orientation of P-axis is found to be dominant in the area, which is consistent with the prevailing north–south compression over the Indian plate.  相似文献   

18.
《Comptes Rendus Geoscience》2019,351(4):303-311
This paper presents the first seismic measurements of serpentinite of Bou Azzer ophiolite, central Anti-Atlas of Morocco, including seismic velocities and anisotropy. Two serpentinite samples collected from the Ait Ahmane fault zone were analyzed in order to define the mineralogical and seismic features of the natural serpentinite of the Bou Azzer ophiolite. The mineralogical features were investigated using microscopic observation and Raman spectroscopy, while the seismic features were performed using an Electron Backscatter Diffraction (EBSD) instrument. Microscopy and spectroscopy analyses confirmed that the investigated serpentinite suffers from a variable degree of serpentinization, and the antigorite is the dominant variety of serpentine minerals in the study area. The crystal preferred orientation (CPO) results show that the axis [001] of the antigorite deformation is aligned subnormal to the foliation, while the axis [010] is mostly aligned subparallel to the lineation. The seismic anisotropy results are depending on serpentine amount in the rock samples. The sample with a low serpentine amount (30%) shows lowest P- and S-wave anisotropy (Vp = 7.2% and AVs = 6.55%), while the sample with a high amount of serpentine (85%) presents highest P-wave and S-wave anisotropy (Vp = 8.6% and AVs = 11.06%). Consequently, the results indicate that seismic anisotropy increases when increasing the antigorite amount.  相似文献   

19.
深地震测深是探测壳幔岩石圈精细速度结构、探讨岩石圈变形和演化过程的一种有效方法,在青藏高原隆升、克拉通裂解等大陆动力学研究中已发挥了重要的作用.然而,地震测深方法与深部动力学研究的结合尚处于现象描述为主的状态.因此,本文对前人利用深地震测深资料进行深部动力学研究的相关方法进行了回顾与总结:宽角反射/折射地震震相特征具有...  相似文献   

20.
Earthquake source parameters and crustal \(Q_{0}\) values for the 138 selected local events of (\(\hbox {M}_{\mathrm{w}}{:}2.5{-}4.4\)) the 2001 Bhuj earthquake sequence have been computed through inversion modelling of S-waves from three-component broadband seismometer data. SEISAN software has been used to locate the identified local earthquakes, which were recorded at least three or more stations of the Kachchh seismological network. Three component spectra of S-wave are being inverted by using the Levenberg–Marquardt non-linear inversion technique, wherein the inversion scheme is formulated based on \(\omega ^{2}\) source model. SAC Software (seismic analysis code) is being utilized for calculating three-component displacement and velocity spectra of S-wave. The displacement spectra are used for estimating corner frequency (in Hz) and long period spectral level (in nm-s). These two parameters play a key role in estimating earthquake source parameters. The crustal \({Q}_{0}\) values have been computed simultaneously for each component of three-component broadband seismograph. The estimated seismic moment (\(M_{0}\)) and source radius (r) using S-wave spectra range from 7.03E+12 to 5.36E+15 N-m and 178.56 to 565.21 m, respectively. The corner frequencies for S-wave vary from 3.025 to 7.425 Hz. We also estimated the radiated energy (\(E_{S}\)) using velocity spectra, which is varying from 2.76E+06 to 4.07E+11 Joules. The estimated apparent stress drop and static stress drop values range from 0.01 to 2.56 and 0.53 to 36.79 MPa, respectively. Our study also reveals that estimated \(Q_{0}\) values vary from 119.0 to 7229.5, with an average \(Q_{0}\) value of 701. Another important parameter, by which the earthquake rupture process can be recognized, is Zuniga parameter. It suggests that most of the Kachchh events follow the frictional overshoot model. Our estimated static stress drop values are higher than the apparent stress drop values. And the stress drop values are quite larger for intraplate earthquakes than the interplate earthquakes.  相似文献   

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

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