Effect of frictional heating and thermal advection on pre-seismic sliding: a numerical simulation using a rate-, state- and temperature-dependent friction law

Laboratory experiments on simulated faults in rocks clearly show the temperature dependence of dynamic rock friction. Since rocks surrounding faults are permeable, we have developed a numerical method to describe the thermo-mechanical evolution of the pre-seismic sliding phase which takes into account both the rate-, state- and temperature-dependent friction law and the heat advection term in the energy equation. We consider a laminar fluid motion perpendicular to a vertical fault plane and assume that fluids move away from the fault plane. A semi-analytical temperature solution which accounts for the variability of slip velocity and stress on the fault has been found. This solution has been generalized to the case of a time varying fluid velocity and then was used to include the thermal pressurization effect. After discretizing the temperature solution, the evolution of the system is obtained by the solution of a system of first order differential equations which allows us to determine the evolution of slip, slip rate, friction coefficient, effective normal stress, temperature and fluid velocity. The numerical solutions are found using a Runge-Kutta method with an adaptative stepsize control in time. When the thermal pressurization effects can be neglected, the heat advection effect gives rise to a delay, with respect to the purely conductive case, of the earthquake occurrence time. This delay increases with increasing permeability H of the system. When the thermal pressurization effects are taken into account the situation is opposite, i.e. the onset of instability tends to precede that of the purely conductive case. The advance in the time of occurrence of instability increases with increasing coefficient of thermal pressurization. In the small permeability range (H ≤ 10^?18 m²), the seismic moment and nucleation length of the pre-seismic phase are significantly smaller than those predicted by the purely conductive model.     

Three‐dimensional texture of natural pseudotachylyte: Pseudotachylyte formation mechanism in hydrous accretionary complex

Melt‐origin pseudotachylyte is the most reliable seismogenic fault rock. It is commonly believed that pseudotachylyte generation is rare in the plate subduction zone where interstitial fluids are abundant and can trigger dynamic fault‐weakening mechanisms such as thermal pressurization. Some recent studies, however, have discovered pseudotachylyte‐bearing faults in exhumed ancient accretionary complexes, indicating that frictional melting also occurrs during earthquakes in subduction zones. To clarify the pseudotachylyte generation mechanism and the variation of slip behavior in the plate subduction zone, a pseudotachylyte found in the exhumed fossil accretionary complex (the Shimanto Belt, Nobeoka, Japan) was re‐focused and microscopic and three‐dimensional observations of the pseudotachylyte‐bearing fault were performed based on optical, electron, and X‐ray microscope images. Based on the patterns contained in the fragment, the pseudotachylyte is divided into four domains, although no clear domain boundaries or layering structures are not found. Three‐dimensional observation also suggests that the pseudotachylyte were fragmented or isolated by cataclasite or carbonate breccia. The pseudotachylyte was rather injected into the surrounding carbonate breccia, which is composed of angular fragments of the host rock and a matrix of tiny crystalline carbonate. The pseudotachylyte volume was extracted from the X‐ray microscope image and the heat abundance consumed by the pseudotachylyte generation was estimated at 2.18 MJ/m², which can be supplied during a slip of approximately 0.5 m. These observations and calculations, together with the results of the previous investigations, suggest hydrofracturing and rapid carbonate precipitation that preceded or accompanied the frictional melting. Dynamic hydrofracturing during a slip can be caused by rapid fluid pressurization, and can induce abrupt decrease in fluid pressure while drastically enhancing the shear strength of the shear zone. Consequently, frictional heating would be reactivated and generate the pseudotachylyte. These deformation processes can explain pseudotachylyte generation in hydrous faults with the impermeable wall rock.     

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.     

Fluid involvement in normal faulting

