Shear zones in porous sand: Insights from ring-shear experiments and naturally deformed sandstones

In a high-resolution small-scale seismic experiment we investigated the shallow structure of the Wadi Araba fault (WAF), the principal fault strand of the Dead Sea Transform System between the Gulf of Aqaba/Eilat and the Dead Sea. The experiment consisted of 8 sub-parallel 1 km long seismic lines crossing the WAF. The recording station spacing was 5 m and the source point distance was 20 m. The first break tomography yields insight into the fault structure down to a depth of about 200 m. The velocity structure varies from one section to the other which were 1 to 2 km apart, but destinct velocity variations along the fault are visible between several profiles. The reflection seismic images show positive flower structures and indications for different sedimentary layers at the two sides of the main fault. Often the superficial sedimentary layers are bent upward close to the WAF. Our results indicate that this section of the fault (at shallow depths) is characterized by a transpressional regime. We detected a 100 to 300 m wide heterogeneous zone of deformed and displaced material which, however, is not characterized by low seismic velocities at a larger scale. At greater depth the geophysical images indicate a blocked cross-fault structure. The structure revealed, fault cores not wider than 10 m, are consistent with scaling from wear mechanics and with the low loading to healing ratio anticipated for the fault.     

NEHRP soil classification and estimation of 1-D site effect of Dehradun fan deposits using shear wave velocity

Shear wave velocity of the near surface soil at nearly 50 sites in the sub Himalayan mountain exit covering Doon fan deposits, was determined using Multi-channel Analysis of Surface Waves (MASW), a seismic reflection technique. Based on the average shear wave velocity of the upper 30 m soil column, sites in the Dehradun fan are predominantly classified as class ‘D’ (180–360 m/s). Similarly, sites located on the northwestern, eastern and southeastern sides of the fan deposit have shear wave velocities (in the upper 30 m soil) greater than 360 m/s, thereby classifying them as class ‘C’ (360–760 m/s) in accordance with NEHRP provisions. Some of the sites towards the southwestern side of the fan deposits had average shear wave velocities less than 180 m/s and could be classified as soil class ‘E’. One dimensional site effects, including amplification and dynamic period were calculated for the majority of the sites. However, some of the representative suite of sites across the north–south profile of Dehradun fan has been discussed here. Although the attenuation is greater on the southwestern side of the Dehradun fan deposits (i.e. thicker, low velocity sediments) and the sites had been classified as class ‘D’ and ‘E’ but the site amplification tends to be greater in the northern and northwestern part of the city due to large impedance contrast with in the near surface soils.     

A GPR study of subsidence-creep-fault processes in Morelia, Michoacán, Mexico

In 1983, inhabitants of the City of Morelia, Michoacán, Mexico, began to observe a series of differential settlements causing damages to constructions along linear trends parallel to a system of regional faults. The same phenomenon occurs in others cities of the Mexican Volcanic Belt (MVB), such as Celaya, Aguascalientes, and Querétaro, and is linked to a structurally controlled subsidence, caused by groundwater withdrawal, and the presence of geological faults. We define this subsidence type as Subsidence-Creep-Fault Processes (SCFP), based on the necessary elements for their generation, and we studied them through geophysical and geotechnical techniques. In Morelia, the geophysical investigations have been carried out using ground-penetrating radar (GPR). GPR profiles, perpendicular to the axis of the surface fault generated by the SCFP were carried out. The common-offset single-fold profiling was used, with a central frequency of 50 MHz. In all cases it has been possible to visualize a fault plane dividing two blocks, the presence of synthetic and antithetic faults, influence zones from 20 m to 40 m, and a maximum “net throw” of 4 m. Exploration trenches followed the same direction of the profiles obtained with GPR (perpendicular to the axis of the surface fault). These trenches exposed a fault plane dividing two blocks with different lithology, generating a maximum “net throw” of 4.40 m; as well they help in the determination of influence zones that varied from 14 m to 40 m.     

Imaging granitic plutons along the IBERSEIS profile

