Streaming potential observations, using geothermal wells and in situ electrokinetic coupling coefficients under high temperature

In an attempt to detect streaming potentials induced by subsurface water flows, we have observed the horizontal electric field (self-potential) variations across stationary electric dipoles near geothermal wells in the Takinoue geothermal area, Japan. We observed variations of self-potential which seem to be associated with the water flows in the aquifer, induced by turning on and off the flow of the wells. Amplitudes of the variations are 3–5 mV across 60–200 m dipoles, and can be explained well with a proposed electrokinetic model: the streaming potential coefficient of − 15 mV/bar and/or the ζ-potential of −50 to −100 mV in the aquifer are appropriate to explain the observed data by the model. The obtained electrokinetic coupling coefficients are in situ ones and determined for crustal rock-water system under high temperature (˜200°C) condition. The present results, together with a laboratory study by Ishido and Mizutani (1981), give fundamental information on electrokinetic coupling coefficients in the earth's interior, and are very important when we make quantitative interpretations of self-potentials generated by geothermal activity on the basis of electrokinetic effects.     

Borehole water and hydrologic model around the Nojima fault, SW Japan

The active fault drilling at Nojima Hirabayashi after the 1995 Hyogoken-nanbu (Kobe) earthquake (M_JMA = 7.2) provides us with a unique opportunity to investigate subsurface fault structure and the in-situ properties of fault and fluid. The borehole intersected the fault gouge of the Nojima fault at a depth interval of 623 m to 625 m. The lithology is mostly Cretaceous granodiorite with some porphyry dikes.The fault core is highly permeable due to fracturing. The borehole water was sampled in 1996 and 2000 from the depth interval between 630 and 650 m, just below the fault core. The chemical and isotopic compositions were analyzed. Carbon and oxygen isotope ratios of carbonates from the fault core were analyzed to estimate the origin of fluid.The following conclusions were obtained. (1) The ionic and isotopic compositions of borehole water did not change from 1996 to 2000. They are mostly derived from local ground water as mentioned by Sato and Takahashi [Sato, T., Takahashi, M., 2000. Chemical and isotopic compositions of groundwater obtained from the Hirabayashi well. Geological Survey of Japan Interim Report No. EQ/00/1, 187–192.]. (2) Geochemical speciation revealed that the borehole water was derived from a relatively deep reservoir, which may be situated at a depth of 3 to 4 km where the temperature is about 80–90 °C. (3) The shallower part of the Nojima fault (shallower than the reservoir depth) has not been healed from the hydrological viewpoints 5 years after the event, in contrast to the rapid healing detected by S wave splitting [Tadokoro, K., Ando, M., 2002. Evidence for rapid fault healing derived from temporal changes in S wave splitting, Geophys. Res. Lett., 29, 10.1029/2001GL013644.]. (4) Precipitation of calcite from the present borehole water since drilling supports the idea of precipitation of some calcite in coseismic hydraulic fractures in the fault core [Boullier, A-M., Fujimoto, K., Ohtani, T., Roman-Ross, G., Lewin, E., Ito, H., Pezard, P., Ildefonse, B., 2004. Textural evidence for recent co-seismic circulation of fluids in the Nojima fault zone, Awaji Island, Japan., Tectonophysics, 378, 165–181.]. (5) Carbon and oxygen isotope ratios of calcite indicated that the meteoric water flux had been localized at the fault core. (6) A difference in the carbon isotope ratio between the footwall and the hanging wall suggests that the fault has been acted as a hydrologic barrier, although the permeability along the fault is still high.     

Temporal change in permeability of the Nojima fault zone by repeated water injection experiments

Three boreholes were drilled near the Nojima fault, which the 1995 Hyogoken–Nanbu earthquake occurred on. In order to research the properties and the healing process of the fault, water injection experiments were conducted every 3 years. In this report, we researched the permeability of the fault as a measurement of crack density or porosity of the fault zone. Pore water pressure changes in rock due to the water injections at one borehole were observed as discharge changes or groundwater level changes at the other borehole. Using numerical calculations, the permeability of the fault fracture zone was estimated for each experiment. The permeability has been decreasing as time passed, which is thought to show the fault healing process of the Nojima fault after the 1995 Hyogoken–Nanbu earthquake.     

