The influence of water on the stress supported by experimentally faulted Westerly granite

Summary. Fault zones in wet Westerly granite deformed at temperatures of 300° and 400°C require markedly lower shear stresses for sliding than when dry, even when the effective confining pressure is held constant between the wet and dry tests, provided that the strain rate is lower than 10⁻⁷s⁻¹. The rate of strength reduction is enhanced by increasing the pore water pressure. The deformation rate is a power function of the applied stress where the stress exponent is approximately 7 for pore water pressure of 100 MPa and 21 for pore water pressure of 20 MPa.
The experimental results are extrapolated to conditions believed to occur at depths of 10 km along the San Andreas Fault Zone. It is suggested that for slow tectonic deformation at strain rates of 10⁻¹¹ and 10⁻¹⁴s⁻¹ the shear stress for sliding on faults in granite would be approximately 60 and 20 MPa, respectively, at pore water pressures equal to the hydrostatic head. Fluid overpressures of c. 0.8 lithostatic pressure are required to lower the shear stress for sliding into the 10 MPa range at the slower strain rate.     

Shear-wave anisotropy: spatial and temporal variations in time delays at Parkfield, Central California

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 M_L = 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 M_L = 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.     

Shear wave anisotropy in the crust around the San Andreas fault near Parkfield: spatial and temporal analysis

We systematically analysed shear wave splitting (SWS) for seismic data observed at a temporary array and two permanent networks around the San Andreas Fault (SAF) Observatory at Depth. The purpose was to investigate the spatial distribution of crustal shear wave anisotropy around the SAF in this segment and its temporal behaviour in relation to the occurrence of the 2004 Parkfield M 6.0 earthquake. The dense coverage of the networks, the accurate locations of earthquakes and the high-resolution velocity model provide a unique opportunity to investigate anisotropy in detail around the SAF zone. The results show that the primary fast polarization directions (PDs) in the region including the SAF zone and the northeast side of the fault are NW–SE, nearly parallel or subparallel to the SAF strike. Some measurements on the southwest side of the fault are oriented to the NNE–SSW direction, approximately parallel to the direction of local maximum horizontal compressive stress. There are also a few areas in which the observed fast PDs do not fit into this general pattern. The strong spatial variations in both the measured fast PDs and time delays reveal the extreme complexity of shear wave anisotropy in the area. The top 2–3 km of the crust appears to contribute the most to the observed time delays; however substantial anisotropy could extend to as deep as 7–8 km in the region. The average time delay in the region is about 0.06 s. We also analysed temporal patterns of SWS parameters in a nearly 4-yr period around the 2004 Parkfield main shock based on similar events. The results show that there are no appreciable precursory, coseismic, or post-seismic temporal changes of SWS in a region near the rupture of an M 6.0 earthquake, about 15 km away from its epicentre.     

P-wave velocities in rocks at high pressure and temperature and the constitution of the central California crust

Summary. Velocities of compressional waves are determined for central California rocks at pressures up to 0.7 GPa (7 kb) and temperatures up to 450°C. These data are used to interpret the seismic velocity structure of the crust in the California Coast Ranges. The seismic data on both sides of the San Andreas fault are consistent with the following model; besides some patches of surface sediments the upper 10—15 km of the crust on the northeast side consists predominantly of sedimentary and metasedimentary rocks of the Franciscan assemblage; the lower crust, of a thickness of 15—20 km, may be composed of gabbroic or other mafic rocks. Across the fault on the south-west side, the entire crustal section is probably a granitic complex similar to that exposed on the surface. The proposed model is shown to be consistent with the observed gravity anomaly.     

Temporal changes in shear-wave splitting at an isolated swarm of small earthquakes in 1992 near Dongfang, Hainan Island, southern China

An isolated swarm of small earthquakes occurred in 1992, near Dongfang on Hainan Island, southern China. The Institute of Geophysics, State Seismological Bureau of China, monitored the swarm with five DCS-302 digital accelerometers for three months from 1992 June 1. 18 earthquakes, with magnitudes M _L ranging from 1.8 to 3.6, were well located by five stations, and shear-wave splitting varying azimuthally was analysed on 27 seismic records from these events. The mean polarization azimuth of the faster shear wave was WNW. Time delays between the split shear waves at two stations varied with time and space. The time delays at one station fell abruptly after earthquakes of magnitudes 3.1 and 3.6, but did not change significantly at the second station. This behaviour is consistent with the delay-time changes being caused by changes in the aspect ratio of vertical liquid-filled (EDA) cracks. Thus, the variation in shear-wave-splitting time delay could be due to changes in crustal stress related to nearby small-magnitude earthquake activity. The connection between earthquake activity and crustal stress variation measured by shear-wave splitting leaves the door open for possible observations of crustal stress transients related to the onset of an earthquake; however, our data cannot be considered as definite evidence for such precursors.     

