Ductile creep and compaction: A mechanism for transiently increasing fluid pressure in mostly sealed fault zones

A simple cyclic process is proposed to explain why major strike-slip fault zones, including the San Andreas, are weak. Field and laboratory studies suggest that the fluid within fault zones is often mostly sealed from that in the surrounding country rock. Ductile creep driven by the difference between fluid pressure and lithostatic pressure within a fault zone leads to compaction that increases fluid pressure. The increased fluid pressure allows frictional failure in earthquakes at shear tractions far below those required when fluid pressure is hydrostatic. The frictional slip associated with earthquakes creates porosity in the fault zone. The cycle adjusts so that no net porosity is created (if the fault zone remains constant width). The fluid pressure within the fault zone reaches long-term dynamic equilibrium with the (hydrostatic) pressure in the country rock. One-dimensional models of this process lead to repeatable and predictable earthquake cycles. However, even modest complexity, such as two parallel fault splays with different pressure histories, will lead to complicated earthquake cycles. Two-dimensional calculations allowed computation of stress and fluid pressure as a function of depth but had complicated behavior with the unacceptable feature that numerical nodes failed one at a time rather than in large earthquakes. A possible way to remove this unphysical feature from the models would be to include a failure law in which the coefficient of friction increases at first with frictional slip, stabilizing the fault, and then decreases with further slip, destabilizing it.     

THE GEOLOGICAL AND ROCK MECHANICAL DISTINCTION EVIDENCE BETWEEN STICK-SLIP AND CREEP IN HOST ROCK SEGMENTS OF FAULT

Paleo-seismic and fault activity are hard to distinguish in host rock areas compared with soft sedimentary segments of fault. However, fault frictional experiments could obtain the conditions of stable and unstable slide, as well as the microstructures of fault gouge, which offer some identification marks between stick-slip and creep of fault. We summarized geological and rock mechanical distinction evidence between stick-slip and creep in host rock segments of fault, and analyzed the physical mechanisms which controlled the behavior of stick-slip and creep. The chemical composition of fault gouge is most important to control stick-slip and creep. Gouge composed by weak minerals, such as clay mineral, has velocity weakening behavior, which causes stable slide of fault. Gouge with rock-forming minerals, such as calcite, quartz, feldspar, pyroxene, has stick-slip behavior under condition of focal depth. To the gouge with same chemical composition, the deformation mechanism controls the frictional slip. It is essential condition to stick slip for brittle fracture companied by dilatation, but creep is controlled by compaction and cataclasis as well as ductile shear with foliation and small fold. However, under fluid conditions, pressure solution which healed the fractures and caused strength recovery of fault, is the original reason of unstable slide, and also resulted in locking of fault with high pore pressure in core of fault zone. Contrast with that, rock-forming minerals altered to phyllosilicates in the gouges by fluid flow through degenerative reaction and hydrolysis reaction, which produced low friction fault and transformations to creep. The creep process progressively developed several wide shear zones including of R, Y, T, P shear plane that comprise gouge zones embedded into wide damage zones, which caused small earthquake distributed along wide fault zones with focal mechanism covered by normal fault, strike-slip fault and reverse fault. However, the stick-slip produced mirror-like slide surfaces with very narrow gouges along R shear plane and Y shear plane, which caused small earthquake distributed along narrow fault zones with single kind of focal mechanism.     

Fault textures in volcanic conduits: evidence for seismic trigger mechanisms during silicic eruptions