Evidence of fluid interaction with normal faults comes from their varied role as flow barriers or conduits in hydrocarbon basins and as hosting structures for hydrothermal mineralisation, and from fault-rock assemblages in exhumed footwalls of steep active normal faults and metamorphic core complexes. These last suggest involvement of predominantly aqueous fluids over a broad depth range, with implications for fault shear resistance and the mechanics of normal fault reactivation. A general downwards progression in fault rock assemblages (high-level breccia-gouge (often clay-rich) → cataclasites → phyllonites → mylonite → mylonitic gneiss with the onset of greenschist phyllonites occurring near the base of the seismogenic crust) is inferred for normal fault zones developed in quartzo-feldspathic continental crust. Fluid inclusion studies in hydrothermal veining from some footwall assemblages suggest a transition from hydrostatic to suprahydrostatic fluid pressures over the depth range 3–5 km, with some evidence for near-lithostatic to hydrostatic pressure cycling towards the base of the seismogenic zone in the phyllonitic assemblages. Development of fault-fracture meshes through mixed-mode brittle failure in rock-masses with strong competence layering is promoted by low effective stress in the absence of thoroughgoing cohesionless faults that are favourably oriented for reactivation. Meshes may develop around normal faults in the near-surface under hydrostatic fluid pressures to depths determined by rock tensile strength, and at greater depths in overpressured portions of normal fault zones and at stress heterogeneities, especially dilational jogs. Overpressures localised within developing normal fault zones also determine the extent to which they may reutilise existing discontinuities (for example, low-angle thrust faults). Brittle failure mode plots demonstrate that reactivation of existing low-angle faults under vertical σ₁ trajectories is only likely if fluid overpressures are localised within the fault zone and the surrounding rock retains significant tensile strength. Migrating pore fluids interact both statically and dynamically with normal faults. Static effects include consideration of the relative permeability of the faults with respect to the country rock, and juxtaposition effects which determine whether a fault is transmissive to flow or acts as an impermeable barrier. Strong directional permeability is expected in the subhorizontal σ₂ direction parallel to intersections between minor faults, extension fractures, and stylolites. Three dynamic mechanisms tied to the seismic stress cycle may contribute to fluid redistribution: (i) cycling of mean stress coupled to shear stress, sometimes leading to postfailure expulsion of fluid from vertical fractures; (ii) suction pump action at dilational fault jogs; and, (iii) fault-valve action when a normal fault transects a seal capping either uniformly overpressured crust or overpressures localised to the immediate vicinity of the fault zone at depth. The combination of σ₂ directional permeability with fluid redistribution from mean stress cycling may lead to hydraulic communication along strike, contributing to the protracted earthquake sequences that characterise normal fault systems.     

龙门山断裂带金河磷矿剖面断层泥的低速至高速摩擦性质研究

本文对龙门山断裂带金河磷矿浅钻岩芯中的三种断层泥开展了低速到高速摩擦滑动的实验研究,并对实验变形样品开展了BET比表面积研究.摩擦实验在干燥和孔隙水压条件下开展,速率范围涵盖20 μm·s^-1~1.4 m·s^-1.实验结果显示,三种断层泥在干燥条件下的摩擦性质差别不大,但在孔隙水压条件下,三者的中低速摩擦强度与层状硅酸盐矿物的种类而非总含量紧密相关,蒙脱石和伊利石相比绿泥石更能有效地弱化断层.三种断层泥在孔隙水压条件下存在中低速率域的速度强化,暗示着对断层的加速滑动存在一定的阻碍作用.孔隙水压下,黄绿色和灰绿色断层泥的初始动态弱化非常迅速并伴随断层泥层的瞬时扩容,凹凸体急剧加热导致的局部热压作用可能是造成这种力学行为的物理机制.在经历高速滑动之后,三种断层泥在干、湿条件下的BET比表面积都显著降低,暗示着可能发生了颗粒烧结.中低速域内,孔隙水的存在使得断层泥呈现分散式的剪切变形,BET比表面积的增加因此比干燥条件下更加明显.对表面能的估算表明,颗粒磨碎所消耗的能量至多不超过摩擦力做功的8%,暗示着断层作用中颗粒磨碎所占的能量比例较低.     

Role of fluids in surface deformation caused by the 1999 Chi‐Chi earthquake in Taiwan

