Source parameters of the ML 4.1 earthquake of November 08, 2006, southeast Beni-Suef, northern Egypt

In the early morning hours on Wednesday November 08, 2006 at 04:32:10(GMT) a small earthquake of M_L 4.1 has occurred at southeast Beni-Suef, approximately 160 km SEE of Cairo, northern Egypt. The quake has been felt as far as Cairo and its surroundings while no casualties were reported. The instrumental epicentre is located at 28.57°N and 31.55°E. Seismic moment is 1.76 E14 Nm, corresponding to a moment magnitude M_w 3.5. Following a Brune model, the source radius is 0.3 km with an average dislocation of 1.8 cm and a 2.4 MPa stress drop. The source mechanism from a first motion fault plane solution shows a left-lateral strike-slip mechanism with a minor dip-slip component along fault NNW striking at 161°, dipping 52° to the west and rake −5°. Trend and plunging of the maximum and minimum principle axes P/T are 125°, 28°, 21°, and 23°, respectively. A comparison with the mechanism of the October, 1999 event shows similarities in faulting type and orientation of nodal planes.Eight small earthquakes (3.0  M_L < 5.0) were also recorded by the Egyptian National Seismological Network (ENSN) from the same region. We estimate the source parameters and fault mechanism solutions (FMS) for these earthquakes using displacement spectra and P-wave polarities, respectively. The obtained source parameters including seismic moments of 4.9 × 10¹²–5.04 × 10¹⁵ Nm, stress drops of 0.2–4.9 MPa and relative displacement of 0.1–9.1 cm. The azimuths of T-axes determined from FMS are oriented in NNE–SSW direction. This direction is consistent with the present-day stress field in Egypt and the last phase of stress field changes in the Late Pleistocene, as well as with recent GPS measurements.     

Static and dynamic fault parameters of the Saitama earthquake of July 1, 1968

The source mechanism of the Saitama earthquake (36.07°N,139.40°E, M_s = 5.4) of July 1, 1968, is studied on the basis of P-wave first motion, aftershock, long-period surface-wave data and low-magnification long-period seismograms recorded in the nearfield. A precise location of the aftershocks is made using P and S—P time data obtained by a micro-earthquake observatory network. The synthetic near-field seismograms based on the Haskell model are directly compared with the observed near-field seismograms for wave form and amplitude to determine the dynamic fault parameters. The results obtained are as follows: source geometry, reverse dip slip with considerable right-lateral strike-slip component; dip direction, N6°E; dip angle 30°; fault dimension, 10 × 6 km²; rupture velocity, 3.4 km/sec in the direction S30°E; average dislocation, 92 cm; average dislocation velocity, 92 cm/sec; seismic moment, 1.9 · 10²⁵ dyn-cm; stress drop, 100 bar. The effective stress is about the same as the stress drop. For major earthquakes in the Japanese Islands, the dislocation velocity, .D, is found to be proportional to the stress drop, σ. This relation can be expressed by .D - (β/μ)σ, where β is the shear velocity and μ is the rigidity. This result has an importance in engineering seismology because the stress drop scales the seismic motion in the vicinity of an earthquake fault.     

The 1980, 1997 and 1998 Azores earthquakes and some seismo-tectonic implications

We have studied the focal mechanisms of the 1980, 1997 and 1998 earthquakes in the Azores region from body-wave inversion of digital GDSN (Global Digital Seismograph Network) and broadband data. For the 1980 and 1998 shocks, we have obtained strike–slip faulting, with the rupture process made up of two sub-events in both shocks, with total scalar seismic moments of 1.9 × 10¹⁹ Nm (M_w = 6.8) and 1.4 × 10¹⁸ Nm (M_w = 6.0), respectively. For the 1997 shock, we have obtained a normal faulting mechanism, with the rupture process made up of three sub-events, with a total scalar seismic moment of 7.7 × 10¹⁷ Nm (M_w = 5.9). A common characteristic of these three earthquakes was the shallow focal depth, less than 10 km, in agreement with the oceanic-type crust. From the directivity function of Rayleigh (LR) waves, we have identified the NW–SE plane as the rupture plane for the 1980 and 1998 earthquakes with the rupture propagating to the SE. Slow rupture velocity, about of 1.5 km/s, has been estimated from directivity function for the 1980 and 1998 earthquakes. From spectral analysis and body-wave inversion, fault dimensions, stress drop and average slip have been estimated. Focal mechanisms of the three earthquakes we have studied, together with focal mechanisms obtained by other authors, have been used in order to obtain a seismotectonic model for the Azores region. We have found different types of behaviour present along the region. It can be divided into two zones: Zone I, from 30°W to 27°W; Zone II, from 27°W to 23°W, with a change in the seismicity and stress direction from Zone I. In Zone I, the total seismic moment tensor obtained corresponded to left-lateral strike–slip faulting with horizontal pressure and tension axes in the E–W and N–S directions, respectively. In Zone II, the total seismic moment tensor corresponded to normal faulting, with a horizontal tension axis trending NE–SW, normal to the Terceira Ridge. The stress pattern for the whole region corresponds to horizontal extension with an average seismic slip rate of 4.4 mm/yr.     