The parameters used for the acquisition of the IBERSEIS deep seismic reflection profile in the SW Iberian Peninsula provide seismic images of the deep crust as well as a high resolution section of the shallow subsurface. A very dense array of sources and receivers allowed high resolution tomographic studies in zones of special interest (granitic plutons). The three dimensional tomographic inversion produced velocity models along a 500 m wide and 1000 m deep strip along the IBERSEIS transect in the areas of La Bazana, La Dehesilla, Feria and Villafranca. In these high resolution velocity models (sampled by 50 × 50 × 50 m cells), high velocity anomalies indicate the geology and extension of the granitic plutons at depth. This directly correlates with the surface outcrops. Moreover, tomographic models provided valuable information for the geometry and characterization of fractured and fresh domains in a rock volume. Furthermore, a piggy back seismic acquisition experiment using additional seismic instrumentation from the University of Paris Sud (40-channel DMT) provided perpendicular, offline recordings of the Vibroseis sources. This additional recording system was deployed perpendicularly to the main IBERSEIS seismic reflection line and provided additional 3D control.     

The Gibraltar Arc seismogenic zone (part 2): Constraints on a shallow east dipping fault plane source for the 1755 Lisbon earthquake provided by tsunami modeling and seismic intensity

The Great Lisbon earthquake has the largest documented felt area of any shallow earthquake and an estimated magnitude of 8.5–9.0. The associated tsunami ravaged the coast of SW Portugal and the Gulf of Cadiz, with run-up heights reported to have reached 5–15 m. While several source regions offshore SW Portugal have been proposed (e.g.— Gorringe Bank, Marquis de Pombal fault), no single source appears to be able to account for the great seismic moment as well as all the historical tsunami amplitude and travel time observations. A shallow east dipping fault plane beneath the Gulf of Cadiz associated with active subduction beneath Gibraltar, represents a candidate source for the Lisbon earthquake of 1755.Here we consider the fault parameters implied by this hypothesis, with respect to total slip, seismic moment, and recurrence interval to test the viability of this source. The geometry of the seismogenic zone is obtained from deep crustal studies and can be represented by an east dipping fault plane with mean dimensions of 180 km (N–S) × 210 km (E–W). For 10 m of co-seismic slip an Mw 8.64 event results and for 20 m of slip an Mw 8.8 earthquake is generated. Thus, for convergence rates of about 1 cm/yr, an event of this magnitude could occur every 1000–2000 years. Available kinematic and sedimentological data are in general agreement with such a recurrence interval. Tsunami wave form modeling indicates a subduction source in the Gulf of Cadiz can partly satisfy the historical observations with respect to wave amplitudes and arrival times, though discrepancies remain for some stations. A macroseismic analysis is performed using site effect functions calculated from isoseismals observed during instrumentally recorded strong earthquakes in the region (M7.9 1969 and M6.8 1964). The resulting synthetic isoseismals for the 1755 event suggest a subduction source, possibly in combination with an additional source at the NW corner of the Gulf of Cadiz can satisfactorily explain the historically observed seismic intensities. Further studies are needed to sample the turbidites in the adjacent abyssal plains to better document the source region and more precisely calibrate the chronology of great earthquakes in this region.     

Refining the image of the San Andreas Fault near Parkfield, California using a finite difference travel time computation technique

The Parkfield Area Seismic Observatory (PASO) was a dense, telemetered seismic array that operated for nearly 2 years in a 15 km aperture centered on the San Andreas Fault Observatory at Depth (SAFOD) drill site. The main objective of this deployment was to refine the locations of earthquakes that will serve as potential targets for SAFOD drilling and in the process develop a high (for passive seismological techniques) resolution image of the fault zone structure. A challenging aspect of the analysis of this data set was the known existence of large (20–25%) contrasts in seismic wavespeed across the San Andreas Fault. The resultant distortion of raypaths could challenge the applicability of approximate ray tracing techniques. In order to test the sensitivity of our hypocenter locations and tomographic image to the particular ray tracing and inversion technique employed, we compare an initial determination of locations and structure developed using a coarse grid and an approximate ray tracer [Thurber, C., Roecker, S., Roberts, K., Gold, M., Powell, M.L. , and Rittger, K., 2003. Earthquake locations and three-dimensional fault zone structure along the creeping section of the San Andreas fault near Parkfield, CA: Preparing for SAFOD, Geophys. Res. Lett., 30 3, 10.1029/2002GL016004.] with one derived from a relatively fine grid and an application of a finite difference algorithm [Hole, J.A., and Zelt, B.C., 1995. 3-D finite-difference reflection traveltimes, Geophys. J. Int., 121, 2, 427–434.]. In both cases, we inverted arrival-time data from about 686 local earthquakes and 23 shots simultaneously for earthquake locations and three-dimensional Vp and Vp/Vs structure. Included are data from an active source seismic experiment around the SAFOD site as well as from a vertical array of geophones installed in the 2-km-deep SAFOD pilot hole, drilled in summer 2002. Our results show that the main features of the original analysis are robust: hypocenters are located beneath the trace of the fault in the vicinity of the drill site and the positions of major contrasts in wavespeed are largely the same. At the same time, we determine that shear wave speeds in the upper 2 km of the fault zone are significantly lower than previously estimated, and our estimate of the depth of the main part of the seismogenic zone decreases in places by  100 m. Tests using “virtual earthquakes” (borehole receiver gathers of picks for surface shots) indicate that our event locations near the borehole currently are accurate to about a few tens of meters horizontally and vertically.     

