Three-dimensional image of the anisotropic bodies beneath central Honshu, Japan

Shear-wave splitting from local deep earthquakes is investigated to clarify the volume and the location of two anisotropic bodies in the mantle wedge beneath central Honshu, Japan. We observe a spatial variation in splitting parameters depending on the combination of sources and receivers, nearly N–S fast in the northern region, nearly E–W fast in the southern region and small time delays in the eastern region. Using forward modelling, two models with 30 and 10 per cent anisotropy are tested by means of a global search for the locations of anisotropic bodies with various volumes. The optimum model is obtained for 30 per cent anisotropy, which means a 5 per cent velocity difference between fast and slow polarized waves. The northern anisotropic body has a volume of 1.00° (longitude) × 0.5° (latitude) × 75 km (depth), with the orientation of the symmetry axis being N20°E. The southern anisotropic body has a volume of 1.25° × 1.25° × 100 km with the symmetry axis along N95°E. Our results show that the anisotropic bodies are located in low-velocity and low- Q regions of the mantle. This, together with petrological data and the location of volcanoes in the arc, suggests that the possible cause of the anisotropy is the preferred alignment of cracks filled with melt.     

Two-dimensional modelling of subduction zone anisotropy with application to southwestern Japan

We present a series of 2-D numerical models of viscous flow in the mantle wedge induced by a subducting lithospheric plate. We use a kinematically defined slab geometry approximating the subduction of the Philippine Sea plate beneath Eurasia. Through finite element modelling we explore the effects of different rheological and thermal constraints (e.g. a low-viscosity region in the wedge corner, power law versus Newtonian rheology, the inclusion of thermal buoyancy forces and a temperature-dependent viscosity law) on the velocity and finite strain field in the mantle wedge. From the numerical flow models we construct models of anisotropy in the wedge by calculating the evolution of the finite strain ellipse and combining its geometry with appropriate elastic constants for effective transversely isotropic mantle material. We then predict shear wave splitting for stations located above the model domain using expressions derived from anisotropic perturbation theory, and compare the predictions to ∼500 previously published shear wave splitting measurements from seventeen stations of the broad-band F-net array located in southwestern Japan. Although the use of different model parameters can have a substantial effect on the character of the finite strain field, the effect on the average predicted splitting parameters is small. However, the variations with backazimuth and ray parameter of individual splitting intensity measurements at a given station for different models are often different, and rigorous analysis of details in the splitting patterns allows us to discriminate among different rheological models for flow in the mantle wedge. The splitting observed in southwestern Japan agrees well with the predictions of trench-perpendicular flow in the mantle wedge along with B-type olivine fabric dominating in a region from the wedge corner to about 125 km from the trench.     

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.     

Distribution of seismic anisotropy in the subduction zone beneath the Wellington region, New Zealand

Shear wave splitting measurements from S arrivals of local earthquakes recorded at the Incorporated Research Institutions for Seismology (IRIS) broadband sensor SNZO are used to determine a basic anisotropic structure for the subduction zone in the Wellington region. With the use of high-frequency filters, fast anisotropic polarization ( φ ) and splitting time ( δt ) measurements typical of crustal anisotropy are evident, but the larger splitting expected from the mantle is often not resolved. The small splitting seen agrees well with the results of previous studies concerning shallow crustal anisotropy. With the use of lower-frequency filters, measurements more consistent with mantle anisotropy are made. Anisotropy of 4.4 ± 0.9 per cent with a fast polarization of 29° ± 38° is calculated for the subducting slab, from 20 to 70  km depth. Using this result in addition to the results of previous studies, a model is proposed. The model requires a frequency-dependent anisotropy of less than 1.4 per cent when measured with a period of ~2  s to be present in the sub-slab mantle.
Separate from this population, a band of events in northern Cook Strait with an 86° ± 10° fast polarization is seen. This is at about 40° from the strike of the Hikurangi margin, and suggests a source of shear strain 40° removed from that found in the majority of the region. The cause of this is probably a deformation in the subducting slab in this region, as it moves towards a greater incline to the south.     

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 and the stress field from borehole recordings at 2.5 km depth at Cajon Pass

53 local earthquakes recorded at 2.5 km depth in the Cajon Pass scientific borehole are analysed for shear-wave splitting. The time delays between the split shear waves can be positively identified for 32 of the events. Modelling these observations of polarizations and time delays using genetic algorithms suggests that the anisotropic structure near Cajon Pass has orthorhombic symmetry. The polarization of the shear waves and the inferred strike of the stress-aligned fluid-filled intergranular microcracks and pores suggests that the maximum horizontal compressional stress direction is approximately N13°W. This is consistent with previous results from earthquake source mechanisms and the right-lateral strike-slip motion on the nearby San Andreas Fault, but not with stresses measured within the uppermost 3 km of the borehole. This study suggests that the San Andreas Fault is driven by deeper tectonic stresses and the present understanding of a weak and frictionless San Andreas Fault may need to be modified. The active secondary faulting and folding close to the fault are probably driven by the relatively shallow stress as measured in the 3.5 km deep borehole.