The 1999 Chi‐Chi earthquake significantly altered the landscape of central Taiwan. Surface deformation produced by the earthquake along the trace of the Chelungpu thrust can be classified into two styles: (1) uplift without significant surface rupture, and (2) uplift accompanied by surface rupture. Here we examine areas that exhibited the first style of deformation (e.g. Wufeng). Seismic stress at the time of the main shock may have been relieved by high pore‐fluid pressure in a 300‐m‐thick sand and gravel aquifer. Along the thrust fault, frictional heating of these sediments resulted in thermal expansion and an increase in pore‐fluid pressure. High pore‐fluid pressure damped seismic‐wave energy and enhanced intergranular slips of unconsolidated sandy and gravel sediments, which were possibly assisted by sulphuric acid corrosion, leading to a high sulphate content in the groundwater (c. 70 mg L^?1). These changes permitted surface folding and terrace‐style uplifting to occur without significant rupture. In contrast, other areas in which the second style of deformation is dominant (e.g. Fengyuen‐Shihkang) have thin (0–10 m) sand and gravel deposits and lower concentrations of sulphate (c. 30 mg L^?1) in groundwater. In these areas, sediments were heated but not sufficiently to produce significant thermal expansion and increase in pore‐fluid pressure; accumulation of stress in these locations led to rupture at the ground surface, with the formation of steep fault scarps. The areas exhibiting the first deformation style are characterized by the presence of high pore‐fluid pressure, frictional heat conduction, and possibly chemical corrosion related to sulphuric acid attack and formation of sulphate, in contrast to those involving significant uplift and surface rupture. The areal distribution of these two surface deformation styles suggests that the aforementioned fluid‐related subsurface processes may have altered the characteristics of sediments and caused diverse responses to the quake. Copyright © 2002 John Wiley & Sons, Ltd.     

汶川地震断层带结构及渗透率

对汶川地震断层带进行了跨断层的渗透率测量.结果显示汶川地震断层由低渗的核部(2.4×10-19～3.8×10-16m2)、高渗的破碎带(3.7×10-16～3.0×10-15m2)以及含裂隙原岩(6.0×10-18～4.3×10-13 m2)组成(有效压力40 MPa),其中新鲜断层泥具有最低的渗透率.断层泥和两侧原岩...     

Earthquake Rupturing in Fluid-Overpressured Crust: How Common?

Whether or not ruptures nucleate in fluid-overpressured crust (λ _v = P _f/σ _v > 0.4) is important because pore-fluids overpressured above hydrostatic lower fault frictional strength and may also vary through the earthquake cycle, acting as an independent variable affecting fault failure. Containment of fluid overpressure is precarious because pressure-dependent activation of faults and fractures allows drainage from overpressured portions of the crust. Discharge of fluids through activated fault-fracture permeability (fault-valve action) decreases overpressure so that subsequent failure depends on the cycling of both overpressure and frictional strength as well as tectonic stress. Geometric and mechanical considerations suggest that fluid overpressures are more likely to develop and be sustained in compressional/transpressional regimes as opposed to extensional/transtensional tectonic settings. On the basis of geophysical observations and force-balance analyses, subduction interface shear zones appear to be strongly but variably overpressured to near-lithostatic levels (λ _v > 0.9) over the full depth range of seismogenic megathrusts. Strong overpressuring at seismogenic depths is also documented in active fold-thrust belts and in areas of ongoing compressional inversion (e.g., northern Honshu) where inherited normal faults are reactivated as steep reverse faults, requiring near-lithostatic overpressures (λ _v → 1.0) at depths of rupture initiation. Evidence for overpressuring around strike-slip faults is less clear but tends to be strongest in areas of transpression. In areas of extensional tectonics coincident with particularly high fluid discharge, there is some evidence of overpressuring concentrated towards the base of the seismogenic zone. In general, because of the limited resolution of geophysical techniques, it is easier to make the case for rupture propagation through overpressured crust than to make a definitive case for the direct involvement of overpressured fluids in rupture nucleation, though in some instances the circumstantial evidence is compelling. An unresolved related issue is the heterogeneity of overpressuring. Do the active fault zones themselves serve as fluid conduits that are locally overpressured with respect to the surrounding crust?     