Rupture length and velocity for earthquakes in the Mid-Atlantic Ridge from directivity effect in body and surface waves

The dimensions and rupture velocities of four earthquakes, two in the Mid-Atlantic Ridge and two in Iceland with strike–slip mechanisms and magnitudes (M_w) between 6.2 and 6.8 were studied using the directivity effects of Rayleigh and body waves. For Rayleigh waves we used the directivity function for different pairs of stations and for body waves the waveforms of P and SH waves corresponding to a simple extended line source. We have found that three have very shallow depths about 3 km and one 8 km, fault lengths between 12 km and 21 km, and a low rupture velocity of about 1.5 km/s to 2.0 km/s which supports the idea of the presence of slow earthquakes in transform faults.     

ESR detection of seismic frictional heating events in the Nojima fault drill core samples, Japan

We have analyzed the Nojima fault NIED 1800 m drill core samples by ESR (Electron Spin Resonance) to detect seismic frictional heating events, especially during the 1995 Kobe Earthquake. Dark gray fault gouge with foliation > 10 cm away from the fault plane at about 1140 m in depth, which was produced by ancient fault movements, has a FMR (ferrimagnetic resonance) signal. Heating experiments show that this FMR signal is derived from ferrimagnetic trivalent ion oxides (γ-Fe₂O₃: maghemite) with imperfect crystallinity, which is produced by thermal dehydration of γ-FeOOH (lepidocrocite) or Fe(OH)₃ (limonite). The existence of the FMR signal means that dry heating such as frictional heating once occurred, and that the frictional heat temperature along the dark gray fault gouge may have risen to over 350 °C during ancient seismic fault slip. In order to detect frictional heating events in fault zones, the increase of the FMR signal and the color change of fault gouge into dark gray or black are important indexes. On the other hand, no FMR signal is detected from the fault gouges just on two fault planes at about 1140 m and 1300 m in depth, which are considered to be possible main fault planes in the 1995 Kobe Earthquake. These two fault planes may not have played an important role of fault slip in the Earthquake.     

Thermal maturity of a fold–thrust belt based on vitrinite reflectance analysis in the Western Foothills complex, western Taiwan

Burial depth, cumulative displacement, and peak temperature of frictional heat of a fault system are estimated by thermal analysis in the fold–thrust belt of the Western Foothills complex, western Taiwan based on the vitrinite reflectance technique. The regional thermal structure across the complex reveals that the rocks were exposed to maximum temperatures ranging from 100 °C to 180 °C, which corresponds to a burial depth of 3.7–6.7 km. A large thermal difference of 90 °C were observed at the Shuilikeng fault which make the eastern boundary of the fold–thrust belt where it is in contact with metamorphic rock of Hsuehshan Range. The large thermal difference corresponds to cumulative displacements on the Shuilikeng fault estimated to be in the range of 5.2–6.9 km. However, thermal differences in across the Shuangtung and Chelungpu faults cannot be determined apparently due to small vertical offsets. The large displacement observed across the Shuilikeng fault is absent at the other faults which are interpreted to be younger faults within the piggyback thrust system. Localized high temperatures adjacent to fault zones were observed in core samples penetrating the Chelungpu fault. Three major fracture zones were observed at core lengths of 225 m, 330 m, and 405 m and the two lower zones which comprise dark gray narrow shear zones. A value of vitrinite reflectance of 1.8%, higher than the background value of 0.8%, is limited at a narrow shear zone of 1 cm thickness at the fracture zone at 330 m. The estimated peak temperature in the range of 550–680 °C in the shear zone is far higher than the background temperature of 130 °C, and it is interpreted as due to frictional heating during seismic faulting.     

