Magnitude of σ1, σ2, and σ3 at mid-crustal levels in an orogenic belt: Microboudin method applied to an impure metachert from Turkey

We revised an equation for estimating palaeostress magnitude using the microboudin technique by incorporating the influence of time on the fracture strength of minerals. The equation was used to estimate triaxial palaeostresses from a rare sample of metachert from Turkey that contains microboudinaged, columnar tourmaline grains in a wide range of orientations within the foliation plane. The estimated principal palaeostresses are σ₁ = 605 MPa, σ₂ = 598 MPa, and σ₃ = 597 MPa. As the microboudinage is considered to have occurred immediately before the rock encountered the brittle-plastic transition during exhumation, these stress values correspond to conditions at approximately 18 km depth and 300 °C within a Cretaceous orogenic belt.     

Talc friction in the temperature range 25°–400 °C: Relevance for Fault-Zone Weakening

Talc is one of the weakest minerals that is associated with fault zones. Triaxial friction experiments conducted on water-saturated talc gouge at room temperature yield values of the coefficient of friction, μ (shear stress, τ/effective normal stress, σ′_N) in the range 0.16–0.23, and μ increases with increasing σ′_N. Talc gouge heated to temperatures of 100°–400 °C is consistently weaker than at room temperature, and μ < 0.1 at slow strain rates in some heated experiments. Talc also is characterized by inherently stable, velocity-strengthening behavior (strength increases with increasing shear rate) at all conditions tested. The low strength of talc is a consequence of its layered crystal structure and, in particular, its very weak interlayer bond. Its hydrophobic character may be responsible for the relatively small increase in μ with increasing σ′_N at room temperature compared to other sheet silicates.Talc has a temperature–pressure range of stability that extends from surficial to eclogite-facies conditions, making it of potential significance in a variety of faulting environments. Talc has been identified in exhumed subduction zone thrusts, in fault gouge collected from oceanic transform and detachment faults associated with rift systems, and recently in serpentinite from the central creeping section of the San Andreas fault. Typically, talc crystallized in the active fault zones as a result of the reaction of ultramafic rocks with silica-saturated hydrothermal fluids. This mode of formation of talc is a prime example of a fault-zone weakening process. Because of its velocity-strengthening behavior, talc may play a role in stabilizing slip at depth in subduction zones and in the creeping faults of central and northern California that are associated with ophiolitic rocks.     

Geometry of the fault zone of the 2003 Ms = 7.5 Chuya earthquake and associated stress fields, Gorny Altai

The co-seismic deformations produced during the September 27, 2003 Chuya earthquake (Ms = 7.5) that affected the Gorny Altai, Russia, are described and discussed along a 30 km long segment. The co-seismic deformations have manifested themselves both in unconsolidated sediments as R- and R′-shears, extension fractures and contraction structures, and in bedrock as the reactivation of preexisting schistosity zones and individual fractures, as well as development of new ruptures and coarse crushing zones. It has been established that the pattern of earthquake ruptures represents a typical fault zone trending NW–SE with a width reaching 4–5 km and a dextral strike–slip kinematics. The initial stress field that produced the whole structural pattern of co-seismic deformations during the Chuya earthquake, is associated with a transcurrent regime with a NNW–SSE, almost N–S, trending of compressional stress axis (σ₁), and a ENE–WSW, almost E–W, trending of tensional stress axis (σ₃). The state of stress in the newly-formed fault zone is relatively uniform. The local stress variations are expressed in insignificant deviation of σ₁ from N–S to NW–SE or NE–SW, in short-term fluctuations of relative stress values in keeping their spatial orientations, or in a local increase of the plunge angle of the σ₁. The geometry of the fault zone associated with the Chuya earthquake has been compared with the mechanical model of fracturing in large continental fault zones with dextral strike–slip kinematics. It is apparent that the observed fracture pattern corresponds to the late disjunctive stage of faulting when the master fault is not fully developed but its segments are already clearly defined. It has been shown that fracturing in widely different rocks follows the common laws of the deformation of solid bodies, even close to the Earth surface, and with high rates of movements.     

A note on fault reactivation

