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.     

The present status of reservoir induced seismicity investigations with special emphasis on Koyna earthquakes

The status of Reservoir Induced Seismicity (RIS) has been reviewed periodically (Rothé, 1968, 1973; Gupta and Rastogi, 1976; Simpson, 1976; Packer et al., 1979). In the present paper, the significant work carried out during the last three years on RIS is reviewed.An earthquake of magnitude occurred on November 14, 1981 in the vicinity of Aswan Lake, Egypt, 17 years after the filling started in 1964. This event occurred 4 days after the seasonal maximum in the reservoir water level and was followed by a long sequence of aftershocks. Another event of magnitude occurred in the vicinity of Aswan Lake on August 20, 1982. Results of preliminary investigations indicate that this seismic activity is reservoir induced. Recent analyses of induced seismic events at Nurek Reservoir U.S.S.R., show that the second stage of filling during August to December 1976, increasing the maximum depth from 120 m to 200 m, was accompanied by an intense burst of shallow seismic activity. An outward migration from the centre of the reservoir, possibly associated with diffusion of pore pressure, is revealed by the temporal distribution of earthquake foci. A variety of investigations including the in situ measurement of tectonic stress, pore pressure, permeability, distribution of faults, etc., in addition to monitoring seismicity, have been undertaken in the vicinity of the Monticello Reservoir, South Carolina. The largest reservoir induced earthquake is predicted not to exceed magnitude 5.The Koyna Reservoir, India, continues to be the most outstanding example of RIS. Three earthquakes of magnitude 5 occurred in September 1980. Earthquakes of magnitude 4 occur frequently in the vicinity of Koyna, the latest being on February 5, 1983. Events that occurred during the period 1967–1973 have been relocated using better procedures and are found to be much shallower and the epicentres less diffused. Location of 12 earthquakes of M_s 4.0, their foreshocks and aftershocks, that occurred during 1973–1976, composite focal mechanism solutions and related studies are consistent with the delineation of a N-S trending fault through the reservoir area. In a couple of interesting studies it has been demonstrated that earthquakes of magnitude 5.0 in the Koyna region are usually preceded by several magnitude 4 earthquakes in the preceding fortnight. Also, a rate of loading of Koyna reservoir of at least 40 ft/week appears to be a necessary, although not sufficient, condition for the occurrence of magnitude 5 earthquakes. Smooth filling/emptying appears to be the key to reduce the hazard of RIS.A map and a table of the reported cases of reservoir induced changes in seismicity through 1982 have been compiled.     

Probabilities for the occurrences of medium to large earthquakes in northeast India and adjoining region

The return periods and occurrence probabilities related to medium and large earthquakes (M _w 4.0–7.0) in four seismic zones in northeast India and adjoining region (20°–32°N and 87°–100°E) have been estimated with the help of well-known extreme value theory using three methods given by Gumbel (1958), Knopoff and Kagan (1977) and Bury (1999). In the present analysis, the return periods, the most probable maximum magnitude in a specified time period and probabilities of occurrences of earthquakes of magnitude M ≥ 4.0 have been computed using a homogeneous and complete earthquake catalogue prepared for the period between 1897 and 2007. The analysis indicates that the most probable largest annual earthquakes are close to 4.6, 5.1, 5.2, 5.5 and 5.8 in the four seismic zones, namely, the Shillong Plateau Zone, the Eastern Syntaxis Zone, the Himalayan Thrusts Zone, the Arakan-Yoma subduction zone and the whole region, respectively. The most probable largest earthquakes that may occur within different time periods have been also estimated and reported. The study reveals that the estimated mean return periods for the earthquake of magnitude M _w 6.5 are about 6–7 years, 9–10 years, 59–78 years, 72–115 years and 88–127 years in the whole region, the Arakan-Yoma subduction zone, the Himalayan Thrusts Zone, the Shillong Plateau Zone and the Eastern Syntaxis Zone, respectively. The study indicates that Arakan-Yoma subduction zone has the lowest mean return periods and high occurrence probability for the same earthquake magnitude in comparison to the other zones. The differences in the hazard parameters from zone to zone reveal the high crustal heterogeneity and seismotectonics complexity in northeast India and adjoining regions.     

Magnitude conversion problem for the Turkish earthquake data

