Modeling of source parameters and moment tensors of local earthquakes occurring in the eastern Indian shield

Earthquake source parameters and crustal Q are being estimated simultaneously through the inversion of S-wave displacement spectra from three-component recordings of ten local cratonic intraplate earthquakes from 3-6 broadband stations in the eastern Indian shield, wherein, an iterative Levenberg-Marquardt inversion technique is used. The estimated seismic moment (M_o) and source radii (r) vary from 7.4 x 10¹² to 7.1 x 10¹⁴ N-m and 144.2 to 211.3 m, respectively, while estimated stress drops (Δσ) and multiplicative factor (E_mo) values range from 0.11 to 4.13 MPa and 1.33 to 2.16, respectively. The corner frequencies range from 6.23 to 8.62 Hz while moment magnitudes vary from 2.44 to 3.57. The radiated seismic energy and apparent stresses range from 8.3 x 10⁶ to 2.0 x 10¹⁰ Joules and 0.06 to 0.94 MPa, respectively, wherein the estimated corner frequencies and seismic moment satisfy the relation Mo ∞ f _c ^–(3+ε) for ε = 12.7. Thus, the source scaling of these events clearly deviates from the self-similarity i.e. f^–3. Estimated Zuniga parameters reveal that all selected events satisfy the partial stress drop model, which is in good agreement with the global observations. Our estimated crustal S-wave quality factors vary from 1091 to 4926 with an average of 3006, suggesting a less heterogeneous crustal structure underlying the study region.We also perform moment tensor inversion of five selected local events using ISOLA software, which reveals that the dominant deformation mode for the eastern Indian shield is left-lateral strike slip motion with minor normal dip-slip component on an almost vertical plane. This observation suggests that neotectonic vertical movements might have played a key role in generating these earthquakes. Our modeling also depicts that the seismically mildly active Singhbhum shear zone and Eastern Ghats mobile belt are characterized by the left-lateral strike motion while two events in the Chotanagpur half graben belt suggest a normal dip-slip motion along a south dipping plane. A north-south orientation of P-axis is found to be dominant in the area, which is consistent with the prevailing north–south compression over the Indian plate.     

Double-planed structure of the deep seismic zone in the northeastern Japan arc

A double-planed structure of deep seismic zone has been found over a wide area more then 300 km × 200 km in the Tohoku District, the northeastern part of Honshu, Japan. This prominent feature of the configuration of the deep seismic zone has been ascertained through a precise determination of the microearthquake hypocenters by using the data from the seismic network of Tohoku University. The two planes are nearly parallel to each other, the distance between the two planes being from 30 to 40 km.Composite focal mechanism solutions are derived from the superposition of the distribution of first motions of P waves, and the different fault types are obtained for the two groups of earthquakes; the earthquakes which occurred in the upper plane are characterized by reverse faulting, some of them by down-dip compressional stresses, and those in the lower plane by down-dip extensional stresses. The evidence obtained here provides valuable information for the definition of the type of mechanism producing the plate motion beneath the island arc.     

Depth variation of seismic source scaling relations: implications for earthquake hazard in southeastern Australia

Advances in earthquake data acquisition and processing techniques have allowed for improved quantification of source parameters for local Australian earthquakes. Until recently, only hypocentral locations and local magnitudes (M_L) had been determined routinely, with little attention given to the inversion of additional source parameters. The present study uses these new source data (e.g. seismic moment, stress drop, source dimensions) to further extend our understanding of seismicity and the continental stress regime of the Australian landmass and its peripheral regions.

Earthquake activity within Australia is typically low, and the proportion of small to large events (i.e. the b value) is also low. It is observed that average stress drops for southeastern Australian earthquakes appear to increase with seismic moment to relatively high levels, up to approximately 10 MPa for M_L 5.0 earthquakes. This is thought to be indicative of high compressive crustal stress, coupled with strong rocks and fault asperities. Furthermore, the data indicates that shallow focus earthquakes (shallower than 6 km) appear to produce lower than average stress drops than deeper earthquakes (between 6 and 20 km) with similar moment.