Tikhonov's regularization for deconvolution in the empirical Green function method and vertical directivity effect

Application of Tikhonov's technique, using input errors for the parameter of regularization estimation, enhances the accuracy and stability of the reconstruction of a source time function (STF) by the empirical Green function (EGF) method that gives us an opportunity to use simultaneously for analysis body and surface waves data, and to estimate the horizontal and vertical directivity effects. Knowledge of the last is particularly useful for the choice of an active nodal plane of earthquakes with the dip slip fault orientation that allows us to classify these earthquakes to the interplate or intraplate types and thereby to reach the better understanding of tectonic processes in the region of interest.By way of illustration, an attempt to estimate average parameters of faulting in a first approximation is made herein for two Russian Far East large events with opposite types of focal mechanism orientation, strike slip and dip slip. The former is not a matter of interest in the context of vertical directivity effect but enables us to test the method.The directivity analysis of pulse durations and inverse amplitudes of the relative source time functions (RSTFs) restored at eight globally distributed stations IRIS indicates that the destruction in the source of the Neftegorsk earthquake (05/27/1995 M_W=7.1) propagated roughly horizontally in the direction 8±11° during 19.2±0.4 s along the rupture extending 35.5±4.9 km. The calculated slip distribution along the rupture coincides within the error with the results of field geological measurements on the causal surface fault that proves that the Neftegorsk earthquake source is well described by the model of the linear unilateral fault and gives a good assessment of the method applied.The average parameters of faulting in the Kamchatka earthquake (03/08/1999 M_W=6.9) have been determined from data of 13 station IRIS. It was shown that the destruction in its source propagated downward at an angle of about 60° with horizon, in the direction about S156° E, during 13.4±0.2 s, along the rupture totaling 25.5±2.3 km in length. Therefore, the nodal plane, steeply dipped to the SE, was active and this event can be regarded as an intraplate type. Two asperities can be selected; the first with the maximum slip 3.3 m located at a distance of about 7 km from the onset of rupture, and the second with the maximum slip about 0.9 m centered at approximately 19 km from that.     

Macroseismic field of the March 4, 1977 Vrancea earthquake

The macroseismic field of the Vrancea earthquake of March 4, 1977, characterized by the following parameters : H = 19h 21m 56s, ø = 45.8°N, γ = 26.8° E, h = 95 Km, M = 7.2, I = VIIIMSK has been analyzed. The following problems were studied : area and shape of the isoseismals of intensity III–VIII ; elongation of the isoseismal ellipses and decrease of intensity with distance. The results confirm our previous studies (Radu and Apopei, 1978) of strong intermediate earthquakes, but render evidence for some peculiarities in the seismic intensity attenuation as well.     

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.     

The earthquake of March 4, 1977, its recording peculiarities and source parameters

A detailed analysis of recording peculiarities at seismic stations of the Uniform System of Seismic Observations (USSO) is presented a complicated nature of the source being shown. Consideration is given to parameters of the earthquake source, including the seismic moment and the length of the rupture.Comparison of magnitudes M_LH and M_PV indicates an anomalous attenuation in surface waves, itis is 3–4 times weaker than it had been noticed in case of other intermediate-depth Carpathian earthquakes.On the basis of comparison of the logarithm of the ratio of P-wave spectra at different epicentral distances (30° –70° ), the fac tor characterizing the absorption of P wave is found to remain practically unchanged.Average value of the seismic moment is estimated to be 2.6 × 10²⁷ dyne × cm, the most reasonable length of the rupture 58 km, and its focus 100 –130 km. The source parameters of the earthquake in question are compared with those of the earthquake of November 10, 1940.     

