Rupture progression along discontinuous oblique fault sets: implications for the Karadere rupture segment of the 1999 Izmit earthquake, and future rupture in the Sea of Marmara

Large earthquakes in strike-slip regimes commonly rupture fault segments that are oblique to each other in both strike and dip. This was the case during the 1999 Izmit earthquake, which mainly ruptured E–W-striking right-lateral faults but also ruptured the N60°E-striking Karadere fault at the eastern end of the main rupture. It will also likely be so for any future large fault rupture in the adjacent Sea of Marmara. Our aim here is to characterize the effects of regional stress direction, stress triggering due to rupture, and mechanical slip interaction on the composite rupture process. We examine the failure tendency and slip mechanism on secondary faults that are oblique in strike and dip to a vertical strike-slip fault or “master” fault. For a regional stress field well-oriented for slip on a vertical right-lateral strike-slip fault, we determine that oblique normal faulting is most favored on dipping faults with two different strikes, both of which are oriented clockwise from the strike-slip fault. The orientation closer in strike to the master fault is predicted to slip with right-lateral oblique normal slip, the other one with left-lateral oblique normal slip. The most favored secondary fault orientations depend on the effective coefficient of friction on the faults and the ratio of the vertical stress to the maximum horizontal stress. If the regional stress instead causes left-lateral slip on the vertical master fault, the most favored secondary faults would be oriented counterclockwise from the master fault. For secondary faults striking ±30° oblique to the master fault, right-lateral slip on the master fault brings both these secondary fault orientations closer to the Coulomb condition for shear failure with oblique right-lateral slip. For a secondary fault striking 30° counterclockwise, the predicted stress change and the component of reverse slip both increase for shallower-angle dips of the secondary fault. For a secondary fault striking 30° clockwise, the predicted stress change decreases but the predicted component of normal slip increases for shallower-angle dips of the secondary fault. When both the vertical master fault and the dipping secondary fault are allowed to slip, mechanical interaction produces sharp gradients or discontinuities in slip across their intersection lines. This can effectively constrain rupture to limited portions of larger faults, depending on the locations of fault intersections. Across the fault intersection line, predicted rakes can vary by >40° and the sense of lateral slip can reverse. Application of these results provides a potential explanation for why only a limited portion of the Karadere fault ruptured during the Izmit earthquake. Our results also suggest that the geometries of fault intersection within the Sea of Marmara favor composite rupture of multiple oblique fault segments.     

Spatial variations in the similarity of earthquake populations: The case of the 1997 Colfiorito–Sellano (northern Apennines, Italy) seismic sequence

The spatial properties of events in the 1997 Colfiorito–Sellano seismic sequence (Northern Apennines, Italy) were investigated using coherence, a parameter derived from seismic moment tensors that quantifies the kinematic similarity between focal mechanisms. The 1997 Colfiorito–Sellano seismic sequence predominantly consists of normal faulting earthquakes, with a few strike-slip and reverse faulting episodes. This kinematic heterogeneity is possibly related to the contemporaneous activity of two different sets of faults: NW–SE normal faults and NNE–SSW sub-vertical faults, the latter inherited from the previous Miocene compressional phase. The study used two independently-derived data sets of the same seismic sequence characterized by a different number of events and by different precision of spatial localisation. Their statistical significances, assessed through a reshuffling procedure, reveal that data sets with at least some hundreds of events and good positional precision are required to obtain significant results through coherence analysis. Results from the better quality data set indicate that this seismic sequence is characterized by a rapid decrease in the kinematic similarity between earthquake pairs within 2 km of separation, particularly along directions sub-perpendicular to the normal fault strike. The decrease rate seems to be controlled by the geometric characteristics of the normal faults, given that the mean along-dip distance between fault segments is 2 km. In proximity to pre-existing tectonic lineaments the relative abundance of strike-slip and reverse faults tends to decrease the kinematic similarity between events but does not influence the coherence decrease rate. The presence of mixed focal mechanisms (normal, reverse and strike-slip) in a single seismic phase implies that mixed fault types are not restricted to polyphase tectonic histories: such heterogeneous kinematics during a single phase may be induced by the presence of inherited discontinuities.     

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.     