Local thickening of the Cascadia forearc crust and the origin of seismic reflectors in the uppermost mantle

Seismic reflection profiles from three different surveys of the Cascadia forearc are interpreted using P wave velocities and relocated hypocentres, which were both derived from the first arrival travel time inversion of wide-angle seismic data and local earthquakes. The subduction decollement, which is characterized beneath the continental shelf by a reflection of 0.5 s duration, can be traced landward into a large duplex structure in the lower forearc crust near southern Vancouver Island. Beneath Vancouver Island, the roof thrust of the duplex is revealed by a 5–12 km thick zone, identified previously as the E reflectors, and the floor thrust is defined by a short duration reflection from a < 2-km-thick interface at the top of the subducting plate. We show that another zone of reflectors exists east of Vancouver Island that is approximately 8 km thick, and identified as the D reflectors. These overlie the E reflectors; together the two zones define the landward part of the duplex. The combined zones reach depths as great as 50 km. The duplex structure extends for more than 120 km perpendicular to the margin, has an along-strike extent of 80 km, and at depths between 30 km and 50 km the duplex structure correlates with a region of anomalously deep seismicity, where velocities are less than 7000 m s^− 1. We suggest that these relatively low velocities indicate the presence of either crustal rocks from the oceanic plate that have been underplated to the continent or crustal rocks from the forearc that have been transported downward by subduction erosion. The absence of seismicity from within the E reflectors implies that they are significantly weaker than the overlying crust, and the reflectors may be a zone of active ductile shear. In contrast, seismicity in parts of the D reflectors can be interpreted to mean that ductile shearing no longer occurs in the landward part of the duplex. Merging of the D and E reflectors at 42–46 km depth creates reflectivity in the uppermost mantle with a vertical thickness of at least 15 km. We suggest that pervasive reflectivity in the upper mantle elsewhere beneath Puget Sound and the Strait of Georgia arises from similar shear zones.     

Geophysical methods applied to fault characterization and earthquake potential assessment in the Lower Tagus Valley, Portugal

The study region is located in the Lower Tagus Valley, central Portugal, and includes a large portion of the densely populated area of Lisbon. It is characterized by a moderate seismicity with a diffuse pattern, with historical earthquakes causing many casualties, serious damage and economic losses. Occurrence of earthquakes in the area indicates the presence of seismogenic structures at depth that are deficiently known due to a thick Cenozoic sedimentary cover. The hidden character of many of the faults in the Lower Tagus Valley requires the use of indirect methodologies for their study. This paper focuses on the application of high-resolution seismic reflection method for the detection of near-surface faulting on two major tectonic structures that are hidden under the recent alluvial cover of the Tagus Valley, and that have been recognized on deep oil-industry seismic reflection profiles and/or inferred from the surface geology. These are a WNW–ESE-trending fault zone located within the Lower Tagus Cenozoic basin, across the Tagus River estuary (Porto Alto fault), and a NNE–SSW-trending reverse fault zone that borders the Cenozoic Basin at the W (Vila Franca de Xira–Lisbon fault). Vertical electrical soundings were also acquired over the seismic profiles and the refraction interpretation of the reflection data was carried out. According to the interpretation of the collected data, a complex fault pattern disrupts the near surface (first 400 m) at Porto Alto, affecting the Upper Neogene and (at least for one fault) the Quaternary, with a normal offset component. The consistency with the previous oil-industry profiles interpretation supports the location and geometry of this fault zone. Concerning the second structure, two major faults were detected north of Vila Franca de Xira, supporting the extension of the Vila Franca de Xira–Lisbon fault zone northwards. One of these faults presents a reverse geometry apparently displacing Holocene alluvium. Vertical offsets of the Holocene sediments detected in the studied geophysical data of Porto Alto and Vila Franca de Xira–Lisbon faults imply minimum slip rates of 0.15–0.30 mm/year, three times larger than previously inferred for active faults in the Lower Tagus Valley and maximum estimates of average return periods of 2000–5000 years for M 6.5–7 co-seismic ruptures.     

