首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
A crack model in antiplane shear configuration is shown representing creep processes interpreted in terms of 'viscous' deformation of a narrow plastic layer, characterized by inhomogeneous rheological properties, embedded within a homogeneous elastic medium. The evolution in time of slip and stress over the crack plane is studied through a truncated expansion in Chebyshev polynomials, and convergence is proved to be fast in the simple examples considered. Finite-stress solutions are found which are compatible with constitutive relations of elasto-plastic materials and furthermore these allow us to simulate creep propagation and stress transfer between locked and unlocked fault segments. This model provides a simple interpretation of the shallow depth of the seismogenic layer observed in several areas of the world and lends itself to modelling creep processes during either post-seismic rebound or pre-seismic stress buildup. Stress transfer is accomplished mostly by the slow extension of the creeping section. During a seismic cycle it is envisaged that different regimes dominate over deep, intermediate and shallow sections of faults: (i) slow pre-seismic stress build-up accompanied by creep and stress migration toward intermediate depths; (ii) brittle fracture over shallow and intermediate sections of faults; (iii) post-seismic rebound over intermediate and deep sections of faults. The present crack model, while providing finite-stress solutions, allows a better understanding of how stress may accommodate at different depths over a fault plane during a seismic cycle.  相似文献   

2.
Summary. A possible mechanism for the occurrence of slow earthquakes is investigated by calculating numerical solutions for the dynamical rupture process on a quasi-three-dimensional fault with heterogeneous frictional strengths. Experimental friction laws for the dependence of sliding frictional stress on slip velocity, which are based on the cohesive properties of fault asperities, are taken into considerations.
It is found that the applied stress does not drop very rapidly with time and the rupture velocity remarkably decreases as the dependence on slip-velocity becomes smaller. These deceleration effects for the rupture propagation are greatly enhanced with increasing heterogeneities in the distribution of frictional strength and as the initial shear stress has lower levels with respect to the average strength. For these cases, the growth of rupture is extremely slow in a nucleus region with the dimension as large as 10 times the initial rupture length, and gains a terminal velocity dependent on the above factors. The displacement-time function becomes noticeably extended in these cases, and indicates a stick—slip-like phenomena in the extended time interval for a strongly heterogeneous fault.
It seems that these results could explain the characteristic features of slow earthquakes.  相似文献   

3.
Geological studies show evidence for temporal clustering of large earthquakes on individual fault systems. Since post-seismic deformation due to the inelastic rheology of the lithosphere may result in a variable loading rate on a fault throughout the interseismic period, it is reasonable to expect that the rheology of the non-seismogenic lower crust and mantle lithosphere may play a role in controlling earthquake recurrence times. We study this phenomenon using a 2-D, finite element method continuum model of the lithosphere containing a single strike-slip fault. This model builds on a previous study using a 1-D spring-dashpot-slider analogue of a single fault system to study the role of Maxwell viscoelastic relaxation in producing non-periodic earthquakes. In our 2-D model, the seismogenic portion of the fault slips when a predetermined yield stress is exceeded; stress accumulated on the seismogenic fault is shed to the viscoelastic layers below and recycled back to the seismogenic fault through viscoelastic relaxation. We find that random variation of the fault yield stress from one earthquake to the next can cause the earthquake sequence to be clustered; the amount of clustering depends on a non-dimensional number, W , called the Wallace number defined as the standard deviation of the randomly varied fault yield stress divided by the effective viscosity of the system times the tectonic loading rate. A new clustering metric based on the bimodal distribution of interseismic intervals allows us to investigate clustering behaviour of systems over a wide range of model parameters and those with multiple viscoelastic layers. For models with   W ≥ 1  clustering increases with increasing W , while those with   W ≤ 1  are unclustered, or quasi-periodic.  相似文献   

