A Probabilistic Approach for Identifying Independent Remote Compressions in an Intraplate Region: The?Iberian Chain (Spain)

The discrimination between distinct remote compressions and multiple local stress deviations within a single compressive stress field has been carried out in the central-eastern Iberian Chain, by using structural criteria, computing palaeostress directions, identifying and ‘filtering’ stress deviations, and analysing time relationships. A probabilistic analysis based on a systematic comparison of real and expected frequencies of coexistence of two compressions is applied by means of the χ ² Test. This allows us to identify those tectonic compressions that behave as independent events from the probabilistic point of view. The results suggest that among five initially defined compression directions only three can be considered as representative of distinct (although partially superposed) externally applied intraplate stress fields: Iberian s.l. (NE-SW), Betic s.l. (NW-SE), and Pyrenean (N-S to NNE-SSW).     

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.     

European continuous active tectonic strain–stress map

This paper shows a new continuous strain–stress map for Europe obtained from the direct inversion of earthquake focal mechanisms calculated from the centroid tensor method. A total of 1608 focal mechanisms have been selected with several quality criteria from different catalogues (CMT Harvard, ETH, Med-Net, I.G.N. and I.A.G.) from 1973 to the present day. Values for the maximum horizontal shortening direction and brittle strain–stress regime defined by the k′ ratio (e_y/e_z, horizontal maximum/vertical strain) have been calculated following in Europe and Pannonian Basin the slip model of tri-axial deformation. The individual results including Dey and the shape of the active brittle strain ellipsoid have been interpolated to a final 15′ regular grid taking into account the relationship between the tectonic horizontal strain–stress value and the vertical load. Both continuous strain regime and maximum horizontal shortening (Dey) maps show a good correlation with the primary tectonic forces generated along the plate boundaries, plate kinematics and also some local perturbations related with main crustal heterogeneities and topography, as well as significant spatial variations in integrated crustal strength.     

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.     

A Probabilistic Approach for Identifying Independent Remote Compressions in an Intraplate Region: The Iberian Chain (Spain)

The discrimination between distinct remote compressions and multiple local stress deviations within a single compressive stress field has been carried out in the central-eastern Iberian Chain, by using structural criteria, computing palaeostress directions, identifying and ‘filtering’ stress deviations, and analysing time relationships. A probabilistic analysis based on a systematic comparison of real and expected frequencies of coexistence of two compressions is applied by means of the χ ² Test. This allows us to identify those tectonic compressions that behave as independent events from the probabilistic point of view. The results suggest that among five initially defined compression directions only three can be considered as representative of distinct (although partially superposed) externally applied intraplate stress fields: Iberian s.l. (NE-SW), Betic s.l. (NW-SE), and Pyrenean (N-S to NNE-SSW).     

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.     

Stress field analysis at the El Teniente Mine: evidence for N–S compression in the modern Andes

There is a large database of triaxial stress measurements at the El Teniente Mine, Central Chile, but the complex geology, severe topography, and proximity of all measurements to extensive mining excavations made interpretation of the stress field difficult. The measurements were analyzed using three-dimensional numerical stress analysis and decomposition of the stress field into gravitational and tectonic components. By removing gravitational stresses plus local effects from the tectonic component of the stress field a calculation of the far-field tectonic stress tensor is made. It is shown that variations in the tectonic component of stress are related to shear zones cutting through the mine. The far-field major principal component of the tectonic stress field was found to be oriented approximately N–S. This is consistent with the most recent direction of local shortening based on kinematic analysis of faults, but is perpendicular to the direction of regional crustal shortening. There appears to be a limiting envelope to the magnitude of the stress field implying that the shear zones are in a state of limiting equilibrium with regional tectonic driving forces.     

Evolution of paleostress fields and brittle deformation of the Tornquist Zone in Scania (Sweden) during Permo-Mesozoic and Cenozoic times