The effect of water on stress relaxation of faulted and unfaulted sandstone

A series of stress relaxation experiments have been carried out on faulted and intact Tennessee sandstone to explore the influence of pore water on strength at different strain rates. Temperatures employed were 20, 300 and 400°C, effective confining pressure was 1.5 kb and strain rates as low as 10^–10 sec^–1 were achieved. Most samples were prefaulted at 2.5 kb confining pressure and room temperature. This is thought to have secured a reproducible initial microstructure.The strength of the dry rock was almost totally insensitive to strain rate in the range 10^–4 to 10^–10 sec^–1. In contrast, the strength of the wet rock decreased rapidly with strain rate at rates less than 10^–6 sec^–1. Brittle fracture of the quartz grains which constitute this rock is the most characteristic mode of failure under the test conditions used.The experimental data are discussed in terms of the possible deformation rate controlling processes, and it is suggested that in the wet experiments at intermediate to high strain rates (10^–7 to 10^–4 sec^–1) the observed deformation rate is controlled by the kinetics of water assisted stress corrosion, whilst deformation at low strain rates (ca. 10^–9 sec^–1) is controlled by a pressure solution process.The results have implications for the rheology of fault rocks at depths of perhaps 10 to 15 km in sialic crust.     

Fluid flow from pore pressure measurements off La Palma, Canary Islands

In situ subseafloor pore pressure results from the western flank of the island of La Palma, Canary Islands, are presented. The data obtained with a Pop Up Pore Pressure Instrument (PUPPI) provide constraints on the fluid circulation and its causes in a very special context: The sediment piles near an intraplate oceanic island built on the continental rise of the Northwest African Margin. The ambient pore pressures estimated from 2 to 4 days long record are negative in almost all cases with values, at depths of a few meters below sea floor, usually on the order of −10 to −70 Pa. Excess pore pressures develop only in the distal most areas. The permeabilities and compressibilities obtained respectively from the decay of the insertion pressures and the amplitude of the tidally induced pore pressure variations range between 2.5×10⁻¹⁸ and 6.6×10⁻¹⁶ m² and, 6.2×10⁻⁹ and 1.5×10⁻⁷ Pa⁻¹. According to these permeabilities fluid flow is estimated to be mostly downward and usually on the range between 0 and −0.3 mm y⁻¹. However, from the excess pore pressure profile a complex pattern of fluid circulation is inferred where horizontal fluid motion cannot be neglected. Horizontal flow is probably controlled by significant contrasts in the permeability of the different layers. The prevailing downward fluid flow is abnormal for a classical passive margin. We thus interpret these results as the superposition to the loss of fluids by sediment compaction (on the continental rise), of a large-scale flow system stimulated by thermal buoyancy (100 km wide) related to the volcanic activity on the island of La Palma.     

Fracturing and hydrothermal alteration in normal fault zones

