10Be surface exposure dating reveals unexpected high deformation rates in the central Andean wedge interior

Understanding deformation associated with active thrust wedges is essential to evaluate seismic hazards. In this study, we investigate the spatial distribution, timing, and rates of deformation across the central Andean backarc of Argentina, where deformation and shortening have been assumed to occur within a narrow wedge‐front zone. The combined results of our geomorphic mapping, topographic surveying, and ¹⁰Be exposure dating demonstrate that fault activity related to the growth of the Andean orogenic wedge is not only limited to a narrow thrust front zone but also occurs in the Andean wedge interior. There, internal structures with deformation rates of ca. 1.3 mm/year have been active during the last ~140 ka. Widely distributed deformation implies that seismic hazards may have been underestimated in the internal part of the Andean orogen.     

Seismic slip rates,potential subsurface rupture areas and seismic potential of the Vienna Basin Transfer Fault

The Vienna Basin Transfer Fault (VBTF) is a slow active fault with moderate seismicity (I _max~8–9, M _max~5.7) passing through the most vulnerable regions of Austria and Slovakia. We use different data to constrain the seismic potential of the VBTF including slip values computed from the seismic energy release during the 20th century, geological data on fault segmentation and a depth-extrapolated 3-D model of a generalized fault surface, which is used to define potential rupture zones. The seismic slip of the VBTF as a whole is in the range of 0.22–0.31 mm/year for a seismogenic fault thickness of 8 km. Seismic slip rates for individual segments vary from 0.00 to 0.77 mm/year. Comparing these data to geologically and GPS-derived slip velocities (>1 mm/year) proofs that the fault yields a significant seismic slip deficit. Segments of the fault with high seismic slip contrast from segments with no slip representing locked segments. Fault surfaces of segments within the seismogenic zone (4–14 km depth) vary from 55 to 400 km². Empirical scaling relations show that these segments are sufficiently large to explain both, earthquakes observed in the last centuries, and the 4th century Carnuntum earthquake, for which archeo-seismological data suggest a magnitude of M ≥ 6. Based on the combination of all data (incomplete earthquake catalog, seismic slip deficits, locked segments, potential rupture areas, indications of strong pre-catalog earthquakes) we argue, that the maximum credible earthquake for the VBTF is in the range M _max = 6.0–6.8, significantly larger than the magnitude of the strongest recorded events (M = 5.7).     

3D modeling of earthquake cycles of the Xianshuihe fault,southwestern China

We perform 3D modeling of earthquake generation of the Xianshuihe fault, southwestern China, which is a highly active strike-slip fault with a length of about 350 km, in order to understand earthquake cycles and segmentations for a long-term forecasting and earthquake nucleation process for a short-term forecasting. Historical earthquake data over the last 300 years indicates repeated periods of seismic activity, and migration of large earthquake along the fault during active seismic periods. To develop the 3D model of earthquake cycles along the Xianshuihe fault, we use a rate- and state-dependent friction law. After analyzing the result, we find that the earthquakes occur in the reoccurrence intervals of 400–500 years. Simulation result of slip velocity distribution along the fault at the depth of 10 km during 2694 years along the Xianshuihe fault indicates that since the third earthquake cycle, the fault has been divided into 3 parts. Some earthquake ruptures terminate at the bending part of the fault line, which may means the shape of the fault line controls how earthquake ruptures. The change of slip velocity and displacement at 10 km depth is more tremendous than the change of the shallow and deep part of the fault and the largest slip velocity occurs at the depth of 10 km which is the exact depth of the seismic zone where fast rupture occurs.     

Salt-influenced normal fault growth and forced folding: The Stavanger Fault System,North Sea

