首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Friction data used in modelling studies of subduction zone seismogenesis are often poorly representative of in situ conditions. We investigated the influence of in situ effective stresses and temperatures on the frictional properties of (simulated) fault gouges, prepared either from Nankai ODP material or illite shale, at sliding velocities approaching those relevant for earthquake nucleation and SSEs. Biaxial (double direct shear) experiments were performed at room temperature, normal stresses of 5–30 MPa, and sliding velocities of 0.16–18 μm/s. All materials exhibited velocity strengthening under these conditions, along with an increase in the friction coefficient and slip hardening rate with increasing normal stress. Illite gouge showed increased velocity strengthening towards higher normal stresses. The effect of temperature was investigated by means of ring shear experiments on illite gouge at 200–300 °C, an effective normal stress of 170 MPa, a pore-fluid pressure of 100 MPa and sliding velocities of 1–100 μm/s. These experiments showed a transition from velocity strengthening to velocity weakening at ∼250 °C. Our results provide a possible explanation for the updip seismogenic limit within subduction zone megathrusts and imply an enhanced tendency for earthquake nucleation and SSEs at low effective normal stresses.  相似文献   

2.
In 2001, a special issue of the Bulletin of the Seismological Society of America (BSSA) featured seismological research for the 1999 Chi–Chi Taiwan earthquake. This study uses source parameters suggested by the first author in this special issue to estimate static stress drop associated with the Chi–Chi earthquake. The waveform simulation method was used to carefully examine these source parameters. The simulation results indicate that source parameters, inferred from near-fault observations, are well determined. According to the rupture area and slip, the static stress drops (Δσs) obtained were distributed between a small value of 47 bars near the epicentral region and a much larger value (>200 bars) to the north. Similar trends in dynamic stress drop (Δσd) were also recognized by the first author in his paper published in 2001 BSSA special issue. Comparing the Δσs with Δσd, satisfies the relation Δσs/Δσd ≈ 1. This relation suggests that fault motion is mostly spent releasing seismic wave energy during the rupture process of the Chi–Chi earthquake. The consistency between static and dynamic stress drops thus provides a measure of energy-moment (Es/M0) ratios, which range from 9.0 × 10−5 to 6.5 × 10−4. The average Es/M0 ratio estimated for the northern portions of the fault is 3.4 × 10−4, which is about 3 times that of the south. Such a high Es/M0 ratio can be interpreted as having low strength in the rupture for the northern portions of the fault, where the fault would release less energy per unit rupture surface to create the new rupture.  相似文献   

3.
Several source parameters (source dimensions, slip, particle velocity, static and dynamic stress drop) are determined for the moderate-size October 27th, 2004 (MW = 5.8), and the large August 30th, 1986 (MW = 7.1) and March 4th, 1977 (MW = 7.4) Vrancea (Romania) intermediate-depth earthquakes. For this purpose, the empirical Green's functions method of Irikura [e.g. Irikura, K. (1983). Semi-Empirical Estimation of Strong Ground Motions during Large Earthquakes. Bull. Dis. Prev. Res. Inst., Kyoto Univ., 33, Part 2, No. 298, 63–104., Irikura, K. (1986). Prediction of strong acceleration motions using empirical Green's function, in Proceedings of the 7th Japan earthquake engineering symposium, 151–156., Irikura, K. (1999). Techniques for the simulation of strong ground motion and deterministic seismic hazard analysis, in Proceedings of the advanced study course seismotectonic and microzonation techniques in earthquake engineering: integrated training in earthquake risk reduction practices, Kefallinia, 453–554.] is used to generate synthetic time series from recordings of smaller events (with 4 ≤ MW ≤ 5) in order to estimate several parameters characterizing the so-called strong motion generation area, which is defined as an extended area with homogeneous slip and rise time and, for crustal earthquakes, corresponds to an asperity of about 100 bar stress release [Miyake, H., T. Iwata and K. Irikura (2003). Source characterization for broadband ground-motion simulation: Kinematic heterogeneous source model and strong motion generation area. Bull. Seism. Soc. Am., 93, 2531–2545.] The parameters are obtained by acceleration envelope and displacement waveform inversion for the 2004 and 1986 events and MSK intensity pattern inversion for the 1977 event using a genetic algorithm. The strong motion recordings of the analyzed Vrancea earthquakes as well as the MSK intensity pattern of the 1977 earthquake can be well reproduced using relatively small strong motion generation areas, which corresponds to small asperities with high stress drops (300–1200 bar) and high particle velocities (3–5 m/s). These results imply a very efficient high-frequency radiation, which has to be taken into account for strong ground motion prediction, and indicate that the intermediate-depth Vrancea earthquakes are inherently different from crustal events.  相似文献   