4.
The coupled plate interface of subduction zones—commonly called the seismogenic zone—has been recognized as the origin of fatal earthquakes. A subset of the after-shock series of the great Antofagasta thrust-type event (1995 July 30; M w = 8.0) has been used to study the extent of the seismogenic zone in northern Chile. To achieve reliable and precise hypocentre locations we applied the concept of the minimum 1-D model, which incorporates iterative simultaneous inversion of velocity and hypocentre parameters. The minimum 1-D model is complemented by station corrections which are influenced by near-surface velocity heterogeneity and by the individual station elevations. By relocating mine blasts, which were not included in the inversion, we obtain absolute location errors of 1  km in epicentre and 2  km in focal depth. A study of the resolution parameters ALE and DSPR documents the importance of offshore stations on location accuracy for offshore events. Based on precisely determined hypo-centres we calculate a depth of 46  km for the lower limit of the seismogenic zone, which is in good agreement with previous studies for this area. For the upper limit we found a depth of 20  km. Our results of an aseismic zone between the upper limit of the seismogenic zone and the surface correlates with a detachment zone proposed by other studies; the results are also in agreement with thermal studies for the Antofagasta forearc region.  相似文献   

5.
We combine Global Positioning System (GPS) measurements with forward modelling of viscoelastic relaxation and after-slip to study the post-seismic deformation of the 1997 Umbria-Marche (Central Apennines) moderate shallow earthquake sequence. Campaign GPS measurements spanning the time period 1999–2003 are depicting a clear post-seismic deformation signal. Our results favour a normal faulting rupture model where most of the slip is located in the lower part of the seismogenic upper crust, consistent with the rupture models obtained from the inversion of strong motion data. The preferred rheological model, obtained from viscoelastic relaxation modelling, consists of an elastic upper crust, underlain by a transition zone with a viscosity of 1018 Pa s, while the rheology of deeper layers is not relevant for the observed time-span. Shallow fault creep and after-slip at the base of the seismogenic upper crust are the first order processes behind the observed post-seismic deformation. The deep after-slip, below the fault zone at about 8 km depth, acting as a basal shear through localized time-dependent deformation, identifies a rheological discontinuity decoupling the seismogenic upper crust from the low-viscosity transition zone.  相似文献   

6.
Source history of the 1905 great Mongolian earthquakes (Tsetserleg, Bolnay)   总被引:1,自引:0,他引:1  
Two great Mongolian earthquakes, Tsetserleg and Bolnay, occurred on 1905 July 9 and 23. We determined the source history of these events using body waveform inversion. The Tsetserleg rupture (azimuth N60°) correspond to a N60° oriented branch of the long EW oriented Bolnay fault.
Historical seismograms recorded by Wiechert instruments are digitized and corrected for the geometrical deformation due to the recording system. We use predictive filters to recover the signals lost at the minute marks.
The total rupture length for the Tsetserleg earthquake may reach up to 190 km, in order to explain the width of the recorded body waves. This implies adding 60 km to the previously mapped fault. The rupture propagation is mainly eastward. It starts at the southwest of the central subsegment, showing a left lateral strike-slip with a reverse component. The total duration of the modelled source function is 65 s. The seismic moment deduced from the inversion is 1021 N m, giving a magnitude   M w = 8  .
The nucleation of the Bolnay earthquake was at the intersection between the main fault (375 km left lateral strike-slip) and the Teregtiin fault (N160°, 80 km long right lateral strike-slip with a vertical component near the main fault). The rupture was bilateral along the main fault: 100 km to the west and 275 km to east. It also propagated 80 km to the southeast along the Teregtiin fault. The source duration was 115 s. The moment magnitude Mw varies between 8.3 and 8.5.
The nucleation and rupture depths remain uncertain. We tested three cases: (1) nucleation and rupture depth limited to the seismogenic zone; (2) nucleation in the seismogenic zone and rupture propagation going to the base of the crust and (3) nucleation within the crust–upper mantle interface and rupture propagation within the upper mantle.  相似文献   