Recurrence estimates were obtained for a discrete seismogenic zone in southeastern Australia. Decreasing b values with increasing focal depth for this zone indicate that larger earthquakes (with high stress drops) tend to occur deeper in the crust. This may offer an explanation for the apparent increase of stress drop with hypocentral depth. Consequently, earthquake hazard estimates that assume a uniform Gutenburg–Richter distribution with depth (i.e. constant b value) may be too conservative and therefore slightly overestimate seismic hazard for surface sites in southeastern Australia.     

Seismic moment tensor for intermediate depth earthquakes at regional distances in Southern Spain

A method has been developed to obtain the seismic moment tensor components by linear inversion of P waves recorded at regional distance for intermediate depth earthquakes. The seismic moment tensor is separated into double couple (DC) and compensated linear vector dipole (CLVD) parts. The method has been applied to four earthquakes (64<h<95 km) which occurred in the Malaga region (southern Spain). Solutions for the 1987 event show a percentage of CLVD of 20% with a short source time function and DC part corresponding to vertical motion. For the 1989, 1990 and 1992 earthquakes, percentages of CLVD between 0% and 6% have been found. Comparison with the results obtained in a previous study [Buforn et al., J. Seismol. 1 (1997) 113] by modelling of P waves using a DC model, shows that the use of a more general representation of the source (seismic moment tensor) gives a fit of data for the 1987 event.     

Source parameters of some recent earthquakes in the Gulf of Aqaba, Egypt

Gulf of Aqaba is recognized as an active seismic zone where many destructive earthquakes have occurred. The estimation of source parameters and coda Q attenuation are the main target of this work. Fifty digital seismic events in eight short-period seismic stations with magnitude 2.5–5.2 are used. Most of these events occurred at hypocentral depths in the range of 7–20 km, indicating that the activity was restricted in the upper crust. Seismic moment, M _o, source radius, r, and stress drop, Δσ, are estimated from P- and S-wave spectra using the Brune’s seismic source model. The average seismic moment generated by the whole sequence of events was estimated to be 4.6E?+?22 dyne/cm. The earthquakes with higher stress drop occur at 10-km depth. The scaling relation between the seismic moment and the stress drop indicates a tendency of increasing seismic moment with stress drop. The seismic moment increases with increasing the source radius. Coda waves are sensitive to changes in the subsurface due to the wide scattering effects generating these waves. Single scattering model of local earthquakes is used to the coda Q calculation. The coda with lapse times 10, 20, and 30 s at six central frequencies 1.5, 3, 6, 12, 18, 24 Hz are calculated. The Q _c values are frequency dependent in the range 1–25 Hz, and are approximated by a least squares fit to the power law [ $ {Q_c}(f) = {Q_o}{(f/{f_o})^\eta } $ ]. The average of Q _c values increases from 53?±?10 at 1.5 Hz to 700?±?120 at 24 Hz. The average of Q _o values ranges from 13?±?1 at 1.5 Hz to 39?±?4 at 24 Hz. The frequency exponent parameter η ranges between 1.3?±?0.008 and 0.9?±?0.001.     

The devastating Muzaffarabad earthquake of 8 October 2005: New insights into Himalayan seismicity and tectonics

