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.     

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.     

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.     

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.     

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 MW = 6.3 2003 Lefkada earthquake (Greece) and induced stress transfer changes

A large earthquake of magnitude M_W = 6.3 occurred on 14 August 2003 NW of the Lefkada Island, which is situated at the Ionian Sea (western Greece). The source parameters of this event are determined using body-wave modeling. The focal depth was found equal to 9 km, the constrained focal mechanism revealed dextral strike–slip motion (φ = 15°, Δ = 80° and λ = 170°), the duration of the source time function was 8 s and the seismic moment 2.9 × 10²⁵ dyn cm. The earthquake occurred close to the northern end of the Kefallinia transform fault, where the 1994 moderate event and its aftershock sequence were also located. The epicentral distribution of the 2003 aftershock sequence revealed the existence of two clusters. The first one is located close to the epicentral area of the mainshock, while the second southern, close to the northwestern coast of the Kefallinia Island. A gap of seismicity is observed between the two clusters. The length of the activated zone is approximately 60 km. The analysis of data revealed that the northern cluster is directly related to the mainshock, while the southern one was triggered by stress transfer caused by the main event.     

Source parameters of the Burma-India border earthquake of July 29, 1970, from body waves

The source parameters are determined for the Burma-India border earthquake of July 29, 1970, from body-wave spectra. We obtain seismic moment [ , ] · 10²⁶ dyne cm, source dimension [ ] km, radiated energy [ , −E_R (S) = 1.35] · 10²⁰ ergs and the stress drop = 11 bars.     

Spectra of the earthquake sequence February-March, 1981, in south-central Sweden

On February 13, 1981 a relatively strong earthquake occurred in the Lake Vänern region in south-central Sweden. The shock had a magnitude ofM_L = 3.3 and was followed within three weeks by three aftershocks, with magnitudes 0.5 ≤ M_L ≤ 1.0. The focal mechanism solution of the main shock indicates reverse faulting with a strike in the N-S or NE-SW direction and a nearly horizontal compressional stress. The aftershocks were too small to yield data for a full mechanism solution, but first motions of P-waves, recorded at two stations, are consistent for the aftershocks. Dynamic source parameters, derived from Pg- and Sg-wave spectra, show similar stress drops for the main shock (2 bar) and the aftershocks (1 bar), while the differences in seismic moment (1.5·10²⁰ resp. 4·10¹⁸dyne cm), fault length (0.7 resp. 0.2 km) and relative displacement (0.15 resp. 0.03 cm) are significant.     

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.     

The Mw = 6.0, 7 September 1999 Athens Earthquake

On 7 September 1999 at 11:56 GMT a destructive earthquake (M_w = 6.0) occurred close to Athens (Greece). The rupture process is examined using data from the Cornet local permanent network, as well as teleseismic recordings. Data recorded by a temporary seismological network were analyzed to study the aftershock sequence. The mainshock was relocated at 38.105°N, 23.565°E, about 20 km northwest of Athens. Four foreshocks were also relocated close to the mainshock. The modeling of teleseismic P and SH waves provides a well-constrained focal mechanism of the mainshock (strike = 105°, dip = 55° and rake = -80°) at a depth of 8 km and a seismic moment M₀ = 1.010²⁵ dyn·cm. The obtained fault plane solution represents normal faulting indicating an almost north-south extension. More than 3500 aftershocks were located, 1813 of which present RMS < 0.1 s and ERH, ERZ < 1.0 km. Two main clusters were distinguished, while the depth distribution is concentrated between 2 and 11 km. Over 1000 fault plane solutions of aftershocks were constrained, the majority of which also correspond to N–S extension. No surface breaks were observed but the fault plane solution of the mainshock is in agreement with the tectonics of the area and with the focal mechanisms obtained by aftershocks. The hypocenter of the mainshock is located on the deep western edge of the fault plane. The relocated epicenter coincides with the fringe that represents the highest deformation observed on the differential interferometric image. The calculated source duration is 5 sec, while the estimated dimensions of the fault are 15 km length and 10 km width. The source process is characterized by unilateral eastward rupture propagation, towards the city of Athens. An evident stop phase observed in the recordings of the Cornet local stations is interpreted as a barrier caused by the Aegaleo Mountain.     

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.     

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.     

Rupture mechanism of the 1993 Killari earthquake, India: constraints from aftershocks and static stress change

The Killari earthquake of September 29, 1993 (Mw=6.2) in peninsular India triggered several aftershocks that were recorded by a network of 21 stations. We computed the change in regional static stress caused by coseismic slip on the earthquake rupture and correlated it with the aftershocks with a view to constrain some of the rupture parameters of this earthquake. We evaluated the six available estimates of fault plane solutions for this earthquake and concluded that reverse slip on a 42° dipping, N112° trending fault, which extends up to the surface from a depth of 7 km, produces maximum correlation between the increased static stress and aftershock distribution. Our analysis suggests that the majority of coseismic slip occurred on the part of the rupture that lies in the depth range of 3–6.5 km.     