Focal mechanism inversion in the Giudicarie–Lessini seismotectonic region (Southern Alps, Italy): Insights on tectonic stress and strain

A set of 41 focal mechanisms (1989–2006) from P-wave first polarities is computed from relocated seismic events in the Giudicarie–Lessini region (Southern Alps). Estimated hypocentral depths vary from 3.1 to 20.8 km, for duration magnitudes (M_D) in the range 2.7–5.1. Stress and strain inversions are performed for two seismotectonic zones, namely G (Giudicarie) and L (Lessini). This subdivision is supported by geological evidence, seismicity distribution, and focal mechanism types. The available number of data (16 in G, 22 in L) does not make possible any further subdivisions. Seismotectonic zones G and L are undergoing different kinematic regimes: thrust with strike-slip component in G, and strike-slip in L. Principal stress and strain axes in each sub-region show similar orientations. The direction of maximum horizontal compressive stress is roughly perpendicular to the thrust fronts along the Giudicarie Belt in zone G, and compatible with right-lateral strike-slip reactivation of the faults belonging to the Schio-Vicenza system in zone L. On the whole, kinematic regimes and horizontal stress orientations show a good fit with other stress data from focal mechanisms and breakouts and with geodetic strain rate axes.     

Induced seismicity in the Castanhão reservoir, NE Brazil — Preliminary results

The Castanhão reservoir was built in the state of Ceará, a dry region in Northeastern Brazil, to regulate the flow of the Jaguaribe River, for irrigation, and for power generation. It is an earth-filled dam, 60 m high, with a water capacity of 4.5 × 10⁹ m³. The seismicity in the area has been monitored since 1998, with a few interruptions, using one analog or one digital station and, during a few periods, a three-station network. The first earthquakes likely to be induced events were detected in 2003, when the water level was about 20 m high. In early 2004 a very heavy rainfall season quickly filled the reservoir. Shortly after, an increase in the seismic activity occurred and many micro-earthquakes were recorded. We suggest that this activity resulted from an increase in pore pressure due to undrained response. Therefore, we may classify this cluster of micro-earthquakes as “initial seismicity.” We deployed a network with four analog stations in the area, following this activity, to determine the epicentral zone. At least three epicentral areas under the reservoir were detected. The spatio-temporal analysis of the available data revealed that the seismicity occurs in clusters and that these were activated at different periods. We identified four sets of faults (N–S-, E–W-, NW–SE-, and NE–SW-oriented), some of which moved in shallow crustal levels and as recently as the Quaternary (1.8 Ma). Under the present-day stress regime, the last two sets moved as strike-slip structures. We suggest a possible correlation between dormant faults and the observed induced seismicity.     

The April 2007 earthquake swarm near Lake Trichonis and implications for active tectonics in western Greece

We investigate the properties of the April 2007 earthquake swarm (Mw 5.2) which occurred at the vicinity of Lake Trichonis (western Greece). First we relocated the earthquakes, using P- and S-wave arrivals to the stations of the Hellenic Unified Seismic Network (HUSN), and then we applied moment tensor inversion to regional broad-band waveforms to obtain the focal mechanisms of the strongest events of the 2007 swarm. The relocated epicentres, cluster along the eastern banks of the lake, and follow a distinct NNW–ESE trend. The previous strong sequence close to Lake Trichonis occurred in June–December 1975. We applied teleseismic body waveform inversion, to obtain the focal mechanism solution of the strongest earthquake of this sequence, i.e. the 31 December 1975 (Mw 6.0) event. Our results indicate that: a) the 31 December 1975 Mw 6.0 event was produced by a NW–SE normal fault, dipping to the NE, with considerable sinistral strike-slip component; we relocated its epicentre: i) using phase data reported to ISC and its coordinates are 38.486°N, 21.661°E; ii) using the available macroseismic data, and the coordinates of the macroseismic epicentre are 38.49°N, 21.63°E, close to the strongly affected village of Kato Makrinou; b) the earthquakes of the 2007 swarm indicate a NNW–SSE strike for the activated main structure, parallel to the eastern banks of Lake Trichonis, dipping to the NE and characterized by mainly normal faulting, occasionally combined with sinistral strike-slip component. The 2007 earthquake swarm did not rupture the well documented E–W striking Trichonis normal fault that bounds the southern flank of the lake, but on the contrary it is due to rupture of a NW–SE normal fault that strikes at a  45° angle to the Trichonis fault. The left-lateral component of faulting is mapped for the first time to the north of the Gulf of Patras which was previously regarded as the boundary for strike-slip motions in western Greece. This result signifies the importance of further investigations to unravel in detail the tectonics of this region.     