The North Maladeta Fault (Spanish Central Pyrenees) as the Vielha 1923 earthquake seismic source: Recent activity revealed by geomorphological and geophysical research

The Spanish Central Pyrenees have been the scenario of at least two damaging earthquakes in the last 800 years. Analysis of macroseismic data of the most recent one, the Vielha earthquake (19 November 1923), has led to the identification of the North Maladeta Fault (NMF) as the seismic source of the event. This E–W trending fault defines the northern boundary of the Maladeta Batholith and corresponds to a segment of the Alpine Gavarnie thrust fault. Our study shows that the NMF offsets a reference Neogene peneplain. The maximum observed vertical displacement is  730 m, with the northern downthrown sector slightly tilting towards the South. This offset provides evidence of normal faulting and together with the presence of tectonic faceted spurs allowed us to geomorphically identify a fault trace of 17.5 km. This length suggests that a maximum earthquake of M_w = 6.5 ± 0.66 could occur in the area. The geomorphological study was improved with a resistivity model obtained at Prüedo, where a unique detritic Late Miocene sequence crops out adjacent to the NMF. The section is made up of 13 audiomagnetotelluric soundings along a 1.5 km transect perpendicular to the fault trace at Prüedo and reveals the structure in depth, allowing us to interpret the Late Miocene deposits as tectonically trapped basin deposits associated with normal faulting of the NMF. The indirect age of these deposits has been constrained between 11.1 and 8.7 Ma, which represents a minimum age for the elevated Pyrenean peneplain in this part of the Pyrenees. Therefore, we propose the maximum vertical dip-slip rate for the NMF to be between 0.06 and 0.08 mm/a. Normal faulting in this area is attributed to the vertical lithospheric stress associated with the thickened Pyrenean crust.     

Rock Damage Structure of the South Longmen-Shan Fault in the 2008 M8 Wenchuan Earthquake Viewed with Fault-Zone Trapped Waves and Scientific Drilling

This article is to review results from scientific drilling and fault-zone trapped waves(FZTWs) at the south Longman-Shan fault(LSF) zone that ruptured in the 2008 May 12 M8 Wenchuan earthquake in Sichuan, China. Immediately after the mainshock, two Wenchuan Fault Scientific Drilling(WFSD) boreholes were drilled at WFSD-1 and WFSD-2 sites approximately 400 m and 1 km west of the surface rupture along the Yinxiu-Beichuan fault(YBF), the middle fault strand of the south LSF zone. Two boreholes met the principal slip of Wenchuan earthquake along the YBF at depths of 589-m and 1230-m, respectively. The slip is accompanied with a 100-200-m-wide zone consisting of fault gouge, breccia, cataclasite and fractures. Close to WFSD-1 site, the nearly-vertical slip of ~4.3-m with a 190-m wide zone of highly fractured rocks restricted to the hanging wall of the YBF was found at the ground surface after the Wenchuan earthquake. A dense linear seismic array was deployed across the surface rupture at this venue to record FZTWs generated by aftershocks. Observations and 3-D finite-difference simulations of FZTWs recorded at this cross-fault array and network stations close to the YBF show a distinct low-velocity zone composed by severely damaged rocks along the south LSF at seismogenic depths. The zone is several hundred meters wide along the principal slip, within which seismic velocities are reduced by ~30–55% from wall-rock velocities and with the maximum velocity reduction in the ~200-m-wide rupture core zone at shallow depth. The FZTW-inferred geometry and physical properties of the south LSF rupture zone at shallow depth are in general consistent with the results from petrological and structural analyses of cores and well log at WFSD boreholes. We interpret this remarkable low-velocity zone as being a break-down zone during dynamic rupture in the 2008 M8 earthquake. We examined the FZTWS generated by similar earthquakes before and after the 2008 mainshock and observed that seismic velocities within fault core zone was reduced by ~10% due to severe damage of fault rocks during the M8 mainshock. Scientific drilling and locations of aftershocks generating prominent FZTWs also indicate rupture bifurcation along the YBF and the Anxian-Guangxian fault(AGF), two strands of the south LSF at shallow depth. A combination of seismic, petrologic and geologic study at the south LSF leads to further understand the relationship between the fault-zone structure and rupture dynamics, and the amplification of ground shaking strength along the low-velocity fault zone due to its waveguide effect.     