Secular change of permeability in the fracture zone near the Nojima fault estimated using strain changes due to water injection experiments

Water injection experiments were performed in 1997, 2000 and 2003 at the 1800 m borehole near the fracture zone of the 1995 Hyogo-ken Nanbu earthquake. During these experiments, a contraction of about 10^− 8–10^− 7 was observed with three-component strainmeters at a bottom of the 800 m borehole, 70 m southwest of the 1800 m borehole. We estimated hydraulic properties of the fracture zone near the Nojima fault by using the strain data to investigate a healing of the fault during the postseismic stage. We calculated pore pressure changes due to the water injection using Darcy's equation and obtained strain changes due to the pore pressure changes as elastic deformations of the crust. The calculated strain changes have a nearly agreement with the observed strain changes. Hydraulic conductivity in 1997, 2000 and 2003 was determined to be 0.9 ± 0.2 × 10^− 6, 0.8 ± 0.2 × 10^− 6 and 0.4 ± 0.1 × 10^− 6 m/s, respectively. The reduced hydraulic conductivities in 2000 and 2003 suggest that the fractures had been healing.     

Probing fault zone heterogeneity on the Nojima fault: Constraints from zircon fission-track analysis of borehole samples

Fission-track (FT) thermochronologic analysis was performed on zircon separates from rocks in and around the Nojima fault, which was activated during the 1995 Kobe earthquake. Samples were collected from the University Group 500 m (UG-500) borehole and nearby outcrops. FT lengths in zircons from localities > 25 m away from the fault plane as well as one 0.1 m away from the fault in the footwall are characterized by concordant mean values of  10–11 μm and unimodal distributions with negative skewness, which showed no signs of appreciable reduction in FT length. In contrast, those adjacent (< 3 m) to the fault at depths on the hanging wall side showed significantly reduced mean track lengths of  6–8 μm and distributions having a peak around 6–7 μm with rather positive skewness. The former pattern is interpreted to reflect cooling through the zircon partial annealing zone (ZPAZ), without later, partial thermal overprints. The latter indicates substantial track shortening due probably to secondary heating by a thermal event(s) that locally perturbed the geothermal structure. Modeled zircon FT length and age data of partially annealed samples from the UG-500 borehole revealed a cooling episode in the ZPAZ that started at  4 Ma within  3 m from the fault plane, whereas those from the Geological Survey of Japan 750 m borehole record cooling started at  31–38 Ma within  25 m from the fault. On the basis of one-dimensional heat conduction modeling as well as the consistency between the degree of FT annealing and the degree of deformation/alteration of borehole rocks, these cooling ages in both boreholes are interpreted as consequences of ancient thermal overprints by heat transfer or dispersion via fluids in the fault zone. Together with the zircon FT data of a pseudotachylyte layer recently analyzed, it is suggested that the present Nojima fault system was reactivated in the Middle Quaternary from an ancient fault initiated at  56 Ma at mid-crustal depths. Also shown is a temporal/spatial variation in terms of the thermal anomalies recorded in the fault rocks, implying heterogeneity of hot fluid flows in the fault zone.     

Tectonic implications of damage zone-related fault-fracture networks revealed in drill core through the Nojima fault, Japan