Brittle–ductile–brittle deformation during cooling of tonalite (Adamello, Southern Italian Alps)

Late- to post-magmatic deformation in slightly diachronous contiguous intrusions of the north-western Adamello batholith (Southern Alps, Italy) is recorded as, from oldest to youngest: (i) joints, (ii) solid-state ductile shear zones, (iii) faults associated with epidote-K-feldspar veins and (iv) zeolite veins and faults. Structures (ii) to (iv) are localized on the pervasive precursory network of joints (i), which developed during the earliest stages of pluton cooling. High temperature ( 500 °C), ductile overprinting of joints produced lineations, defined by aligned biotite and hornblende, on the joint surfaces and highly localized mylonites. The main phase of faulting, producing cataclasites and pseudotachylytes, occurred at  250 °C and was associated with extensive fluid infiltration. Cataclasites and pseudotachylytes are clustered along different E–W-striking dextral strike-slip fault zones correlated with the activity of the Tonale fault, a major tectonic structure that bounds the Adamello batholith to the north. Ductile deformation and cataclastic/veining episodes occurred at P = 0.25–0.3 GPa during rapid cooling of the batholith to the ambient temperatures ( 250 °C) that preceded the exhumation of the batholith. Timing of the sequence of deformation can be constrained by ³⁹Ar–⁴⁰Ar ages of  30 Ma on pseudotachylytes and various existing mineral ages. In the whole composite Adamello batholith, multiple magma pulses were intruded over the time span 42–30 Ma and each intrusive body shows the same ductile-to-brittle structural sequence localized on the early joint sets. This deformation sequence of the Adamello might be typical of intrusions undergoing cooling at depths close to the brittle–ductile transition.     

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.     

Surface kinematics in the Alpine–Carpathian–Dinaric and Balkan region inferred from a new multi-network GPS combination solution