Shear-wave polarizations in the Peter the First Range indicating crack-induced anisotropy in a thrust-fault regime

Summary. Three-component seismograms of small local earthquakes recorded in the Peter the First Range of mountains near Garm, Tadzhikistan SSR, display shear-wave splitting similar to that previously observed near the North Anatolian Fault in Turkey. The Peter the First Range is in a region of compressional tectonics, whereas the North Anatolian Fault is a comparatively simple strike-slip fault. Detailed analysis of the Turkish records suggests that the splitting is diagnostic of crack-induced anisotropy caused by vertical microcracks aligned parallel to the direction of maximum compression. Preliminary examination of paper records from Garm shows that most shear waves arriving within the shear-wave window display shear-wave splitting, and that the polarizations of leading shear-waves are consistently aligned in a NE/SW direction. The area is complicated and the tectonics are not well-understood, but the NE/SW direction is approximately perpendicular to the compressional axis in many of the fault-plane mechanisms of the earthquakes. These earthquakes are usually at depths between 5 and 12 km, although there are some deeper events nearby.
Parallel shear-wave polarizations, such as those observed, are expected to indicate the strike of nearly vertical parallel microcracks, which would be aligned parallel to the direction of maximum compression. Thus the shear-wave polarizations in the Peter the First Range indicate that the directions of principal stress are reversed in the rock above the earthquake foci where thrust faulting is taking place.     

Shallow earthquakes in a viscoelastic shear zone with depth-dependent friction and rheology

Summary. Most crustal earthquakes of the world are observed to occur within a seismogenic layer which extends from the Earth's surface to a depth of a few tens of kilometres at most. A model is proposed in which the shear zone along a transcurrent plate margin is represented as a viscoelastic medium with depth-dependent power-law rheology. A frictional resistance linearly increasing with depth is assumed on a vertical transcurrent fault within the shear zone. Such a model is able to reproduce a continuous transition from the brittle behaviour of the upper crust to the ductile behaviour at depth. Assuming that the shear zone is subjected to a constant strain rate from the opposite motions of the two adjacent plates, it is found that there exists a maximum depth H below which tectonic stress can never reach the frictional threshold: this may be identified as the maximum depth of earthquake nucleation. The value of H is consistent with observations for plausible values of the model parameters. The stress evolution in the shear zone is calculated in the linear approximation of the constitutive equation. A change in rigidity with depth, which is also introduced in the model, may reproduce the high vertical gradient of shear stress, which has been measured across the San Andreas fault, and the fact that most earthquakes are nucleated at some depth in the seismogenic layer. A crack which drops the ambient stress to the dynamic frictional level is then introduced in the model. To this aim, a crack solution is employed without a stress singularity at its edges, which is compatible with a frictional stress threshold criterion for fracture. A constraint on the vertical friction gradient is obtained if such cracks are assumed to be entirely confined within the seismogenic layer.     

Ductile shear zones: some aspects of constant slip velocity and constant shear stress models

Summary The thermomechanical differential equations governing deformation in viscous shear zones have been solved for both constant velocity and constant stress boundary conditions. The solutions show that the inertial term in these equations can be neglected everywhere.
The starting condition of the constant velocity model has been shown to be a constant velocity gradient and not a Heaviside function. The temperature anomaly produced by shear heating at the centre of the shear zone is shown to increase gradually and continuously with time, not reaching an asymptotic value. Conclusions for the constant velocity boundary condition are otherwise generally similar to those presented by Yuen et al , and agree with Fleitout & Froidevaux. The temperatures reached by constant velocity shears are sufficient for partial melting.
Constant stress boundary condition shear zone models show an initially broad shear zone with uniform shear velocity gradient. Depending on the level of applied shear stress and ambient temperature, localized intense shear heating may develop followed by thermal runaway. At lower ambient temperatures relatively high stresses are required to produce thermal runaway.
The broadening of the constant velocity shear zone proceeds more rapidly with increased ambient temperature. This can be used to show that shear zones broaden with depth. The merging of parallel shear zone pairs has been investigated and shear zones separated by distances of less than 10km coalesce to form a single shear zone within 3 Myr. Only shear zones separated by 50km or more remain distinct over periods of tens of millions of years.     