Scaling properties within the Gulf of Corinth, Greece; comparison between offshore and onshore active faults

The Gulf of Corinth, Greece, is a 110-km-long by 30-km-wide active graben displaying strong seismicity hosted both on north and south dipping normal faults. This complex fault pattern consists of two fault populations, offshore and onshore. The offshore fault population is investigated by densely arranged seismic reflection profiles during the last 20 years, whereas the onshore fault population displays spectacular and well exposed faults, delineated by high accuracy mapping. We analyzed fault length and throw, in order to study the scaling properties of 136 well-determined offshore and onshore faults and the comparison between the two datasets. We examined the statistical properties on both fault populations, in order to describe the role of segmentation in the growth of faults and the different stages of the evolution of the fault networks.Our results on power law relationships associated with the scaling properties of the fault zones in the Gulf of Corinth, suggest that both fault populations are bi-fractal, providing the initiation of a sature state in deformation. In addition, the vertical throw of faults shows that both fault populations have similar properties but different distributions below and above 5 km, respectively. Displacement–length ratios, show that faults larger than 9 km appear to accumulate throw without any dramatic change to their length. These observations combined with other geophysical studies within the Gulf, suggest that the characteristic fault lengths of 5 km and 9 km can be correlated to the crustal mechanical structure and the seismicity of the Gulf.     

Insights into Wasatch fault vertical slip rates using the age of sediments in Timpanogos Cave, Utah

Timpanogos Cave, located near the Wasatch fault, is about 357 m above the American Fork River. Fluvial cave sediments and an interbedded carbonate flowstone yield a paleomagnetic and U–Th depositional age of 350 to 780 ka. Fault vertical slip rates, inferred from calculated river downcutting rates, range between 1.02 and 0.46 mm yr^− 1. These slip rates are in the range of the 0–12 Ma Wasatch Range exhumation rate ( 0.5–0.7 mm yr^− 1), suggesting that the long-term vertical slip rate remained stable through mid-Pleistocene time. However, the late Pleistocene (0–250 ka) decelerated slip rate ( 0.2–0.3 mm yr^− 1) and the accelerated Holocene slip rate ( 1.2 mm yr^− 1) are consistent with episodic fault activity. Assuming that the late Pleistocene vertical slip rate represents an episodic slowing of fault movement and the long-term (0–12 Ma) average vertical slip rate, including the late Pleistocene and Holocene, should be  0.6 mm yr^− 1, there is a net late Pleistocene vertical slip deficit of  50–75 m. The Holocene and late Pleistocene slip rates may be typical for episodes of accelerated and slowed fault movement, respectively. The calculated late Pleistocene slip deficit may mean that the current accelerated Wasatch fault slip rate will extend well into the future.     

First seismic record for intra-arc strike-slip tectonics along the Liquiñe-Ofqui fault zone at the obliquely convergent plate margin of the southern Andes

A temporal seismic network recorded local seismicity along a 130 km long segment of the transpressional dextral strike-slip Liquiñe-Ofqui fault zone (LOFZ) in southern Chile. Seventy five shallow crustal events with magnitudes up to M_w 3.8 and depths shallower than 25 km were observed in an 11-month period mainly occurring in different clusters. Those clusters are spatially related to the LOFZ, to the volcanoes Chaitén, Michinmahuida and Corcovado, and to active faulting on secondary faults. Further activity along the LOFZ is indicated by individual events located in direct vicinity of the surface expression of the LOFZ. Focal mechanisms were calculated using deviatoric moment tensor inversion of body wave amplitude spectra which mostly yield strike-slip mechanisms indicating a NE–SW direction of the P-axis for the LOFZ at this latitude. The seismic activity reveals the present-day activity of the fault zone. The recent M_w 6.2 event near Puerto Aysén, Southern Chile at 45.4°S on April 21, 2007 shows that the LOFZ is also capable of producing large magnitude earthquakes and therefore imposing significant seismic hazard to this region.     

Constraints on seismic velocity anomalies beneath the Siberian craton from xenoliths and petrophysics