Seismotectonics of the Nepal Himalaya from a local seismic network

The National Seismological Network of Nepal consists of 17 short period seismic stations operated since 1994. It provides an exceptional view of the microseismic activity over nearly one third of the Himalayan arc, including the only segment, between longitudes 78°E and 85°E, that has not produced any M>8 earthquakes over the last century. It shows a belt of seismicity that follows approximately the front of the Higher Himalaya with most of the seismic moment being released at depths between 10 and 20 km. This belt of seismicity is interpreted to reflect interseismic stress accumulation in the upper crust associated with creep in the lower crust beneath the Higher Himalaya. The seismic activity is more intense around 82°E in Far-Western Nepal and around 87°E in Eastern Nepal. Western Nepal, between 82.5 and 85°E, is characterized by a particularly low level of seismic activity. We propose that these lateral variations are related to segmentation of the Main Himalayan Thrust Fault. The major junctions between the different segments would thus lie at about 87°E and 82°E with possibly an intermediate one at about 85°E. These junctions seem to coincide with some of the active normal faults in Southern Tibet. Lateral variation of seismic activity is also found to correlate with lateral variations of geological structures suggesting that segmentation is a long-lived feature. We infer four 250–400 km long segments that could produce earthquakes comparable to the M=8.4 Bihar–Nepal earthquake that struck eastern Nepal in 1934. Assuming the model of the characteristic earthquake, the recurrence interval between two such earthquakes on a given segment is between 130 and 260 years.     

Source parameters of the 2001 Mw 7.7 Bhuj earthquake,Gujarat, India,aftershock sequence

We present the estimated source parameters from SH-wave spectral modeling of selected 463 aftershocks (2002–06) of the 26 January 2001 Bhuj earthquake, the well-recorded largest continental intraplate earthquake. The estimated seismic moment (M_o), corner frequency (f_c), source radius (r) and stress drop (Δσ) for aftershocks of moment magnitude 1.7 to 5.6 range from 3.55×10¹¹ to 2.84×10¹⁷ N-m, 1.3 to 11.83 Hz, 107 to 1515 m and 0.13 to 26.7 MPa, respectively, while the errors in f_c and Δσ are found to be 1.1 Hz and 1.1 MPa, respectively. We also notice that the near surface attenuation factor (k) values vary from 0.02 to 0.03. Our estimates reveal that the stress drop values show more scatter (Mo^{0.5 to 1} is proportional to Δσ) toward the larger M_o values (≥10^14.5 N-m), while they show a more systematic nature (Mo³ is proportional to Δσ) for smaller M_o values (<10^14.5 N-m), which can be explained as a consequence of a nearly constant rupture radius for smaller aftershocks in the region. The large stress drops (= 10 MPa) associated with events on the north Wagad fault (at 15–30 km depth) and Gedi fault (at 3–15 km depth) can be attributed to the large stress developed at hypocentral depths as a result of high fluid pressure and the presence of mafic intrusive bodies beneath these two fault zones.     

The Gibraltar Arc seismogenic zone (part 2): Constraints on a shallow east dipping fault plane source for the 1755 Lisbon earthquake provided by tsunami modeling and seismic intensity

The Great Lisbon earthquake has the largest documented felt area of any shallow earthquake and an estimated magnitude of 8.5–9.0. The associated tsunami ravaged the coast of SW Portugal and the Gulf of Cadiz, with run-up heights reported to have reached 5–15 m. While several source regions offshore SW Portugal have been proposed (e.g.— Gorringe Bank, Marquis de Pombal fault), no single source appears to be able to account for the great seismic moment as well as all the historical tsunami amplitude and travel time observations. A shallow east dipping fault plane beneath the Gulf of Cadiz associated with active subduction beneath Gibraltar, represents a candidate source for the Lisbon earthquake of 1755.Here we consider the fault parameters implied by this hypothesis, with respect to total slip, seismic moment, and recurrence interval to test the viability of this source. The geometry of the seismogenic zone is obtained from deep crustal studies and can be represented by an east dipping fault plane with mean dimensions of 180 km (N–S) × 210 km (E–W). For 10 m of co-seismic slip an Mw 8.64 event results and for 20 m of slip an Mw 8.8 earthquake is generated. Thus, for convergence rates of about 1 cm/yr, an event of this magnitude could occur every 1000–2000 years. Available kinematic and sedimentological data are in general agreement with such a recurrence interval. Tsunami wave form modeling indicates a subduction source in the Gulf of Cadiz can partly satisfy the historical observations with respect to wave amplitudes and arrival times, though discrepancies remain for some stations. A macroseismic analysis is performed using site effect functions calculated from isoseismals observed during instrumentally recorded strong earthquakes in the region (M7.9 1969 and M6.8 1964). The resulting synthetic isoseismals for the 1755 event suggest a subduction source, possibly in combination with an additional source at the NW corner of the Gulf of Cadiz can satisfactorily explain the historically observed seismic intensities. Further studies are needed to sample the turbidites in the adjacent abyssal plains to better document the source region and more precisely calibrate the chronology of great earthquakes in this region.     