Large normal fault zones are characterized by intense fracturing and hydrothermal alteration. Displacement is localized in a slip zone of cataclasite, breccia and phyllonite surrounding corrugated and striated fault surfaces. Slip zone rock grades into fractured, but less comminuted and hydrothermally altered rock in the transition zone, which in turn grades abruptly into the wall rock. Fracturing and fluid flow is episodic, because permeability generated during earthquakes is destroyed by hydrothermal processes during the time between earthquakes.Fracture networks are described by a fracture fabric tensor (F). The permeability tensor (k) is used to estimate fluid transport properties if the trace of F is sufficiently large. Variations in elastic moduli and seismic velocities between fault zone and wall rock are estimated as a function of fracture density (). Fracturing decreases elastic moduli in the transition zone by 50–100% relative to the country rock, and similar or even greater changes presumably occur in the slip zone.P-andS-wave velocity decrease, andV _p /V _s increases in the fault zone relative to the wall rock. Fracture permeability is highly variable, ranging between 10^–13 m² and 10^–19 m² at depths near 10 km. Changes in permeability arise from variations in effective stress and fracture sealing and healing.Hydrothermal alteration of quartzo-feldspathic rock atT>300°C creates mica, chlorite, epidote and alters the quartz content. Alteration changes elastic moduli, but the changes are much less than those caused by fracturing.P-andS-wave velocities also decrease in the hydrothermally altered fault rock relative to the country rock, and there is a slight decrease inV _p /V _s , which partially offsets the increase inV _p /V _s caused by fracturing.Fracturing and hydrothermal alteration affect fault mechanics. Low modulus rock surrounding fault surfaces increases the probability of exceeding the critical slip distance required for the onset of unstable slip during rupture initiation. Boundaries between low modulus fault rock and higher modulus wall rock also act as rupture guides and enhance rupture acceleration to dynamic velocity. Hydrothermal alteration at temperatures in excess of 300°C weakens the deeper parts of the fault zone by producingphyllitic mineral assemblages. Sealing of fracture in time periods between large earthquakes generates pods of abnormally pressured fluid which may play a fundamental role in the initiation of large earthquakes.     

Effect of displacement rate on the real area of contact and temperatures generated during frictional sliding of Tennessee sandstone

Summary The real area of contact has been determined, and measurements of the maximum and average surface temperatures generated during frictional sliding along precut surfaces in Tennessee sand-stone have been made, through the use of thermodyes. Triaxial tests have been made at 50 MPa confining pressure and constant displacement rates of 10^–2 to 10^–6 cm/sec, and displacements up to 0.4 om. At 0.2 cm of stable sliding, the maximum temperature decreases with decreasing nominal displacement rate from between 1150° to 1175°C at 10^–2 cm/sec to between 75° to 115°C at 10^–3 cm/sec. The average temperature of the surface is between 75 and 115°C at 10^–2 cm/sec, but shows no rise from room temperature at 10^–3 cm/sec. At 0.4 cm displacement, and in the stick-slip mode, as the nominal displacement rate decreases from 10^–3 to 10^–6 cm/sec, the maximum temperature decreases from between 1120° to 1150°C to between 1040° to 1065°C. The average surface temperature is 115° to 135°C at displacement rates from 2.6×10^–3 to 10^–4 cm/sec.With a decrease in the displacement rate from 10^–2 to 10^–6 cm/sec, the real area of contact increases from about 5 to 14 percent of the apparent area; the avergge area of asperity contact increases from 2.5 to 7.5×10^–4 cm². Although fracture is the dominate mechanism during stick-up thermal softening and creep may also contribute to the unstable sliding process.     

Pore Creation due to Fault Slip in a Fluid-permeated Fault Zone and its Effect on Seismicity: Generation Mechanism of Earthquake Swarm

—Spatio-temporal variation of rupture activity is modeled assuming fluid migration in a narrow porous fault zone formed along a vertical strike-slip fault in a semi-infinite elastic medium. Pores are assumed to be created in the fault zone by fault slip. The effective stress principle coupled to the Coulomb failure criterion introduces mechanical coupling between fault slip and pore fluid. The fluid is assumed to flow out of a localized high-pressure fluid compartment in the fault with the onset of earthquake rupture. The duration of the earthquake sequence is assumed to be considerably shorter than the recurrence period of characteristic events on the fault. The rupture process is shown to be significantly dependent on the rate of pore creation. If the rate is large enough, a foreshock–mainshock sequence is never observed. When an inhomogeneity is introduced in the spatial distribution of permeability, high complexity is observed in the spatio-temporal variation of rupture activity. For example, frequency-magnitude statistics of intermediate-size events are shown to obey the Gutenberg–Richter relation. Rupture sequences with features of earthquake swarms can be simulated when the rate of pore creation is relatively large. Such sequences generally start and end gradually with no single event dominating in the sequence. In addition, the b values are shown to be unusually large. These are consistent with seismological observations on earthquake swarms.     