Displacement ratio (D_r) is the ratio between salt thickness (T_v) and sub-salt normal fault displacement (D) (D_r = T_v/D), and it is typically used to predict the degree of geometric and kinematic linkage between sub- and supra-salt fault populations, and the overall structural style in salt-influenced extensional settings. However, we currently lack natural examples of how D_r and the underlying geological controls vary, and how these may control the three-dimensional geometry and evolution of salt-influenced normal fault systems. Furthermore, it is currently unknown if kinematic coherence in salt-influenced extensional settings can be maintained over relatively long length-scales (10¹–10³ m) and for relatively long timeframes, and how this may impact the growth and geometry of large-throw (>500 m), salt-influenced normal fault systems. In this paper we use a 3600 km², high-quality 3D seismic reflection dataset and borehole data from the Stavanger Fault System (SFS), Egersund Basin, eastern North Sea Basin to investigate; (i) how pre-rift salt thickness (T_v) and sub-salt fault throw (T) control the structural style and evolution of a basin-bounding, salt-influenced normal fault system; and (ii) the role salt plays in maintaining kinematic coherence in normal fault systems. We demonstrate that; (i) pre-rift salt distribution (T_v), specifically its presence in the proto-footwall (i.e., when T_v > 0), is the primary control on partitioning of faulting and (forced) folding along the fault system, and the style of linkage (i.e., hard- or soft-linkage) between sub- and supra-salt fault populations; and (ii) sub- and supra-salt fault populations represent brittle elements of a single, geometrically and kinematically coherent structure, the geometry and evolution of which is related to the ductile translation of strain on a scale (up to 8 km) and duration (c. 65 Myr) that believe is significantly greater and longer than previously documented.     

Do faults preserve a record of seismic slip: A second opinion

Exhumed fault zones offer insights into deformation processes associated with earthquakes in unparalleled spatial resolution; however it can be difficult to differentiate seismic slip from slow or aseismic slip based on evidence in the rock record. Fifteen years ago, Cowan (1999) defined the attributes of earthquake slip that might be preserved in the rock record, and he identified pseudotachylyte as the only reliable indicator of past earthquakes found in ancient faults. This assertion was based on models of frictional heat production (Sibson, 1975, 1986) providing evidence for fast slip. Significant progress in fault rock studies has revealed a range of reaction products which can be used to detect frictional heating at peak temperatures less than the melt temperature of the rock. In addition, features formed under extreme transient stress conditions associated with the propagating tip of an earthquake rupture can now be recognized in the rock record, and are also uniquely seismic. Thus, pseudotachylyte is no longer the only indicator of fossilized earthquake ruptures.We review the criteria for seismic slip defined by Cowan (1999), and we determine that they are too narrow. Fault slip at rates in the range 10⁻⁴−10¹ m/s is almost certainly dynamic. This implies that features reproduced in experiments at rates as low as 10⁻⁴ m/s may be indicators of seismic slip. We conclude with a summary of the rock record of seismic slip, and lay out the current challenges in the field of earthquake geology.     

Geogas transport in fractured hard rock – Correlations with mining seismicity at 3.54 km depth, TauTona gold mine, South Africa

An on-site gas monitoring study has been conducted in the framework of an earthquake laboratory (The International NELSAM–DAFGAS projects) at the TauTona gold mine, South Africa. Five boreholes up to 60 m long were drilled at 3.54 km depth into the highly fractured Pretorius Fault Zone and instruments for chemical and seismic monitoring installed therein. Over the span of 4 years sensitive gas monitoring devices were continuously improved to enable the direct observation of geogas concentration variations in the DAFGAS borehole. The major gas concentrations are constant and air-like with about 78% N₂, 21% O₂, 1% Ar. The geogas components CO₂, CH₄, He and H₂ show the most interesting trends and variations on the minute-by-minute basis and significantly correlate with seismic data, while the ²²²Rn activity remains constant. Time series and cross correlation analysis allow the identification of different gas components (geogas and tunnel air) and the identification of two processes influencing the borehole gas composition: (1) pumping-induced tunnel air breakthrough through networks of initially water-saturated fault fractures; and (2) seismicity induced permeability enhancement of fault fractures to above ∼5 × 10^-10 m². The current set-up of the gas monitoring system is sensitive enough to quantify the resulting geogas transport during periods of intense blasting activities (including recorded blasts with seismic moment ?1 × 10⁹ Nm, located within 1000 m of the cubby) and, it is suggested, also during induced earthquakes, a final goal of the project.     

Seismotectonics of the 2008 and 2009 Qaidam Earthquakes and its Implication for Regional Tectonics