The recent earthquake of 8 October 2005 in the Muzaffarabad region in western Himalaya destroyed several parts of Pakistan and the north Indian state of Jammu and Kashmir. The earthquake of magnitude 7.6 claimed more than 80,000 lives, clearly exposing the poor standards of building construction — a major challenge facing the highly populated, earthquake prone, third world nations today. In this paper, we examine variations in the stress field, seismicity patterns, seismic source character, tectonic setting, plate motion velocities, GPS results, and the geodynamic factors relating to the geometry of the underlying subsurface structure and its role in generation of very large earthquakes. Focal mechanism solutions of the Muzaffarabad earthquake and its aftershocks are found to have steep dip angles comparable to the Indian intra-plate shield earthquakes rather than the typical Himalayan earthquakes that are characterized by shallow angle northward dips. A low p-value of 0.9 is obtained for this earthquake from the decay pattern of 110 aftershocks, which is comparable to that of the 1993 Latur earthquake in the Indian shield — the deadliest Stable Continental Region (SCR) earthquake till date. Inversion of focal mechanisms of the Harvard CMT catalogue indicates distinct stress patterns in the Muzaffarabad region, seemingly governed by an overturned Himalayan thrust belt configuration that envelops this region, adjoined by the Pamir and Hindukush regions. Recent developments in application of seismological tools like the receiver function technique have enabled accurate mapping of the dipping trends of the Moho and Lithosphere–Asthenosphere Boundary (LAB) of Indian lithosphere beneath southern Tibet. These have significantly improved our understanding of the collision process, the mechanism of Himalayan orogeny and uplift of the Tibetan plateau, besides providing vital constraints on the seismic hazard threat posed by the Himalaya. New ideas have also emerged through GPS, macroseismic investigations, paleoseismology and numerical modeling approaches. While many researchers suggest that the Himalayan front is already overdue for several 8.0 magnitude earthquakes, some opine that most of the front may not really be capable of sustaining the stress accumulation required for generation of great earthquakes. We propose that the occurrence of great earthquakes like those of 1897 in Shillong and 1950 in Assam have a strong correlation with their proximity to multiple plate junctions conducive for enormous stress build up, like the eastern Himalayan syntaxis comprising the junction of the India, Eurasia plates, and the Burma, Sunda micro-plates.     

Seismic moment release data in earthquake catalogue: Application of Hurst statistics in delineating temporal clustering and seismic vulnerability

Sequential cumulative moment release data of macroearthquakes (Mw≥4.3) of seventeen seismic zones (A to Q) belonging to NE-Himalaya, Burmese-Andaman arc and West- Sunda arc are analysed by Hurst analysis, a non-parametric statistical procedure to identify clustering of low and high values in a time series. The moment release in a zone occurs in alternate positive, negative and positive sloping segments forming a wave like pattern with intervening small horizontal segment. The negative sloping segments indicate decelerated moment release pattern or temporal slackening of elastic strain release with high b–value (>0.95). The horizontal segment indicates temporal clustering of moderate magnitude events/seismic moments with moderate b-values (0.8–0.95). The positive segment is characterised by accelerated moment release within a short span of time indicating temporal clustering of larger magnitude earthquakes/seismic moments and exhibit lowest b–value (<0.7). All zones attest moderate to high Hurst K values, range 0.7-0.86. The pattern in Hurst plots, specially a reversal of trend after prolong negative slope is used for earthquake prognostication in the seismic zones. Our analysis shows that most of the zones register a notable reversal of Hurst clustering trend after a prolonged negative slope which is accompanied by a major earthquake near its end. However, South Burma region (Zone-I) and Tripura fold belt and Bangladesh Plain (Zone-K) do not show any moderate or large shock around the end of the negative sloping trend in Hurst plot. Hence, these two zones can be considered more prone to produce moderate to larger earthquakes in future.     

Tectonic implications of the large shioya-oki earthquakes of 1938

The tectonic processes taking place along the southern part of the Japan trench are discussed on the basis of the focal mechanism of the 1938 Shioya-Oki event which consists of the five large earthquakes of M_s = 7.4, 7.7, 7.8, 7.7 and 7.1. Detailed analyses of seismic waves and tsunamis are made for each of these earthquakes, and the dislocation parameters are obtained. The total seismic moment amounts to 2.3 · 10²⁸ dyn.cm. The five earthquakes are grouped into either a low-angle thrust type or a nearly vertical normal-fault type. These mechanisms are common with other great earthquakes of the northwestern Pacific belt, and can be explained in terms of the interaction between the oceanic and continental plates. The vertical displacement inferred from the seismic results is in approximate agreement with the precise level data over the period from 1939 and 1897. This agreement suggests that the rate of the strain accumulation at the preseismic time is very small in the epicentral area. Repeated levelings at the postseismic time reveal a large-scale recovery of the coseismic subsidence. The postseismic deformation is one-third to one-half of the coseismic displacement. The time constant of the recovery is estimated to be 5 years or less. This type of deformation may be a manifestation of viscoelasticity of a weak zone underlying the continent. The amount of dislocation, together with the longterm seismicity, suggests a seismic slip rate of about 0.4 cm/year, which is one order of magnitude smaller than that for the adjacent regions. This suggests that a large part of the plate motion is taking place aseismically in this region. The tectonic process now taking place in the southern Japan trench can be considered to represent a stage just prior to a complete detachment of the sinking portion of the oceanic plate.     