4.
The Mw 9.0 Tohoku-Oki earthquake that occurred off the Pacific coast of Japan on March 11, 2011, was followed by thousands of aftershocks, both near the plate interface and in the crust of inland eastern Japan. In this paper, we report on two large, shallow crustal earthquakes that occurred near the Ibaraki-Fukushima prefecture border, where the background seismicity was low prior to the 2011 Tohoku-Oki earthquake. Using densely spaced geodetic observations (GPS and InSAR datasets), we found that two large aftershocks in the Iwaki and Kita-Ibarake regions (hereafter referred to as the Iwaki earthquake and the Kita-Ibarake earthquake) produced 2.1 m and 0.44 m of motion in the line-of-sight (LOS), respectively. The azimuth-offset method was used to obtain the preliminary location of the fault traces. The InSAR-based maximum offset and trace of the faults that produced the Iwaki earthquake are consistent with field observations. The fault location and geometry of these two earthquakes are constrained by a rectangular dislocation model in a multilayered elastic half-space, which indicates that the maximum slips for the two earthquakes are 3.28 m and 0.98 m, respectively. The Coulomb stress changes were calculated for the faults following the 2011 Mw 9.0 Tohoku-Oki earthquake based on the modeled slip along the fault planes. The resulting Coulomb stress changes indicate that the stresses on the faults increased by up to 1.1 MPa and 0.7 MPa in the Iwaki and Kita-Ibarake regions, respectively, suggesting that the Tohoku-Oki earthquake triggered the two aftershocks, supporting the results of seismic tomography.  相似文献   

5.
Pseudotachylytes are generally considered to be remnent products of past earthquakes, conveying significant information, which provide improved insight into the fault behavior and their mechanical properties. Pseudotachylytes and cataclastic rocks are exposed in the Pengguan complex along the Yingxiu-Beichuan fault zone, which ruptured during the 2008 Wenchuan earthquake, at the eastern margin of the Tibetan Plateau in Sichuan province, China. From outcrop investigations, pseudotachylyte veins are gray, dark-gray, brown and black, with thicknesses ranging from several mm to ∼10 cm. Microstructural observations, via optical microscope and SEM, show that distinct textures, such as spherulites, microlites, honeycomb-like vesicle structures, flow textures, irregularly shaped and deep embayed clasts, are present in the pseudotachylytes, suggesting that the pseudotachylytes are indeed melt-origin. Different colored veins and their overprinting relations observed at both macroscopic and microscopic scales, demonstrate that seismic faulting events and the associated generations of pseudotachylyte occurred repeatedly along the Yingxiu-Beichuan fault zone. Friction melt lubrication, thermal pressurization and mechanical lubrication, all can exist during one fault slip motion, indicating that multiple fault weakening mechanisms may work during a single coseismic event.  相似文献   

6.
Two and a half years prior to China’s M7.9 Wenchuan earthquake of May 2008, at least 300 million metric tons of water accumulated with additional seasonal water level changes in the Minjiang River Valley at the eastern margin of the Longmen Shan. This article shows that static surface loading in the Zipingpu water reservoir induced Coulomb failure stresses on the nearby Beichuan thrust fault system at <17 km depth. Triggering stresses exceeded levels of daily lunar and solar tides and perturbed a fault area measuring 416 ± 96 km2. These stress perturbations, in turn, likely advanced the clock of the mainshock and directed the initial rupture propagation upward towards the reservoir on the "Coulomb-like" Beichuan fault with rate- and state-dependent frictional behavior. Static triggering perturbations produced up to 60 years (0.6%) of equivalent tectonic loading, and show strong correlations to the coseismic slip. Moreover, correlations between clock advancement and coseismic slip, observed during the mainshock beneath the reservoir, are strongest for a longer seismic cycle (10kyr) of M > 7 earthquakes. Finally, the daily event rate of the micro-seismicity (M ≥ 0.5) correlates well with the static stress perturbations, indicating destabilization.  相似文献   