Numerical simulations of tectonic processes in Southern California

Summary. A relaxation model for crustal tectonics, using the Finite Element method, was developed and applied to the region of Southern California. Both the existence of the Transverse Ranges and the geometry of the San Andreas fault in the region of Southern California, imply a stress pattern deviating from the simple horizontal shear which parallels the spreading between the Pacific and the North American plates. A number of possible mechanisms responsible for this anomalous stress field were examined quantitatively in the light of seismicity and other tectonic observations, and in particular to the Palmdale uplift which was reported to have occurred between the years 1960–65.     

Detailed S-wave structure in the Dublin Basin and its northern margin

Summary. The seismic structure has been measured to a depth of about 3 km along a 30 km seismic profile in east central Ireland. This profile is unusual in that it is the S -wave velocity—depth structure that has been measured to a degree of precision more normally associated with P -wave results. One reason for this is that the sources used were quarry blasts which generated strong S -waves and short-period surface waves but rather weak P -waves.
The results show a layer of Carboniferous limestone with shear velocity 2.65 km⁻¹ s overlying a layer with a velocity of 3.06 km s⁻¹. This second layer was interpreted as Lower Palaeozoic strata (Silurian/Ordovician) since this velocity was evident in an inlier seen at the surface at the northern end of the line. A third refraction horizon, shear velocity 3.45 km s⁻¹ and displaying a basinal structure, was also recognized. This may be Cambrian or Precambrian basement.     

Finite-frequency traveltime tomography of San Francisco Bay region crustal velocity structure

Seismic velocity structure of the San Francisco Bay region crust is derived using measurements of finite-frequency traveltimes. A total of 57 801 relative traveltimes are measured by cross-correlation over the frequency range 0.5–1.5 Hz. From these are derived 4862 'summary' traveltimes, which are used to derive 3-D P -wave velocity structure over a 341 × 140 km² area from the surface to 25 km depth. The seismic tomography is based on sensitivity kernels calculated on a spherically symmetric reference model. Robust elements of the derived P -wave velocity structure are: a pronounced velocity contrast across the San Andreas fault in the south Bay region (west side faster); a moderate velocity contrast across the Hayward fault (west side faster); moderately low velocity crust around the Quien Sabe volcanic field and the Sacramento River delta; very low velocity crust around Lake Berryessa. These features are generally explicable with surface rock types being extrapolated to depth ∼10 km in the upper crust. Generally high mid-lower crust velocity and high inferred Poisson's ratio suggest a mafic lower crust.     

Implications of stress concentration on a strike-slip fault in an elastic plate subject to basal shear stress

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.     

Small-scale mantle flows and induced lithospheric stress near island arcs

Summary. This paper explores the middle ground between complex thermally-coupled viscous flow models and simple corner flow models of island arc environments. The calculation retains the density-driven nature of convection and relaxes the geometrical constraints of corner flow, yet still provides semianalytical solutions for velocity and stress. A novel aspect of the procedure is its allowance for a coupled elastic lithosphere on top of a Newtonian viscous mantle. Initially, simple box-like density drivers illustrate how vertical and horizontal forces are transmitted through the mantle and how the lithosphere responds by trench formation. The flexural strength of the lithosphere spatially broadens the surface topography and gravity anomalies relative to the functional form of the vertical flow stresses applied to the plate base. I find that drivers in the form of inclined subducting slabs cannot induce self-driven parallel flow; however, the necessary flow can be provided by supplying a basal drag of 1–5 MPa to the mantle from the oceanic lithosphere. These basal drag forces create regional lithospheric stress and they should be quantifiable through seismic observations of the neutral surface. The existence of a shallow elevated phase transition is suggested in two slab models of 300 km length where a maximum excess density of 0.2 g cm⁻³ was needed to generate an acceptable mantle flow. A North New Hebrides subduction model which satisfies flow requirements and reproduces general features of topography and gravity contains a high shear stress zone (75 MPa) around the upper slab surface to a depth of 150 km and a deviatoric tensional stress in the back arc to a depth of 70 km. The lithospheric stress state of this model suggests that slab detachment is possible through whole plate fracture.     