It is proposed that fault textures in two dissected rhyolitic conduits in Iceland preserve evidence for shallow seismogenic faulting within rising magma during the emplacement of highly viscous lava flows. Detailed field and petrographic analysis of such textures may shed light on the origin of long-period and hybrid volcanic earthquakes at active volcanoes. There is evidence at each conduit investigated for multiple seismogenic cycles, each of which involved four distinct evolutionary phases. In phase 1, shear fracture of unrelaxed magma was triggered by shear stress accumulation during viscous flow, forming the angular fracture networks that initiated faulting cycles. Transient pressure gradients were generated as the fractures opened, which led to fluidisation and clastic deposition of fine-grained particles that were derived from the fracture walls by abrasion. Fracture networks then progressively coalesced and rotated during subsequent slip (phase 2), developing into cataclasite zones with evidence for multiple localised slip events, fluidisation and grain size reduction. Phase 2 textures closely resemble those formed on seismogenic tectonic faults characterised by friction-controlled stick-slip behaviour. Increasing cohesion of cataclasites then led to aseismic, distributed ductile deformation (phase 3) and generated deformed cataclasite zones, which are enriched in metallic oxide microlites and resemble glassy pseudotachylite. Continued annealing and deformation eventually erased all structures in the cataclasite and formed microlite-rich flow bands in obsidian (phase 4). Overall, the mixed brittle–ductile textures formed in the magma appear similar to those formed in lower crustal rocks close to the brittle–ductile transition, with the rheological response mediated by strain-rate variations and frictional heating. Fault processes in highly viscous magma are compared with those elsewhere in the crust, and this comparison is used to appraise existing models of volcano seismic activity. Based on the textures observed, it is suggested that patterns of long-period and hybrid earthquakes at silicic lava domes reflect friction-controlled stick-slip movement and eventual healing of fault zones in magma, which are an accelerated and smaller-scale analogue of tectonic faults.Editorial responsibility: J. Stix     

Laboratory Generated M -6 Earthquakes

We consider whether mm-scale earthquake-like seismic events generated in laboratory experiments are consistent with our understanding of the physics of larger earthquakes. This work focuses on a population of 48 very small shocks that are foreshocks and aftershocks of stick–slip events occurring on a 2.0 m by 0.4 m simulated strike-slip fault cut through a large granite sample. Unlike the larger stick–slip events that rupture the entirety of the simulated fault, the small foreshocks and aftershocks are contained events whose properties are controlled by the rigidity of the surrounding granite blocks rather than characteristics of the experimental apparatus. The large size of the experimental apparatus, high fidelity sensors, rigorous treatment of wave propagation effects, and in situ system calibration separates this study from traditional acoustic emission analyses and allows these sources to be studied with as much rigor as larger natural earthquakes. The tiny events have short (3–6 μs) rise times and are well modeled by simple double couple focal mechanisms that are consistent with left-lateral slip occurring on a mm-scale patch of the precut fault surface. The repeatability of the experiments indicates that they are the result of frictional processes on the simulated fault surface rather than grain crushing or fracture of fresh rock. Our waveform analysis shows no significant differences (other than size) between the M -7 to M -5.5 earthquakes reported here and larger natural earthquakes. Their source characteristics such as stress drop (1–10 MPa) appear to be entirely consistent with earthquake scaling laws derived for larger earthquakes.     

Breaking Up: Comminution Mechanisms in Sheared Simulated Fault Gouge

The microstructural state and evolution of fault gouge has important implications for the mechanical behaviour, and hence the seismic slip potential of faults. We use 3D discrete element (DEM) simulations to investigate the fragmentation processes operating in fault gouge during shear. Our granular fault gouge models consist of aggregate grains, each composed of several thousand spherical particles stuck together with breakable elastic bonds. The aggregate grains are confined between two blocks of solid material and sheared under a given normal stress. During shear, the grains can fragment in a somewhat realistic way leading to an evolution of grain size, grain shape and overall texture. The ‘breaking up’ of the fault gouge is driven by two distinct comminution mechanisms: grain abrasion and grain splitting. The relative importance of the two mechanisms depends on applied normal stress, boundary wall roughness and accumulated shear strain. If normal stress is sufficiently high, grain splitting contributes significantly to comminution, particularly in the initial stages of the simulations. In contrast, grain abrasion is the dominant mechanism operating in simulations carried out at lower normal stress and is also the main fragmentation mechanism during the later stages of all simulations. Rough boundaries promote relatively more grain splitting whereas smooth boundaries favor grain abrasion. Grain splitting (plus accompanying abrasion) appears to be an efficient mechanism for reducing the mean grain size of the gouge debris and leads rapidly to a power law size distribution with an exponent that increases with strain. Grain abrasion (acting alone) is an effective way to generate excess fine grains and leads to a bimodal distribution of grain sizes. We suggest that these two distinct mechanisms would operate at different stages of a fault’s history. The resulting distributions in grain size and grain shape may significantly affect frictional strength and stability. Our results therefore have implications for the earthquake potential of seismically active faults with accumulations of gouge. They may also be relevant to the susceptibility of rockslides since non-cohesive basal shear zones will evolve in a similar way and potentially control the dynamics of the slide.     