1800 m of drill core through the Nojima fault zone, Japan, reveals subsidiary fault and fracture networks that developed in the fault zone that triggered the 1995 Ms 7.2 Kobe earthquake. The subsidiary fault zones contain a fault gouge of < 1 cm bounded by thin zones of foliated cataclasite or breccia. Fractures are filled with calcite veins, calcite-cemented breccias, clay, and iron-oxide and carbonate alternation of the granitic host rock. These features are typical of extensional fractures that form the conduit network for fluid flux close to a major fault zone. The zone of distributed deformation surrounding the main fault is 50 m in width, and the dip of the Nojima fault at > 1 km depth is 75°. The fault-fracture networks associated with the Nojima fault zone are coseismic and were filled with carbonate and fine-grained material during repeated seismic-related infiltration of the fault zone by carbonate-bearing subsurface water. This study shows that fault-related fracture networks plays an important role as fluid flow conduits within seismically active faults, and can change in character from zones of high permeability to low permeability due to cementation and/or pore collapse.     

Textural evidence for recent co-seismic circulation of fluids in the Nojima fault zone, Awaji island, Japan

The Hirabayashi borehole (Awaji Island, Japan) was drilled by the Geological Survey of Japan (GSJ) 1 year after the Hyogo-ken Nanbu (Kobe) earthquake (1995, M_JMA=7.2). This has enabled scientists to study the complete sequence of deformation across the active Nojima fault, from undeformed granodiorite to the fault core. In the fault core, different types of gouge and fractures have been observed and can be interpreted in terms of a complex history of faulting and fluid circulation. Above the fault core and within the hanging wall, compacted cataclasites and gouge are cut by fractures which show high apparent porosity and are filled by 5–50 μm euhedral and zoned siderite and ankerite crystals. These carbonate-filled fractures have been observed within a 5.5-m-wide zone above the fault, but are especially abundant in the vicinity (1 m) of the fault. The log-normal crystal size distributions of the siderite and ankerite suggest that they originated by decaying-rate nucleation accompanied by surface-controlled growth in a fluid saturated with respect to these carbonates. These carbonate-filled fractures are interpreted as the result of co-seismic hydraulic fracturing and upward circulation of fluids in the hanging wall of the fault, with the fast nucleation of carbonates attributed to a sudden fluid or CO₂ partial pressure drop due to fracturing. The fractures cut almost all visible structures at a thin section scale, although in some places, the original idiomorphic shape of carbonates is modified by a pressure-solution mechanism or the carbonate-filled fractures are cut and brecciated by very thin gouge zones; these features are attributed to low and high strain-rate mechanisms, respectively. The composition of the present-day groundwater is at near equilibrium or slightly oversaturated with respect to the siderite, calcite, dolomite and rhodochrosite. Taken together, this suggests that these fractures formed very late in the evolution of the fault zone, and may be induced by co-seismic hydraulic fracturing and circulation of a fluid with a similar composition to the present-day groundwater. They are therefore potentially related to recent earthquake activity (<1.2 Ma) on the Nojima fault.     

Apatite and zircon fission-track dating from the Hirabayashi-NIED borehole, Nojima Fault, Japan: Evidence for anomalous heating in fracture zones

Fission-track (FT) analysis using apatite and zircon was performed on samples from two fracture zones (FZ) at the depths of 1140 and 1310 m within the 1838 m borehole core penetrating the Ryoke Granitic Rocks in the Nojima Fault at Nojima-Hirabayashi, Awaji Island, Japan, drilled just after the 1995 Hyogo-ken Nanbu earthquake. Clear discordance in apatite and zircon FT age was found for two samples located at  2 m below the central part of each FZs where the presence of pseudotachylyte and/or fault gouge would predict the largest amount of slip. Asymmetric distribution was identified by discordant ages with respect to the central part of FZs. These very local discordant ages in the fault reflect thermal anomalies caused by secondary heating with an inferred maximum temperature in the region between apatite and zircon closure temperatures at a time post-48 Ma. As a source of the secondary heating, heat transfer or dispersion via geothermal fluids caused presumably the observed similarity in asymmetric distribution of discordant FT ages at two different FZs. Other samples yield concordant FT zircon and apatite ages and these indicate rapid cooling within the bounds of two closure temperatures of these minerals at  60 Ma of the Ryoke Granitic Rocks.     

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.     

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.     

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.     