Earthquake catalogues which form the main input in seismic hazard analysis generally report earthquake magnitudes in different scales. Magnitudes reported in different scales have to be converted to a common scale while compiling a seismic data base to be utilized in seismic hazard analysis. This study aims at developing empirical relationships to convert earthquake magnitudes reported in different scales, namely, surface wave magnitude, M _S, local magnitude, M _L, body wave magnitude, m _b and duration magnitude, M _d, to the moment magnitude (M _w). For this purpose, an earthquake data catalogue is compiled from domestic and international data bases for the earthquakes occurred in Turkey. The earthquake reporting differences of various data sources are assessed. Conversion relationships are established between the same earthquake magnitude scale of different data sources and different earthquake magnitude scales. Appropriate statistical methods are employed iteratively, considering the random errors both in the independent and dependent variables. The results are found to be sensitive to the choice of the analysis methods.     

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.     

Association of subsurface radon changes in Alaska and the northeastern United States with earthquakes

Changes of radon concentration in the ground have been reported to correlate with seismic activity in many places over widely varying distances from the earthquake epicenters, sometimes over 100's and occasionally over 1000's of km depending on the magnitude M of the earthquakes. For proper use of such correlations in forecasting the location and size of impending earthquakes, it would be useful to establish a relation between M and the maximum distance x_M at which radon can be significantly altered by preseismic effects.We have monitored radon in the ground continuously at Blue Mountain Lake, NY starting Dec. 3. 1975 and at three locations in Alaska—Icy Bay, Yakataga, and Sand Point (starting Nov. 12, 1979; March 5, 1980; and June 29, 1980 respectively). Sample correlations from visual examination of the recent radon record and tentative use of a proposed scaling relation show a number of possible earthquake-related signals at these locations. The data are compatible with scaling relations that were derived from two separate models of premonitory elastic strains. In this work x_M = 10^0.48M, where x_M is in units of km and M ≥ 3. Since tilt and strain signals correlate similarly with magnitude and distance, it is likely that most earthquake-related radon signals are mechanically induced.     

Seismicity of Gujarat

Paper describes tectonics, earthquake monitoring, past and present seismicity, catalogue of earthquakes and estimated return periods of large earthquakes in Gujarat state, western India. The Gujarat region has three failed Mesozoic rifts of Kachchh, Cambay, and Narmada, with several active faults. Kachchh district of Gujarat is the only region outside Himalaya-Andaman belt that has high seismic hazard of magnitude 8 corresponding to zone V in the seismic zoning map of India. The other parts of Gujarat have seismic hazard of magnitude 6 or less. Kachchh region is considered seismically one of the most active intraplate regions of the World. It is known to have low seismicity but high hazard in view of occurrence of fewer smaller earthquakes of M????6 in a region having three devastating earthquakes that occurred during 1819 (M _w7.8), 1956 (M _w6.0) and 2001 (M _w7.7). The second in order of seismic status is Narmada rift zone that experienced a severely damaging 1970 Bharuch earthquake of M5.4 at its western end and M????6 earthquakes further east in 1927 (Son earthquake), 1938 (Satpura earthquake) and 1997 (Jabalpur earthquake). The Saurashtra Peninsula south of Kachchh has experienced seismicity of magnitude less than 6.     

Continued seismic activity at the Koyna reservoir site, India

Following the impounding of the Shivaji Sager Lake in 1962, tremors became prevalent in the Koyna region, considered previously to be aseismic. During ensuing years the tremor frequency appears to have been dependent on the rate of increase of water level, maximum water level reached, and the period for which high levels were retained. This culminated in a burst of seismic activity from September 1967 to January 1968 following the record water levels in the reservoir and included the earthquake of September 13, 1967 with magnitude 5.5 and the damaging December 10, 1967 earthquake of magnitude 6.0. During the next five years water levels were kept low and no significant earthquakes occurred subsequent to the October 29, 1968 earthquake of magnitude 5.

The reservoir was filled to maximum capacity during September 1973 and this was followed by a conspicuous increase in seismic activity which included an earthquake of magnitude 5.1 on October 17, 1973. However, seismic activity during 1973 was much less severe than that of 1967. This relative decrease in seismicity may indicate that (a) the “threshhold level” for relatively large magnitude earthquakes had increased; (b) a major portion of the accumulated strains had been released; and/or (c) the importance of the longer period of high loading in 1967. Similar observations have been made at other seismically active reservoir sites.     

原位地应力测量与实时监测在地震预报研究中的作用和意义