Monte Carlo Inversion of DInSAR Data for Dislocation Modeling: Application to the 1997 Umbria-Marche Seismic Sequence (Central Italy)

— The study of surface deformation due to seismic activity is often made using dislocations with uniform slip and simple geometries. A better modeling of coseismic and postseismic surface displacements can be obtained by using dislocations with variable slip and nonregular shapes. This is consistent with the asperity model of fault surfaces, assuming a friction distribution on faults made of locked zones with much higher friction than surrounding zones. In this paper we consider the 1997–1998 Colfiorito seismic sequence. The coseismic surface displacements in the Colfiorito zone are used in order to infer the slip distribution on the fault surface at different stages of the sequence. The displacement field has been modeled varying the slip distribution on the fault, and comparing the deformation observed by SAR and GPS techniques with model results. The slip distribution is calculated by Monte Carlo simulations on a normal fault with the dip angle equal to 40°. A good approximation is obtained by using square asperity units of 1.5×1.5 km². In the first stage, we employed a simplified model with uniform slip, in which each asperity unit is allowed to slip a constant amount or not to slip at all, and in the second stage, we evaluate the slip distribution in the dislocation area determined by the Monte Carlo inversion: in this case we allow unit cells to undergo different values of slip in order to refine the initial dislocation model. The results show that the 1997 seismic events of the sequence can be modeled by irregular dislocations, obtaining a good fit to the DInSAR and GPS observations. The model also confirms the results of previous studies by a different methodology, defining the distribution of asperities on the fault plane using the fault geometry, the geodetic data and the seismic moment of the 1997–1998 Colfiorito seismic sequence. Furthermore, the analysis of 1997 aftershocks in the seismogenic region shows a strong correlation between most events and the asperity distribution, which can be considered as an independent test of the validity of the model.     

Vent fluid chemistry in Bahía Concepción coastal submarine hydrothermal system, Baja California Sur, Mexico

Shallow submarine hydrothermal activity has been observed in the Bahía Concepción bay, located at the Gulf coast of the Baja California Peninsula, along faults probably related to the extensional tectonics of the Gulf of California region. Diffuse and focused venting of hydrothermal water and gas occurs in the intertidal and shallow subtidal areas down to 15 m along a NW–SE-trending onshore–offshore fault. Temperatures in the fluid discharge area vary from 50 °C at the sea bottom up to 87 °C at a depth of 10 cm in the sediments.Chemical analyses revealed that thermal water is enriched in Ca, As, Hg, Mn, Ba, HCO₃, Li, Sr, B, I, Cs, Fe and Si, and it has lower concentrations of Cl, Na, SO₄ and Br than seawater. The chemical characteristics of the water samples indicate the occurrence of mixing between seawater and a thermal end-member. Stable isotopic oxygen and hydrogen composition of thermal samples plot close to the Local Meteoric Water Line on a mixing trend between a thermal end-member and seawater. The composition of the thermal end-member was calculated from the chemistry of the submarine samples data by assuming a negligible amount of Mg for the thermal end-member. The results of the mixing model based on the chemical and isotopic composition indicate a maximum of 40% of the thermal end-member in the submarine vent fluid.Chemical geothermometers (Na/Li, Na–K–Ca and Si) were applied to the thermal end-member concentration and indicate a reservoir temperature of approximately 200 °C. The application of K–Mg and Na/Li geothermometers for vent fluids points to a shallow equilibrium temperature of about 120 °C.Results were integrated in a hydrogeological conceptual model that describes formation of thermal fluids by infiltration and subsequent heating of meteoric water. Vent fluid is generated by further mixing with seawater.     

Fault dimensions,displacements and growth