Aftershocks of the june 20, 1978, Greece earthquake: A multimode faulting sequence

A 10-station portable seismograph network was deployed in northern Greece to study aftershocks of the magnitude (m_b) 6.4 earthquake of June 20, 1978. The main shock occurred (in a graben) about 25 km northeast of the city of Thessaloniki and caused an east-west zone of surface rupturing 14 km long that splayed to 7 km wide at the west end. The hypocenters for 116 aftershocks in the magnitude range from 2.5 to 4.5 were determined. The epicenters for these events cover an area 30 km (east-west) by 18 km (north-south), and focal depths ranges from 4 to 12 km. Most of the aftershocks in the east half of the aftershock zone are north of the surface rupture and north of the graben. Those in the west half are located within the boundaries of the graben. Composite focalmechanism solutions for selected aftershocks indicate reactivation of geologically mapped normal faults in the area. Also, strike-slip and dip-slip faults that splay off the western end of the zone of surface ruptures may have been activated.The epicenters for four large (M 4.8) foreshocks and the main shock were relocated using the method of joint epicenter determination. Collectively, those five epicenters form an arcuate pattern convex southward, that is north of and 5 km distant from the surface rupturing. The 5-km separation, along with a focal depth of 8 km (average aftershock depth) or 16 km (NEIS main-shock depth), implies that the fault plane dips northward 58° or 73°, respectively. A preferred nodal-plane dip of 36° was determined by B.C. Papazachos and his colleagues in 1979 from a focal-mechanism solution for the main shock. If this dip is valid for the causal fault and that fault projects to the zone of surface rupturing, a decrease of dip with depth is required.     

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.     

Source parameters of intermediate-depth Vrancea (Romania) earthquakes from empirical Green's functions modeling

Several source parameters (source dimensions, slip, particle velocity, static and dynamic stress drop) are determined for the moderate-size October 27th, 2004 (M_W = 5.8), and the large August 30th, 1986 (M_W = 7.1) and March 4th, 1977 (M_W = 7.4) Vrancea (Romania) intermediate-depth earthquakes. For this purpose, the empirical Green's functions method of Irikura [e.g. Irikura, K. (1983). Semi-Empirical Estimation of Strong Ground Motions during Large Earthquakes. Bull. Dis. Prev. Res. Inst., Kyoto Univ., 33, Part 2, No. 298, 63–104., Irikura, K. (1986). Prediction of strong acceleration motions using empirical Green's function, in Proceedings of the 7th Japan earthquake engineering symposium, 151–156., Irikura, K. (1999). Techniques for the simulation of strong ground motion and deterministic seismic hazard analysis, in Proceedings of the advanced study course seismotectonic and microzonation techniques in earthquake engineering: integrated training in earthquake risk reduction practices, Kefallinia, 453–554.] is used to generate synthetic time series from recordings of smaller events (with 4 ≤ M_W ≤ 5) in order to estimate several parameters characterizing the so-called strong motion generation area, which is defined as an extended area with homogeneous slip and rise time and, for crustal earthquakes, corresponds to an asperity of about 100 bar stress release [Miyake, H., T. Iwata and K. Irikura (2003). Source characterization for broadband ground-motion simulation: Kinematic heterogeneous source model and strong motion generation area. Bull. Seism. Soc. Am., 93, 2531–2545.] The parameters are obtained by acceleration envelope and displacement waveform inversion for the 2004 and 1986 events and MSK intensity pattern inversion for the 1977 event using a genetic algorithm. The strong motion recordings of the analyzed Vrancea earthquakes as well as the MSK intensity pattern of the 1977 earthquake can be well reproduced using relatively small strong motion generation areas, which corresponds to small asperities with high stress drops (300–1200 bar) and high particle velocities (3–5 m/s). These results imply a very efficient high-frequency radiation, which has to be taken into account for strong ground motion prediction, and indicate that the intermediate-depth Vrancea earthquakes are inherently different from crustal events.     