当2013年4月20日四川芦山Ms 7.0级大地震发生时,中国科学家已不再像2008年汶川Ms 8.0大地震发生时那样茫然和不知所措。其根本原因在于,2008年汶川大地震发生后,龙门山地区开展了大量的科学研究工作,已超前初步认知龙门山断裂带西南端具有潜在地震危险性,特别是原位地应力测量和实时监测已发现绝对地应力大小高值异常和相对地应力大小临震异常变化。论文简要介绍了地震预报国际主流观点与认识,梳理了地应力在地震预报研究中的作用和认识,探讨了2004 Parkfield earthquake钻孔应变监测结果给予的启示,详细介绍了原位地应力测量与实时监测在地震预报研究中应用的实践与探索。实践证明:地震预报是值得探索和研究的,原位地应力测量与实时监测是地震预报的有效方法之一。      

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.     

Process Analysis of In-situ Strain during the Ms8.0 Wenchuan Earthquake—Data from the Stress Monitoring Station at Shandan

Abstract: There were huge life and property losses during the Ms8.0 Wenchuan earthquake on May 12, 2008. Strain fluctuation curves were completely recorded at stress observatory stations in the Qinghai-Tibet plateau and its surroundings in the process of the earthquake. This paper introduces the geological background of the Wenchuan earthquake and the profile of in-situ stress monitoring stations. In particular, data of 174 earthquakes (Ms4.0-Ms8.5) were processed and analyzed with various methods, which were recorded at the Shandan station from August 2007 to December 2008. The results were compared with other seismic data, and further analyses were done for the recoded strain seismic waves, co-seismic strain stepovers, pre-earthquake strain valleys, Earth’s free oscillations before and after the earthquake and their physical implications. During the Wenchuan earthquake, the strainmeter recorded a huge extensional strain of 70 seconds, which shows that the Wenchuan earthquake is a rupture process predominated by thrusting. Significant precursory strain anomalies were detected 48 hours, 30 hours, 8 hours and 37 minutes before the earthquake. The anomalies are very high and their forms are very similar to that of the main shock. Similar anomalies can also be found in strain curves of other shocks greater than Ms7.0, indicating that such anomalies are prevalent before a great earthquake. In this paper, it is shown that medium aftershocks (Ms5.5-6.0) can also cause Earth’s free oscillations. Study of free oscillations is of great significance to understand the internal structure of the Earth and focal mechanisms of earthquakes and to recognize slow shocks, thus providing a scientific basis for the prevention and treatment of geological disasters and the prediction of future earthquakes.     

Late-Quaternary Slip Rate and Seismic Activity of the Xianshuihe Fault Zone in Southwest China

The Xianshuihe fault zone is a seismo-genetic fault zone of left-lateral slip in Southwest China. Since 1725, a total of 59 Ms ≥ 5.0 earthquakes have occurred along this fault zone, including 18 Ms 6.0–6.9 and eight Ms ≥ 7.0 earthquakes. The seismic risk of the Xianshuihe fault zone is a large and realistic threat to the western Sichuan economic corridor. Based on previous studies, we carried out field geological survey and remote sensing interpretation in the fault zone. In addition, geophysical surveys, trenching and age-dating were conducted in the key parts to better understand the geometry, spatial distribution and activity of the fault zone. We infer to divide the fault zone into two parts: the northwest part and the southeast part, with total eight segments. Their Late Quaternary slip rates vary in a range of 11.5 mm/a –(3±1) mm/a. The seismic activities of the Xianshuihe fault zone are frequent and strong, periodical, and reoccurred. Combining the spatial and temporal distribution of the historical earthquakes, the seismic hazard of the Xianshuihe fault zone has been predicted by using the relationship between magnitude and frequency of earthquakes caused by different fault segments. The prediction results show that the segment between Daofu and Qianning has a possibility of Ms ≥ 7.0 earthquakes, while the segment between Shimian and Luding is likely to have earthquakes of about Ms 7.0. It is suggested to establish a GPS or In SAR-based real-time monitoring network of surface displacement to cover the Xianshuihe fault zone, and an early warning system of earthquakes and post seismic geohazards to cover the major residential areas.     

Characteristics of ring seismicity in different depth ranges before large and great earthquakes in the Sumatra region