The understanding of the intraplate tectonics of Central Europe requires a detailed picture of how stress is transferred from the interaction of the Eurasian, Nubian and Anatolian plates to the Alpine, Carpathian, Pannonian and Dinaric regions. Recent strain distribution is controlled by the Adria horizontal push, by the Vrancea vertical slab pull and associated horizontal displacements, and by the Aegean/Anatolia extension and slab-roll back. We present a horizontal velocity field for the Alpine-Carpathian-Pannonic-Dinaric and Balkan regions resulting from a new combination of seven different GPS networks formed from permanent and campaign stations. Dedicated velocity profiles in two specific regions are studied in detail. One is the Alpine Pannonian region, with a detailed picture of the NS indentation of the Adria microplate into the Southern Alps, in NE Italy, the deformation in the Tauern Window and the eastwards kinematics of a Pannonian plate fragment. The second study region includes Transylvania, the Southern Carpathians up to the Aegean sea and Albania, where a major right lateral shear deformation exists as a consequence of the NE convergence of the Apulia platform towards the Dinarids, and the SSW motion of Macedonia, Western Bulgaria and Rumania, related to the Hellenic arc dynamics in the Eastern Mediterranean. The profiles in the Alpine–Pannonian area indicate that a velocity drop of 2.5 +/− 0.4 mm/yr associated with the Adria indentation concentrates on a segment of some 50 km south of the Periadriatic fault. The deformation becomes extensional by a similar amount just north of the Periadriatic fault, in the Tauern Window, where the updoming of the Tauern Window implies vertical motion which could well be associated with surface extension. In the EW profile, we observe a sudden velocity change of 1.5 +/− 0.2 mm/yr in 20 km, corresponding to the right lateral Lavant fault, which seems to mark the border between dominant indentation kinematics to the West and dominant extrusion kinematics to the East.Three profiles are considered in Southern and Eastern Europe: one across the lower Adriatic sea from Apulia in Italy to the southern Dinarides, which enables it to constrain the velocity drop associated with the subduction of the Adria microplate into the Dinarides to 3.2 +/− 0.5 mm/yr in 140 km. The second profile is longitudinal and constrains the velocity inversion of 7.4 +/− 1.0 mm/yr in 350 km associated with right lateral shear faults in Albania. The third profile crosses the Transylvania region with a shortening of 2.3 +/− 1.0 mm/yr in 220 km, and the Wallachian–Moesian region up to the Chalcidic peninsula in N Greece. This lower part of the profile implies an extensional stretch of the upper crust of 3.2 +/− 0.9 mm/yr in 440 km, culminating in the Hellenic arc. Strain rate maps are presented in this regional scale, showing the excellent agreement between fault plane solutions of crustal earthquakes and the eigenvectors of the GPS derived two dimensional strain rate tensor.Three profiles are considered in the Balkan and SE Carpathians: one across the lower Adriatic sea from Apulia in Italy to the southern Dinarides, which enables to constrain the velocity drop associated to the subduction of the Adria microplate into the Dinarides to 3.2 +/− 0.5 mm/yr in 140 km. The second profile is longitudinal and constrains the velocity inversion of 7.4 +/− 1.0 mm/yr in 350 km associated to right lateral shear faults in Macedonia, a highly seismic region. The third profile crosses the Transylvania with a shortening2.3 +/− 1.0 mm/yr in 220 km, and the Wallachian–Moesian region up to the Chalcidic peninsula in N Greece. This lower part of the profile implies an extensional stretch of the upper crust of 3.2 +/− 0.9 mm/yr in 440 km, culminating in the Hellenic arc.     

Local seismic array observations at north Evoikos, central Greece, delineate crustal deformation between the North Aegean Trough and Corinthiakos Rift

In order to better constrain and define the microseismic activity at the north Evoikos Gulf and its surrounding area we deployed an onshore/offshore seismic array consisting of 31 three-component seismic digital stations. The array was active from 30 June to 24 October 2003, and covered an area of 2500 km². We located more than 2000 seismic events ranging from 0.7 to 4.5 M_L by using six stations as a minimum in order to define the foci parameters. Recorded seismicity delineated three major zones of deformation: from south to north, the Eretria–Parnis–eastern Corinthiakos zone, the Psachna–Viotia zone, and the Northern Sporades–North Evia–Bralos zone. Alignments of the recorded seismicity follow the tectonic trends and their orientation in the above zones. The whole area accommodates the stress field between the North Aegean Trough and the Corinthiakos Gulf. Rate of deformation intensifies from north to south, as revealed also by historical and instrumental seismicity. The successive change of orientation between the two stress fields fragments the crust in relatively small units and the fault systems developed do not permit the generation of major earthquakes in the north Evoikos area and its immediate vicinity. This is also supported by the instrumental seismicity of the last century. Larger events reported in historical times are probably overestimated.Most seismic activity is crustal. Subcrustal events were recorded mainly below the Lichades area and are interpreted as the consequence of the subduction of the Ionian oceanic lithosphere below the Hellenides. The Lichades volcano is the most northern end of the Hellenic volcanic arc.At present the highest seismic activity is associated with the Psachna region of north Evia that has been continuously active since 2001. Considering, however, the development of the seismic activity during the last decade, there has been a sequence of large events, i.e., Parnis in 1999, Skyros in 2001 and Psachna in 2001–2003. This demonstrates the fact that the tectonic deformation in all this area is intense and important for the accommodation of the stress field of the North Aegean Trough to that of the Corinthiakos Rift.     

The role of the Kazakhstan orocline in the late Paleozoic amalgamation of Eurasia