Reactivation of existing faults whose normal lies in the σ₁^′σ₃^′ plane of a stress field with effective principal compressive stresses σ₁^′ >σ₂^′ >σ₃^′ is considered for the simplest frictional failure criterion, τ = μσ_n^′ = μ(σ_n − P), where τ and σ_n are respectively the shear and normal stresses to the existing fault, P is the fluid pressure and μ is the static friction. For a plane oriented at θ to σ₁^′, the stress ratio for reactivation is (σ₁^′/σ₃^′) = (1 + μ cot θ)/(1 − μ tan θ). This ratio has a minimum positive value at the optimum angle for reactivation given by (1/μ) but reaches infinity when θ = 2θ^*, beyond which σ₃^′ < 0 is a necessary condition for reactivation. An important consequence is that for typical rock friction coefficients, it is unlikely that normal faults will be reactivated as high-angle reverse faults or thrusts as low-angle normal faults, unless the effective least principal stress is tensile.     

A statistical analysis applied to the Rf/φ method

A statistical analysis for the R_f/φ method is developed in the paper. The technique provides estimates of R_s, R_i and θ according to R_f and φ. Moreover, we can restore the pattern of initial orientations of marker objects by means of their frequency histogram, provided the strain is homogeneous. It is unnecessary to postulate random initial orientations if the principal directions of strain ellipsoid are known.     

A new method of estimating the ratio between in situ rock stresses and tectonics based on empirical and probabilistic analyses

This paper describes a new procedure for assessing the ratio between in situ stresses in rock masses by means of K (K = σ_H / σ_v, being σ_H and σ_v principal stress) and tectonics for purposes of engineering geology and rock mechanics. The method combines the use of the logic decision tree and the empirical relationship between the Tectonic Stress Index, TSI, and a series of K in situ values obtained from an extensive database. The decision tree considers geological and geophysical factors affecting stress magnitudes both on the regional and local scale. The TSI index is defined by geological and geomechanical parameters. The method proposed provides an assessment of the magnitude of horizontal stresses of tectonic origin. Results for several regions of Europe are presented and the possible applications of the procedure are discussed.     

Neotectonic evolution of the Central Betic Cordilleras (Southern Spain)

Paleostress orientations were calculated from fault-slip data of 36 sites located along a traverse through the Central Betic Cordilleras (southern Spain). Heterogeneous fault sets, which are frequent in the area, have been divided into homogeneous subsets by cross-cutting relationships observed in the field and by a paleostress stratigraphy approach applied on each individual fault population. The state of stress was sorted according to main tectonic events and a new chronology is presented of the Miocene to Recent deformation in the central part of the Betic Cordilleras. The deviatoric stress tensors fall into four distinct groups that are regionally consistent and correlate with three Late Oligocene–Aquitanian to Recent major tectonic events in the Betic Cordilleras. The new chronology of the neotectonic evolution includes, from oldest to youngest, the following main tectonic phases:

(1) Late Oligocene–Aquitanian to Early Tortonian: σ₁ subhorizontal N–S, partly E–W directed, σ₃ subvertical; compressional structures (thrusting of nappes, large-scale folding) and strike-slip faulting in the Alborán Domain and the External Zone of the Betic Cordilleras;

(2) Early Tortonian to Pliocene–Pleistocene: σ₁ subvertical, σ₃ subhorizontal NW–SE, partly N–S directed or E–W-directed (radial extension); large-scale normal faulting in the Central Betic Cordilleras and in the oldest Neogene formations of the Granada Basin related to the gravitational collapse of the Betic Cordilleras and the exhumation of the intensely metamorphosed rock series of the Internal Zones, at the same time formation of the Alborán Basin and intramontane basins such as the Granada Basin;

(3) Pleistocene to Recent: (3a) σ₁ subvertical, σ₃ subhorizontal NE–SW with prominent normal faulting, but coevally; (3b) σ₁ subhorizontal NW directed, σ₃ NE–SW subhorizontal with strike-slip faulting. Extensional structures and strike-slip faulting are related to the ongoing convergence of the Eurasian and African Plates and coeval uplift of the Betic Cordilleras. Reactivation of pre-existing fractures and faults was frequently observed. Phase 3 is interpreted as periodic strike-slip and normal faulting events due to a permutation of the principal stress axes, mainly σ₁ and σ₂.

Determination of Shear Strength and Three-dimensional Yield Strength for the Hoek-Brown Criterion