红河和曲江断裂带断层泥的特征及其地震地质意义

活断层产物的研究近年来受到国内外地震地质学者的关注。笔者对地质、地震证据充分的活断层——红河断裂带和曲江断裂带的断层泥进行多种手段的分析研究。它们的显微构造、粒度分布,组构特征以及石英碎砾的SEM显微形貌等显示了活断层中的长期蠕滑研磨作用、地震过程的快速剪切运动以及活断层特定条件下的退变质作用和其他化学作用等成因特征。综合分析指出,红河断裂带南段可能是蠕滑运动为主的活断层,其北段和曲江断裂带则为粘滑运动为主的活断层     

Constitutive behavior and stability of frictional sliding of granite

An understanding of the frictional sliding on faults that can lead to earthquakes requires a knowledge of both constitutive behavior of the sliding surfaces and its mechanical interaction with the loading system. We have determined the constitutive parameters for frictional sliding of initially bare surfaces of Westerly granite, using a recently developed high pressure rotary shear apparatus that allows long distances of sliding and therefore a greater assurance of attaining steady state behavior. From experiments conducted at room temperature and normal stresses of 27–84 MPa several important results have been found. (1) A gouge layer 100 to 200 m thick was developed from the initially bare rock surfaces after 18 to 70 mm of sliding. (2) The steady state frictional resistance, attained after about 10 mm of sliding, is proportional to the negative of the logarithm of the sliding velocity. (3) Abrup changes in the velocity of sliding result in initial changes in the frictional resistance, which have the same sign as the velocity change, and are followed by a gradual decay to a new steady state value over a characteristic distance of sliding. This velocity weakening behavior is essentially identical with that found by several previous workers on the same material at lower normal stress. (4) Our results are well described by a two state variable constitutive law. The values of the constitutive parameters are quite similar to those found previously at low normal stress, but the characteristic distance is about an order of magnitude smaller than that found at 10 MPa normal stress with thicker layers of coarser gouge. (5) We have approximated our results with a one state variable constitutive law and compared the results with the predictions of existing nonlinear stability analysis; in addition, we have extended the stability analysis to systems possessing two state variables. With such formulations good agreement is found between the experimentally observed and theoretically predicted transitions between stable and unstable sliding. These results allow a better understanding of the instabilities that lead to earthquakes.     

利用滑动速率、古地震资料估算走滑断裂的粘滑及蠕滑量的方法

断层走滑包括粘滑和蠕滑,前者是伴有强震发生的快速运动,后者则是一种缓慢的无震稳滑运动。这２ 种运动常随时间交替出现,共同构成了断层的基本运动方式。在利用断层滑动速率讨论大震重复率的问题中,人们最关心的问题则是如何从总滑动量中分辨和划分出其中所包含的粘滑和蠕滑量,特别是对蠕滑量的划分,因为它直接影响着大震重复率的正确性。笔者在野外考察的基础上,对昌马活动断裂的位移量进行了分级,确定出昌马断裂带全新世以来的水平滑动量大致可划分为５个级别：３１～４１ ｍ ;２５～３１ ｍ ; １５～２２ ｍ ; ８～１３ ｍ ;１～５．５ ｍ ,同时依据古地震学方法并结合１４Ｃ断代法及断层崖形成年代的数学模拟计算,求得全新世以来在昌马断裂带上共发生５ 次古地震事件。在上述２ 项资料确定的基础上,进一步对昌马断裂带的粘滑及蠕滑量进行了划分,并给出了它们随时间的变化情况     

Evolution of Stress Fields and Faulting in Seismic Zones