Three magnitude >6 earthquakes struck Qaidam, Qinghai province, China, in November 10th 2008, August 28th and 31st 2009 respectively. The Zongwulongshan fault has often been designated as the active seismogenic structure, although it is at odd with the data. Our continuous GPS station (CGPS), the Xiao Qaidam station, located in the north of the Qaidam basin, is less than 30 km to the southwest of the 2008 earthquake. This CGPS station recorded the near field co-seismic deformation. Here we analyzed the co-seismic dislocation based on the GPS time series and the rupture processes from focal mechanism for the three earthquakes. The aftershocks were relocated to constrain the spatial characteristics of the 2008 and 2009 Qaidam earthquakes. Field geological and geomorphological investigation and interpretation of satellite images show that the Xitieshan fault and Zongwulongshan fault were activated as left lateral thrust during the late Quaternary. Evidence of folding can also be identified. Integrated analyses based on our data and the regional tectonic environment show that the Xitieshan fault is the fault responsible for the 2008 Qaidam earthquake, which is a low dip angle thrust with left lateral strike slip. The Zongwulongshan fault is the seismogenic fault of the 2009 earthquakes, which is a south dipping back thrust of the northern marginal thrust system of the Qaidam basin. Folding takes a significant part of the deformation in the northern marginal thrust system of the Qaidam basin, dominating the contemporary structure style of the northern margin of the Qaidam basin and Qilianshan tectonic system. In this region, this fault and fold system dominates the earthquake activities with frequent small magnitude earthquakes.     

Focal parameters of seismic sources during the 1981 and 1983 eruption at Mt. Etna volcano (Sicily,Italy)

In this study 50 seismic events, preceding and accompanying the eruptions occurring in 1981 and 1983, have been considered. Seismic moments, fault radii, stress drops and seismic energies have been calculated using Brune’s model (J Geophys Res 75:4997–5009, 1970; J Geophys Res 76:5002, 1971); site, anelastic attenuation along the propagation path, geometrical spreading and interaction with the free surface effects are taken into account. For each event we have also estimated the equivalent Wood–Anderson magnitude (MWAeq) (Scherbaum and Stoll in Bull Seism Soc Am 73:1321–1343, 1983); relations among all these source parameters have been determined. Furthermore, the hypothesis of self-similarity (Aki in J Geophys Res 72:1217–1231, 1967) is not verified for events with seismic moments <10¹² N-m: in fact the relationship between log-stress drop and log-moment is linear up to a moment of 10¹² N-m (events of 1981 eruption), while for higher moments (events of 1983 eruption) the slope of the regression line is not significantly different from zero. We suppose that such a behaviour is related to a heterogeneous medium with barriers on the faults. Finally, the main conclusion is that eruptions of 1981 and 1983 differ from one another both in eruptive and seismic aspects; analysis of seismic energies indicates an increase in Mt. Etna’s activity, confirmed by studies performed on the following lateral eruption of 1991–1993 (Patanè et al. in Bull Volcanol 47:941–952, 1995), occurring on the same structural trend.     

The Mogangling giant landslide triggered by the 1786 Moxi M 7.75 earthquake,China

A good understanding of seismic giant landslides could provide favourable guidance for seismic stability evaluation of nearby slopes. Here, an excellent example of a catastrophic seismic landslide named the Mogangling giant landslide (MGL), located upstream along the Dadu River and triggered by the 1786 Moxi M 7.75 earthquake, is analysed for its deposit characteristics, failure mechanism and dammed lake. The MGL, with a volume of approximately 4500?×?10⁴ m³, 450 m long and 1000 m wide, blocked the Dadu River completely and caused over 100 000 deaths when the landslide dam broke. The MGL occurred on the upper part of a narrow granite ridge; a potentially unstable wedge-shaped rock mass was separated from the remaining massif by unloading fissures and an active fault (Detuo fault) that just crossed the slope foot. The Moxi earthquake coupled with strong site amplification triggered the MGL, which blocked the Dadu River; the elevation of the dam crest was approximately 130 m higher than the present river level. The dammed lake had a volume of approximately 9.504?×?10⁸ m³, an area of 19.91 km² and a length of approximately 31 km; the peak flow of the outburst flood was larger than 7100 m³/s. After hundreds of years of concave bank erosion, the deposit is divided into the right bank deposit (main deposit) and left bank deposit (residual deposit).

     