Summary. Most currently used techniques for analysing the stability of near surface structures, such as rock slopes, are based on the application of the effective Coulomb shear strength parameters cohesion c′, and the angle of friction φ′ on some known or anticipated shear surface subjected to an effective normal stress σ′_n. The most widely used of these techniques are the variants of the method of slices and related upper bound techniques. If the Hoek-Brown criterion is to be used to model the strength of near surface fractured rocks, it is necessary to determine equivalent Coulomb shear strength parameters for the specified level of effective normal stress. Calculation of the equivalent Coulomb parameters for the Hoek-Brown criterion for cases when a ≠ 0.5 is not a straightforward matter. A simple procedure for calculating instantaneous values of c′_i and φ′_i has been developed based on spreadsheet calculations and the application of a numerical optimisation routine. This procedure can also be applied to calculating the Hoek-Brown envelope plotted in shear stress/normal stress space. A simple closed form solution for c′_i and tan φ′_i has also been developed for the special case when a = 1. A three-dimensional version of the Hoek-Brown criterion has been developed by combining it with the Drucker-Prager criterion. This new yield criterion has been implemented by numerical solution of the governing equations. A simplification of this three-dimensional yield criterion has been developed by introducing an intermediate principal stress weighting factor. Comparison with published results demonstrates that this simplified criterion has the capacity to model the results of true triaxial tests for a range of different rock types over a wide range of stress levels. The new three-dimensional yield criterion has the advantage that its input parameters can be determined from routine uniaxial compression tests and mineralogical examination.     

CO2 and sulfuric acid controls of weathering and river water composition

Reactions of CO₂ with carbonate and silicate minerals in continental sediments and upper part of the crystalline crust produce HCO₃⁻ in river and ground waters. H₂SO₄ formed by the oxidation of pyrite and reacting with carbonates may produce CO₂ or HCO₃⁻. The ratio, ψ, of atmospheric or soil CO₂ consumed in weathering to HCO₃⁻ produced depends on the mix of CO₂ and H₂SO₄, and the proportions of the carbonates and silicates in the source rock. An average sediment has a CO₂ uptake potential of ψ = 0.61. The potential increases by inclusion of the crystalline crust in the weathering source rock. A mineral dissolution model for an average river gives ψ = 0.68 to 0.72 that is within the range of ψ = 0.63 to 0.75, reported by other investigators using other methods. These results translate into the CO₂ weathering flux of 20 to 24 × 10¹²mol/yr.     

Fault slip analysis of Quaternary faults in southeastern Korea

The Quaternary stress field has been reconstructed for southeast Korea using sets of fault data. The subhorizontal direction of the maximum principal stress (σ₁) trended ENE and the direction of the minimum principal stress (σ₃) was nearly vertical. The stress ratio (Φ = (σ₂ − σ₃) / (σ₁ − σ₃)) value was 0.65. Two possible interpretations for the stress field can be made in the framework of eastern Asian tectonics; (1) The σ_Hmax trajectory for southeast Korea fits well with the fan-shaped radial pattern of maximum principal stress induced by the India–Eurasia collision. Thus, we suggest that the main source for this recent stress field in southeast Korea is related to the remote India–Eurasia continental collision. (2) The stress field in Korea shows a pattern similar to that in southwestern Japan. The origin for the E–W trending σ_Hmax in Japan is known to be related to the mantle upwelling in the East China Sea. Thus, it is possible that Quaternary stress field in Korea has evolved synchronously with that in Japan. We suggest further studies (GPS and in situ stress measurement) to test these hypotheses.     

Stress fields and fault reactivation angles of the 2000 western Tottori aftershocks and the 2001 northern Hyogo swarm in southwest Japan

We estimated the stress fields of the aftershocks of the 2000 western Tottori earthquake (M_w 6.6) and the northern Hyogo swarm (max M_w 5.2) by a stress tensor inversion of moment tensor solutions reported from the National Research Institute for Earth Science and Disaster Prevention (Japan). The maximum principal stress direction of the western Tottori sequence was estimated as N107°E with a strike–slip regime. In the northern Hyogo swarm, the orientations of the principal stress directions could not be well constrained by the observed data, but after examining the detailed characteristics of the solution, we obtained a most probable solution of N113°E for the σ₁ direction. These solutions are consistent with the maximum horizontal directions roughly estimated from the strike directions of large earthquakes occurring geographically between these two seismic activities. We measured the angle between each fault–slip direction and maximum principal stress direction to investigate the frictional properties of earthquakes. The distribution of the angles was forward modeled to estimate the coefficient of friction and the stress ratio, assuming uniformly distributed fault orientations. For the western Tottori sequence, a homogeneous stress field with a coefficient of friction less than 0.4 was estimated. A high stress level was also suggested because very little change occurred in the stress field during the mainshock. For the northern Hyogo sequence, the coefficient of friction was estimated to be between 0.5 and 1.0.     