Maximum total displacement (D) is plotted against fault or thrust width(W) for 65 faults, thrusts, and groups of faults from a variety of geological environments. Displacements range from 0.4 m to 40 km and widths from 150 m to 630 km, and there is a near linear relationship betweenD andW ². The required compatibility strains (e _s) in rocks adjacent to these faults increases linearly withW and with and ranges frome _s=2×10^–4 toe _s=3×10^–1. These are permanent ductile strains, which compare with values ofe _s=2×10^–5 for the elastic strains imposed during single slip earthquake events, which are characterised by a linear relationship between slip (u) andW.The data are consisten with a simple growth model for faults and thrusts, in which the slip in successive events increases by increments of constant size, and which predicts a relationship between displacement and width of the formD=cW ². Incorporation of constant ductile strain rate into the model shows that the repreat time for slip events remains constant throughout the life of a fault, while the displacement rate increases with time. An internally consistent model withe _s=2×10^–5, giving repeat times of 160 years and instantaneous displacement rates of 0.02 cm/yr, 0.2 cm/yr, and 2.0 cm/yr when total displacement is 1 m, 100 m, and 10 km, and slip increasing by 0.5 mm with each event, gives a good approximation of the data. The model is also applicable to stable sliding, the slip rate varying with ductile strain rate and withW ².     

Visco-elastic stress triggering model of Tangshan earthquake sequence

We calculated the Coulomb failure stress change generated by the 1976 Tangshan earthquake that is projected onto the fault planes and slip directions of large subsequent aftershocks.Results of previous studies on the seismic fail-ure distribution,crustal velocity and viscosity structures of the Tangshan earthquake are used as model constraints.Effects of the local pore fluid pressure and impact of soft medium near the fault are also considered.Our result shows that the subsequent Luanxian and Ninghe earthquakes occurred in the regions with a positive Coulomb fail-ure stress produced by the Tangshan earthquake.To study the triggering effect of the Tangshan,Luanxian,and Ninghe earthquakes on the follow-up small earthquakes,we first evaluate the possible focal mechanisms of small earthquakes according to the regional stress field and co-seismic slip distributions derived from previous studies,assuming the amplitude of regional tectonic stress as 10 MPa.By projecting the stress changes generated by the above three earthquakes onto the possible fault planes and slip directions of small earthquakes,we find that the "butterfly" distribution pattern of increased Coulomb failure stress is consistent with the spatial distribution of follow-up earthquakes,and 95% of the aftershocks occurred in regions where Coulomb failure stresses increase,indicating that the former large earthquakes modulated occurrences of follow-up earthquakes in the Tangshan earthquake sequence.This result has some significance in rapid assessment of aftershock hazard after a large earthquake.If detailed failure distribution,seismogenic fault in the focal area and their slip features can be rapidly determined after a large earthquake,our algorithm can be used to predict the locations of large aftershocks.     

Nanoscale porosity in SAFOD core samples (San Andreas Fault)

With transmission electron microscopy (TEM) we observed nanometer-sized pores in four ultracataclastic and fractured core samples recovered from different depths of the main bore hole of the San Andreas Fault Observatory at Depth (SAFOD). Cutting of foils with a focused ion beam technique (FIB) allowed identifying porosity down to the nm scale. Between 40 and 50% of all pores could be identified as in-situ pores without any damage related to sample preparation. The total porosity estimated from TEM micrographs (1–5%) is comparable to the connected fault rock porosity (2.8–6.7%) estimated by pressure-induced injection of mercury. Permeability estimates for cataclastic fault rocks are 10^? 21–10^? 19 m² and 10^? 17 m² for the fractured fault rock. Porosity and permeability are independent of sample depth. TEM images reveal that the porosity is intimately linked to fault rock composition and associated with deformation. The TEM-estimated porosity of the samples increases with increasing clay content. The highest porosity was estimated in the vicinity of an active fault trace. The largest pores with an equivalent radius > 200 nm occur around large quartz and feldspar grains or grain-fragments while the smallest pores (equivalent radius < 50 nm) are typically observed in the extremely fine-grained matrix (grain size < 1 μm). Based on pore morphology we distinguish different pore types varying with fault rock fabric and alteration. The pores were probably filled with formation water and/or hydrothermal fluids at elevated pore fluid pressure, preventing pore collapse. The pore geometry derived from TEM observations and BET (Brunauer, Emmett and Teller) gas adsorption/desorption hysteresis curves indicates pore blocking effects in the fine-grained matrix. Observations of isolated pores in TEM micrographs and high pore body to pore throat ratios inferred from mercury injection suggest elevated pore fluid pressure in the low permeability cataclasites, reducing shear strength of the fault.     