7.
Summary We consider a long strike-slip fault in a lithosphere modelled as an elastic slab. To the base of the slab a shear stress distribution is applied which simulates the viscous drag exerted by the asthenosphere. The resulant stress on the fault plane may directly fracture the lithosphere in its brittle upper portion; alternatively it may give rise at first to a stable aseismic sliding in the lower portion. In the latter case, stress concentration due to the deep aseismic slip is the relevant feature of the pre-seismic stress acting on the upper section of the lithosphere. The two cases are examined by use of dislocation theory and their observable effects compared. Different depths of the aseismic slip zone and the presence or absence of a uniform friction on the seismic fault are allowed for. If the model is applied to the San Andreas fault region, where a steady sliding condition actually seems to be present at shallow depth, it turns out that the slip amplitudes commonly associated with large earthquakes are consistent with average basal stress values which can be substantially lower than a few bars, a value often quoted as the steady state basal stress due to a velocity gradient in the upper asthenosphere.  相似文献   

8.
Shear-wave splitting is analysed on data recorded by the High Resolution Seismic Network (HRSN) at Parkfield on the San Andreas fault, Central California, during the three-year period 1988-1990. Shear-wave polarizations either side of the fault are generally aligned in directions consistent with the regional horizontal maximum compressive stress, at some 70° to the fault strike, whereas at station MM in the immediate fault zone, shear-wave polarizations are aligned approximately parallel to the fault. Normalized time delays at this station are found to be about twice as large as those in the rock mass either side. This suggests that fluid-filled cracks and fractures within the fault zone are elastically or seismically different from those in the surrounding rocks, and that the alignment of fault-parallel shear-wave polarizations are associated with some fault-specific phenomenon.
Temporal variations in time delays between the two split shear-waves before and after a ML = 4 earthquake can be identified at two stations with sufficient data: MM within the fault zone and VC outside the immediate fault zone. Time delays between faster and slower split shear waves increase before the ML = 4 earthquake and decrease near the time of the event. The temporal variations are statistically significant at 68 per cent confidence levels. Earthquake doublets and multiplets also show similar temporal variations, consistent with those predicted by anisotropic poroelasticity theory for stress modifications to the microcrack geometry pervading the rock mass. This study is broadly consistent with the behaviour observed before three other earthquakes, suggesting that the build-up of stress before earthquakes may be monitored and interpreted by the analysis of shear-wave splitting.  相似文献   

9.
The Benevento region is part of the southern Apennines seismogenic belt, which experienced large destructive seismic events both in historical and in recent times. The study area lies at the northern end of the Irpinia fault, which ruptured in 1980 with a Ms = 6.9 normal faulting event, which caused about 3000 casualties. The aims of this paper are to image lateral heterogeneities in the upper crust of the Benevento region, and to try to identify the fault segments that are expected to generate such large earthquakes. This work is motivated by the recognition that lithological heterogeneities along major fault zones, inferred from velocity anomalies, reflect the presence of fault patches that behave differently during large rupture episodes. In this paper, we define the crustal structure of the Benevento region by using the background seismicity recorded during 1991 and 1992 by a local seismic array. These data offer a unique opportunity to investigate the presence of structural discontinuities of a major seismogenic zone before the occurrence of the next large earthquake. The main result that we obtained is the delineation of two NW-trending high-velocity zones (HVZs) in the upper crust beneath the Matese limestone massif. These high velocities are interpreted as high-strength regions that extend for 30-40 km down to at least 12 km depth. The correspondence of these HVZs with the maximum intensity regions of historical earthquakes (1688 AD, 1805 AD) suggests that these anomalies delineate the extent of two fault segments of the southern Apenninic belt capable of generating M = 6.5−7 earthquakes. The lateral offset observed between the two segments from tomographic results and isoseismal areas is possibly related to transverse right-lateral faults.  相似文献   