Characteristics of the seismicity in depth ranges 0–33 and 34–70 km before ten large and great (M _w = 7.0−9.0) earthquakes of 2000–2008 in the Sumatra region are studied, as are those in the seismic gap zones where no large earthquakes have occurred since at least 1935. Ring seismicity structures are revealed in both depth ranges. It is shown that the epicenters of the main seismic events lie, as a rule, close to regions of overlap or in close proximity to “shallow” and “deep” rings. Correlation dependences of ring sizes and threshold earthquakes magnitudes on energy of the main seismic event in the ring seismicity regions are obtained. Identification of ring structures in the seismic gap zones (in the regions of Central and South Sumatra) suggests active processes of large earthquake preparation proceed in the region. The probable magnitudes of imminent seismic events are estimated from the data on the seismicity ring sizes.     

Tsunami magnitudes in Taiwan, the Philippines and Indonesia

The regional characteristics of tsunami magnitudes in the SE Asia region are discussed in relation to earthquake magnitudes during the period from 1960 to 1994. Tsunami magnitudes on the Imamura-Iida scale are investigated by the author's method (Hatori 1979, 1986) using the data of inundation heights near the source area and tide-gauge records observed in Japan. The magnitude values of the Taiwan tsunamis showed relatively to be small. On the contrary, the magnitudes of tsunamis in the vicinities of the Philippines and Indonesia exceed more than 1–2 grade (tsunami heights: 2–5 times) compared to earthquakes with similar size on the circum-Pacific zone. The relation between tsunami magnitude, m, and earthquake magnitude, M _s, is expressed as m = 2.66 M _s– 17.5 for these regions. For example, the magnitudes for the 1976 Mindanao tsunami (M _s= 7.8, 3702 deaths) and the 1992 Flores tsunami (M _s= 7.5, 1713 deaths) were determined to be m = 3 and m = 2.5, respectively. The focal depth of tsunamigenic earthquakes is shallower thand< 36 km, and the detectively of tsunamis is small for deep earthquakes being d > 40 km. For future tsunamis, it is indispensable to take precautions against shallow earthquakes having the magnitudes M _s> 6.5.     

Seismic hazard across Bulgaria and neighbouring areas: regional and site-specific maximum credible magnitudes and earthquake perceptibility

A probabilistic seismic hazard assessment is developed here using maximum credible earthquake magnitude statistics and earthquake perceptibility hazard. Earthquake perceptibility hazard is defined as the probability a site perceives ground shaking equal to or greater than a selected ground motion level X, resulting from an earthquake of magnitude M, and develops estimates for the most perceptible earthquake magnitude, M _P(max). Realistic and usable maximum magnitude statistics are obtained from both whole process and part process statistical recurrence models. These approaches are extended to develop relationships between perceptible earthquake magnitude hazard and maximum magnitude recurrence models that are governed by asymptotic and finite return period properties, respectively. Integrated perceptibility curves illustrating the probability of a specific level of perceptible ground motion due to all earthquakes over the magnitude range extending from ?∞ to a magnitude M _i are then developed from reviewing site-specific magnitude perceptibility. These lead on to achieving site-specific annual probability of exceedance hazard curves for the example cities of Sofia and Thessaloniki for both horizontal ground acceleration and ground velocity. Both the maximum credible earthquake magnitude M ₃ and the most perceptible earthquake magnitude M _P(max) are of importance to the earthquake engineer when approaching anti-seismic building design. Both forms of hazard are illustrated using contoured hazard maps for the region bounded by 39°–45°N, 19°–29°E. Patterns are observed for these magnitude hazard estimates—especially M _P(max) specific to horizontal ground acceleration and horizontal ground velocity—and compared to inferred patterns of crustal deformation across the region. The full geographic region considered is estimated to be subject to a maximum credible earthquake magnitude M ₃—estimated using cumulative seismic moment release statistics—of 7.53 M _w, calculated from the full content of the adopted earthquake catalogue, while Bulgaria’s capital, Sofia, is estimated a comparable value of 7.36 M _w. Sofia is also forecast most perceptible earthquake magnitudes for the lowest levels considered for horizontal ground acceleration of M _PA(50) = 7.20 M _w and horizontal ground velocity of M _PV(5) = 7.23 M _w for a specimen focal depth of 15 km.     

Three-dimensional intraplate stress distributions associated with topography and crustal density inhomogeneities beneath the Deccan Volcanic Province