Kamchatka earthquake of 17 August 1983: A large intermediate event at the junction of the Aleutian and Kuril-Kamchatka arcs

This paper gives macroseismic and instrumental data on the 17 August 1983 Kamchatka earthquake which occurred in the center of the Kamchatka Gulf bend at a depth of 98 km with an epicentral intensity of VI–VII (MSK-64 scale), energy class 15, and magnitude MLH = 6.9. The focal mechanism represents a thrust along the inclined surface across the strike of the Kamchatka Gulf coastal line. The Primary-wave seismic moment M ₀ is 6.3 × 10¹⁹ Nm, the Rayleigh wave M ₀ is 1.6 × 10¹⁹ Nm, and the stress drop is 2.5 MPa. Copies of displacement and acceleration records are presented and the temporal and spatial distribution of the aftershocks is analyzed.     

The earthquake in the Swabian Jura of 16 November 1911 and present concepts of seismotectonics

Beginning with the Swabian Jura earthquake in 1911 the seismic activity in Central Europe is concentrated to this area. A comparison with other events of the same epicentral region shows that the largest earthquake in Germany has the character of a left-lateral horizontal strike slip striking N — NNE. The focal parameters can be assumed within the following intervalls: Seismic moment M_o = 1…8·10¹⁷ Nm; focal area F_o = 18…53 km²; average dislocation d_o = 20…53 cm and stress drop Δp_o = 13…19 bar.     

Seismic hazard map of Coimbatore using subsurface fault rupture

This study presents the future seismic hazard map of Coimbatore city, India, by considering rupture phenomenon. Seismotectonic map for Coimbatore has been generated using past earthquakes and seismic sources within 300 km radius around the city. The region experienced a largest earthquake of moment magnitude 6.3 in 1900. Available earthquakes are divided into two categories: one includes events having moment magnitude of 5.0 and above, i.e., damaging earthquakes in the region and the other includes the remaining, i.e., minor earthquakes. Subsurface rupture character of the region has been established by considering the damaging earthquakes and total length of seismic source. Magnitudes of each source are estimated by assuming the subsurface rupture length in terms of percentage of total length of sources and matched with reported earthquake. Estimated magnitudes match well with the reported earthquakes for a RLD of 5.2% of the total length of source. Zone of influence circles is also marked in the seismotectonic map by considering subsurface rupture length of fault associated with these earthquakes. As earthquakes relive strain energy that builds up on faults, it is assumed that all the earthquakes close to damaging earthquake have released the entire strain energy and it would take some time for the rebuilding of strain energy to cause a similar earthquake in the same location/fault. Area free from influence circles has potential for future earthquake, if there is seismogenic source and minor earthquake in the last 20 years. Based on this rupture phenomenon, eight probable locations have been identified and these locations might have the potential for the future earthquakes. Characteristic earthquake moment magnitude (M _w ) of 6.4 is estimated for the seismic study area considering seismic sources close to probable zones and 15% increased regional rupture character. The city is divided into several grid points at spacing of 0.01° and the peak ground acceleration (PGA) due to each probable earthquake is calculated at every grid point in city by using the regional attenuation model. The maximum of all these eight PGAs is taken for each grid point and the final PGA map is arrived. This map is compared to the PGA map developed based on the conventional deterministic seismic hazard analysis (DSHA) approach. The probable future rupture earthquakes gave less PGA than that of DSHA approach. The occurrence of any earthquake may be expected in near future in these eight zones, as these eight places have been experiencing minor earthquakes and are located in well-defined seismogenic sources.     

Source parameters of earthquakes in the reservoir-triggered seismic (RTS) zone of Koyna–Warna, Western India

New empirical relations are derived for source parameters of the Koyna–Warna reservoir-triggered seismic zone in Western India using spectral analysis of 38 local earthquakes in the magnitude range M _L 3.5–5.2. The data come from a seismic network operated by the CSIR-National Geophysical Research Institute, India, during March 2005 to April 2012 in this region. The source parameters viz. seismic moment, source radius, corner frequency and stress drop for the various events lie in the range of 10¹³–10¹⁶ Nm, 0.1–0.4 km, 2.9–9.4 Hz and 3–26 MPa, respectively. Linear relationships are obtained among the seismic moment (M ₀), local magnitude (M _L), moment magnitude (M _w), corner frequency (fc) and stress drop (?σ). The stress drops in the Koyna–Warna region are found to increase with magnitude as well as focal depths of earthquakes. Interestingly, accurate depths derived from moment tensor inversion of earthquake waveforms show a strong correlation with the stress drops, seemingly characteristic of the Koyna–Warna region.