After the 2005 Kashmir earthquake, we mapped surface ground fractures in Tangdhar, Uri, Rajouri and Punch sectors and liquefaction features in Jammu area lying close to the eastern side of the Line of Control (LOC) in Kashmir, India. The NW trending ground fractures occurred largely in the hanging wall zone of the southeastern extension of the causative fault in Tangdhar and Uri sectors. The principal compressive stress deduced from the earthquake induced ground fractures is oriented at N10°, whereas the causative Balakot–Bagh fault strikes 330°. The fault-plane solution indicates primarily SW thrusting of the causative fault with a component of strike–slip motion. The ground fractures reflect pronounced strike–slip together with some tensile component. The Tangdhar area showing left-lateral strike–slip motion lies on the hanging wall, and the Uri region showing right-lateral strike–slip movement is located towards the southeastern extension of the causative fault zone. The shear fractures are related to static stress that was responsible for the failure of causative fault. The tensile fractures with offsets are attributed to combination of both static and dynamic stresses, and the fractures and openings without offsets owe their origin due to dynamic stress. In Punch–Rajouri and Jammu area, which lies on the footwall, the fractures and liquefactions were generated by dynamic stress. The occurrence of liquefaction features in the out board part of the Himalayan range front near Jammu is suggestive of stress transfer  230 km southeast of the epicenter. The Balakot–Bagh Fault (BBF), the Muzaffarabad anticline, the rupture zone of causative fault and the zone of aftershocks — all are aligned in a  25 km wide belt along the NW–SE trending regional Himalayan strike of Kashmir region and lying between the MBT and the Riasi Thrust (Murree Thrust), suggesting a seismogenic zone that may propagate towards the southeast to trigger an earthquake in the eastern part of the Kashmir region.     

Active faulting and seismicity along the Siculo–Calabrian Rift Zone (Southern Italy)

Southern Italy is dominated by extensional tectonics that in the Calabrian arc and Eastern Sicily produced the development of the Siculo–Calabrian Rift Zone (SCRZ). This zone is represented by a ≈ 370 km-long fault belt consisting of 10 to 50 km long distinct fault segments which extend both offshore and on land being also responsible of the crustal seismicity of this region. The geological and morphological observations indicate that the active normal faults of the SCRZ are characterized by throw-rates ranging from 0.7 to 3.1 mm/a. They accommodate an almost uniform horizontal extension-rate of about 3.0 mm/a along a WNW–ESE regional extension direction. Based on our field observations and following empirical relationships between magnitude and surface rupture length connections between large crustal earthquakes and distinct fault segments of the SCRZ have been also tentatively tested. Our data indicate moreover that the magnitudes (M) of the historical and instrumental earthquakes are consistent with the estimated values and that the geometry and kinematics of the fault segments and the related different crustal features of the SCRZ control the different seismic behaviours of adjacent portions of the active rift zone.     

Evidence for coupled reverse and normal active faulting in W Iberia: The Vidigueira–Moura and Alqueva faults (SE Portugal)

The Vidigueira–Moura fault (VMF) is a 65 km long, E–W trending, N dipping reverse left-lateral late Variscan structure located in SE Portugal (W Iberia), which has been reactivated during the Cenozoic with reverse right-lateral slip. It is intersected by, and interferes with the NE–SW trending Alentejo–Plasencia fault. East of this intersection, for a length of 40 km the VMF borders an intracratonic tectonic basin on its northern side, thrusting Paleozoic schists, meta-volcanics and granites, on the north, over Cenozoic continental sediments preserved in the basin, on the south. West of the faults intersection, evidence of Cenozoic reactivation is scarce. In the eastern sector, Plio-Quaternary VMF reactivation is indicated by geomorphologic, stratigraphic, and structural data, showing reverse movement with a right-lateral strike-slip component, in response to a NW–SE trending compressive stress. An average vertical displacement rate of 0.06 to 0.08 mm/yr since late Pliocene (roughly the last 2.5 Ma) is estimated. The Alqueva fault (AF) is a subparallel, northward dipping, 7.5 km long anastomosing fault zone that affects Palaeozoic basement rocks, and is located 2.5 km north and on the hanging block of the VMF. The AF is also a reverse left-lateral late Variscan structure, which has been reactivated during the Tertiary with reverse right-lateral slip; however, Plio-Quaternary reactivation was normal left-lateral, as shown by abundant kinematical criteria (slickensides) and geomorphic evidence. It shows an average displacement rate of 0.02 mm/yr for the vertical component of movement in the approximately last 2.5 Ma. It is proposed that the normal displacements on the AF result from tangential longitudinal strain on the upthrown block of the VMF above a convex ramp of this main reverse structure. According to this model of faults interaction, the AF is interpreted to work as a bending-moment fault sited above the VMF thrust ramp. Consequently, it is expected that the displacements on the AF increase towards the topographic surface with the increase in the imposed extension, declining downwards until they vanish above or at the VMF ramp. In order to constrain the proposed scheme, numerical modeling was performed, aiming at the reproduction of the present topography across the faults using different geodynamic models and fault geometries and displacements.     