—Measurements indicate that stress magnitudes in the crust are normally limited by the frictional equilibrium on pre-existing, optimally oriented faults. Fault zones where these limitations are frequently reached are referred to as seismic zones. Fault zones in the crust concentrate stresses because their material properties are different from those of the host rock. Most fault zones are spatially relatively stable structures, however the associated seismicity in these zones is quite variable in space and time. Here we propose that this variability is attributable to stress-concentration zones that migrate and expand through the fault zone. We suggest that following a large earthquake and the associated stress relaxation, shear stresses of a magnitude sufficient to produce earthquakes occur only in those small parts of the seismic zone that, because of material properties and boundary conditions, encourage concentration of shear stress. During the earthquake cycle, the conditions for seismogenic fault slip migrate from these stress-concentration regions throughout the entire seismic zone. Thus, while the stress-concentration regions continue to produce small slips and small earthquakes throughout the seismic cycle, the conditions for slip and earthquakes are gradually reached in larger parts of, and eventually the whole, seismogenic layer of the seismic zone. Prior to the propagation of an earthquake fracture that gives rise to a large earthquake, the stress conditions in the zone along the whole potential rupture plane must be essentially similar. This follows because if they were not, then, on entering crustal parts where the state of stress was unfavourable to this type of faulting, the fault propagation would be arrested. The proposed necessary homogenisation of the stress field in a seismic zone as a precursor to large earthquakes implies that by monitoring the state of stress in a seismic zone, its large earthquakes may possibly be forecasted. We test the model on data from Iceland and demonstrate that it broadly explains the historical, as well as the current, patterns of seismogenic faulting in the South Iceland Seismic Zone.     

Frictional resistance to accelerating slip

General analytical expressions for the friction stress and state variable, based on a rate and state-dependent constitutive friction law proposed by Dieterich and Ruina, have been obtained as an explicit function of slip rateV or slip timet or slip displacement δ under the assumption that slip accelerationa is constant or piecewise constant. Properties of the solutions have been discussed, and reviewed, for uniformly accelerating (or decelerating) slip, the following.

1. Frictional stress increases (or decreases) with increasing time, or slip rate, or slippage at the beginning of motion, until a maximum (or a minimum) value (when it exists) has been reached, then decreases (or increases), and finally approaches a special frictional state, namely a steady state, for which stress depends on instantaneous slip rate.
2. The maximal value of frictional stress is dependent on accelerationa; the larger thea, the larger the magnitude of the maximum.

The energy expense ε needed to overcome the frictional stress has been estimated. And therefore the optimal value of acceleration for saving energy for a slipping mechanical system has been obtained. The energy release rateG for an abrupt rupture process of a crack or fault has been estimated.     

Intermediate-term,pre-earthquake phenomena in California, 1975–1986, and preliminary forecast of seismicity for the next decade

Intermediate-term observations preceding earthquakes of magnitude 5.7 or greater in California from 1975 through 1986 suggest that: (1) The sudden appearance of earthquakes in a previously inactive area indicates an increased likelihood of a significant earthquake in that area for a period from days to years; (2) these larger earthquakes tend to occur towards the ends of creeping fault segments; (3) one large earthquake in a region increases the likelihood of a subsequent significant event in the adjacent area; and (4) marginal evidence for the occurrence of a regional deformation event suggests that such events increase the probability of earthquake occurrence throughout the entire area. A common element in many of these observed patterns appears to be the transmission and amplification of tectonic stress changes by the mechanism of fault creep, and suggests that surface fault creep is a sensitive indicator of changes in stress. The preceding critieria are used to construct a preliminary forecast of the likely locations of significant earthquakes over the next decade.     

Principal Slip Zones in Limestone: Microstructural Characterization and Implications for the Seismic Cycle (Tre Monti Fault, Central Apennines, Italy)