Magnetic susceptibility anisotropy in the Cambrian slate belt of North Wales and correlation with strain

The magnetic susceptibility anisotropy of 275 specimens comprising 38 sites from the Cambrian slate belt in North Wales was measured to determine the magnetic fabric of the slates. The susceptibility ellipsoid is oblate for all sites, and the maximum/intermediate susceptibility plane always coincides with the cleavage plane of the slates which has a Caledonian strike and is nearly vertical. The maximum axes align sub-vertically and the intermediate axes sub-horizontally, trending NE-SW. The minimum susceptibility axes are normal to this foliation plane and coincide with the poles to the slaty cleavage. The orientations of the principal susceptibility axes are found to be in excellent agreement with the orientations of the principal strain directions, determined by X-ray goniometry on one of the samples from almost all of the sites. Correlation of the magnetic susceptibility anisotropy with predicted March strains (March, 1932) shows that the principal magnitudes of susceptibility can be related to those of the strain by: (for i = 1, 2, 3. The orthogonal principal axes), where χ_f and χ₀ are the final and initial susceptibilities along a given axis i and l_f and l_i are final and initial axial dimensions in the same direction i of a principal strain axis. The exponent a for the North Wales slates was found to be 0.145 ± 0.005. Knowledge of such a relationship may permit rapid approximate determinations of a petrofabric in similar rocks from their magnetic fabrics. However, the exponent a will probably have to be recalibrated for each rock type.     

Determination of stress directions from plagioclase fabrics in high grade deformed rocks (Além Paraíba shear zone, Ribeira fold belt, southeastern Brazil)

A new method to determine stress directions using the preferential orientation of plagioclase mechanical twins has been applied to high-temperature mylonitic rocks from the Além Paraíba shear zone, Ribeira fold belt, southeastern Brazil. We have measured the lattice-preferred orientation of plagioclase grains and calculated the orientation of the stress axes possible for the observed twin orientations. The maximum compressive stress direction (σ₁), determined for all studied samples, is a function of the mechanical twin orientations of a number of distinct plagioclase populations. The σ₁ direction is generally subperpendicular to the (010) plane. The statistical treatment for most of the plagioclase grains examined for each sample shows that σ₁ is almost perpendicular to the foliation plane, suggesting a significant coaxial component in the deformation process of these rocks.     

Constraining stress magnitudes using petroleum exploration data in the Cooper–Eromanga Basins, Australia

The magnitude of the in situ stresses in the Cooper–Eromanga Basins have been determined using an extensive petroleum exploration database from over 40 years of drilling. The magnitude of the vertical stress (S_v) was calculated based on density and velocity checkshot data in 24 wells. Upper and lower bound values of the vertical stress magnitude are approximated by S_v = (14.39 × Z)^1.12 and S_v = (11.67 × Z)^1.15 functions respectively (where Z is depth in km and S_v is in MPa). Leak-off test data from the two basins constrain the lower bound estimate for the minimum horizontal stress (S_hmin) magnitude to 15.5 MPa/km. Closure pressures from a large number of minifrac tests indicate considerable scatter in the minimum horizontal stress magnitude, with values approaching the magnitude of the vertical stress in some areas. The magnitude of the maximum horizontal stress (S_Hmax) was constrained by the frictional limits to stress beyond which faulting occurs and by the presence of drilling-induced tensile fractures in some wells. The maximum horizontal stress magnitude can only be loosely constrained regionally using frictional limits, due to the variability of both the minimum horizontal stress and vertical stress estimates. However, the maximum horizontal stress and thus the full stress tensor can be better constrained at individual well locations, as demonstrated in Bulyeroo-1 and Dullingari North-8, where the necessary data (i.e. image logs, minifrac tests and density logs) are available. The stress magnitudes determined indicate a predominantly strike-slip fault stress regime (S_Hmax > S_v > S_hmin) at a depth of between 1 and 3 km in the Cooper–Eromanga Basins. However, some areas of the basin are transitional between strike-slip and reverse fault stress regimes (S_Hmax > S_v ≈ S_hmin). Large differential stresses in the Cooper–Eromanga Basins indicate a high upper crustal strength for the region, consistent with other intraplate regions. We propose that the in situ stress field in the Cooper–Eromanga Basins is a direct result of the complex interaction of tectonic stresses from the convergent plate boundaries surrounding the Indo-Australian plate that are transmitted into the center of the plate through a high-strength upper crust.     