Aseismio fault slip and block deformation in North China

In North China, the tectonic fault-block system enables us to use the Discontinuous Deformation Analysis (DDA) method to simulate the long-term cross-fault survey and other geodetic data related to aseismic tectonic deformation. By the simulation we have found that: (1) Slips on faults with different orientation are generally in agreement with the ENE-WSW tectonic stress field, but the slip pattern of faulting can vary from nearly orthogonal, to pure shear along the strike of the faults, this pattern cannot be explained by simple geometric relation between the strike of the fault and the direction of the tectonic shortening. This phenomenon has been observed at many sites of cross-fault geodetic surveys, and might be caused by the interactions between different blocks and faults. (2) According to the DDA model, if the average aseismic slip rate along major active faults is at the order of several tenths of millimeter per year as observed by the cross-fault geodetic surveys, the typical strain rate inside a block is at the order of 10^–8 year^–1 or less, so that the rate of 10^–6 year^–1, as reported by observations in smaller areas, cannot be the representative deformation rate in this region. (3) Between the slips caused by regional compression and block rotation, there is a possibility that the sense of slip caused by rigid body rotation in two adjacent blocks is opposite to the slip caused by the tectonic compression. But the magnitude of slip resulting from the tectonic compression is much larger than that due to the block rotation. Thus, in general, the slip pattern on faults as a whole agrees with the sense of tectonic compression in this region. That is to say, the slip caused by regional compression dominates the entire slip budget. (4) Based on (3), some observed slips in contradiction to ENE tectonic stress field may be caused by more localized sources, and have no tectonic significance.     

Postglacial fault movement and palaeoseismicity in western Scotland: A reappraisal of the Kinloch Hourn fault, Kintail

The Kinloch Hourn fault is the most prominent of a number of suspectedpostglacial faults in the western Scottish Highlands. These faults areinterpreted to have been reactivated by repeated large (M > 6)palaeoseismic events following deglaciation 10,000–13,000 years ago.Based on inferred deflections of drainage courses, previous studies of thefault have estimated 160 ± 40 m cumulative left-lateral displacementalong a 14 km long active segment during postglacial times. Reportedsoft-sediment deformation phenomena imply that activity on the KinlochHourn fault has persisted into the late Holocene, with the most recentmovement having been associated with a magnitude 5.5–6.0 surface-faultingevent between 3500 and 2400 years ago. The marked contrast betweensuch palaeoseismic activity and the present-day seismic quiescence ofwestern Scotland has stimulated this critical reappraisal of the KinlochHourn fault.This paper reassesses the key lines of evidence for postglacial fault activityand palaeoseismicty on the Kinloch Hourn fault, combining the analysis of1:15,000-scale air photos, field-based geomorphic mapping andpalaeoenvironmental investigations. Our reappraisal of inferred drainagedeflections across the fault contends that previous reports of significant(10² m) left-lateral slip on the fault during the Holocene arespurious. Instead, incidences of Holocene channel abandonment along thefault line are non-synchronous and probably reflect non-tectonic drainagechanges. The timing of soft-sediment deformation in the vicinity of the faultis revised to an early Holocene date (8990–8580 calendar years BP), whichis in accord with both the palaeoenvironmental history of the site andconsistent with published ages of earthquake-induced liquefactionphenomena documented elsewhere in western Scotland. An alleged recent(post-2400 radiocarbon years BP) ground rupture on the fault isquestioned in the light of uncertainty about both the nature of the faultedsoil deposit and the late Holocene age attributed to it.The study concludes that there is no convincing evidence for postglacialsurface rupture on the Kinloch Hourn fault and speculates that the casefor significant (10¹–10² m) postglacial movement on otherfaults in western Scotland may be similarly `unproven'.