Earthquakes in central Italy, and in other areas worldwide, often nucleate within and rupture through carbonates in the upper crust. During individual earthquake ruptures, most fault displacement is thought to be accommodated by thin principal slip zones. This study presents detailed microstructural observations of the slip zones of the seismically active Tre Monti normal fault zone. All of the slip zones cut limestone, and geological constraints indicate exhumation from <2?km depth, where ambient temperatures are ?100°C. Scanning electron microscope observations suggest that the slip zones are composed of 100% calcite. The slip zones of secondary faults in the damage zone contain protocataclastic and cataclastic fabrics that are cross-cut by systematic fracture networks and stylolite dissolution surfaces. The slip zone of the principal fault has much more microstructural complexity, and contains a 2?C10?mm thick ultracataclasite that lies immediately beneath the principal slip surface. The ultracataclasite itself is internally zoned; 200?C300???m-thick ultracataclastic sub-layers record extreme localization of slip. Syn-tectonic calcite vein networks spatially associated with the sub-layers suggest fluid involvement in faulting. The ultracataclastic sub-layers preserve compelling microstructural evidence of fluidization, and also contain peculiar rounded grains consisting of a central (often angular) clast wrapped by a laminated outer cortex of ultra-fine-grained calcite. These ??clast-cortex grains?? closely resemble those produced during layer fluidization in other settings, including the basal detachments of catastrophic landslides and saturated high-velocity friction experiments on clay-bearing gouges. An overprinting foliation is present in the slip zone of the principal fault, and electron backscatter diffraction analyses indicate the presence of a weak calcite crystallographic preferred orientation (CPO) in the fine-grained matrix. The calcite c-axes are systematically inclined in the direction of shear. We suggest that fluidization of ultracataclastic sub-layers and formation of clast-cortex grains within the principal slip zone occurred at high strain rates during propagation of seismic ruptures whereas development of an overprinting CPO occurred by intergranular pressure solution during post-seismic creep. Further work is required to document the range of microstructures in localized slip zones that cross-cut different lithologies, and to compare natural slip zone microstructures with those produced in controlled deformation experiments.     

Nonuniform friction as a physical basis for earthquake mechanics

A review of simple models and observations suggests that the main first-order features of active faulting-mechanical instability, the frequency-magnitude relations, seismic and aseismie slip, seismic radiation, incoherency and rupture stoppage — may be explained by a single characteristic of crustal faults: the spatial variation of the effective frictional stress, which resists slippage on faults. Faultoffset data suggest that rupture propagation ceases in regions of high resistance which act, as barriers. In these regions slippage is associated with negative stress drop. The spacing and the amplitudeA() of the barriers, as inferred from the frequency-magnitude and moment relation for earthquakes, obeys a simple statistical relationA()^p. On the scale of particle motion, this variability of frictional stress provides a mechanical instability which may be associated with the concept of dynamic friction. Invariably, the rapid particle motion in the model is always preceded by accelerated creep. The particle acceleration is highly irregular, giving rise to an almost random acceleration record on the fault. The particle displacement is relatively smooth, giving rise to simple displacement time function in the far field. Rupture propagation time is approximately proportional to the gradient of frictional stress along the fault. Consequently sharp changes of this stress may cause multiple events and other long period irregularities in the fault motion.The power density spectrum associated with the frictional stress implies that stress may be related to a Poisson distribution of lengths. The autocorrelation of such type of distribution yields a correlation lengthk _L ^–1 , similar perhaps toHaskell's (1964) andAki's (1967) correlation lengths inferred from spectral analysis of seismic waves. The partial incoherency of faulting implies that preseismic deformation may be significantly incoherent, consequently the prediction of small moderate earthquakes may be subject to inherent uncertainties. We conclude that frictional stress heterogeneities may be necessary and sufficient to explain active faulting associated with small and moderate earthquakes.     

Earthquake mechanism and predictability shown by a laboratory fault

Slip events generated in a laboratory fault model consisting of a circulinear chain of eight spring-connected blocks of approximately equal weight elastically driven to slide on a frictional surface are studied. It is found that most of the input strain energy is released by a relatively few large events, which are approximately time predictable. A large event tends to roughen stress distribution along the fault, whereas the subsequent smaller events tend to smooth the stress distribution and prepare a condition of simultaneous criticality for the occurrence of the next large event. The frequency-size distribution resembles the Gutenberg-Richter relation for earthquakes, except for a falloff for the largest events due to the finite energy-storage capacity of the fault system. Slip distributions, in different events are commonly dissimilar. Stress drop, slip velocity, and rupture velocity all tend to increase with event size. Rupture-initiation locations are usually not close to the maximum-slip locations.     

Fault rock analysis of the northern part of the Chelungpu Fault and its relation to earthquake faulting of the 1999 Chi-Chi earthquake, Taiwan