Source parameters and ground motion of the Suez-Cairo shear zone earthquakes,Eastern Desert,Egypt

Three felt earthquakes with local magnitudes 4.0 (June 29th, 2000), 4.2 (July 07th, 2005) and 3.7 (October 30th, 2007) occurred to the southeast of Cairo along the Suez-Cairo shear zone. Being the most well recorded events by the Egyptian National Seismic Network (ENSN) in this area, they provide us an excellent opportunity to study the tectonics, the stress field, the source parameters, and the expected ground motion at Cairo. The focal mechanisms of the three events based on the first motion P-wave polarities indicate mainly normal faulting with a slight strike-slip component along nodal planes trending nearly EW and NW–SE. The three focal solutions imply a rejuvenation of the pre-existing EW and NW–SE faults due to a partly transfer of rifting deformation from the Red Sea-Gulf of Suez along these trends. Dynamic source parameters of these events are estimated from the P-wave spectra of the closest ENSN stations. The average values of the seismic moment, stress drop, rupture radius, and fault dislocation are estimated from six stations and interpreted in the context of the tectonic setting. These parameters in addition to the effects of site and path are used to synthesize the ground motion seismograms using stochastic modeling technique at the recorded stations and at some strategic sites.     

Mechanism of the reservoir impounding earthquakes at Hsinfengkiang and a preliminary endeavour to discuss their cause

After the filling up of the Hsinfengkiang Reservoir Kwangtung Province, seismicity was greatly increased. The majority of earthquakes occurred in the deep water gorge close to the dam, concentrated within a northwest belt. They are usually of shallow focal depths. A strong earthquake with magnitude 6.1 took place on March 19, 1962, about two and a half years since the impounding of the reservoir.

According to the results of analysis of data from geodetic leveling and the spectra of seismic waves, the fault parameters of the main shock were determined. The fault plane solutions of 150 small earthquakes, occurring within a period of 18 months before and after the main shock were determined from the amplitudes of the first motion of P wave. The directions of the earthquake generating stress of about 2000 small earthquakes were obtained by smoothing the first motion patterns. Displacement field and stress field in the rock bodies underneath the reservoir caused by the loading of the reservoir water were calculated. Variations of the velocity ratio of the P and S waves prior to the main shock and several strong aftershocks were analysed.

In consideration of the seismicity as well as the geological background, we endeavour to discuss the cause of reservoir impounding earthquakes at Hsinfengkiang. We have the opinion that the penetration of water along fissures becomes the most important cause of the main shock of March 19, 1962 at Hsinfengkiang.     

Source parameters of selected earthquakes on the central and western margin of Afar

Spectral analysis of 196 short-period and long-period vertical- and horizontal-component seismograms from ten earthquakes on the central and western margin of Afar is performed to determine source parameters and discuss their tectonic implications. For the earthquakes in the regions under study, the stress drop varies from 2 to 31 bar while the seismic moment varies from 2 × 10²⁴ to 154 × 10²⁴ dyn cm. In general, low stress-drop values are obtained indicating the presence of softer material (especially for central Afar) at a shallower depth. It is observed that there is an increase in stress drop with the increase in moment-magnitude which in turn is obtained from the calculated average seismic moment. Energy estimates show that the mode of energy release is different in the two regions indicating that different tectonic processes are involved in the two regions. The slip rate obtained for the Serdo area is of the order of 1.6 cm/yr and is in close agreement with the spreading rate already obtained for central Afar. Spreading rates obtained earlier and that of the present study show a low spreading rate for Afar and neighbouring regions as compared to those of the other regions of the world.     

Crustal deformation in northern Central America