Evidence for Holocene activity of the Yilan-Yitong fault,northeastern section of the Tan-Lu fault zone in Northeast China

The Tan-Lu fault zone (TLFZ) is the largest of the major faults in eastern China. Many strong earthquakes have occurred on its section in North China, but no quake greater than M ⩾ 6 has been documented in history at its northeastern section, the Yilan-Yitong fault (YYF) in Northeast China. It is usually considered that this fault has been inactive since late Quaternary and incapable of generating moderate-sized quakes. This conclusion is, however, questioned by our recent work based on high-resolution satellite image interpretation and field investigation. We found a 70-km-long surface scarp near Fangzheng county in Heilongjiang province (HLJP) and a 20-km-long scarp near Shulan county in Jilin province (JLP), and both are associated with the YYF. The trenches across these two scarps reveal a ¹⁴C displacement date of 1730 ± 40 years BP at Fangzheng and of 4410 ± 30 years BP at Shulan. The dextral offsets of the Songhua River and Second Songhua River and nearly horizontal fault striations indicate that the new activity of the YYF has been dominated by dextral strike slipping with a normal component. These new data suggest that, at least for partial sections, the YYF has been active since the Holocene, implying a potential seismic hazard. However, current quake-protection standards in this region are very low due to the previous view that the YYF fault has not been active since the late Quaternary. If an M ⩾ 7 quake takes place on this fault, it will be a devastating event. Therefore, it is necessary to conduct a detailed study on the whole YYF and to reassess its future seismic risk.     

Co-seismic displacements associated to the Molise (Southern Italy) earthquake sequence of October–November 2002 inferred from GPS measurements

The 2002 earthquake sequence of October 31 and November 1 (main shocks Mw = 5.7) struck an area of the Molise region in Southern Italy. In this paper we analyzed the co-seismic deformation related to the Molise seismic sequence, inferred from GPS data collected before and after the earthquake, that ruptured a rather deep portion of crust releasing a moderate amount of seismic energy with no surface rupture. The GPS data have been reduced using two different processing strategies and softwares (Bernese and GIPSY) to have an increased control over the result accuracy, since the expected surface displacements induced by the Molise earthquake are in the order of the GPS reliability. The surface deformations obtained from the two approaches are statistically equivalent and show a displacement field consistent with the expected deformation mechanism and with no rupture at the surface. In order to relate this observation with the seismic source, an elastic modeling of fault dislocation rupture has been performed using seismological parameters as constraints to the model input and comparing calculated surface displacements with the observed ones. The sum of the seismic moments (8.9 × 10¹⁷ Nm) of the two main events have been used as a constraint for the size and amount of slip on the model fault while its geometry has been constrained using the focal mechanisms and aftershocks locations. Since the two main shocks exhibit the same fault parameters (strike of the plane, dip and co-seismic slip), we modelled a single square fault, size of 15 km × 15 km, assumed to accommodate the whole rupture of both events of the seismic sequence. A vertical E–W trending fault (strike = 266°) has been modeled, with a horizontal slip of 120 mm. Sensitivity tests have been performed to infer the slip distribution at depth. The comparison between GPS observations and displacement vectors predicted by the dislocation model is consistent with a source fault placed between 5 and 20 km of depth with a constant pure right-lateral strike-slip in agreement with fault slip distribution analyses using seismological information. The GPS strain field obtained doesn't require a geodetic moment release larger than the one inferred from the seismological information ruling out significant post-seismic deformation or geodetic deformation released at frequencies not detectable by seismic instruments. The Molise sequence has a critical seismotectonic significance because it occurred in an area where no historical seismicity or seismogenic faults are reported. The focal location of the sequence and the strike-slip kinematics of main shocks allow to distinguish it from the shallower and extensional seismicity of the southern Apennines being more likely related to the decoupling of the southern Adriatic block from the northern one.     

Background seismicity in the Central Apennines of Italy: The Abruzzo region case study