10.
We use teleseismic waveform analysis and locally recorded aftershock data to investigate the source processes of the 2004 Baladeh earthquake, which is the only substantial earthquake to have occurred in the central Alborz mountains of Iran in the modern instrumental era. The earthquake involved slip at 10–30 km depth, with a south-dipping aftershock zone also restricted to the range 10–30 km, which is unusually deep for Iran. These observations are consistent with co-seismic slip on a south-dipping thrust that projects to the surface at the sharp topographic front on the north side of the Alborz. This line is often called the Khazar Fault, and is assumed to be a south-dipping thrust which bounds the north side of the Alborz range and the south side of the South Caspian Basin, though its actual structure and significance are not well understood. The lack of shallower aftershocks may be due to the thick pile of saturated, overpressured sediments in the South Caspian basin that are being overthrust by the Alborz. A well-determined earthquake slip vector, in a direction different from the overall shortening direction across the range determined by GPS, confirms a spatial separation ('partitioning') of left-lateral strike-slip and thrust faulting in the Alborz. These strike-slip and thrust fault systems do not intersect within the seismogenic layer on the north side, though they may do so on the south. The earthquake affected the capital, Tehran, and reveals a seismic threat posed by earthquakes north of the Alborz, located on south-dipping thrusts, as well as by earthquakes on the south side of the range, closer to the city.  相似文献   

11.
In a tectonically active setting large earthquakes are always threats; however, they may also be useful in elucidating the subsurface geology. Instrumentally recorded seismicity is, therefore, widely utilized to extend our knowledge into the deeper crust, especially where basement is involved. It is because the earthquakes are triggered by underground stress changes that usually corresponding to the framework of geological structures. Hidden faults, therefore, can be recognized and their extension as well as orientation can be estimated. Both above are of relevance for assessment on seismic hazard of a region, since the active faults are supposed to be re-activated and cause large earthquakes. In this study, we analysed the 1999 October 22 earthquake sequence that occurred in southwestern Taiwan. Two major seismicity clusters were identified with spatial distribution between depths of 10 and 16 km. One cluster is nearly vertical and striking 032°, corresponding to the strike-slip Meishan fault (MSF) that generated the 1906 surface rupture. Another cluster strikes 190° and dips 64° to the west, which is interpreted as west-vergent reverse fault, in contrast to previous expectation of east vergence. Our analysis of the focal solutions of all the larger earthquakes in the 1999 sequence with the 3-D distribution of all the earthquakes over the period 1990–2004 allows us reinterpret the structural framework and suggest previously unreognized seismogenic sources in this area. We accordingly suggest: (1) multiple detachment faults are present in southwestern Taiwan coastal plain and (2) additional seismogenic sources consist of tear faults and backthrust faults in addition to sources associated with west-vergent fold-and-thrust belt.  相似文献   

12.
Scaling relationships between seismic moment, rupture length, and rupture width have been examined. For this purpose, the data from several previous studies have been merged into a database containing more than 550 events. For large earthquakes, a dependence of scaling on faulting mechanism has been found. Whereas small and large dip-slip earthquakes scale in the same way, the self-similarity of earthquakes breaks down for large strike-slip events. Furthermore, no significant differences in scaling could be found between normal and reverse earthquakes and between earthquakes from different regions. Since the thickness of the seismogenic layer limits fault widths, most strike-slip earthquakes are limited to rupture widths of between 15 and 30 km while the rupture length is not limited. The aspect ratio of dip-slip earthquakes is similar for all earthquake sizes. Hence, the limitation in rupture width seems to control the maximum possible rupture length for these events. The different behaviour of strike-slip and dip-slip earthquakes can be explained by rupture dynamics and geological fault growth. If faults are segmented, with the thickness of the seismogenic layer controlling the length of each segment, strike-slip earthquakes might rupture connected segments more easily than dip-slip events, and thus could produce longer ruptures than dip-slip events of the same width  相似文献   

13.
Destructive upper crustal earthquakes in Central America are often located between active volcanic centres—a geometric relationship that we study using finite element Coulomb failure stress (CFS) models that incorporate the rheologically heterogeneous nature of the volcanic arc. Volcanoes are simulated as mechanically weak zones within a stronger crust. We find that deformation of the volcanic centres within a regional stress field dominated by dextral shear causes stress increases in surrounding crust, with a maximum CFS change between neighbouring volcanoes. This increase in CFS enhances the probability of fault slip on arc-normal faults that are located between volcanic centres; for example, the Tiscapa fault, which ruptured during the 1972 December 13,   M s   6.2 Managua earthquake. The amount of stress increase due to long-term (100 yr) volcano shearing is on the order of 0.1–0.6 bars, similar to values estimated for subduction zone earthquakes.  相似文献   