Evaluation of the seismic moment tensor for earthquakes on plate boundary is a standard procedure to determine the relative velocity of plates, which controls the seismic deformation rate predicted from the slip on a single fault. The moment tensor is also decomposed into an isotropic and a deviatoric part to discover the relationship between the average strain rate and the relative velocity between two plates. We utilize this procedure to estimate the rates of deformation in northern Central America where plate boundaries are seismically well defined. Four different tectonic environments are considered for modelling of the plate motions. The deformation rates obtained here compare well with those predicted from the plate motions models and are in good agreement with actual observations. Deformation rates on faults are increasingly being used to estimate earthquake recurrence from information on fault slip rate and more on how we can incorporate our current understanding into seismic hazard analyses.     

The present-day active deformation in the central and northern parts of the Gulf of Suez area,Egypt, from earthquake focal mechanism data

The average seismic strain rate is estimated for the seismotectonic zone of the northern/central parts of the Gulf of Suez. The principal strain rate tensor and velocity tensor were derived from a combination of earthquake focal mechanisms data and seismic moment of small-sized earthquakes covering a time span of 13 years (1992–2004). A total of 17 focal mechanism solutions have been used in the calculation of the moment tensor summation. The local magnitudes (MLs) of these events range from 2.8 to 4.7. The analysis indicates that the dominant mode of deformation in the central and northern parts of the Gulf of Suez is extension at a rate of 0.008 mm/year in N28°E direction and a small crustal thinning of 0.0034 mm/year. This low level of strain means that this zone experienced a little seismic deformation. There is also a right lateral shear motion along the ESE–WNW direction. This strain pattern is consistent with the predominant NW–SE normal faulting and ESE–WNW dextral transtensive faults in this zone. Comparing the results obtained from both stress and strain tensors, we find that the orientations of the principal axes of both tensors have the same direction with a small difference between them. Both tensors show a predominantly extensional domain. The nearly good correspondence between principal stress and strain orientations in the area suggests that the tectonic strength is relatively uniform for this crustal volume.     

The large-scale distribution of great shallow intraplate earthquakes in mainland China and adjacent areas and seismic coupling along plate margins

In mainland China and its surroundings the large-scale distribution of great, shallow, intraplate earthquakes shows that there are four main areas of high intraplate seismicity which are (a) the Northern China Seismic Area in east China (30°–42°N); (b) the Southeast-coast Seismic Area in eastern China (19°–25°N); (c) the North-south trending Seismic Area in western China and its surroundings (Burma–China–Mongolia); (d) the Central Asian Seismic Area in west China and its surroundings (Pamirs–Tianshan Mts–Baikal). These four intraplate seismic areas are approximately perpendicular to those sections of the Eurasian plate boundary where the Eurasian plate has a strong seismic coupling with the North American–Pacific Ocean–Philippine Sea plates, and with the Indian plate. The large-scale uneven distribution of intraplate seismicity in China and its surroundings may be controlled by heterogeneity in the stress state on different sections of the plate boundary.     

Double-sided subduction with contrasting polarities beneath the Pamir-Hindu Kush: Evidence from focal mechanism solutions and stress field inversion

The Pamir-Hindu Kush region at the western end of the Himalayan-Tibet orogen is one of the most active regions on the globe with strong seismicity and deformation and provides a window to evaluate continental collision linked to two intra-continental subduction zones with different polarities. The seismicity and seismic tomography data show a steep northward subducting slab beneath the Hindu Kush and southward subducting slab under the Pamir. Here, we collect seismic catalogue with 3988 earthquake events to compute seismicity images and waveform data from 926 earthquake events to invert focal mechanism solutions and stress field with a view to characterize the subducting slabs under the Pamir-Hindu Kush region. Our results define two distinct seismic zones: a steep one beneath the Hindu Kush and a broad one beneath the Pamir. Deep and intermediate-depth earthquakes are mainly distributed in the Hindu Kush region which is controlled by thrust faulting, whereas the Pamir is dominated by strike-slip stress regime with shallow and intermediate-depth earthquakes. The area where the maximum principal stress axis is vertical in the southern Pamir corresponds to the location of a high-conductivity low-velocity region that contributes to the seismogenic processes in this region. We interpret the two distinct seismic zones to represent a double-sided subduction system where the Hindu Kush zone represents the northward subduction of the Indian plate, and the Pamir zone shows southward subduction of the Eurasian plate. A transition fault is inferred in the region between the Hindu Kush and the Pamir which regulates the opposing directions of motion of the Indian and Eurasian plates.     