We investigate background seismic activity of the Abruzzo region, a 5000 km² area located within the Central Apennines of Italy, where in the past 600 years at least 5 large earthquakes (I = XI–X) have occurred.Between April 2003 and September 2004, a dense temporary seismic network composed of 30 digital three-component seismic stations recorded 850 earthquakes with 0.9 < M_L < 3.7. We present earthquake locations and focal mechanisms obtained by standard procedures and an optimized velocity model computed with a search technique based on genetic algorithms.The seismicity occurs at a low and constant rate of  2.6 e− 04 events/daykm² and is sparsely distributed within the first 15 km of the crust. Minor increases in the seismicity rate are related to the occurrence of small and localised seismic sequences that occur at the tip of major active normal faults along secondary structures.We observe that during the 16 months of study period, the Fucino fault system responsible for the 1915 Fucino earthquake (M_S = 7.0), and the major normal faults of the area, did not produce significant seismic activity.Fault plane solutions evaluated using P-wave polarity data show the predominance of normal faulting mechanisms ( 55%) with NE-trending direction of extension coherent with the regional stress field active in this sector of the Apennines. Around 27% of the focal solutions have pure strike–slip mechanisms and the rest shows transtensional faulting mechanisms that mainly characterise the kinematics of the secondary structures activated by the small sequences.We hypothesize that the largest known NW-trending normal faults are presently locked and we propose that in the case of activation, the secondary structures located at their tips may act as transfer faults accommodating a minor part of the extensional deformation with strike–slip motion.     

Fault growth and interactions in a multiphase rift fault network: Horda Platform,Norwegian North Sea

Physical models predict that multiphase rifts that experience a change in extension direction between stretching phases will typically develop non-colinear normal fault sets. Furthermore, multiphase rifts will display a greater frequency and range of styles of fault interactions than single-phase rifts. Although these physical models have yielded useful information on the evolution of fault networks in map view, the true 3D geometry of the faults and associated interactions are poorly understood. Here, we use an integrated 3D seismic reflection and borehole dataset to examine a range of fault interactions that occur in a natural multiphase fault network in the northern Horda Platform, northern North Sea. In particular we aim to: i) determine the range of styles of fault interaction that occur between non-colinear faults; ii) examine the typical geometries and throw patterns associated with each of these different styles; and iii) highlight the differences between single-phase and multiphase rift fault networks. Our study focuses on a ca. 350 km² region around the >60 km long, N–S-striking Tusse Fault, a normal fault system that was active in the Permian–Triassic and again in the Late Jurassic-to-Early Cretaceous. The Tusse Fault is one of a series of large (>1500 m throw) N–S-striking faults forming part of the northern Horda Platform fault network, which includes numerous smaller (2–10 km long), lower throw (<100 m), predominantly NW–SE-striking faults that were only active during the Late Jurassic to Early Cretaceous. We examine how the 2^nd-stage NW–SE-striking faults grew, interacted and linked with the N–S-striking Tusse Fault, documenting a range of interaction styles including mechanical and kinematic isolation, abutment, retardation and reactivated relays. Our results demonstrate that: i) isolated, and abutting interactions are the most common fault interaction styles in the northern Horda Platform; ii) pre-existing faults can act as sites of nucleation for 2^nd-stage faults or may form mechanical barriers to propagation; iii) the throw distribution on reactivated 1^st-stage faults will be modified in a predictable manner if they are intersected or influenced by 2^nd-stage faults; iv) sites of fault linkage and relay-breaching associated with the first phase of extension can act as preferential nucleation sites for 2^nd-stage faults; and v) the development of fault intersections is a dynamic process, involving the gradual transition from one style to another.     

The last major earthquakes along the Magallanes–Fagnano fault system recorded by disturbed trees (Tierra del Fuego,South America)