14.
We perform a systematic parameter space study of the seismic response of a large fault with different levels of heterogeneity, using a 3-D elastic framework within the continuum limit. The fault is governed by rate-and-state friction and simulations are performed for model realizations with frictional and large scale properties characterized by different ranges of size scales. We use a number of seismicity and stress functions to characterize different types of seismic responses and test the correlation between hypocenter locations and the employed distributions of model parameters. The simulated hypocenters are found to correlate significantly with small L values of the rate-and-state friction. The final sizes of earthquakes are correlated with physical properties at their nucleation sites. The obtained stacked scaling relations are overall self-similar and have good correspondence with properties of natural earthquakes.  相似文献   

15.
Summary. The 1973 Hawaii earthquake occurred north of Hilo, at a depth of 40 to 50km. The location was beneath the east flank of Mauna Kea, a volcano dormant historically, but active within the last 4000 yr. Aftershocks were restricted to a depth of 55–35km. The event and its aftershock sequence are located in an area not normally associated with the seismicity of the Mauna Loa and Kilauea calderas. The earthquake was a double event, the epicentres trending NE-SW. The events were of similar size and faulting mechanism. The fault plane solutions obtained by seismic waveform analysis are a strike-slip fault striking EW and dipping 55° S, the auxiliary plane a NS vertical plane with a faulting plunge of 35°. The axis of maximum compressive stress is aligned with the direction of the gravity gradient associated with the island of Hawaii. The fault plane striking EW parallels a surface feature, the Mauna Kea east rift zone. The earthquakes were clearly not associated with volcanic activity normally associated with Mauna Loa and Kilauea and may indicate a deep seated prelude to a resumption of activity at Mauna Kea.  相似文献   

16.
Two distinct phases are commonly observed at the initial part of seismograms of large shallow earthquakes: low-frequency and low-amplitude waves following the onset of a P wave ( P 1) are interrupted by the arrival of the second impulsive phase P2 enriched with high-frequency components. This observation suggests that a large shallow earthquake involves two qualitatively different stages of rupture at its nucleation.
We propose a theoretical model that can naturally explain the above nucleation behaviour. The model is 2-D and the deformation is assumed to be anti-plane. A key clement in our model is the assumption of a zone in which numbers of pre-existing cracks are densely distributed; this cracked zone is a model for the fault zone. Dynamic crack growth nucleated in such a zone is intensely affected by the crack interactions, which exert two conflicting effects: one tends to accelerate the crack growth, and the other tends to decelerate it. The accelerating and decelerating effects are generally ascribable to coplanar and non-coplanar crack interactions, respectively. We rigorously treat the multiple interactions among the cracks, using the boundary integral equation method (BIEM), and assume the critical stress fracture criterion for the analysis of spontaneous crack propagation.
Our analysis shows that a dynamic rupture nucleated in the cracked zone begins to grow slowly due to the relative predominance of non-coplanar interactions. This process radiates the P1 phase. If the crack continues to grow, coalescence with adjacent coplanar cracks occurs after a short time. Then, coplanar interactions suddenly begin to prevail and crack growth is accelerated; the P2 phase is emitted in this process. It is interpreted that the two distinct phases appear in the process of the transition from non-coplanar to coplanar interaction predominance.  相似文献   