The NW-SE oriented Sorgenfrei–Tornquist Zone (STZ) has been thoroughly studied during the last 25 years, especially by means of well data and seismic profiles. We present the results of a first brittle tectonic analysis based on about 850 dykes, veins and minor fault-slip data measured in the field in Scania, including paleostress reconstruction. We discuss the relationships between normal and strike-slip faulting in Scania since the Permian extension to the Late Cretaceous–Tertiary structural inversions. Our paleostress determinations reveal six successive or coeval main stress states in the evolution of Scania since the Permian. Two stress states correspond to normal faulting with NE-SW and NW-SE extensions, one stress state is mainly of reverse type with NE-SW compression, and three stress states are strike-slip in type with NNW-SSE, WNW-ESE and NNE-SSW directions of compression.The NE-SW extension partly corresponds to the Late Carboniferous–Permian important extensional period, dated by dykes and fault mineralisations. However extension existed along a similar direction during the Mesozoic. It has been locally observed until within the Danian. A perpendicular NW-SE extension reveals the occurrence of stress permutations. The NNW-SSE strike-slip episode is also expected to belong to the Late Carboniferous–Permian episode and is interpreted in terms of right-lateral wrench faulting along STZ-oriented faults. The inversion process has been characterised by reverse and strike-slip faulting related to the NE-SW compressional stress state.This study highlights the importance of extensional tectonics in northwest Europe since the end of the Palaeozoic until the end of the Cretaceous. The importance and role of wrench faulting in the tectonic evolution of the Sorgenfrei–Tornquist Zone are discussed.     

Immature and mature transform zones near a hot spot: The South Iceland Seismic Zone and the Tjörnes Fracture Zone (Iceland)

We describe and compare the two transform zones that connect the Icelandic rift segments and the mid-Atlantic Ridge close to the Icelandic hot spot, in terms of geometry of faulting and stress fields. The E–W trending South Iceland Seismic Zone is a diffuse shear zone with a Riedel fault pattern including N0°–N20°E trending right-lateral and N60°–N70°E trending left-lateral faults. The dominant stress field in this zone is characterised by NW–SE extension, in general agreement with left-lateral transform motion. The Tjörnes Fracture Zone includes three major lineaments at different stages of development. The most developed, the Húsavík–Flatey Fault, presents a relatively simple geometry with a major fault that trends ESE–WNW. The stress pattern is however complex, with two dominant directions of extension, E–W and NE–SW on average. Both these extensions are compatible with the right-lateral transform motion and reveal different behaviours in terms of coupling. Transform motion has unambiguous fault expression along a mature zone, a situation close to that of the Tjörnes Fracture Zone. In contrast, transform motion along the immature South Iceland Seismic Zone is expressed through a more complicate structural pattern. At the early stage of the transform process, relatively simple stress patterns prevail, with a single dominant stress field, whereas, when the transform zone is mature, moderate and low coupling situations may alternate, as a function of volcanic–tectonic crises and induce changes in stress orientation.     

From continental margin extension to collision orogen: structural development and tectonic rotation of the Hengchun peninsula, southern Taiwan

As a result of oblique collision, the Taiwan orogen propagates southward. The Hengchun peninsula in the southern tip of the Taiwan Central Range, preserving the youngest, the least deformed and the most complete accretionary prism sequences, allows therefore better understanding of the tectonic evolution of Taiwan orogen. On the Hengchun peninsula, four main stages of paleostress can be recognized by the analysis of brittle tectonics. After recording the first two stages of paleostress, rocks of the Hengchun peninsula (the Hengchun block) have undergone both tilting and counterclockwise rotation of about 90°. The structural boundaries of this rotated Hengchun block are: the Kenting Mélange zone in the southwest, the Fongkang Fault in the north, and a submarine backthrust in the east. The angle of this rotation is principally calculated by the paleomagnetic analysis data and a physical model experiment. Through a systematic back-tilting and back-rotating restoration, the original orientations of the four paleostress stages of Hengchun peninsula are recognized. They are, from the ancient to the recent, a NW–SE extension, a combination of NW–SE transtension and NE–SW transpression, a NE–SW compression, and finally a combination of NE–SW transtension and NW–SE transpression. This result can be explained by a phenomenon of stress axes permutation, instead of a complex polyphase tectonism. This stress axes permutation is caused by the horizontal compression increase accompanying the propagation of the accretionary prism. Combining the tectonic and paleomagnetic data with paleocurrent and stratigraphic data enables us to reconstruct the tectonic evolution of the Hengchun peninsula. This reconstruction corresponds to the deformation history of a continental margin basin, from its opening to its intense deformation in the accretionary prism.     

Paleo-deviatoric stress magnitudes from calcite twins and related structural permeability evolution in minor faults: Example from the toarcian shale of the French Causses Basin, Aveyron, France