The generation of intraplate earthquakes has been attributed to perturbations in the stress regime, either due to surface and sub-surface loading or strength weakening of the rock mass. The present work aims at estimating the intraplate stresses associated with topography and crustal density inhomogeneities beneath the Deccan Volcanic Province (DVP). A layered crustal model with irregular interfaces of small amplitude has been used for elastostatic stress calculations. The computed principal stress differences show a significant concentration at 5–20 km depths beneath the western side of the region. The maximum magnitude of principal stress difference occurs beneath the Karad at a depth of 10 km with a value of 60 MPa. The deviatoric stress estimates are further superposed on inferred stresses due to the regional plate tectonic forces. These results show principal stress difference concentrations beneath the Koyna, Poona and Karad regions which may thus be more vulnerable to brittle failure. It is also seen that the principal total stress directions point to the strike slip motion at Koyna, similar to that which is associated with the 1967 Koyna earthquake.     

Crustal heterogeneity and seismotectonics of the region around Beijing, China

A detailed three-dimensional (3-D) P-wave velocity model of the crust and uppermost mantle under the Chinese capital (Beijing) region is determined with a spatial resolution of 25 km in the horizontal direction and 4–17 km in depth. We used 48,750 precise P-wave arrival times from 2973 events of local crustal earthquakes, controlled seismic explosions and quarry blasts. These events were recorded by a new digital seismic network consisting of 101 seismic stations equipped with high-sensitivity seismometers. The data are analyzed by using a 3-D seismic tomography method. Our tomographic model provides new insights into the geological structure and tectonics of the region, such as the lithological variations and large fault zones across the major geological terranes like the North China Basin, the Taihangshan and the Yanshan mountainous areas. The velocity images of the upper crust reflect well the surface geological and topographic features. In the North China Basin, the depression and uplift areas are imaged as slow and fast velocities, respectively. The Taihangshan and Yanshan mountainous regions are generally imaged as broad high-velocity zones, while the Quaternary intermountain basins show up as small low-velocity anomalies. Velocity changes are visible across some of the large fault zones. Large crustal earthquakes, such as the 1976 Tangshan earthquake (M=7.8) and the 1679 Sanhe earthquake (M=8.0), generally occurred in high-velocity areas in the upper to middle crust. In the lower crust to the uppermost mantle under the source zones of the large earthquakes, however, low-velocity and high-conductivity anomalies exist, which are considered to be associated with fluids. The fluids in the lower crust may cause the weakening of the seismogenic layer in the upper and middle crust and thus contribute to the initiation of the large crustal earthquakes.     

Prediction of strong earthquake within the southwestern shelf of Sakhalin Island and its realization during the August 2, 2007, Nevelsk earthquake

Materials of the long- and short-term predictions of the destructive earthquake with the magnitude M _LH = 6.6 ± 0.6 within the southwestern shelf of Sakhalin Island are described. The long-term prediction was issued in December 2005 and was affirmed by the Russian Council of Experts on Earthquake Forecasting and Seismic Hazard Assessment in August 2006. The August 17(18), 2006, Gornozavodsk earthquake with a magnitude of M _w = 5.6 was the beginning of the realization of this prediction. Six days after its occurrence, the short-term prediction of a much more serious seismic event in the alarm region was prepared. One year later, the prediction of the August 2, 2007, Nevelsk earthquake with a magnitude of M _w = 6.2 (M _LH = 6.2) proved to be correct.     

Seismic hazard monitoring in Kamchatka and its applications to the M = 6.6 earthquake of 6 October 1987

The earthquake of 6 October 1987 (M = 6.6), which occurred near the Shipunsky Cape, Kamchatka, was the largest crustal event in the vicinity of the main city of Kamchatka — Petropavlovsk-Kamchatsky — during the last three decades. It was followed by numerous aftershocks. This earthquake allowed us to test the effectiveness of the seismic hazard monitoring in Kamchatka, including the seismological, geodetic and hydrogeochemical surveys. The seismic survey provided the location and source nature of the main shock and aftershocks and the seismic environment of the main shock. The geodetic and hydrogeochemical surveys have yielded data on the response to earthquakes of the Earth's surface deformations, water level, and chemical elements concentration in the underground water. As a result, the following data were obtained:

u

The earthquake of 6 October had a seismic moment 4–10 E18 Nm, thrust type of faulting and the source volume of 20 × 20 × 10 km³. The maximum intensity was VI–VII (MSK-64 scale) and maximum acceleration 88 cm/s².

Before this event, a relative increase in the number of the upper mantle (depth more than 100 km) moderate magnitude earthquakes during 5 years and a one-year period of seismic quiescence for small shallow earthquakes, were recognized. Significant anomalies in HCO₃ and H₃BO₃ concentrations in the underground waters were observed in the wells a week before the main shock.