7.
We localized crustal earthquakes in the Andean arc, between 35°S and 36°S, from December 2009 to May 2010. This research shows a seismicity increase, in a narrow longitudinal area, of more than nine times after the great Mw 8.8 Maule earthquake.The localized seismicity defines an area of ∼80 km long and ∼18 km wide and NNW to NNE trend. The Md magnitudes varied from 0.7 to 3.1 except for two earthquakes with Mw of 3.9 and 4.5, located in the northern end of the area. The focal mechanisms for these two last events were normal/strike-slip and strike-slip respectively.During 2011, a network of 13 temporary stations was installed in the trasarc region in Malargüe, Argentina. Sixty earthquakes were localized in the study region during an 8 month period.We explored how changes in Coulomb conditions associated with the mega-thrust earthquake triggered subsequent upper-plate events in the arc region. We assumed the major proposed structures as receiver faults and used previously published earthquake source parameters and slip distribution for the Maule quake. The largest contribution to static stress change, up to 5 bars, derives from unclamping resulting consistent with co-seismic dilatational deformation inferred from GPS observations in the region and subsidence in nearby volcanoes caused by magma migration.Three different Quaternary tectonic settings–extensional, strike-slip and compressional-have been proposed for the arc region at these latitudes. We found that the unclamping produced by the Maule quake could temporarily change the local regime to normal/strike-slip, or at least it would favor the activation of Quaternary NNE to N-trending dextral strike-slip faults with dextral transtensional movement.  相似文献   

8.
A damaging and widely felt moderate (Mw 5.0) earthquake occurred in the Talala region of Saurashtra, Gujarat (western India) on November 6, 2007. The highly productive sequence comprised about 1300 micro earthquakes (M > 0.5) out of which 325 of M ? 1.5 that occurred during November 6, 2007–January 10, 2008 were precisely located. The spatial aftershock distribution revealed a NE–SW striking fault in accordance with the centroid moment tensor solution, which in turn implies left-lateral motion. The orientation and sense of shear are consistent with similarly orientated geological fault identified in the area from satellite imagery and field investigation.The aftershocks temporal decay, b-value of frequency–magnitude distribution, spatial fractal dimension, D, and slip ratio (ratio of the slip occurred on the primary fault to the total slip) were examined with the purpose to identify the properties of the sequence. The high b-value (1.18 ± 0.01) may be attributed to the paucity of the larger (M ? 4.0) aftershocks and reveals crustal heterogeneity and low stress regime. The high p-value (1.10 ± 0.39), implying fast decay rate of aftershocks, evidences high surface heat flux. A value of the spatial fractal dimension (D) equal to 2.21 ± 0.02 indicates random spatial distribution and source in a two-dimensional plane that is being filled-up by fractures. A slip ratio of 0.42 reveals that more slip occurred on secondary fault systems.The static Coulomb stress changes due to the coseismic slip of the main shock, enhanced off fault aftershock occurrence. The occurrence of a moderate earthquake (Mw 4.3) on October 5, 2008 inside a region of positive Coulomb stress changes supports the postulation on aftershock triggering. When the stress changes were resolved on a cross section including the stronger (M4.8) foreshock plane that is positioned adjacent to the main fault, it became evident that the activity continued there due to stress transfer from the main rupture.  相似文献   

9.
The generation of intraplate earthquakes has been attributed to perturbations in the stress regime, either due to surface and sub-surface loading or strength weakening of the rock mass. The present work aims at estimating the intraplate stresses associated with topography and crustal density inhomogeneities beneath the Deccan Volcanic Province (DVP). A layered crustal model with irregular interfaces of small amplitude has been used for elastostatic stress calculations. The computed principal stress differences show a significant concentration at 5–20 km depths beneath the western side of the region. The maximum magnitude of principal stress difference occurs beneath the Karad at a depth of 10 km with a value of 60 MPa. The deviatoric stress estimates are further superposed on inferred stresses due to the regional plate tectonic forces. These results show principal stress difference concentrations beneath the Koyna, Poona and Karad regions which may thus be more vulnerable to brittle failure. It is also seen that the principal total stress directions point to the strike slip motion at Koyna, similar to that which is associated with the 1967 Koyna earthquake.  相似文献   