The relationship between deformation and so-called fluid paleotransfers in minor faults has been analysed in an argillaceous formation located in the Causses Basin in France. The fluid paleotransfers are related to the fault activity to a large extent. We attempt to estimate the intensity of paleo-deviatoric stress magnitudes under which the fault activity may have occurred and consequently, the change in the structural fault permeability. The paleo-deviatoric stress magnitudes were calculated with the inverse method of Etchecopar applied to calcite twinning. The measured crystals are contained within the core zone of minor faults and this study is based on a previous complete microtectonic and microstructural analysis of the faults. In this paper, analysis of calcite twinning has been applied for the first time to vein fillings associated small faults in a context of relatively weak deformation, a condition ensured by the tectonic and structural analysis. Calculation and discussion of the paleo-deviatoric stress tensors in relation to the evolution of the structural fault permeability and to the hydraulic behaviour of the faults are the aim of this paper. The analysed faults, created and active during the same tectonic event, were permeable under a (σ₁–σ₃) mean value of 40–50 MPa. On the other hand, the reactivation of faults during a second tectonic event implies mean (σ₁–σ₃) value higher than 40–50 MPa, especially for the faults that are poorly oriented with respect to the principal tectonic stress directions. The core zone of these faults remained sealed and impermeable or became permeable by development of microcracks inside the pre-existing fillings.     

Variation with depth of the stress tensor anisotropy inferred from microfault analysis

In order to study the Pyrenean tectonic phase, a quantitative method of stress analysis using microfault measurements is used on a calcareous plateau located in southern France. The method developed here allows the determination of several tectonic events and the evaluation of (with σ σ₂ σ₃).The Pyrenean compression is seen to occur in two stages, confirming previous geological studies.A recent canyon allows the study of the variation with depth of the R ratio on a vertical cross-section (300 m). With a simple model of gravity and tectonic stresses, the vertical variation of R can be used to estimate quantitatively the Pyrenean tectonic stresses. For realistic values of the parameters in this model, the horizontal tectonic stresses are obtained in the following range: 50–200 bar for the maximum horizontal principal stress, 10–25 bar for the other horizontal stress.These results seem to be consistent with in situ stress measurements, but they are much lower than those predicted by experimental rock mechanics.     

Striated and pitted pebbles as paleostress markers: an example from the central transect of the Betic Cordillera (SE Spain)

Striated and pitted pebbles provide scarce structures that preserve information on the stresses that their host rocks have undergone. This information can be obtained by the measurement of a large number of microfaults with striae and solution marks within a small rock volume. For non-rotational deformation, the statistical procedures for microfault analysis provide a valid tool for determining the overprinting of successive stress ellipsoids, including their axial ratios and the orientations of the main axes. The trends of compressions obtained from striae can be compared with the determinations from the pole of pebble solution pits. However, in complex tectonics settings, the solution pits of several deformation phases are mixed and only striae analysis allows overprinted paleostresses to be accurately distinguished. The analysis of several pebbles from the same outcrop, including five from moderately complex settings, allows determination of the homogeneity of the paleostresses at outcrop scale, the detection of redeposited pebbles, and supports the results of microtectonic analysis for large areas. Solution mark distributions on pebbles depend on the burial and tectonic stresses. Conglomerates from shallow levels, such as those from Quaternary fluvial terraces, only record horizontal compressional solution marks because the minimum vertical stress needed to develop these structures are not reached by burial.In the central Betic Cordillera, striated and pitted pebbles are composed of carbonate surrounded by a matrix containing siliciclastic elements. The study of several outcrops located across a transect of the Cordillera shows a change in the recent stress field. While conglomerates near the Internal–External zone boundary show extensional stresses that may be related to the uplift of the Cordillera since Tortonian times, the outcrops located in the External Zone and up to the mountain front indicate the existence of horizontal NW–SE and NE–SW compressions related to prolate ellipsoids. These two compression directions, which affect conglomerates up to the Quaternary in the same outcrop, may be produced by a local permutation of stress axes, which in general indicates NW–SE compression related to the Eurasia–Africa plate boundary convergence, but which locally may switch to an orthogonal compression.     

Stress field in the Ryukyu–Kyushu Wadati–Benioff zone by inversion of earthquake focal mechanisms