Seismic activity associated with the subducting motion of the Philippine Sea plate beneath the Kanto district, Japan

The Pacific plate and the Philippine Sea plate overlap and subduct underneath the Kanto region, central Japan, causing complex seismic activities in the upper mantle. In this research, we used a map selection tool with a graphic display to create a data set for earthquakes caused by the subducting motion of the Philippine Sea plate that are easily determined. As a result, we determined that there are at least four earthquake groups present in the upper mantle above the Pacific plate. Major seismic activity (Group 1) has been observed throughout the Kanto region and is considered to originate in the uppermost part of mantle in the subducted Philippine Sea plate, judging from the formation of the focal region and comparison with the 3D structure of seismic velocity. The focal mechanism of these earthquakes is characterized by the down-dip compression. A second earthquake layer characterized by down-dip extension (Group 2), below the earthquakes in this group, is also noted. The focal region for those earthquakes is considered to be located at the lower part of the slab mantle, and the Pacific plate located directly below is considered to influence the activity. Earthquakes located at the shallowest part (Group 3) form a few clusters distributed directly above the Group 1 focal region. Judging from the characteristics of later phases in these earthquakes and comparing against the 3D structure of seismic velocity, the focal regions for the earthquakes are considered to be located near the upper surface of the slab. Another earthquake group (Group 4) originates further below Group 2; it is difficult to consider these earthquakes within a single slab. The seismic activities representing the upper area of the Philippine Sea plate are Group 3. This paper proposes a slab geometry model that is substantially different from conventional models by strictly differentiating the groups.     

Tectonic stress,seismicity, and seismic hazard in the southeastern Carpathians

Intermediate-depth earthquakes in the Vrancea region occur in response to stress generation due to descending lithosphere beneath the southeastern Carpathians. In this article, tectonic stress and seismicity are analyzed in the region on the basis of a vast body of observations. We show a correlation between the location of intermediate-depth earthquakes and the predicted localization of maximum shear stress in the lithosphere. A probabilistic seismic hazard assessment (PSHA) for the region is presented in terms of various ground motion parameters on the utilization of Fourier amplitude spectra used in engineering practice and risk assessment (peak ground acceleration, response spectra amplitude, and seismic intensity). We review the PSHA carried out in the region, and present new PSHA results for the eastern and southern parts of Romania. Our seismic hazard assessment is based on the information about the features of earthquake ground motion excitation, seismic wave propagation (attenuation), and site effect in the region. Spectral models and characteristics of site-response on earthquake ground motions are obtained from the regional ground motion data including several hundred records of small and large earthquakes. Results of the probabilistic seismic hazard assessment are consistent with the features of observed earthquake effects in the southeastern Carpathians and show that geological factors play an important part in the distribution of the earthquake ground motion parameters.     

Seismic source parameters of large earthquake in Vrancea

The seismic source parameters for five large Vrancea earthquakes are calculated : corner frequency, focal radius, seismic moment and stress drop.The multiplicity character of some events is put in evidence, using spectral and time domain analysis.     

近500年中国及邻区强震空间格局及迁移模式初步分析

介绍了中国及邻区地震目录数据库的开发状况，从历史强震、近代强震及强震震源深度等几个方面，结合我国及邻区板块构造以及应力状态、壳幔结构的东西部差异对我国强震的空间格局进行了较为系统的分析．公元1900年以前我国华北地区的强震呈现出大梯形格局，而公元1900年之后我国强震则以西部的巨大扇形为主要特征，其中震源深度最深的区域分别位于大扇形的2个顶点处，即兴都库什-帕米尔以及缅印交界地区．就我国强震所特有的空间分布格局对我国强震的线性以及区域性迁移模式进行了总结和分析，其中线性迁移模式又可细分为前进跳跃式和钟摆式．对文中涉及的强震迁移机制分别利用断层破裂、弹簧-滑块、壳幔的结构性差异以及它们之间的相对运动等模型进行了初步的解释．