On the assessment of maximum earthquakes in the alps and adjacent areas

Strong tectonic earthquakes within the crust always occur on already existing faults, and they have the property of a shear rupture. Such earthquakes with surface-wave magnitudes M < 7 obviously have a geometric similarity. Because of this similarity and the validity of the Gutenberg and Richter's energy—magnitude relation, the expression M = 2 log₁₀ L + const., with L = focal length, is valid.The expression L_maxL^* for the maximum focal length, is also valid if L^* is the length of the rectilinear extent of the seismic line on which the maximum earthquake occurs. The bounds of L^* may be given by sharp bends and/or by traversing deep faults. Thus the maximum imaginable earthquake on a seismic line with the length L^* has the magnitude M_max = 2 log₁₀ L^* + const.For the investigated region — the Alps and adjacent areas — from the data of recent and historical strong earthquakes, it follows that M_max = 2 log₁₀ L^* + 1.7, if L^* is measured in kilometres. These limiting values lie in the centre-field of the magnitude range for maximum earthquakes, published by Shebalin in 1970. By the aid of this equation it is also possible to assess the upper limiting value of the accompanying maximum scale intensity.     

Recent plate tectonics of the Arctic Basin and of northeastern Asia

The recent tectonics of the Arctic Basin and northeastern Asia are considered as a result of interaction between three lithospheric plates: North-America, Eurasia and Spitsbergen. Seismic zones (coinciding in the Norway-Greenland basin with the Kolbeinsey, Mohns and Knipovich ridges, and in the Arctic Ocean with the Gakkel Ridge) clearly mark the boundaries between them. In southernmost Svalbard (Spitsbergen), the secondary seismic belt deviates from the major seismic zone. This belt continues into the seismic zone of the Franz Josef Land and then merges into the seismic zone of the Gakkel Ridge at 70°–90°E. The smaller Spitsbergen plate is located between the major seismic zone and its secondary branch.Within northeastern Asia, earthquake epicenters with magnitude over 4.5 are concentrated within a 300-km wide belt crossing the Eurasian continent over a distance of 3000 km from the Lena estuary to the Komandorskye Islands. A single seismic belt crosses the northern sections of the Verkhoyansky Ridge and runs along the Chersky Ridge to the Kolymo-Okhotsk Divide.To compute the poles of relative rotation of the Eurasian, North-American and Spitsbergen plates we use 23 new determinations of focal-mechanism solutions for earthquakes, and 38 azimuths of slip vectors obtained by matching of symmetric mountain pairs on both sides of the Knipovich and Gakkel ridges; we also use 14 azimuths of strike-slip faults within the Chersky Ridge determined by satellite images. The following parameters of plate displacement were obtained: Eurasia/North America: 62.2°N, 140.2°E (from the Knipovich Ridge section south of the triple junction); 61.9°N, 143.1°E (from fault strikes in the Chersky Ridge); 60.42°N, 141.56°C (from the Knipovich section and from fault strikes in the Chersky Ridge); 59.48°N, 140.83°E, α = 1.89 · 10⁻⁷ deg/year (from the Knipovich section, from fault strikes in the Chersky Ridge and from the Gakkel Ridge section east of the triple junction). The rate was calculated by fitting the 2′ magnetic lineations within the Gakkel Ridge).North-America/Spitsbergen: 70.96°N, 121.18°E, α = −2.7 · 10⁻⁷ deg/year from the Knipovich Ridge section north of the triple junction, from earthquakes in the Spitsbergen fracture zone and from the Gakkel Ridge section west of the triple junction). Eurasia/Spitsbergen: 70.7°N, 25.49°E, α = −0.99 · 10⁻⁷ deg/year (from closure of vector triangles).     

Focal mechanisms of three earthquakes in Finland and their relation to surface faults