The study addresses the space distribution of the stress field in the Kyushu–Ryukyu Wadati–Benioff zone based on homogeneous data of earthquake focal mechanisms and the inverse technique by Gephart and Forsyth [J. Geophys. Res. 89 (1984) 9305]. The used data set consists of 148 Harvard CMT solutions and 22 earthquake focal mechanisms listed in previous studies. The stress field parameters are determined for 0–40, 41–100 and h>100 km depth ranges. The top 100-km layer of the Wadati–Benioff zone (WBZ) is characterized by strike normal maximum compression σ₁ and steeper than the slab minimum compression σ₃, the last indicating for unbalanced slab pull force. The Tokara channel ‘divides’ the subduction into two parts of different stress regime at depth greater than 100 km. To the south of the channel the slab is under slab parallel σ₁ and slab normal σ₃ while its northern part, beneath Kyushu, is under slab parallel extension and slab normal compression. The results of recent studies on the regional velocity structure and geochemistry of the volcanic lava indicate that the most plausible reason for the observed stress field difference below 100 km in the northern and rest part of the arc is the presence of hot low viscosity upper mantle west of Kyushu.The results of this study indicate that the forces involved in the contemporary subduction dynamics in the Ryukyu–Kyushu Wadati–Benioff zone are related to the convergence between the Philippine Sea Plate and the Eurasian plate, the trench suction force, slab pull, the slab anchor force and, in the southern-central part of the arc, mantle resistance.     

Paleostress analysis of the Cretaceous rocks in the eastern margin of the Dead Sea transform, Jordan

This paper presents the first paleostress results from fault-slip data on Cretaceous limestone at the eastern rim of the Dead Sea transform (DST) in Jordan. Stress inversion of fault-slip data is performed using an improved right dieder method, followed by rotational optimization (Delvaux, TENSOR Program). The orientation of the principal stress axes (σ₁, σ₂ and σ₃) and the ratio of the principal stress differences ( ) show two main paleostress fields marking two main stress regimes, strike-slip and extensional. The first is characterized by NNW–SSE compression and ENE–WSW extension and related to Middle Miocene-Recent sinistral movement along the Dead Sea transform and the opening of the Red Sea. The second paleostress field is a WNW–ESE compression and NNE–SSW extension restricted to the northern part of the investigated area. This stress field could be associated with the development of the Syrian Arc fold belt which started during the Turonian, or it may be due to an anticlockwise rotation of the first stress field.     

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.     

Photoelastic and numerical investigation of stress distributions around fault models under biaxial compressive loading conditions

This paper presents the results of a photoelastic and numerical study of stress distributions (contours and trajectories) around fault models of various geometries, submitted to a biaxial compressive load. It aims to describe typical biaxial stress behaviours and emphasize the existing differences with the well-known uniaxial compressive load case. Stress trajectories are sometimes shown by joint sets acting as markers of a paleostress field, and they can be interpreted by particular shallow tension–compression situations. At depth fractures can be reactivated, or can dilate under conditions of triaxial compression, and behaviour is essentially controlled by a high stress ratio (high σ₃/σ₁ ratio). In spite of the potential importance of such stress states on fracture permeability and fluid flow, and although they are frequently found at depth in a reservoir context, such stress conditions are poorly investigated, particularly in terms of stress perturbations.The presented analogue experiments consisted of compression tests done on polymethylmethacrylate (PMMA) plates; this material has mechanical properties comparable to those of brittle rocks in the upper crust, and presents birefringence. The samples contained open defects acting as faults, and the stress trajectories around these faults were investigated using a photoelastic device. Comparable numerical experiments were realised with a finite-element code (Franc 2D), using frictionless fault models.First, the effect of an increasing biaxial compressive load ratio σ₃/σ₁ on stress trajectories around an isolated open defect was explored. It was shown that the stress trajectories were drastically modified when σ₃/σ₁ increased from 0.2 to 0.4, this result being consistent with previous studies. In particular, when σ₃/σ₁ was superior or equal to 0.4, external isotropic points around where trajectories diverged (called repulsive isotropic points) appeared near the tips, but away from the defects. They tended to move away from the defects towards the main load direction when σ₃/σ₁ increased. The described isotropic points were points of stress trajectory divergence, i.e. points where stress decreased, implying that zones around them were strongly unfavourable to shear reactivation.Second, stress trajectories around fault models of various geometries (oblique isolated defect, dilational jog, compressive jog, and complex patterns) were studied, the applied biaxial compressive load ratio being 0.7. These biaxial stress trajectories were compared with similar uniaxial stress trajectories in order to explore the existing differences between the two regimes. It enabled new stress trajectory geometries to be described and interpreted. In particular, numerous external repulsive isotropic points were observed, and defect tips were shown to be zones of high convergence of stress trajectories. Furthermore, in contrast to the uniaxial compressive load case, stress trajectories and stress contours were geometrically similar within dilational and compressive jogs under biaxial compressive load. For both jog types, the centre of the overlapping zone, and the areas along the wall of a jog-defining fault and facing the overlapping tip of the other fault, proved to be zones of low mean stress, implying that fluids may migrate towards them from the tips of the faults, in response to mean stress gradients. Furthermore, the centre of the jogs exhibiting a high differential stress was a favourable area in terms of fracture reactivation, which may facilitate fluids transfer and storage.It was also observed that for both the uniaxial and biaxial compressive loads, isotropic zones were localised at the acute angle between branching defects.     