Triaxial stress state deep in orogenic belts: an example from Turkey

A rare metachert pebble containing amphibole grains with microboudin structures in a wide range of orientations provides an opportunity to perform stress analysis in two orthogonal orientations on the foliation surface. The sample was analysed by the microboudin method to infer the triaxial stress state during microboudinage. Stress parameters proportional to the far-field differential stress were determined for sodic amphibole grains in the two orientations. The ratio of the stresses in the two orthogonal orientations (σ₁−σ₂)/(σ₁−σ₃) was calculated to be 0.64, indicating that σ₂ lies closer to the midpoint between σ₁ and σ₃ than to σ₃.     

3-D velocity structure of the 2003 Bam earthquake area (SE Iran): Existence of a low-Poisson's ratio layer and its relation to heavy damage

To understand the generation mechanism of the Bam earthquake (Mw 6.6), we studied three-dimensional V_P, V_S and Poisson's ratio (σ) structures in the Bam area by using the seismic tomography method. We inverted accurate arrival times of 19490 P waves and 19015 S waves from 2396 aftershocks recorded by a temporal high-sensitivity seismic network. The 3-D velocity structure of the seismogenic region was well resolved to a depth of 14 km with significant velocity variations of up to 5%. The general pattern of aftershock distribution was relocated by using the 3-D structure to delineate a source fault for a length of approximately 20 km along a line 4.5 km west of the known geological Bam fault; this source fault dips steeply westward and strikes a nearly north–south line. The main shallow cluster of aftershocks south of the city of Bam is distributed just under the minor surface ruptures in the desert. The 3-D velocity structure shows a thick layer of high V_S and low σ (minimum: 0.20) at a depth range of 2–6 km. The deeper layer, with a thickness of about 2 km, appears to have a low V_S and high σ (maximum: 0.28) from 6 km depth beneath Bam to a depth of 9 km south of the city. The inferred increase of Poisson's ratio from 2 to 10 km in depth may be associated with a change from rigid and SiO₂-rich rock to more mafic rock, including the probable existence of fluids. The main seismic gap of aftershock distribution at the depth range of 2 to 7 km coincides well with the large slip zone in the shallow thick layer of high V_S and low σ. The large slip propagating mainly in the shallow rigid layer may be one of the main reasons why the Bam area suffered heavy damage.     

Cenozoic post-collisional brittle tectonic history and stress reorientation in the High Zagros Belt (Iran, Fars Province)

The Zagros fold-and-thrust belt of SW-Iran is among the youngest continental collision zones on Earth. Collision is thought to have occurred in the late Oligocene–early Miocene, followed by continental shortening. The High Zagros Belt (HZB) presents a Neogene imbricate structure that has affected the thick sedimentary cover of the former Arabian continental passive margin. The HZB of interior Fars marks the innermost part of SE-Zagros, trending NW–SE, that is characterised by higher elevation, lack of seismicity, and no evident active crustal shortening with respect to the outer (SW) parts. This study examines the brittle structures that developed during the mountain building process to decipher the history of polyphase deformation and variations in compressive tectonic fields since the onset of collision. Analytic inversion techniques enabled us to determine and separate different brittle tectonic regimes in terms of stress tensors. Various strike–slip, compressional, and tensional stress regimes are thus identified with different stress fields. Brittle tectonic analyses were carried out to reconstruct possible geometrical relationships between different structures and to establish relative chronologies of corresponding stress fields, considering the folding process. Results indicate that in the studied area, the main fold and thrust structure developed in a general compressional stress regime with an average N032° direction of σ₁ stress axis during the Miocene. Strike–slip structures were generated under three successive strike–slip stress regimes with different σ₁ directions in the early Miocene (N053°), late Miocene–early Pliocene (N026°), and post-Pliocene (N002°), evolving from pre-fold to post-fold faulting. Tensional structures also developed as a function of the evolving stress regimes. Our reconstruction of stress fields suggests an anticlockwise reorientation of the horizontal σ₁ axis since the onset of collision and a significant change in vertical stress from σ₃ to σ₂ since the late stage of folding and thrusting. A late right-lateral reactivation was also observed on some pre-existing belt-parallel brittle structures, especially along the reverse fault systems, consistent with the recent N–S plate convergence. However, this feature was not reflected by large structures in the HZB of interior Fars. The results should not be extrapolated to the entire Zagros belt, where the deformation front has propagated from inner to outer zones during the younger events.     