Focal mechanisms for three recent earthquakes in Finland are determined using P-wave polarities together with SV/P and SH/P phase amplitude ratios. The events occurred on May 11, 2000 in Toivakka, Central Finland (M_L=2.4), on September 15, 2000 in Kuusamo, northeastern Finland (M_L=3.5), and on May 2, 2001 in Kolari, western Finnish Lapland (M_L=2.9).In order to obtain reliable estimates of the source parameters, one-dimensional crust and upper mantle velocity models are derived for the epicenter areas from deep-seismic sounding results. The starting models are modified by one-dimensional ray tracing using the earthquake observations. The events are relocated by employing P- and S-phase arrival times from the nearest seismic stations and the final velocity models. Synthetic waveforms, calculated with the reflectivity method, are used to further constrain and verify the source and structural parameters.The Toivakka earthquake indicates thrust- or reverse-faulting mechanism at a depth of 5 km. After comparison with aeromagnetic and topographic data we suggest the eastward dipping nodal plane (358°/42°) was the fault plane. The best-fitting fault plane solution of the Kolari earthquake suggests pure thrust-faulting at a depth of 5 km. The nodal plane striking 035°/30° correlates well with surface observations of the postglacial, possibly listric fault systems in the source area. The Kuusamo earthquake (focal depth 14 km) has a normal-faulting mechanism with the nodal planes trending 133°/47° or 284°/47°. Preference is given to the SW-dipping nodal plane, as it seems to coincide with topographic and magnetic lineament directions that have been active after the last ice age.The three earthquakes have occurred in old Precambrian faults and shear zones, which have been reactivated. The reactivated faults are favourably oriented in the local stress field.     

The rupture process during the 1999 Düzce, Turkey, earthquake from joint inversion of teleseismic and strong-motion data

Teleseismic and strong-motion data are inverted to determine the rupture process during the November 1999 Düzce earthquake in NW Turkey. The fault geometry, rise time and rupture velocity are determined from the aftershock distribution and preliminary inversions of the teleseismic data. Joint inversion of the teleseismic and strong-motion data is then carried out for the slip distribution. We obtain the strike 264°, dip 64°, rake −172°, seismic moment 5.0×10¹⁹ N m (M_w 7.1), and average stress drop 7 MPa. This earthquake was characterized by bilateral fault rupture and asymmetric slip distribution. Two asperities (areas of large slip) are identified, the eastern one being 1.5 times larger than the western one. The derived slip distribution is consistent with the aftershock distribution, surface rupture and damage. The point of rupture initiation in this Düzce earthquake coincided with the eastern tip of the aftershock distribution of the August 1999 Izmit earthquake.     

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

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

Statistical occurrence analysis and spatio-temporal distribution of earthquakes in the Apennines (Italy)

We present two examples of statistical analysis of seismicity conducted by integrating geological, geophysical and seismological data with the aim to characterize the active stress field and to define the spatio-temporal distribution of large earthquakes. Moreover, our data will help to improve the knowledge of the “seismogenic behavior” of the areas and to provide useful information for seismic hazard evaluation.The earthquakes are described by two non-parametric statistical procedures integrating also tectonic-physical parameters to study the spatio-temporal variability.The results show that the areas are characterized by: 1) a stress regime with mainly extensional kinematics; 2) tectonic structures mainly oriented with the active stress field (S_hmin = N44° ± 18° in the southern Apennines and S_hmin = N50° ± 17° in the central Apennines); 3) cluster distribution of seismicity and 4) a high probability of earthquake occurrence (M > 5.5) in the next 10 years.     

Shear wave splitting beneath South Kamchatka during the 3-year period associated with the 1997 Kronotsky Earthquake

Shear wave splitting measurements in South Kamchatka during the 3-year period (1996–1998) in which the Kronotsky Earthquake (M=7.7, December 5, 1997) occurred are used to determine anisotropic parameters of the subduction zone and shear wave splitting variations with time. The local small seismic events recorded at the Petropavlovskaya IRIS station (PET) were analyzed. The dominant azimuths of the fast shear wave polarizations for the 3-year period are defined within N95±15°E, which are consistent with the general Pacific Plate motion direction. Modeling of fast shear wave polarizations shows that HTI model with the symmetry axis oriented along N15°E±10° fit well the observed data for events the focal depths of which are less than 80 km. For the greater depths, the orthorhombic symmetry of medium is not excluded. The anisotropy coefficient increases generally with depth from 1–2% in the crust to 4–7.5% in the subducting plate. Variations in time delays show a general increase up to 10–15 ms/km during 1996–1997 before the large crustal earthquake series (M≈5.5–7) in the Avacha Bay and before the Kronotsky Earthquake. Analysis of fast S-wave azimuths of mantle events reveals a temporal cyclic variation. The most regular variations are observed for fast azimuths of deep events with a period of about 172 days over the 3-year period. The fast polarizations of crustal events behave comparatively stable. It is assumed that the major instabilities in stress state are localized in the descending slab and influenced the upper mantle and comparatively stable crust.     