Stress tensor analysis of the 1998–1999 tectonic swarm of northern Quito related to the volcanic swarm of Guagua Pichincha volcano, Ecuador

The phreatic activity and the subsequent dacitic dome growth in 1998–1999 at Guagua Pichincha volcano, Ecuador, were associated with two seismic swarms: one located in the northern part of Quito (population: 1,500,000) and another one, just below the active volcano, about 15–20 km SW from the first one. Quito swarm tectonic events have high frequencies (from 1 to 10–15 Hz). We registered more than 3200 events (among which 2354 events of 1.4≤M_L≤4.2) between June 1998 and December 1999 at the −2- and −17-km depth. The volcanic events below the Guagua Pichincha caldera have high (from 1 to 10–15 Hz) and low (less than 3 Hz) frequencies. Approximately, 130,000 events were registered between September 1998 and December 1999 at the +2.4- and −3.5-km depth. Here, we study the stress tensors of these two swarms deduced from the polarities of P first motions and compare them to the regional stress tensor deduced from CMT Harvard focal mechanisms. The Quito swarm stress tensor is relatively close to the regional stress tensor (the σ₁ axis was oriented N117°E close to the N102°E direction of the plate motion found by the GPS measurement, and σ₃ is nearly vertical). The difference may be due to the action of the closely active Guagua Pichincha volcano. The Guagua Pichincha stress tensor is very different from the regional tectonic one. The σ₁ axis of the volcano is oriented N214°E, almost perpendicular to the σ₁ of the swarm of Quito and σ₃ is almost horizontal. Even if these two tensors are different, they can be explained in a more general tectonic scheme. The almost horizontal direction of σ₃ just below the volcano is compatible with an extensional horizontal direction that may be expected in the shallow extrados part of a compressional region and consistent with an opening of the top of the Guagua Pichincha volcano. The movement of the fluids (magma, gas and/or groundwater) produced by the closely active Guagua Pichincha volcano seems to have an influence in the acceleration of the generation of seismic events.     

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.     

Palaeomagnetism and magnetostratigraphy of the Permian–Triassic northwest central Siberian Trap Basalts

A detailed palaeomagnetic and magnetostratigraphic study of the Permian–Triassic Siberian Trap Basalts (STB) in the Noril'sk and Abagalakh regions in northwest Central Siberia is presented. Thermal (TH) and alternating field (AF) demagnetisation techniques have been used and yielded characteristic magnetisation directions. The natural remanent magnetisation of both surface and subsurface samples is characterised by a single component in most cases. Occasionally, a viscous overprint can be identified which is easily removed by TH or AF demagnetisation.The resulting average mean direction after tectonic correction for the 95 flows sampled in outcrops is D=93.7°, I=74.7° with k=19 and α₉₅=3.3°. The corresponding pole position is 56.2°N, 146.0°E.Unoriented samples from four boreholes cores in the same regions have also been studied. They confirm the reversed–normal succession found in outcrops. The fact that only one reversal of the Earth's magnetic field has been recorded in the traps can be taken as evidence for a rather short time span for the major eruptive episode in this region. However, there is evidence elsewhere that the whole volcanic activity associated with the emplacement of the STB was much longer and lasted several million years.