Detection and characterization of hydraulically active fractures in a carbonate aquifer: results from self-potential, temperature and fluid electrical conductivity logging in the Combioula hydrothermal system in the southwestern Swiss Alps

A geophysical and geochemical study has been conducted in a fractured carbonate aquifer located at Combioula in the southwestern Swiss Alps with the objective to detect and characterize hydraulically active fractures along a 260-m-deep borehole. Hydrochemical analyses, borehole diameter, temperature and fluid electrical conductivity logging data were integrated in order to relate electrokinetic self-potential signals to groundwater flow inside the fracture network. The results show a generally good, albeit locally variable correlation of variations of the self-potential signals with variations in temperature, fluid electrical conductivity and borehole diameter. Together with the hydrochemical evidence, which was found to be critical for the interpretation of the self-potential data, these measurements not only made it possible to detect the hydraulically active fractures but also to characterize them as zones of fluid gain or fluid loss. The results complement the available information from the corresponding litholog and illustrate the potential of electrokinetic self-potential signals in conjunction with temperature, fluid electrical conductivity and hydrochemical analyses for the characterization of fractured aquifers, and thus may offer a perspective for an effective quantitative characterization of this increasingly important class of aquifers and geothermal reservoirs.     

Late Quaternary activity and dextral strike-slip movement on the Karakax Fault Zone, northwest Tibet

Field observations and interpretations of satellite images reveal that the westernmost segment of the Altyn Tagh Fault (called Karakax Fault Zone) striking WNW located in the northwestern margin of the Tibetan Plateau has distinctive geomorphic and tectonic features indicative of right-lateral strike-slip fault in the Late Quaternary. South-flowing gullies and N–S-trending ridges are systematically deflected and offset by up to ~ 1250 m, and Late Pleistocene–Holocene alluvial fans and small gullies that incise south-sloping fans record dextral offset up to ~ 150 m along the fault zone. Fault scarps developed on alluvial fans vary in height from 1 to 24 m. Riedel composite fabrics of foliated cataclastic rocks including cataclasite and fault gouge developed in the shear zone indicate a principal right-lateral shear sense with a thrust component. Based on offset Late Quaternary alluvial fans, ¹⁴C ages and composite fabrics of cataclastic fault rocks, it is inferred that the average right-lateral strike-slip rate along the Karakax Fault Zone is ~ 9 mm/a in the Late Quaternary, with a vertical component of ~ 2 mm/a, and that a M 7.5 morphogenic earthquake occurred along this fault in 1902. We suggest that right-lateral slip in the Late Quaternary along the WNW-trending Karakax Fault Zone is caused by escape tectonics that accommodate north–south shortening of the western Tibetan Plateau due to ongoing northward penetration of the Indian plate into the Eurasian plate.     

Ultrashallow seismic experiment on a trenched section of the Chelunpu fault zone, Taiwan

Ultrashallow P-wave seismic reflection experiments were conducted at a model test site and in a trenched shallow fault zone along the Chelunpu fault line. The field layout was designed to have the shallowest undistorted reflection from about 1 m depth with 0.5 m vertical resolution. The smallest group interval tested in this study was 0.05 m with a 0.25 ms sample interval, which can avoid spatial aliasing of ground roll if the target is very shallow and the velocities are low. Data processing was designed to be simple but consistent. As the ultrashallow reflections may be contaminated with high-amplitude coherent noise in many aspects, first break muting and surgical muting were performed on each file as detailed as possible, and f–k filtering was applied mainly for the purpose of attenuating the aliasing energy and back-scattered noise. Data acquired in this study show that the low P-wave velocities (< 200 m/s) and high dominant frequencies (120–200 Hz) of near-surface layers may have a potential vertical resolution of 0.4 m or even better.Comparing the test profile with the ground-penetrating radar (GPR) control profile of the same test site and correlating the results obtained from the study site with those of the geologic cross-section of the trench, this experiment demonstrates the possibility of using seismic methods in investigating shallow structures at depths of less than a few meters with vertical resolution comparable to the GPR technique.     

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.     