Forests situated above active fault zones may record hillslope evolution, thus holding information about recent seismic events. Lenga trees (Nothofagus pumilio) extend across the Magallanes–Fagnano fault system (MFFS), the active transform boundary between the South American and Scotia plates. Coseismic surface ruptures along the fault scarp tilt trees located uphill. During the interseismic period, tree growth curves the trunks. Annual tree rings from the study area show abrupt changes from concentric to asymmetric, allowing the timing of major historical earthquakes to be established. In this case, tree‐ring analysis suggests rupture on the MFFS fault scarp in 1883 ± 5 and 1941 ± 10, coinciding with the February 1, 1879 (Modified Mercalli Scale, VI) and the December 17, 1949 (Ms 7.8) earthquakes in Tierra del Fuego. Our results provide evidence that this fault system was the source of these earthquakes, which has implications for seismic hazard in the study region.     

Frictional properties of saponite-rich gouge from a serpentinite-bearing fault zone along the Gokasho-Arashima Tectonic Line,central Japan

We studied a serpentinite-bearing fault zone in Gokasho-Arashima Tectonic Line, Mie Prefecture, central Japan, characterizing its internal structures, mineral assemblage, permeability, and frictional properties. The fault core situated between the serpentinite breccia and the adjacent sedimentary rocks is characterized by a zone locally altered to saponite. The clayey gouge layer separates fault rocks of serpentinite origin containing talc and tremolite from fault rocks of sedimentary origin containing chlorite but no quartz. The minerals that formed within the fault are the products of metasomatic reaction between the serpentinite and the siliceous rocks. Permeability measurements show that serpentinite breccia and fault gouge have permeability of 10⁻¹⁴–10⁻¹⁷ m² and 10⁻¹⁵–10⁻¹⁸ m², respectively, at 5–120 MPa confining pressure. Frictional coefficient of the saponite-rich clayey fault gouge ranged between 0.20 and 0.35 under room-dry condition, but was reduced to 0.06–0.12 when saturated with water. The velocity dependence of friction was strongly positive, mostly ranging between 0.005 and 0.006 in terms of a–b values. The governing friction law is not constrained yet, but we find that the saponite-rich gouge possesses an evolutional behavior in the opposite direction to that suggested by the rate and state friction law, in addition to its direct velocity dependence.     

Palaeopermeability structure within fault-damage zones: A snap-shot from microfracture analyses in a strike-slip system

Understanding fault zone permeability and its spatial distribution allows the assessment of fluid-migration leading to precipitation of hydrothermal minerals. This work is aimed at unraveling the conditions and distribution of fluid transport properties in fault zones based on hydrothermally filled microfractures, which reflect the ‘‘frozen-in’’ instantaneous advective hydrothermal activity and record palaeopermeability conditions of the fault-fracture system. We studied the Jorgillo Fault, an exposed 20 km long, left-lateral strike-slip fault, which juxtaposes Jurassic gabbro against metadiorite belonging to the Atacama Fault System in northern Chile. Tracings of microfracture networks of 19 oriented thin sections from a 400 m long transect across the main fault trace was carried out to estimate the hydraulic properties of the low-strain fault damagezone, adjacent to the high-strain fault core, by assuming penny-shaped microfractures of constant radius and aperture within an anisotropic fracture system. Palaeopermeability values of 9.1*10⁻¹¹ to 3.2*10⁻¹³ m² in the gabbro and of 5.0*10⁻¹⁰ to 1.2*10⁻¹³ m² in the metadiorite were determined, both decreasing perpendicularly away from the fault core. Fracture porosity values range from 40.00% to 0.28%. The Jorgillo Fault has acted as a left-lateral dilational fault-bend, generating large-scale dilation sites north of the JF during co-seismic activity.     

Fault zone architecture and fluid flow in interlayered basaltic volcaniclastic-crystalline sequences

Faults in continental flood basalt sequences potentially control subsurface fluid flow. We present field and microstructural observations from fault zones cutting interlayered basaltic volcaniclastic-crystalline sequences within the North Atlantic Igneous Province. Fractures likely initiate within lava units, before linking through the volcaniclastic units. Through-going faults show refraction, with subvertical faults in the lavas joined to variably inclined faults in the volcaniclastic layers. At >1.0 m displacement, volcaniclastic units are progressively dragged into the fault plane forming a smear. Volcaniclastic sandstones deform by flow. Claystones fracture, and are incorporated into smears as breccia. Experimentally measured host and fault rock sample permeabilities, at aquifer to reservoir pressures (i.e., 10–90 MPa; ∼0.3–3.0 km depth) show fault rocks from low displacement faults have relatively low permeability (10⁻¹⁷–10⁻²⁰ m²); fault rocks from higher displacement structures have comparatively high permeability (10⁻¹⁵–10⁻¹⁷ m²). Our observations suggest that permeability is determined by the opposing influences of clay mineralization, which decreases permeability, versus the development of interconnected, higher permeability zeolite veins. Brecciation and the formation of zeolite vein networks within claystone smears results in high permeability. Zeolite veins in volcaniclastic units form poorly-connected, spaced sets, parallel to the slip plane, hence sequence permeability remains low.     