High seismic Vp velocity anomalies (8.7–9.0 km s^− 1) have long been known about in regions of the uppermost mantle of the Siberian craton, often in association with kimberlite fields. Laboratory measurement of seismic properties of five xenoliths, three peridotites and two eclogites, from the Udachnaya kimberlite under confining pressures up to 600 MPa were extrapolated to uppermost mantle P–T conditions of 1500 MPa and 500 °C, however none of the velocities are high enough to explain the observations. Eclogites or peridotites are commonly considered to be the source of anomalous high velocities. We prefer a peridotitic source to an eclogitic source due to the unusual chemistry and regional uniformity of eclogitic garnets required, maximum velocity limitations on laboratory measurements of seismic properties of natural eclogites, and purported abundance of eclogites in the lithosphere. Alternatively, a highly depleted peridotite, such as dunite or harzburgite, can produce velocities high enough to match observations. Olivine petrofabrics in most peridotites, including the three peridotites used in this study, are great enough to produce the observed high velocities provided olivine petrofabrics are continuous enough and correctly oriented to be seismically detectable and the modal proportion of olivine is high. There have been suggestions by other authors that the Siberian upper mantle is highly depleted and that a lithosphere-scale shear zone exists, which may have acted to organize fabrics into segments large enough for detection. Anomalously high Vp–Vs velocity ratios of greater than 1.8 are expected parallel to the olivine [100] maxima required to be present in a high-velocity olivine-dominated upper mantle. Vp–Vs velocity ratios can serve as a means of inferring large-scale anisotropy when limited seismic data are available, as in Siberia.     

Tomographic image of crustal structures across the Chelungpu fault: Is the seismogenic layer structure- or depth-dependent?

In 1999, a large earthquake (M_w = 7.6) occurred along the Chelungpu fault in the fold-and-thrust belt of western Taiwan. To shed more light on the subsurface structures and the seismogenic layers, three-dimensional velocity structures were inverted by using the travel times of both P- and S-waves from 2391 aftershocks recorded by the Central Weather Bureau during the 15 months that followed. From tomography, a typical image of the large-scale thrusting structures in the upper crust across the Chelungpu fault was obtained. In general, high velocities beneath the Western Foothills and Central Ranges are separated from low velocities beneath the Coastal Plain by an east-dipping boundary that is roughly consistent with the Chelungpu fault on the surface. The contrast in velocity on either side of the Chelungpu fault is indicative of about a 7- to 9-km vertical offset in the upper crust. The relocated hypocenter for the Chi-Chi earthquake shifts by 2.2 km toward the northwest, and its focal depth decreases by 0.7 km. A plot of focal depths versus rock velocities where the aftershocks occurred shows earthquakes are more inclined to occur in rock with a velocity of around 5.6 km/s. This strongly suggests the seismogenic layer in the fold-and-thrust belt of Taiwan is more structure-dependent than depth-dependent.     

Time evolution of hydraulic and electrokinetic parameters around the Nojima fault, Japan, estimated by an electrokinetic method

The Nojima Fault Zone Probe was designed to study the properties and healing processes of the Nojima fault, which is the surface fault rupture of the Hyogo-ken Nanbu earthquake (M7.2) of 1995 (1995 Kobe earthquake). In this project, water injection experiments were conducted in a borehole of 1800 m depth at the Nojima fault. We set up electrodes around the borehole and observed self-potential variations to investigate the magnitude of electrokinetic and hydraulic parameters around the Nojima fault zone. In the 1997 experiment, self-potential variations were in the range of a few to about 20 mV across 320–450 m electrode dipoles with hydraulic pressure variations from 3.5 to 4 MPa. In the 2000 experiment, self-potential variations were in the range of a few to about 85 mV across 160–260 m electrode dipoles with the hydraulic pressure variations from 3 to 4.5 MPa. In the 2003 experiment, self-potential variations were in the range of a few to about 30 mV across 20–80 m electrode dipoles with hydraulic pressure of 4 MPa. These observed self-potential variations were explained well with an electrokinetic effect due to the underground flow of the injected water. From the observed results, we estimated that the ratio of hydraulic parameters (permeability, porosity, and tortuosity) to electrokinetic parameters (zeta potential and dielectric constant) decreased approximately 40% during eight years after the earthquake. This result suggests that the healing process around the fault zone progress.