Fracture-zone conditions on a recently active fault: insights from mineralogical and geochemical analyses of the Hirabayashi NIED drill core on the Nojima fault, southwest Japan, which ruptured in the 1995 Kobe earthquake

An 1800-m-deep borehole into the Nojima fault zone was drilled at Nojima-Hirabayashi, Japan, after the 1995 Hyogo-ken Nanbu (Kobe) earthquake. Three possible fracture zones were detected at depths of about 1140, 1300, and 1800 m. To assess these fracture zones in this recently active fault, we analyzed the distributions of fault rocks, minerals, and chemical elements in these zones. The central fault plane in the shallowest fracture zone was identified by foliated blue-gray gouge at a depth of 1140 m. The degree of fracturing was evidently greater in the hanging wall than in the footwall. Minerals detected in this zone were quartz, orthoclase, plagioclase, and biotite, as in the parent rock (granodiorite), and also kaolinite, smectite, laumontite, stilbite, calcite, ankerite, and siderite, which are related to hydrothermal alteration. Biotite was absent in both the hanging wall and footwall across the central fault plane, but it was absent over a greater distance from the central fault plane in the hanging wall than in the footwall. Major element compositions across this zone suggested that hydrothermal alteration minerals such as kaolinite and smectite occurred across the central fault plane for a greater distance in the hanging wall than in the footwall. Similarly, H₂O+ and CO₂ had higher concentrations in the hanging wall than in the footwall. This asymmetrical distribution pattern is probably due to the greater degree of wall–rock fracturing and associated alteration in the hanging wall. We attributed the characteristics of this zone to fault activity and fluid–rock interactions. We analyzed the other fracture zones along this fault in the same way. In the fracture zone at about 1300 m depth, we detected the same kinds of hydrothermal alteration minerals as in the shallower zone, but they were in fewer samples. We detected relatively little H₂O+ and CO₂, and little evidence for movement of the major chemical elements, indicating little past fluid–rock interaction. In the fracture zone at about 1800 m depth, H₂O+ and CO₂ were very enriched throughout the interval, as in the fracture zone at about 1140 m depth. However, smectite was absent and chlorite was present, indicating the occurrence of chloritization, which requires a temperature of more than 200 °C. Only smectite can form under the present conditions in these fracture zones. The chloritization probably occurred in the past when the fracture zone was deeper than it is now. These observations suggest that among the three fracture zones, that at about 1140 m depth was the most activated at the time of the 1995 Hyogo-ken Nanbu (Kobe) earthquake.     

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.     

An empirical model for estimating hydraulic conductivity of highly disturbed clastic sedimentary rocks in Taiwan

This paper presents the establishment of an empirical HC model for estimating rock mass hydraulic conductivity of highly disturbed clastic sedimentary rocks in Taiwan using high-resolution borehole acoustic televiewer and double packer hydraulic tests. Four geological parameters including rock quality designation (RQD), depth index (DI), gouge content designation (GCD), and lithology permeability index (LPI) were adopted for establishing the empirical HC model. To verify rationality of the proposed HC model, 22 in-situ hydraulic tests were carried out to measure the hydraulic conductivity of the highly disturbed clastic sedimentary rocks in three boreholes at two different locations in Taiwan. Besides, the model verification using another borehole data with four additional in-situ hydraulic tests from similar clastic sedimentary rocks was also conducted to further verify the feasibility of the proposed empirical HC model. The field results indicated that the rock mass in the study area has a conductivity between the order of 10^− 10 m/s and 10^− 6 m/s at the depth between 34 m and 275 m below ground surface. Results demonstrate that the empirical HC model may provide a useful tool to predict hydraulic conductivity of the highly disturbed clastic sedimentary rocks in Taiwan based on measured HC-values.     

Touhuanping Fault, an active wrench fault within fold-and-thrust belt in northwestern Taiwan, documented by spatial analysis of fluvial terraces