Precursory accelerating seismic crustal deformation in the Northwestern Anatolia Fault Zone

We present the results of a systematic search for the identification of accelerating seismic crustal deformation in the broader northern Aegean area and in northwestern Turkey. We found that accelerating seismic deformation release, expressed by the generation of intermediate magnitude earthquakes, is currently observed in NW Turkey. On the basis of the critical earthquake model and by applying certain constraints which hold between the basic quantities involved in this phenomenon, it can be expected that this accelerating seismic activity may culminate in the generation of two strong earthquakes in this area during the next few years.The estimated epicenter coordinates of the larger of these probably impending earthquakes are 39.7°N–28.8°E, its magnitude is 7.0 and its occurrence time t_c=2003.5. The second strong event is expected to occur at t_c=2002.5 with a magnitude equal to 6.4 and epicenter coordinates 40.0°N–27.4°E. The uncertainties in the calculated focal parameters for these expected events are of the order of 100 km for the epicenter, ±0.5 for their magnitude and ±1.5 years for their occurrence time.     

Recent activity and seismotectonics of the Eastern Pyrenees

Results from a recent earthquake in the Eastern Pyrenees are presented and the seismotectonics of the region is analyzed from the presently available data. On 26 September 1984 an earthquake (M_L = 4.4) took place in the area of the historical destructive earthquake of 1428. Several portable stations installed in the epicentral area to record aftershocks permitted of defining a precise location at 42°19.2′N, 2°10.2′E and 5 km depth. A maximum felt intensity of V (MSK) is obtained from macroseismic data. The epicentral location lies within a block bounded by E-W-trending structures and the focal solution shows right-lateral shearing with a NW-SE pressure axis.The seismicity in the Eastern Pyrenees shows a complex pattern which can be associated with both E-W fractures and NE-SW fault systems. Focal solutions of another two recent earthquakes of M_L ~ 4, with differences in horizontal pressure axis, are also discussed.     

Complex geometry and segmentation of the surface rupture associated with the 14 November 2001 great Kunlun earthquake, northern Tibet, China

The 14 November 2001 Kunlun, China, earthquake with a moment magnitude (M_w) 7.8 occurred along the Kusai Lake–Kunlun Pass fault of the Kunlun fault system. We document the spatial distribution and geometry of surface rupture zone produced by this earthquake, based on high-resolution satellite (Landsat ETM, ASTER, SPOT and IKONOS) images combined with field measurements. Our results show that the surface rupture zone can be divided into five segments according to the geometry of surface rupture, including the Sun Lake, Buka Daban–Hongshui River, Kusai Lake, Hubei Peak and Kunlun Pass segments from west to east. These segments, each 55 to 130 km long, are separated by step-overs. The Sun Lake segment extends about 65 km with a strike of N45° 75°W (between 90°05′E 90°50′E) along the previously unrecognized West Sun Lake fault. A gap of about 30 km long exists between the Sun Lake and Buka Daban Peak where no obvious surface ruptures can be observed either from the satellite images or field observations. The Buka Daban–Hongshui River, Kusai Lake, Hubei Peak and Kunlun Pass segments run about 365 km striking N75° 85°W along the southern slope of the Kunlun Mountains (between 91°07′E 94°58′E). This segmentation of the surface rupture is well correlated with the pattern of slip distribution measured in the field. Detailed mapping suggest that these five first-order segments can be further separated into over 20 second-order segments with a length of 10–30 km, linked by smaller scale step-overs or bends.Our result also shows that the total coseismic surface rupture length produced by the 2001 Kunlun earthquake is about 430 km (excluding the 30-km-long gap), which is the longest coseismic surface rupture for an intracontinental earthquake ever recorded.Finally, we suggest a multiple bilateral rupture propagation model that shows the rupture process of the 2001 M_w 7.8 earthquake is complex. It consists of westward and eastward rupture propagations and interaction of these bilateral rupture processes.