Thessaloniki–Gerakarou Fault Zone (TGFZ): the western extension of the 1978 Thessaloniki earthquake fault (Northern Greece) and seismic hazard assessment

Active faulting and seismic properties are re-investigated in the eastern precinct of the city of Thessaloniki (Northern Greece), which was seriously affected by two large earthquakes during the 20th century and severe damage was done by the 1759 event. It is suggested that the earthquake fault associated with the occurrence of the latest destructive 1978 Thessaloniki earthquake continues westwards to the 20-km-long Thessaloniki–Gerakarou Fault Zone (TGFZ), which extends from the Gerakarou village to the city of Thessaloniki. This fault zone exhibits a constant dip to the N and is characterised by a complicated geometry comprised of inherited 100°-trending faults that form multi-level branching (tree-like fault geometry) along with NNE- to NE-trending faults. The TGFZ is compatible with the contemporary regional N–S extensional stress field that tends to modify the pre-existing NW–SE tectonic fabric prevailing in the mountainous region of Thessaloniki. Both the 1978 earthquake fault and TGFZ belong to a ca. 65-km-long E–W-trending rupture fault system that runs through the southern part of the Mygdonia graben from the Strymonikos gulf to Thessaloniki. This fault system, here called Thessaloniki–Rentina Fault System (TRFS), consists of two 17–20-km-long left-stepping 100°-trending main fault strands that form underlapping steps bridged by 8–10-km-long ENE–WSW faults. The occurrence of large (M6.0) historical earthquakes (in 620, 677 and 700 A.D.) demonstrates repeated activation, and therefore the possible reactivation of the westernmost segment, the TGFZ, could be a major threat to the city of Thessaloniki. Changes in the Coulomb failure function (ΔCFF) due to the occurrence of the 1978 earthquake calculated out in this paper indicate that the TGFZ has been brought closer to failure, a convincing argument for future seismic hazard along the TGFZ.     

Seismic faults in the Aegean area

Reliable fault plane solutions of shallow earthquakes and information on surface fault traces in combination with other seismic, geomorphological and geological information have been used to determine the orientation and other properties of the seismic faults in the Aegean and surrounding area.Thrust faults having an about NW-SE strike occur in the outer seismic zone along western Albania-westernmost part of mainland of Greece-Ionian Sea-south of Crete-south of Rhodes.The inner part of the area is dominated by strike-slip and normal faulting. Strike-slip with an about NE-SW slip direction occurs in the inner part of the Hellenic arc along the line Peloponnesus-Cyclades-Dodecanese-southwest Turkey as well as along a zone which is associated with the northern Aegean trough and the northwesternmost part of Anatolia. All other regions in the inner part of the area are characterized by normal faulting. The slip direction of the normal faults has an about SW-NE direction in Crete (N38°E) and an about E-W direction (N81°E) in a zone which trends N-S in eastern Albania and its extension to western mainland of Greece. In all other regions (central Greece-southern Yugoslavia and Bulgaria, western Turkey) the slip of the normal faults has an about N-S direction.     

The 1984 Abruzzo earthquake (Italy): an example of seismogenic process controlled by interaction between differently oriented synkinematic faults