Amplification of ground motion and waveform complexity in fault zones: examples from the San Andreas and Calaveras Faults

Summary. P -wave seismograms at ranges less than 10 km are synthesized by asymptotic ray theory and by summation of Gaussian beams for point sources located in a low-velocity wedge surrounding a fault. The computations are performed using models of the wedge inferred from the analysis of reflection and refraction experiments across the San Andreas and Hayward-Calaveras faults. Calculations in these models show that the 10–20Hz vertical displacements of earthquakes located at 3–10km depth are amplified by up to an order of magnitude in a 1–2km wide region centred on the fault trace compared to displacements predicted by laterally homogeneous models of the crust. This amplification is not cancelled by high attentuation in the fault zone and compensates for the reduction in amplitudes directly above the source predicted from the radiation pattern of a strike-slip earthquake. Depending on the source depth of the earthquake and the structure and velocity contrast of the wedge, multiple triplications in the travel-time curve of direct P - and S -waves will occur at stations in the fault zone. A wedge model successfully predicts the triplications observed in the P waveforms of aftershocks of the Coyote Lake earthquake recorded in the fault zone, showing that body waves from microearthquakes can be used to determine the three-dimensional velocity structure of the fault zone. The amplification, waveform complexity, and distortion of ray paths introduced by the low- velocity wedge suggest that its effects should be included in the interpretation of strong ground motions and travel times observed in the fault zone. For realistic models of the wedge, asymptotically approximate methods of calculating the body waveforms are strictly valid for frequencies greater than 20Hz. Numerical methods may be necessary to calculate accurately the wavefield at lower frequencies.     

Preliminary Results on the 3-dimensional Seismic Structure of the Lithosphere under the USGS Central California Seismic Array

About 1500 readings of teleseismic P -time residuals obtained from the US Geological Survey seismograph network in central California have been used to obtain a three-dimensional image of seismic velocity anomalies for this area by the method of Aki, Christoffersson & Husebye We found that the California network is less suitable than the LASA and NORSAR arrays for this kind of studies because of its greater proportion of peripheral blocks in which the resolution is very poor for the stochastic inverse solution and the random error effect is severe for the generalized inverse solution. Nevertheless, the resultant velocity anomalies show a remarkable correlation with the San Andreas fault zone to a depth of 75 km. The anomaly pattern changes drastically as the depth exceeds 75 km, suggesting that the asthenosphere has been reached.     

On the stability of frictionally heated shear flows in the asthenosphere

Summary . Frictional heating in upper mantle shear flows may lead to localized thermal runaway and partial melting in the asthenosphere, but only as the result of a finite-amplitude disturbance. A rigorous two-dimensional stability analysis shows that asthenospheric shear flows are stable to small-amplitude perturbations whether such flows are supercritical (shear stress decreases with increasing plate velocity) or subcritical (shear stress increases with increasing plate velocity). Disturbances which maintain a shear stress larger than the critical value for sufficiently long will lead to runaway. The response of the asthenosphere to events which do not satisfy this criterion must be determined by a non-linear analysis. Reasonable models of flow in the asthenosphere could be driven to runaway, at a superexponential growth rate, by sudden increases in shear stress of less than 10 bar. Disturbances resulting from plate collisions may maintain large enough stresses for sufficiently long times to initiate runaways, while stress changes associated with large earthquakes probably occur too rapidly to do so.     

Present‐day stress orientation in the Clarence‐Moreton Basin of New South Wales,Australia: a new high density dataset reveals local stress rotations