This study aims at the recent activity and development of an active wrench fault, the Touhuanping Fault in northwestern Taiwan. Northwestern Taiwan has been proposed in a current situation between the mature to waning collision in terms of tectonic evolution. The main drainage in this area, the Chungkang River, flows close to the trace of the fault mentioned above. We examined various types of deformation of fluvial terraces along the Chungkang River as a key to understanding the nature and rate of the late Quaternary tectonics. The E–W trending Touhuanping Fault has long been mapped as a geological boundary fault, but its recent activity was suspected. Field survey revealed that its late Quaternary activity is recorded in the offset fluvial terraces. Our result shows dextral slip and vertical offset with upthrown side on the south, and activated at least twice since the emergence of terrace 4 (older terrace 3 with OSL date of ca. 80 ka). Total amount of offset recorded in the Touhuanping terrace sequence is 15 m for dextral and 10 m for vertical offset. Estimated recurrence time of earthquake rupture may be a few tens of thousand years. Uplift on the upthrown side of the Touhuanping Fault also resulted in the formation of drowned valleys which were graded to terrace 4. Other deformation features, such as back-tilting, westward warping, and a range-facing straight scarp, were also identified. A second-order anticline roughly parallel to the Touhuanping Fault is suggested to be the origin of the northward tilting on terrace 3; it could have resulted from a flower structure on the Touhuanping Fault at shallow depth. This may demonstrate that the buried segment of the Touhuanping Fault has also been active since 80 ka. In the northern study area, the westward warping at terrace 2 probably represents late Quaternary activity of another NE–SW trending Hsincheng Fault.     

Seismicity, deformation and seismic hazard in the western rift of Corinth: New insights from the Corinth Rift Laboratory (CRL)

This paper presents the main recent results obtained by the seismological and geophysical monitoring arrays in operation in the rift of Corinth, Greece. The Corinth Rift Laboratory (CRL) is set up near the western end of the rift, where instrumental seismicity and strain rate is highest. The seismicity is clustered between 5 and 10 km, defining an active layer, gently dipping north, on which the main normal faults, mostly dipping north, are rooting. It may be interpreted as a detachment zone, possibly related to the Phyllade thrust nappe. Young, active normal faults connecting the Aigion to the Psathopyrgos faults seem to control the spatial distribution of the microseismicity. This seismic activity is interpreted as a seismic creep from GPS measurements, which shows evidence for fast continuous slip on the deepest part on the detachment zone. Offshore, either the shallowest part of the faults is creeping, or the strain is relaxed in the shallow sediments, as inferred from the large NS strain gradient reported by GPS. The predicted subsidence of the central part of the rift is well fitted by the new continuous GPS measurements. The location of shallow earthquakes (between 5 and 3.5 km in depth) recorded on the on-shore Helike and Aigion faults are compatible with 50° and 60° mean dip angles, respectively. The offshore faults also show indirect evidence for high dip angles. This strongly differs from the low dip values reported for active faults more to the east of the rift, suggesting a significant structural or rheological change, possibly related to the hypothetical presence of the Phyllade nappe. Large seismic swarms, lasting weeks to months, seem to activate recent synrift as well as pre-rift faults. Most of the faults of the investigated area are in their latest part of cycle, so that the probability of at least one moderate to large earthquake (M = 6 to 6.7) is very high within a few decades. Furthermore, the region west to Aigion is likely to be in an accelerated state of extension, possibly 2 to 3 times its mean interseismic value. High resolution strain measurement, with a borehole dilatometer and long base hydrostatic tiltmeters, started end of 2002. A transient strain has been recorded by the dilatometer, lasting one hour, coincident with a local magnitude 3.7 earthquake. It is most probably associated with a slow slip event of magnitude around 5 ± 0.5. The pore pressure data from the 1 km deep AIG10 borehole, crossing the Aigion fault at depth, shows a 1 MPa overpressure and a large sensitivity to crustal strain changes.