Physical properties of fault zones within a granite body: Example of the Soultz-sous-Forêts geothermal site

In EGS projects, fault zones are considered as the structures controlling deep flow at the reservoir scale. Using a large set of petrophysical properties (porosity, density, permeability, thermal conductivity [TC]) measured on cores collected along the EPS-1 borehole, a model of fault zone is proposed to describe them. A fault zone is a complex structure, showing different parts with different kinds of deformations and/or materials that could explain chemical and physical processes observed during fluid-rock interactions. The different parts composing the fault zone are: (1) the fault core or gauge zone; (2) the damage zone; (3) and the protolith. They are usually heterogeneous and show different physical properties. The damage zone is a potential high permeability channel and could become the main pathway for fluids if secondary minerals seal the fault core. Porosity is the lowest within the protolith, between 0.5 and 1%, but can go up to 15% in the fault zone. Permeability ranges from 10^?20 m² in the fresh granite to, at least, 10^?15 m² in the fault core, and TC ranges from 2.5 W K^?1m^?1 to 3.7 W K^?1m^?1. Finally, variations in specific surface are set over two orders of magnitude. If the lowest values usually characterize the fresh granite far from fault zones, physical properties could show variations spread over their whole respective ranges within these fault zones.     

Clay–clast aggregates in fault gouge: An unequivocal indicator of seismic faulting at shallow depths?

A common problem encountered in studies of gouge-bearing natural faults is the difficulty of ascertaining whether the observed gouge was sheared seismically or aseismically; this problem arises because of the scarcity of indicators of fault slip rates for gouge. Recently, clay–clast aggregates (CCAs; a CCA comprises a clastic core mantled by a rim of ultrafine particles) were proposed as a possible indicator of seismic slip in gouge, on the basis of shear experiments on gouge at seismic slip rates. To examine the processes and conditions of CCA formation, we conducted rotary shear experiments on quartz and quartz–bentonite gouges under normal stresses (0.3–3.0 MPa) and slip rates (0.0005–1.3 m s⁻¹), and in both room-humidity (room-dry) and water-saturated (wet) conditions. We found that CCAs could be produced in room-dry gouges even at the lowest slip rates, which are considerably slower than actual seismic slip rates. This finding demonstrates that thermal pressurization and fluidization at elevated temperature during seismic slip are not necessarily needed for the formation of CCAs, contrary to previous views. Given the occurrence of CCAs over a wide range of slip rates, we suggest that the presence of CCAs is not an unequivocal indicator of fault slip at seismic slip rates.     

Earthquake generation in Cyprus revealed by the evolving stress field

Strong earthquake occurrence (M ≥ 6.0) onshore and offshore the Cyprus Island constitutes significant seismic hazard because they occur close to populated areas. Seismicity is weak south of the Island along the Cyprean Arc and strong events are aligned along the Paphos transform fault and Larnaka thrust fault zone that were already known and the Lemessos thrust fault zone that defined in the present study. By combining the past history of strong (M ≥ 6.0) events and the long-term tectonic loading on these major fault zones, the evolution of the stress field from 1896 until the present is derived. Although uncertainties exist in the location, magnitude and fault geometries of the early earthquakes included in our stress evolutionary model, the resulting stress field provides an explanation of later earthquake triggering. It was evidenced that the locations of all the strong events were preceded by a static stress change that encouraged failure. The current state of the evolved stress field may provide evidence for the future seismic hazard. Areas of positive static stress changes were identified in the southwestern offshore area that can be considered as possible sites of future seismic activity.