Early phases of the Australian Stress Map project revealed that plate boundary forces acting on the Indo‐Australian Plate control the long wavelength of the maximum horizontal present‐day stress orientation in the Australian continent. However, all numerical models of the stress field to date are unable to predict the observed orientation of maximum horizontal stress in the northeast of New South Wales, Australia. Recent coal seam gas exploration in the Clarence‐Moreton Basin, eastern Australia, provides an opportunity to better evaluate the state of crustal stress in this part of the continent where only limited information was available prior to this study. Herein, we conduct the first analysis of the present‐day tectonic stress in the Clarence‐Moreton Basin, from drilling‐induced tensile fractures and borehole breakouts interpreted using 11.3 km of acoustic image logs in 27 vertical wells. A total of 2822 drilling‐induced stress indicators suggest a mean orientation of N069°E (±23°) for the maximum horizontal present‐day stress in the basin which is different from that predicted by published geomechanical‐numerical models. In addition, we find significant localised perturbations of borehole breakouts, both spatially and with depth, that are consistent with stress variations near faults, fractures and lithological contrasts, indicating that local structures are an important source of stress in the basin. The observation that structures can have a major control on the stresses in the basin suggests that, while gravity and plate boundary forces have the major role in the long wavelength (first‐order) stress pattern of the continent, local perturbations are significant and can lead to substantial changes in the orientation of the maximum horizontal present‐day stress, particularly at the basin scale. These local perturbations of stress as a result of faults and fractures have important implications in borehole stability and permeability of coal seam gas reservoirs for safe and sustainable extraction of methane in this area.     

A geomagnetic study of the Great Glen Fault

Summary. An attempt is made to determine the range of two-dimensional current models consistent with the measured magnetovariational response, for periods from 5–30 min, near the Great Glen Fault in northern Scotland. All current models must be symmetric about the fault line but, because of uncertainty about the magnitude of the ocean effect, models ranging from a line current at 80 km depth to a uniform current sheet, 60 km wide, at 10 km depth are equally acceptable. Comparison with other geophysical studies of the same area suggests that a suitable conducting zone is unlikely to be present at shallow depths and interpretation in terms of a conducting zone in the 20–80 km depth range is favoured, although no such zone has been resolved by the other studies.     

Experimental study of the anisotropy of longitudinal and transverse waves from local earthquake records

Summary. The paper gives the results of a study of the anisotropy of seismic wave velocities within the Ashkhabad test field in Central Asia. The anisotropy was studied by analysing variations in the values of apparent velocities of first arrivals for epicentral distances ranging from 30 to 130 km and by analysing the delays (Δ t_s1-_s2 ) between the arrival times of shear waves with different polarizations.
The velocities of P -waves vary with azimuth from 5.3 to 6.27 km s^-1 and the velocities of S -waves vary from 3.15 to 3.5 km s^-1.
The delay times Δ t_S1 - _S2 depend on the direction of the propagation. The character of the variation of the propagation velocity of the longitudinal wave, the presence of two differently polarized shear waves S 1 and S 2 propagating at different velocities, and the character of the distribution of Δ t_S1 - _S2 on the stereogram suggest that the symmetry of the anisotropic medium is close to hexagonal with a nearly horizontal symmetry axis coinciding with the direction of maximal velocity. The azimuth of the symmetry axis of the medium is 140° and coincides with the direction of geological faults.     

Sensitivity and inversion of borehole flexural dispersions for formation parameters

Propagation characteristics of borehole flexural waves are functions of five borehole, fluid and formation parameters. Uncertainties in any of these model parameters may introduce errors in the estimate of formation shear-wave speed. Other sources of error in the estimate may be caused by deviations from the assumed circular borehole cross-section or heterogeneity in the material properties of formation in the vicinity of a borehole from the usually assumed homogeneous properties. The influence of uncertainties in model parameters on borehole flexural dispersion has been calculated from a general perturbation model based on Hamilton's principle. A sensitivity analysis of the flexural dispersion to small variations in the model parameters shows that the formation shear speed has by far the most dominant influence in a slow formation. In contrast, the flexural dispersion in a fast formation is significantly influenced by three of the five model parameters: the formation shear speed, the borehole fluid compressional speed and the borehole diameter. The frequency dependence of these sensitivity functions indicates that the inversion of flexural dispersion for formation shear speed is optimal in the range 2–4 kHz for a borehole of diameter 25.4 cm. The range of validity of the perturbation model has been estimated by comparing results of concentric annuli of different thicknesses, and shear-wave speeds different from that of the formation with those from exact numerical solutions from a modal search program. Generally, the perturbation results for altered annuli of thicknesses up to 15 cm are accurate to within 1 per cent for a shear-wave speed 10 per cent lower than that of the formation. This difference between the perturbation and exact results increases to approximately 1 to 3 per cent for a shear-wave speed 20 per cent lower than that of the formation.