17.
A hypothesis for the seismogenesis of a double seismic zone   总被引:1,自引:0,他引:1  
The seismogenesis of a double seismic zone, in particular the lower layer of a double seismic zone, has not been adequately explained in the literature. On the basis of seismic data and geothermal structures along three well-studied cross-sections in the Kuril-Kamchatka and Japan subduction zones, we investigate the temperature/pressure conditions associated with seismogenic structures of the double seismic zones. the corresponding T/P loci seem to suggest that earthquakes observed in the lower layer and in the lower part (below approximately 130 ± 20 km) of the top layer of a double seismic zone were caused by metastable phase transition-a mechanism similar to that responsible for deep-focus earthquakes only at lower temperature/pressure conditions. Under this hypothesis, the wedge-shaped configuration of a double seismic zone is interpreted to represent the loci of the kinetic boundary of the phase transition. According to theoretical/experimental studies and the constraints imposed by our observations, a likely candidate for such a phase transition is the metastable Al-rich enstatite decomposing into the assemblage of Al-poor enstatite plus garnet. Earthquakes in the upper part of the top layer were most probably due to conventional mechanisms such as dehydration of subducted materials and/or facies change from basalt to eclogite. That the top layer involves more than one seismogenic mechanism is also implied by the distinct behaviour of seismicity in the vicinity of 130 ± 20 km. Because the presence of deviatoric stress is critical to the reaction rate of a metastable phase transition, it is inferred that single seismic zones are also caused by the same mechanisms, except that the implicit layer of a supposed double seismic zone is missing, due to the insufficient amount of appropriate metastable minerals or to the lack of appropriate deviatoric stresses in the source region.  相似文献   

18.
Summary. Current models of the structure of an active fault zone recognize two important subdivisions – an upper zone, extending to mid-crustal depths, in which processes associated with brittle fracture and friction dominate the fault behaviour, and a lower zone, extending into the mantle, within which stresses may be relieved by ductile flow. Anisotropy directly or indirectly induced by stress might occur throughout the fault zone, especially if caused by some form of stress-induced crack alignment. Dilatancy associated with high stresses is likely to be a very localized phenomenon in the vicinity of high strength regions (asperities), but alignments caused by subcritical crack growth at low stress and strain rate ( extensive-dilatancy anisotropy ) could give rise to anisotropy throughout the fault region.  相似文献   

19.
High-frequency body waves recorded by a temporary seismic array across the surface rupture trace of the 1992 Landers, California, earthquake were used to determine fault-zone structures down to the seismogenic depth. We first developed a technique to use generalized ray theory to compute synthetic seismograms for arbitrarily oriented tabular low-velocity fault-zone models. We then generated synthetic waveform record sections of a linear array across a vertical fault zone. They show that both arrival times and waveforms of P and S waves vary systematically across the fault due to transmissions and reflections from boundaries of the low-velocity fault zone. The waveform characteristics and arrival-time patterns in the record sections allow us to locate the boundaries of the fault zone and to determine its P - and S -wave velocities independently as well as its depth extent. Therefore, the trade-off between the fault-zone width and velocities can be avoided. Applying the method to the Landers waveform data reveals a low-velocity zone with a width of 270–360 m and a 35–60 per cent reduction in P and S velocities relative to the host rock. The analysis suggests that the low-velocity zone extends to a depth of ∼7 km. The western boundary of the low-velocity zone coincides with the observed main surface rupture trace.  相似文献   

20.
Large historical earthquakes in Italy define a prominent gap in the Pollino region of the southern Apennines. Geomorphic and palaeoseismological investigations in this region show that the Castrovillari fault (CF) is a major seismogenic source that could potentially fill the southern part of this gap. The surface expression of the CF is a complex, 10–13 km long set of prominent scarps. Trenches across one scarp indicate that at least four surface-faulting earthquakes have occurred along the CF since Late Pleistocene time, each producing at least 1 m of vertical displacement. The length of the fault and the slip per event suggest M =6.5-7.0 for the palaeoearthquakes. Preliminary radiocarbon dating coupled with historical considerations imply that the most recent of these earthquakes occurred between 380 BC and 1200 AD, and probably soon after 760 AD; no evidence for this event has been found in the historical record. We estimate a minimum recurrence interval of 1170 years and a vertical slip rate of 0.2-0.5 mm yr-1 for the CF, which indicates that the seismic behaviour of this fault is comparable to other major seismogenic faults of the central-southern Apennines. The lack of mention or the mislocation of the most recent event in the historical seismic memory of the Pollino region clearly shows that even in Italy, which has one of the longest historical records of seismicity, a seismic hazard assessment based solely on the historical record may not be completely reliable, and shows that geological investigations are critical for filling possible information gaps.  相似文献   

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

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