Stress field in the Northern Rhine area, Central Europe, from earthquake fault plane solutions

Fault plane solutions (FPS) from 110 earthquakes in the northern Rhine area with local magnitudes, ranging from 1.0 to 6.1, and occurring between 1976 and 2002 are determined. FPS are retrieved from P-wave first motions using a grid search approach allowing a detailed exploration of the parameter space. The influence of the 1D velocity model on take-off angles and resulting FPS is examined. All events were relocated with a recently developed minimum 1D model of the velocity structure [J. Geophys. Res. (2003)]. Rose diagrams of the orientation of P, T and B axes show a clear preference of trends of P and T axes at N292°E and N27°E, respectively. The majority of B axes trend in northerly directions. Plunges of P and T axes are mostly around 45° while most B axes are subhorizontal. The main direction of the maximum horizontal stress directly inferred from the fault plane solutions is N118°E.To calculate the orientations of the principal stress axes and the shape of the stress tensor, the inversion method of Gephard and Forsyth [J. Geophys. Res. 89 (1984) 9305] was applied to the whole data set and to several subsets of data. The subsets were formed by grouping events from various geological and tectonic areas and by grouping events into different depth ranges. The subset areas include the Lower Rhine Embayment, the Rhenish Massif, the middle Rhine area, the Neuwied Basin and the area known as the Stavelot–Venn Massif. Inversion of the entire data set shows some ambiguity between a strike-slip and extensional stress regime, with a vertical axis for the medium principal stress and a trend of N305°E and N35°E for the σ₁ and σ₃ axis, respectively, as the best fitting tensor. Earthquakes from the Lower Rhine Embayment and, to some degree, from the middle Rhine area indicate an extensional stress regime. In the Lower Rhine Embayment, plunge and trend of the σ₁ axis are 76° and N162°E and for the σ₃ axis 7° and N42°E. The best fitting solution for the area of the Stavelot–Venn Massif is a strike-slip regime with subhorizontal σ₁ and σ₃ axes with a trend of N316°E and N225°E, respectively. Stress orientations found here agree overall with the results from earlier studies based on smaller data sets. The directions of the maximum and minimum horizontal stresses inverted from focal mechanisms agree well with the stress field predicted by the European Stress Map. This confirms earlier interpretations that the stress field of the Rhine Graben system is controlled by plate driving forces acting on the plate boundaries. However, amplitudes of the stresses change on a local scale and with depth. Estimates of the absolute magnitude of principal stresses favor a normal faulting regime in the shallow crust (above 12-km depth) and a strike-slip regime in the lower crust.     

Tectonic analysis of an oceanic transform fault zone based on fault-slip data and earthquake focal mechanisms: the Húsavík–Flatey Fault zone, Iceland