The evolution of the seismogenic process associated with the M_s 5.8 Sangro Valley earthquake of May 1984 (Abruzzo, central Italy) is closely controlled by the Quaternary extensional tectonic pattern of the area. This pattern is characterised by normal faults mainly NNW striking, whose length is controlled by pre-existing Mio–Pliocene N100±10° left-lateral strike-slip fault zones. These are partly re-activated as right-lateral normal-oblique faults under the Quaternary extensional regime and behave as transfer faults.Integration of re-located aftershocks, focal mechanisms and structural features are used to explain the divergence between the alignment of aftershocks (WSW–ENE) and the direction of seismogenic fault planes defined by the focal mechanisms (NNW–SSE) of the main shock and of the largest aftershock (M_s=5.3).The faults that appear to be involved in the seismogenic process are the NNW–SSE Barrea fault and the E–W M. Greco fault. There is field evidence of finite Quaternary deformation indicating that the normal Barrea fault re-activates the M. Greco fault as right-lateral transfer fault. No surface faulting was observed during the seismic sequence. The apparently incongruent divergence between aftershocks and nodal planes may be explained by interpreting the M. Greco fault as a barrier to the propagation of earthquake rupturing. The rupture would have nucleated on the Barrea fault, migrating along-strike towards NNW. The sharp variation in direction from the Barrea to the M. Greco fault segments would have represented a structural complexity sufficient to halt the rupture and subsequent concentration of post-seismic deformation as aftershocks around the line of intersection between the two fault planes.Fault complexities, similar to those observed in the Sangro Valley, are common features of the seismic zone of the Apennines. We suggest that the zones of interaction between NW–SE and NNW–SSE Plio-Quaternary faults and nearly E–W transfer faults, extending for several kilometres in the same way as M. Greco does, might act as barriers to the along-strike propagation of rupture processes during normal faulting earthquakes. This might have strong implications on seismic hazard, especially for the extent of the maximum magnitude expected on active faults during single rupture episodes.     

Relocation of aftershocks, focal mechanisms and stress inversion: Implications toward the seismo-tectonics of the causative fault zone of Mw7.6 2001 Bhuj earthquake (India)

The HYPODD relocation of 1172 aftershocks, recorded on 8–17 three-component digital seismographs, delineate a distinct south dipping E–W trending aftershock zone extending up to 35 km depth, which involves a crustal volume of 40 km × 60 km × 35 km. The relocated focal depths delineate the presence of three fault segments and variation in the brittle–ductile transition depths amongst the individual faults as the earthquake foci in the both western and eastern ends are confined up to 28 km depth whilst in the central aftershock zone they are limited up to 35 km depth. The FPFIT focal mechanism solutions of 444 aftershocks (using 8–12 first motions) suggest that the focal mechanisms ranged between pure reverse and pure strike slip except some pure dip slip solutions. Stress inversion performed using the P and T axes of the selected focal mechanisms reveals an N181°E oriented maximum principal stress with a very shallow dip (= 14°). The stress inversions of different depth bins of the P and T axes of selected aftershocks suggest a heterogeneous stress regime at 0–30 km depth range with a dominant consistent N–S orientation of the P-axes over the aftershock zone, which could be attributed to the existence of varied nature and orientation of fractures and faults as revealed by the relocated aftershocks.     

Crustal deformation in the Central Ionian Islands (Greece): Results from DGPS and DInSAR analyses (1995–2006)