Abstract   The 1999 Chi-Chi earthquake in Taiwan ( M _w = 7.6) produced a surface rupture along the north–south-striking Chelungpu thrust fault with pure dip-slip (east side up) and left lateral strike-slip displacements. Near-field strong-motion data for the northern part of the fault illustrate a distinct lack of the high-frequency seismic radiation associated with a large slip (10–15 m) and a rapid slip velocity (2–4 m/s), suggesting a smooth seismic slip associated with low dynamic frictional resistance on the fault. A drillhole was constructed at shallow depths in the possible fault zones of the northern part of the Chelungpu Fault, which may have slipped during the 1999 earthquake. One of the zones consists of a 20-cm-thick, unconsolidated fault breccia with a chaotic texture lacking both discrete slip surfaces (e.g. Riedel shears) and grain crushing. Other possible fault zones are marked by the narrow (less than a few centimeters) gouge zone in which clayey material intrudes into the damaged zone outside of the gouge zone. These characteristic fault rock textures suggest that the slip mechanisms at shallow levels during the earthquake involved either granular flow of initially unconsolidated material or slip localization under elevated pore pressure along the narrow clayey gouge zone. Because both mechanisms lead to low dynamic frictional resistance on the fault, the rapid seismic slip in the deep portions of the fault (i.e. the source region of strong-motion radiation) could have been accommodated by frictionless slip on the shallow portions of the fault. The combination of strong-motion data and fault rock analysis suggests that smooth slip associated with low dynamic friction occurred on both the deep and shallow portions of the fault, resulting in a large slip between the source region and the surface in the northern region.     

Rupture process of large earthquakes in the northern Mexico subduction zone

The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M _s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M _s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M _s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M _s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M _s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.     

The Fluid Dynamics of Solid Mechanical Shear Zones

Shear zones in outcrops and core drillings on active faults commonly reveal two scales of localization, with centimeter to tens of meters thick deformation zones embedding much narrower zones of mm-scale to cm-scale. The narrow zones are often attributed to some form of fast instability such as earthquakes or slow slip events. Surprisingly, the double localisation phenomenon seem to be independent of the mode of failure, as it is observed in brittle cataclastic fault zones as well as ductile mylonitic shear zones. In both, a very thin layer of chemically altered, ultra fine grained ultracataclasite or ultramylonite is noted. We present an extension to the classical solid mechanical theory where both length scales emerge as part of the same evolutionary process of shearing the host rock. We highlight the important role of any type of solid-fluid phase transitions that govern the second degree localisation process in the core of the shear zone. In both brittle and ductile shear zones, chemistry stops the localisation process caused by a multiphysics feedback loop leading to an unstable slip. The microstructural evolutionary processes govern the time-scale of the transition between slow background shear and fast, intermittent instabilities in the fault zone core. The fast cataclastic fragmentation processes are limiting the rates of forming the ultracataclasites in the brittle domain, while the slow dynamic recrystallisation prolongs the transition to ultramylonites into a slow slip instability in the ductile realm.     

Retardations in fault creep rates before local moderate earthquakes along the San Andreas fault system,central California