The Húsavík–Flatey Fault (HFF) is an oblique dextral transform fault, part of the Tjörnes Fracture Zone (TFZ), that connects the North Volcanic Zone of Iceland and the Kolbeinsey Ridge. We carry out stress inversion to reconstruct the paleostress fields and present-day stress fields along the Húsavík–Flatey Fault, analysing 2700 brittle tectonic data measured on the field and about 700 earthquake focal mechanisms calculated by the Icelandic Meteorological Office. This allows us to discuss the Latest Cenozoic finite deformations (from the tectonic data) as well as the present-day deformations (from the earthquake mechanisms). In both these cases, different tectonic groups are reconstructed and each of them includes several distinct stress states characterised by normal or strike-slip faulting. The stress states of a same tectonic group are related through stress permutations (σ₁−σ₂ and σ₂−σ₃ permutations as well as σ₁−σ₃ reversals). They do not reflect separate tectonic episodes. The tectonic groups derived from the geological data and the earthquake data have striking similarity and are considered to be related. The obliquity of the Húsavík–Flatey Fault implies geometric accommodation in the transform zone, resulting mainly from a dextral transtension along an ENE–WSW trend. This overall mechanism is subject to slip partitioning into two stress states: a Húsavík–Flatey Fault-perpendicular, NE–SW trending extension and a Húsavík–Flatey Fault-parallel, NW–SE trending extension. These three regimes occur in various local tectonic successions and not as a regional definite succession of tectonic events. The largest magnitude earthquakes reveal a regional stress field tightly related to the transform motion, whereas the lowest magnitude earthquakes depend on the local stress fields. The field data also reveal an early extension trending similar to the spreading vector. The focal mechanism data do not reflect this extension, which occurred earlier in the evolution of the HFF and is interpreted as a stage of structural development dominated by the rifting process.     

Seismotectonics of strike–slip earthquakes within the deep crust of southern Italy: Geometry, kinematics, stress field and crustal rheology of the Potenza 1990–1991 seismic sequences (Mmax 5.7)

We present a revision and a seismotectonic interpretation of deep crust strike–slip earthquake sequences that occurred in 1990–1991 in the Southern Apennines (Potenza area). The revision is motivated by: i) the striking similarity to a seismic sequence that occurred in 2002  140 km NNW, in an analogous tectonic context (Molise area), suggesting a common seismotectonic environment of regional importance; ii) the close proximity of such deep strike–slip seismicity with shallow extensional seismicity (Apennine area); and iii) the lack of knowledge about the mechanical properties of the crust that might justify the observed crustal seismicity. A comparison between the revised 1990–1991 earthquakes and the 2002 earthquakes, as well as the integration of seismological data with a rheological analysis offer new constraints on the regional seismotectonic context of crustal seismicity in the Southern Apennines. The seismological revision consists of a relocation of the aftershock sequences based on newly constrained velocity models. New focal mechanisms of the aftershocks are computed and the active state of stress is constrained via the use of a stress inversion technique. The relationships among the observed seismicity, the crustal structure of the Southern Apennines, and the rheological layering are analysed along a crustal section crossing southern Italy, by computing geotherms and two-mechanism (brittle frictional vs. ductile plastic strength) rheological profiles. The 1990–1991 seismicity is concentrated in a well-defined depth range (mostly between 15 and 23 km depths). This depth range corresponds to the upper pat of the middle crust underlying the Apulian sedimentary cover, in the footwall of the easternmost Apennine thrust system. The 3D distribution of the aftershocks, the fault kinematics, and the stress inversion indicate the activation of a right-lateral strike–slip fault striking N100°E under a stress field characterized by a sub-horizontal N142°-trending σ1 and a sub-horizontal N232°-trending σ3, very similar to the known stress field of the Gargano seismic zone in the Apulian foreland. The apparent anomalous depths of the earthquakes (> 15 km) and the confinement within a relatively narrow depth range are explained by the crustal rheology, which consists of a strong brittle layer at mid crustal depths sandwiched between two plastic horizons. This articulated rheological stratification is typical of the central part of the Southern Apennine crust, where the Apulian crust is overthrusted by Apennine units. Both the Potenza 1990–1991 and the Molise 2002 seismic sequences can be interpreted to be due to crustal E–W fault zones within the Apulian crust inherited from previous tectonic phases and overthrusted by Apennine units during the Late Pliocene–Middle Pleistocene. The present strike–slip tectonic regime reactivated these fault zones and caused them to move with an uneven mechanical behaviour; brittle seismogenic faulting is confined to the strong brittle part of the middle crust. This strong brittle layer might also act as a stress guide able to laterally transmit the deviatoric stresses responsible for the strike–slip regime in the Apulian crust and may explain the close proximity (nearly overlapping) of the strike–slip and normal faulting regimes in the Southern Apennines. From a methodological point of view, it seems that rather simple two-mechanism rheological profiles, though affected by uncertainties, are still a useful tool for estimating the rheological properties and likely seismogenic behaviour of the crust.