10.
To characterize the fault-related rocks within the Chelungpu fault, we performed X-ray computed tomography (CT) image analyses and microstructural observations of Hole B core samples from the Taiwan Chelungpu-fault Drilling Project. We identified the slip zone associated with the 1999 Chi-Chi earthquake, within the black gouge zone in the shallowest major fault zone, by comparison with previous reports. The slip zone was characterized by low CT number, cataclastic (or ultracataclastic) texture, and high possibility to have experienced a mechanically fluidized state. Taking these characteristics and previous reports of frictional heating in the slip zone into consideration, we suggested that thermal pressurization was the most likely dynamic weakening mechanism during the earthquake.  相似文献   

11.
The Great Lisbon earthquake has the largest documented felt area of any shallow earthquake and an estimated magnitude of 8.5–9.0. The associated tsunami ravaged the coast of SW Portugal and the Gulf of Cadiz, with run-up heights reported to have reached 5–15 m. While several source regions offshore SW Portugal have been proposed (e.g.— Gorringe Bank, Marquis de Pombal fault), no single source appears to be able to account for the great seismic moment as well as all the historical tsunami amplitude and travel time observations. A shallow east dipping fault plane beneath the Gulf of Cadiz associated with active subduction beneath Gibraltar, represents a candidate source for the Lisbon earthquake of 1755.Here we consider the fault parameters implied by this hypothesis, with respect to total slip, seismic moment, and recurrence interval to test the viability of this source. The geometry of the seismogenic zone is obtained from deep crustal studies and can be represented by an east dipping fault plane with mean dimensions of 180 km (N–S) × 210 km (E–W). For 10 m of co-seismic slip an Mw 8.64 event results and for 20 m of slip an Mw 8.8 earthquake is generated. Thus, for convergence rates of about 1 cm/yr, an event of this magnitude could occur every 1000–2000 years. Available kinematic and sedimentological data are in general agreement with such a recurrence interval. Tsunami wave form modeling indicates a subduction source in the Gulf of Cadiz can partly satisfy the historical observations with respect to wave amplitudes and arrival times, though discrepancies remain for some stations. A macroseismic analysis is performed using site effect functions calculated from isoseismals observed during instrumentally recorded strong earthquakes in the region (M7.9 1969 and M6.8 1964). The resulting synthetic isoseismals for the 1755 event suggest a subduction source, possibly in combination with an additional source at the NW corner of the Gulf of Cadiz can satisfactorily explain the historically observed seismic intensities. Further studies are needed to sample the turbidites in the adjacent abyssal plains to better document the source region and more precisely calibrate the chronology of great earthquakes in this region.  相似文献   

12.
After the 2005 Kashmir earthquake, we mapped surface ground fractures in Tangdhar, Uri, Rajouri and Punch sectors and liquefaction features in Jammu area lying close to the eastern side of the Line of Control (LOC) in Kashmir, India. The NW trending ground fractures occurred largely in the hanging wall zone of the southeastern extension of the causative fault in Tangdhar and Uri sectors. The principal compressive stress deduced from the earthquake induced ground fractures is oriented at N10°, whereas the causative Balakot–Bagh fault strikes 330°. The fault-plane solution indicates primarily SW thrusting of the causative fault with a component of strike–slip motion. The ground fractures reflect pronounced strike–slip together with some tensile component. The Tangdhar area showing left-lateral strike–slip motion lies on the hanging wall, and the Uri region showing right-lateral strike–slip movement is located towards the southeastern extension of the causative fault zone. The shear fractures are related to static stress that was responsible for the failure of causative fault. The tensile fractures with offsets are attributed to combination of both static and dynamic stresses, and the fractures and openings without offsets owe their origin due to dynamic stress. In Punch–Rajouri and Jammu area, which lies on the footwall, the fractures and liquefactions were generated by dynamic stress. The occurrence of liquefaction features in the out board part of the Himalayan range front near Jammu is suggestive of stress transfer  230 km southeast of the epicenter. The Balakot–Bagh Fault (BBF), the Muzaffarabad anticline, the rupture zone of causative fault and the zone of aftershocks — all are aligned in a  25 km wide belt along the NW–SE trending regional Himalayan strike of Kashmir region and lying between the MBT and the Riasi Thrust (Murree Thrust), suggesting a seismogenic zone that may propagate towards the southeast to trigger an earthquake in the eastern part of the Kashmir region.  相似文献   