Records of shallow aseismic slip (fault creep) obtained along parts of the San Andreas and Calaveras faults in central California demonstrate that significant changes in creep rates often have been associated with local moderate earthquakes. An immediate postearthquake increase followed by gradual, long-term decay back to a previous background rate is generally the most obvious earthquake effect on fault creep. This phenomenon, identified as aseismic afterslip, usually is characterized by above-average creep rates for several months to a few years. In several cases, minor step-like movements, called coseismic slip events, have occurred at or near the times of mainshocks. One extreme case of coseismic slip, recorded at Cienega Winery on the San Andreas fault 17.5 km southeast of San Juan Bautista, consisted of 11 mm of sudden displacement coincident with earthquakes ofM _L =5.3 andM _L =5.2 that occurred 2.5 minutes apart on 9 April 1961. At least one of these shocks originated on the main fault beneath the winery. Creep activity subsequently stopped at the winery for 19 months, then gradually returned to a nearly steady rate slightly below the previous long-term average.The phenomena mentioned above can be explained in terms of simple models consisting of relatively weak material along shallow reaches of the fault responding to changes in load imposed by sudden slip within the underlying seismogenic zone. In addition to coseismic slip and afterslip phenomena, however, pre-earthquakeretardations in creep rates also have been observed. Onsets of significant, persistent decreases in creep rates have occurred at several sites 12 months or more before the times of moderate earthquakes. A 44-month retardation before the 1979M _L =5.9 Coyote Lake earthquake on the Calaveras fault was recorded at the Shore Road creepmeter site 10 km northwest of Hollister. Creep retardation on the San Andreas fault near San Juan Bautista has been evident in records from one creepmeter site for the past 5 years. Retardations with durations of 21 and 19 months also occurred at Shore Road before the 1974 and 1984 earthquakes ofM _L =5.2 andM _L =6.2, respectively.Although creep retardation remains poorly understood, several possible explanations have been discussed previously. (1) Certain onsets of apparent creep retardation may be explained as abrupt terminations of afterslip generated from previous moderate-mainshock sequences. (2) Retardations may be related to significant decreases in the rate of seismic and/or aseismic slip occurring within or beneath the underlying seismogenic zone. Such decreases may be caused by changes in local conditions related to growth of asperities, strain hardening, or dilatancy, or perhaps by passage of stress-waves or other fluctuations in driving stresses. (3) Finally, creep rates may be lowered (or increased) by stresses imposed on the fault by seismic or aseismic slip on neighboring faults. In addition to causing creep-rate increases or retardations, such fault interactions occasionally may trigger earthquakes.Regardless of the actual mechanisms involved and the current lack of understanding of creep retardation, it appears that shallow fault creep is sensitive to local and regional effects that promote or accompany intermediate-term preparation stages leading to moderate earthquakes. A strategy for more complete monitoring of fault creep, wherever it is known to occur, therefore should be assigned a higher priority in our continuing efforts to test various hypotheses concerning the mechanical relations between seismic and aseismic slip.     

Propagation of slip along frictional surfaces

The propagation of slip along a pre-existing frictional plane is formulated for the faults containing interstitial fluid. When normal and shear stresses satisfy the effective frictional law, a frictional sliding occurs stably or unstably, depending on the inhomogeneities of the surface. For a slip to rapidly sweep the surface, the pore pressure of fluid must exceed the critical value that is related to the physical or geometrical irregularities of the surface. If this condition fails, stable sliding is expected, analogous to a seismic fault creep. This prediction makes the role of water in seismic faulting more clear.     

A discussion on corioli force effect and aftershock activity tendency of the <Emphasis Type="Italic">M</Emphasis>=8.1 Kunlun Mountain Pass earthquake on Nov. 14, 2001

Following the theory and definition of the Corioli force in physics, the Corioli force at the site of the M=8.1 Kunlun Mountain Pass earthquake on November 14, 2001, is examined in this paper on the basis of a statistical research on relationship between the Corioli force effect and the maximum aftershock magnitude of 20 earthquakes with M≥7.5 in Chinese mainland, and then the variation tendency of aftershock activity of the M=8.1 earthquake is discussed. The result shows: a) Analyzing the Corioli force effect is an effective method to predict maximum aftershock magnitude of large earthquakes in Chinese mainland. For the sinistral slip fault and the reverse fault with its hanging wall moving toward the right side of the cross-focus meridian plane, their Corioli force pulls the two fault walls apart, decreasing frictional resistance on fault plane during the fault movement and releasing elastic energy of the mainshock fully, so the maximum magnitude of aftershocks would be low. For the dextral slip fault, its Corioli force presses the two walls against each other and increases the frictional resistance on fault plane, prohibiting energy release of the mainshock, so the maximum magnitude of aftershocks would be high. b) The fault of the M=8.1 Kunlun Mountain earthquake on Nov. 14, 2001 is essentially a sinistral strike-slip fault, and the Corioli force pulled the two fault walls apart. Magnitude of the induced stress is about 0.06 MPa. After a comparison analysis, we suggest that the aftershock activity level will not be high in the late period of this earthquake sequence, and the maximum magnitude of the whole aftershocks sequence is estimated to be about 6.0.