Ground deformation studies based on Differential GPS (DGPS) and Differential Interferometric SAR (DInSAR) analyses have been conducted in the seismically active area of the Central Ionian Islands. Local GPS networks were installed in Cephallonia (2001) and Zakynthos (2005). The Cephallonian network has been remeasured five times and Zakynthos' once as of July 2006. The studies have yielded detailed information regarding both local and regional deformations that are occurring in the area.For Lefkas Island, DInSAR analysis (March to September 2003) revealed 56 mm of uplift in the central and western parts and is attributed to the August 2003 earthquake (M_w = 6.3) that occurred offshore to the west. Synthetic DInSAR modelling of the magnitude and extent of deformation is consistent with the seismologically deduced parameters for the ruptured segment along the Lefkas Transform Fault. Subsidence (< 28 mm) along the northern part of the island is attributed to local conditions unrelated to the earthquake. For Zakynthos Island, large-magnitude earthquakes that occurred offshore to the south in October 2005 and April 2006 most likely contributed to the observed deformation as deduced from DGPS measurements for an encompassing period (August 2005 to July 2006). The largest amount of horizontal deformation occurred in the south, where its western part moved in a W–NW direction, while the eastern part moved towards the NE, with magnitudes ranging from 15 to 26 mm. The southern part of the island uplifted a maximum of 65 mm whereas the north subsided from 12 to 28 mm.For Cephallonia Island, DInSAR analysis (1995 to 1998) indicated ground deformation up to 28 mm located in small sections of the island. Further interferometric analysis for the period 2003 to 2004, encompassing the occurrence of the Lefkas earthquake in August 2003, indicated 28 mm of uplift in the northern part, while during the next two years (2004 to 2005), further uplift of at least 56 mm had taken place at the western and northern part of the island.DGPS measurements for the period 2001 to 2006 revealed a clockwise rotation of the island with respect to a centrally located station on Aenos Mt. The horizontal component of deformation generally ranged from 6 to 34 mm, with the largest values at the western and northern parts of the island. Considering the vertical deformation, two periods are distinguished. The first one (2001 to 2003) is consistent with anticipated motions associated with the main geological and tectonic features of the island. The second one (2003 to 2006) has been tentatively attributed to dilatancy in which relatively small uplift (20–40 mm) occurred along the southern and southeastern parts of the island, while larger magnitudes (> 50 mm) happened at the western part (Paliki Peninsula). These large magnitudes of uplift over an extended area (> 50 km), in conjunction with an accelerated Benioff strain determined from the analysis of the seismicity in the broader region, are consistent with dilatancy. This effect commenced some time after 2003 and is probably centered in the area between Zakynthos and Cephallonia. If this interpretation is correct, it may foreshadow the occurrence of a very strong earthquake(s) sometime during 2007 to 2008 in the above designated region.     

Kinematics of the Iberia–Maghreb plate contact from seismic moment tensors and GPS observations

The Iberian Peninsula and the Maghreb experience moderate earthquake activity and oblique,  NW–SE convergence between Africa and Eurasia at a rate of  5 mm/yr. Coeval extension in the Alboran Basin and a N35°E trending band of active, left-lateral shear deformation in the Alboran–Betic region are not straightforward to understand in the context of regional shortening, and evidence complexity of deformation at the plate contact. We estimate 86 seismic moment tensors (M_W 3.3 to 6.9) from time domain inversion of near-regional waveforms in an intermediate period band. Those and previous moment tensors are used to describe regional faulting style and calculate average stress tensors. The solutions associated to the Trans-Alboran shear zone show predominantly strike-slip faulting, and indicate a clockwise rotation of the largest principal stress orientation compared to the regional convergence direction (σ₁ at N350°E). At the N-Algerian and SW-Iberian margins, reverse faulting solutions dominate, corresponding to N350°E and N310°E compression, respectively. Over most of the Betic range and intraplate Iberia, we observe predominately normal faulting, and WSW–ENE extension (σ₃ at N240°E). From GPS observations we estimate that more than 3 mm/yr of African (Nubian)–Eurasian plate convergence are currently accommodated at the N-Algerian margin,  2 mm/yr in the Moroccan Atlas, and  2 mm/yr at the SW-Iberian margin. 2 mm/yr is a reasonable estimate for convergence within the Alboran region, while Alboran extension can be quantified as  2.5 mm/yr along the stretching direction (N240°E). Superposition of both motions explains the observed left-lateral transtensional regime in the Trans-Alboran shear zone. Two potential driving mechanisms of differential motion of the Alboran–Betic–Gibraltar domain may coexist in the region: a secondary stress source other than plate convergence, related to regional-scale dynamic processes in the upper mantle of the Alboran region, as well as drag from the continental-scale motion of the Nubian plate along the southern limit of the region. In the Atlantic Ocean, the  3.5 mm/yr, westward motion of the Gibraltar Arc relative to intraplate Iberia can be accommodated at the transpressive SW-Iberian margin, while available GPS observations do not support an active subduction process in this area.