13.
The Killari earthquake of September 29, 1993 (Mw=6.2) in peninsular India triggered several aftershocks that were recorded by a network of 21 stations. We computed the change in regional static stress caused by coseismic slip on the earthquake rupture and correlated it with the aftershocks with a view to constrain some of the rupture parameters of this earthquake. We evaluated the six available estimates of fault plane solutions for this earthquake and concluded that reverse slip on a 42° dipping, N112° trending fault, which extends up to the surface from a depth of 7 km, produces maximum correlation between the increased static stress and aftershock distribution. Our analysis suggests that the majority of coseismic slip occurred on the part of the rupture that lies in the depth range of 3–6.5 km.  相似文献   

14.
A common problem encountered in studies of gouge-bearing natural faults is the difficulty of ascertaining whether the observed gouge was sheared seismically or aseismically; this problem arises because of the scarcity of indicators of fault slip rates for gouge. Recently, clay–clast aggregates (CCAs; a CCA comprises a clastic core mantled by a rim of ultrafine particles) were proposed as a possible indicator of seismic slip in gouge, on the basis of shear experiments on gouge at seismic slip rates. To examine the processes and conditions of CCA formation, we conducted rotary shear experiments on quartz and quartz–bentonite gouges under normal stresses (0.3–3.0 MPa) and slip rates (0.0005–1.3 m s−1), and in both room-humidity (room-dry) and water-saturated (wet) conditions. We found that CCAs could be produced in room-dry gouges even at the lowest slip rates, which are considerably slower than actual seismic slip rates. This finding demonstrates that thermal pressurization and fluidization at elevated temperature during seismic slip are not necessarily needed for the formation of CCAs, contrary to previous views. Given the occurrence of CCAs over a wide range of slip rates, we suggest that the presence of CCAs is not an unequivocal indicator of fault slip at seismic slip rates.  相似文献   

15.
剪切带的流体压力因关系着地震断层的滑动强度以及应力释放过程一直备受重视。断层带的渗透率结构对认识断层带内流体活动状况、断层强度、摩擦稳定性以及同震过程等至关重要。文中介绍了断层带渗透率研究的原理、方法,总结了国内外最新的研究成果。脆性域断层的渗透率表现为核部低、两侧高的典型结构,断层泥通常具有最低的渗透率和强各向异性。在地震周期中,断层带的渗透性表现为突然增强-逐渐降低的过程。胶结、水岩反应及溶解-沉淀是间震期断层愈合的主要方式。热压作用是一种有效的断层弱化机制。当断层泥的渗透率小于10-18 m2量级时,地震会伴随强烈的同震热压效应。  相似文献   

16.
Analysis of vertical crustal deformation data in the southwestern part of Shikoku, southwest Japan, suggests that the Nankaido earthquake of 1946 (Mw = 8.1), which is a principal interplate thrust earthquake, was accompanied by subsidiary faulting on a splay fault adjacent to the coast of Shikoku. Discarding crustal movement resulting from the main thrusting of the Nankaido earthquake, local leveling data are explained by slip on a simple rectangular thrust fault located just offshore of Shikoku. Although it is difficult to constrain the fault location, a possible result is a high-angle thrust dipping landward at an angle of about 70°, with a dislocation of about 1.5 m, and source dimensions of 30 × 13 km along strike and dip. respectively. This result indicates that the fault may be one of the steeply dipping subsidiary faults branching from the main low-angle thrust, as was the case in the Alaska earthquake of 1964. Although several lines of evidence suggest that this faulting occurred as slow aseismic slip, its discrimination from the main seismic event is extremely difficult. This kind of high-angle thrusting just offshore of the coast would play an important role for the formation of the marine terraces during the late Quaternary period.  相似文献   

17.
Seismic imaging together with global positioning system (GPS) and crustal stress data analyses show that the Mw7.9 2008 Wenchuan earthquake occurred within a distinct area of high crustal stress (～ 17.5 MPa) and high Poisson's ratio (7–10%) anomalies centered on the Longmen-Shan (Shan means Mountain in Chinese) tectonic fault belt. Low P-wave and S-wave velocities in the southwest (SW) segment contrast with high-velocity anomalies in the central portion (CP) and northeast (NE) segment within the uppermost ～ 15 km depths along the tectonic fault belt, though a presumably ductile zone with low-velocity anomalies separates the CP and NE segment. The rupture initiated near the southwestern end of the CP at a zone of high Poisson's ratio (σ) which extends down into the lower crust. These low-velocity and high-σ anomalies immediately below the source hypocenter, together with the high crustal stress, indicate the presence of high-pressure fluids from the lower crust, which might have reduced the mechanical strength of the fractured rock matrix and triggered the earthquake. Our study suggests that the structural heterogeneity and high crustal stress played an important role in the nucleation of the Wenchuan earthquake and its rupture process.  相似文献   

18.
Teleseismic and strong-motion data are inverted to determine the rupture process during the November 1999 Düzce earthquake in NW Turkey. The fault geometry, rise time and rupture velocity are determined from the aftershock distribution and preliminary inversions of the teleseismic data. Joint inversion of the teleseismic and strong-motion data is then carried out for the slip distribution. We obtain the strike 264°, dip 64°, rake −172°, seismic moment 5.0×1019 N m (Mw 7.1), and average stress drop 7 MPa. This earthquake was characterized by bilateral fault rupture and asymmetric slip distribution. Two asperities (areas of large slip) are identified, the eastern one being 1.5 times larger than the western one. The derived slip distribution is consistent with the aftershock distribution, surface rupture and damage. The point of rupture initiation in this Düzce earthquake coincided with the eastern tip of the aftershock distribution of the August 1999 Izmit earthquake.  相似文献   

19.
The 2009 L'Aquila sequence activated a normal fault system 50 km long in the Central Apennines, composed of two main NW-trending faults 12–16 km long: the main high angle L'Aquila segment and the Campotosto listric fault.The MW 6.1 L'Aquila mainshock nucleated on the Paganica fault at a depth of ∼8.6 km and cut through the upper crust producing coseismic surface slip of up to 10 cm observed along a strike length of ∼13 km. Analysis of historical seismicity and data collected in paleo-seismological trenches suggest that this event filled a >500-year gap. In contrast, the blind Campotosto listric fault is composed of different fault segments displaying abrupt changes in dip at a depth where major events nucleate suggesting a rheological and geometrical control on stress concentration.A foreshock sequence that started around 4 months before the L'Aquila mainshock activated the deepest portion of the Paganica fault and marked the onset of large variations in elastic properties of the crustal volume. The variations have been modelled in terms of dilatancy and diffusion processes, corroborating the hypothesis that fluids play a key role in the nucleation process of extensional faults in the crust.  相似文献   

20.
To investigate the strength of frictional sliding and stability of mafic lower crust, we conducted experiments on oven-dried gabbro gouge of 1 mm thick sandwiched between country rock pieces (with gouge inclined 35° to the sample axis) at slip rates of 1.22 × 10− 3 mm/s and 1.22 × 10− 4 mm/s and elevated temperatures up to 615 °C. Special attention has been paid to whether transition from velocity weakening to velocity strengthening occurs due to the elevation of temperature.Two series of experiments were conducted with normal stresses of 200 MPa and 300 MPa, respectively. For both normal stresses, the friction strengths are comparable at least up to 510 °C, with no significant weakening effect of increasing temperature. Comparison of our results with Byerlee's rule on a strike slip fault with a specific temperature profile in the Zhangbei region of North China shows that the strength given by experiments are around that given by Byerlee's rule and a little greater in the high temperature range.At 200 MPa normal stress, the steady-state rate dependence a − b shows only positive values, probably still in the “run-in” process where velocity strengthening is a common feature. With a normal stress of 300 MPa, the values of steady-state rate dependence decreases systematically with increasing temperature, and stick-slip occurred at 615 °C. Considering the limited displacement, limited normal stress applied and the effect of normal stress for the temperatures above 420 °C, it is inferred here that velocity weakening may be the typical behaviour at higher normal stress for temperature above 420 °C and at least up to 615 °C, which covers most of the temperature range in the lower crust of geologically stable continental interior. For a dry mafic lower crust in cool continental interiors where frictional sliding prevails over plastic flow, unstable slip nucleation may occur to generate earthquakes.  相似文献   

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

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