Seismicity assessment in and around Syria based on instrumental data: application of Gumbel distributions and Gutenberg-Richter relationship

An instrumental earthquake catalog covering the time span between 1903 and 2007 and for the area bounded by 32°N–38°N and 35°E–43°E has been compiled in this research. The catalog has a magnitude of completeness (M _c ) with 3.5. Least squares and statistical probability Gumbel’s techniques with different approaches have been applied on the instrumental events in order to assess the average recurrence time periods for different earthquake magnitudes. The constants a and b of Gutenberg-Richter and the average recurrence times have been computed firstly for the study area and secondly for the central and northern parts of Dead Sea fault system. The different statistical computations using Knopoff and Kagan formalism are generally in agreement and suggest an average recurrence time of 203 years for an earthquake of magnitude 7 for the region. The occurrence of large well-documented historical earthquakes in Lebanon and western Syria, the existence of active fault segments, the absence of large earthquakes during the study period, the increasing number of the low-magnitude earthquakes, and the continued accumulation of the strain since 1900 indicate therefore the probability of an earthquake occurrence of a large magnitude. This should be permanently taken into consideration in seismic hazard assessment on the local and regional scales.     

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.     

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.     

Aftershock sequences of two great Sumatran earthquakes of 2004 and 2005 and simulation of the minor tsunami generated on September 12, 2007 in the Indian Ocean and its effect

We present the seismic energy, strain energy, frequency–magnitude relation (b-value) and decay rate of aftershocks (p-value) for the aftershock sequences of the Andaman–Sumatra earthquakes of December 26, 2004 (M _w 9.3) and March 28, 2005 (M _w 8.7). The energy released in aftershocks of 2004 and 2005 earthquake was 0.135 and 0.365% of the energy of the respective mainshocks, while the strain release in aftershocks was 39 and 71% for the two earthquakes, respectively. The b-value and p-value indicate normal value of about 1. All these parameters are in normal range and indicate normal stress patterns and mechanical properties of the medium. Only the strain release in aftershocks was considerable. The fourth largest earthquake in this region since 2004 occurred in September 2007 off the southern coast of Island of Sumatra, generating a relatively minor tsunami as indicated by sea level gauges. The maximum wave amplitude as registered by the Padang, tide gauge, north of the earthquake epicenter was about 60 cm. TUNAMI-N2 model was used to investigate ability of the model to capture the minor tsunami and its effect on the eastern Indian Coast. A close comparison of the observed and simulated tsunami generation, propagation and wave height at tide gauge locations showed that the model was able to capture the minor tsunami phases. The directivity map shows that the maximum tsunami energy was in the southwest direction from the strike of the fault. Since the path of the tsunami for Indian coastlines is oblique, there were no impacts along the Indian coastlines except near the coast of epicentral region.     

Estimation of earthquake hazard parameters for incomplete and uncertain data files

The maximum likelihood estimation of earthquake hazard parameters (maximum regional magnitudem _max, activity rate λ, and theb parameter in the Gutenberg-Richter distribution) is extended to the cases of incomplete and uncertain data. The method accepts mixed data containing only large (extreme) events and a variable quality of complete data with different threshold magnitude values. Uncertainty of earthquake magnitude is specified by two values, the lower and upper magnitude limits. It is assumed that such an interval contains the real unknown magnitude. The proposed approach allows the combination of different quality catalog parts, e.g. those where the assignment of magnitude is questionable and those with magnitudes precisely determined. As an illustration of the method, the seismic hazard analysis for western Norway and adjacent sea area (4–8°E, 58–64°N) is presented on the basis of the strongest earthquakes felt during the period 1831–1889 and three complete catalog parts, covering the period 1890–1987.     

Spatial-temporal variability of seismic hazard in peninsular India

This paper examines the variability of seismic activity observed in the case of different geological zones of peninsular India (10°N–26°N; 68°E–90°E) based on earthquake catalog between the period 1842 and 2002 and estimates earthquake hazard for the region. With compilation of earthquake catalog in terms of moment magnitude and establishing broad completeness criteria, we derive the seismicity parameters for each geologic zone of peninsular India using maximum likelihood procedure. The estimated parameters provide the basis for understanding the historical seismicity associated with different geological zones of peninsular India and also provide important inputs for future seismic hazard estimation studies in the region. Based on present investigation, it is clear that earthquake recurrence activity in various geologic zones of peninsular India is distinct and varies considerably between its cratonic and rifting zones. The study identifies the likely hazards due to the possibility of moderate to large earthquakes in peninsular India and also presents the influence of spatial rate variation in the seismic activity of this region. This paper presents the influence of source zone characterization and recurrence rate variation pattern on the maximum earthquake magnitude estimation. The results presented in the paper provide a useful basis for probabilistic seismic hazard studies and microzonation studies in peninsular India.     

Effect of Aftershocks on Earthquake Hazard Estimation: An Example from the North Anatolian Fault Zone

In order to investigate the effect of aftershocks on earthquake hazard estimation, earthquake hazard parameters (_m, b and M_max) have been estimated by the maximum likelihood method from the main shocks catalogue and the raw earthquakes catalogue for the North Anatolian Fault Zone (NAFZ). The main shocks catalogue has been compiled from the raw earthquake catalogue by eliminating the aftershocks using the window method. The raw earthquake catalogue consisted of instrumentally detected earthquakes between 1900 and 1992, and historical earthquakes that occurred between 1000–1900. For the events of the mainshock catalogue the Poisson process is valid and for the raw earthquake catalogue it does not fit. The paper demonstrates differences in the hazard outputs if on one hand the main catalogues and on the other hand the raw catalogue is used. The maximum likelihood method which allows the use of the mixed earthquake catalogue containing incomplete (historical) and complete (instrumental) earthquake data is used to determine the earthquake hazard parameters. The maximum regional magnitude (M_max, the seismic activity rate (_m), the mean return period (R) and the b value of the magnitude-frequency relation have been estimated for the 24°–31° E, 31°–41° E, 41°–45° E sections of the North Anatolian Fault Zone from the raw earthquake catalogue and the main shocks catalogue. Our results indicate that inclusion of aftershocks changes the b value and the seismic activity rate _m depending on the proportion of aftershocks in a region while it does not significantly effect the value of the maximum regional magnitude since it is related to the maximum observed magnitude. These changes in the earthquake hazard parameters caused the return periods to be over- and underestimated for smaller and larger events, respectively.     

Seismic hazard in mega city Kolkata, India

The damages caused by recent earthquakes in India have been a wake up call for people to take proper mitigation measures, especially the major cities that lie in the high seismic hazard zones. Kolkata City, with thick sediment deposit (∼12 km), one of the earliest cities of India, is an area of great concern as it lies over the Bengal Basin and lies at the boundary of the seismic zones III and IV of the zonation map of India. Kolkata has been affected by the 1897 Shillong earthquake, the 1906 Calcutta earthquake, and the 1964 Calcutta earthquake. An analysis on the maximum magnitude and b-value for Kolkata City region is carried out after the preparation of earthquake catalog from various sources. Based on the tectonic set-up and seismicity of the region, five seismic zones are delineated, which can pose a threat to Kolkata in the event of an earthquake. They are broadly classified as Zone 1: Arakan-Yoma Zone (AYZ), Zone 2: Himalayan Zone (HZ), Zone 3: Shillong Plateau Zone (SPZ), Zone 4: Bay of Bengal Zone (BBZ), and Zone 5: Shield Zone (SZ). The maximum magnitude (m _max) for Zones 1, 2, 3, 4, and 5 are 8.30 ± 0.51, 9.09 ± 0.58, 9.20 ± 0.51, 6.62 ± 0.43 and 6.61 ± 0.43, respectively. A probability of 10% exceedance value in 50 years is used for each zone. The probabilities of occurrences of earthquakes of different magnitudes for return periods of 50 and 100 years are computed for the five seismic zones. The Peak Ground Acceleration (PGA) obtained for Kolkata City varies from 0.34 to 0.10 g.     

Seismogenic Tectonics and Dynamics of the 2011 Ms5.9 Yingjiang Earthquake in Yunnan,China

In the southern South–North Seismic Zone, China, seismic activity in the Yingjiang area of western Yunnan increased from December 2010, and eventually a destructive earthquake of Ms5.9 occurred near Yingjiang town on 10 March 2011. The focal mechanism and hypocenter location of the mainshock suggest that the Dayingjiang Fault was the site of the mainshock rupture. However, most of foreshocks and all aftershocks recorded by a portable seismic array located close to the mainshock occurred along the N–S-striking Sudian Fault, indicating that this fault had an important influence on these shocks. Coulomb stress calculations show that three strong(magnitude ≥5.0) earthquakes that occurred in the study region in 2008 increased the coulomb stress along the plane parallel to the Dayingjiang Fault. This supports the Dayingjiang Fault, and not the Sudian Fault, as the seismogenic fault of the 2011 Ms5.9 Yingjiang earthquake. The strong earthquakes in 2008 also increased the Coulomb stress at depths of ≤5 km along the entire Sudian Fault, and by doing so increased the shallow seismic activity along the fault. This explains why the foreshocks and aftershocks of the 2011 Yingjiang earthquake were located mostly on the Sudian Fault where it cuts the shallow crust. The earthquakes at the intersection of the Sudian and Dayingjiang faults are distributed mainly along a belt that dips to the southeast at ~40°, suggesting that the Dayingjiang Fault in the mainshock area also dips to the southeast at ~40°.     

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.     

Maximum Likelihood Estimation of Seismic Hazard for Sweden

The maximum magnitude, the activity rate, and the Gutenberg-Richterb parameter as earthquake hazard parameters, have been evaluated for Sweden. The maximum likelihood method permits the combination of historical and instrumental data. The catalog used consists of 1100 earthquakes in the time interval 1375–1989. The extreme part of the catalog contains only the strongest historical earthquakes, whereas the complete part is divided into several subcatalogs, each assumed complete above a specified threshold magnitude. The uncertainty in magnitude determination was taken into account. For southern Sweden, the calculations giveb-values of 1.04 (0.05) for the whole area south of 60° N and 0.98 (0.06) for a subregion of enhanced seismicity in the Lake Vänern area. For the whole area north of 60° N, theb-value is 1.35 (0.06) and for the seismicity zone along the Gulf of Bothnia 1.26 (0.06). The number of annually expected earthquakes with magnitude equal to or larger than 2.4 [M_L(UPP) or M_M(UPP)] is 1.8 for the whole southern Sweden, 1.3 for the Lake Vänern region, 3.7 for northern Sweden, and 2.4 for the region along the Gulf of Bothnia. The maximum expected regional magnitude is calculated to 4.9 (0.5) for a time span of 615 years for southern Sweden and the Lake Vänern subregion, and 4.3 (0.5) for a time span of 331 years for northern Sweden and the Gulf of Bothnia subregion. However, several historical earthquakes with magnitude above 5 in nearby areas of Norway indicate that the seismic potential may be higher.     

Seismicity of the Baikal rift system from regional network observations

In the paper we report the state-of-the-art of seismicity study in the Baikal rift system and the general results obtained. At present, the regional earthquake catalog for fifty years of the permanent instrumental observations consists of over 185,000 events. The spatial distribution of the epicenters, which either gather along well-delineated belts or in discrete swarms is considered in detail for different areas of the rift system. At the same time, the hypocenters are poorly constrained making it difficult to identify the fault geometry. Clustered events like aftershock sequences or earthquake swarms are typical patterns in the region; moreover, aftershocks of M ⩾ 4.7 earthquakes make up a quarter of the whole catalog. The maximum magnitude of earthquakes recorded instrumentally is M_LH7.6 for a strike-slip event in the NE part of the Baikal rift system and M_LH6.8 for a normal fault earthquake in the central part of the rift system (Lake Baikal basin). Predominant movement type is normal faulting on NE striking faults with a left lateral strike-slip component on W–E planes. In conclusion, some shortcomings of the seismic network and data processing are pointed out.     

Low b-value prior to the Indo-Myanmar subduction zone earthquakes and precursory swarm before the May 1995 M 6.3 earthquake

Some 455 events (m_b ⩾ 4.5) in the Indo-Myanmar subduction zone are compiled using the ISC/EHB/NEIC catalogues (1964–2011) for a systematic study of seismic precursors, b-value and swarm activity. Temporal variation of b-value is studied using the maximum likelihood method beside CUSUM algorithm. The b-values vary from 0.95 to 1.4 for the deeper (depth ⩾60 km) earthquakes, and from 0.85 to 1.3 for the shallower (depth <60 km) earthquakes. A sudden drop in the b-value, from 1.4 to 0.9, prior to the occurrence of larger earthquake(s) at the deeper depth is observed. It is also noted that the CUSUM gradient reversed before the occurrence of larger earthquakes. We further examined the seismicity pattern for the period 1988–1995 within a radius of 150 km around the epicentre (latitude: 24.96°N; longitude: 95.30°E) of a deeper event M 6.3 of May 6, 1995 in this subduction zone. A precursory swarm during January 1989 to July 1992 and quiescence during August 1992 to April 1995 are identified before this large earthquake. These observations are encouraging to monitor seismic precursors for the deeper events in this subduction zone.     

Recent plumbing system of the Krakatau volcano revealed by teleseismic earthquake distribution

Spatial and temporal analysis of global seismological data 1964–2005 reveals a distinct teleseismic earthquake activity producing a columnar-like formation in the continental wedge between the Krakatau volcano at the surface and the subducting slab of the Indo-Australian plate. These earthquakes occur continuously in time, are in the body-wave (m _b) magnitude range 4.5–5.3 and in the depth range 1–100 km. The Krakatau earthquake cluster is vertical and elongated in the azimuth N30°E, suggesting existence of a deep-rooted fault zone cutting the Sunda Strait in the SSW-NNE direction. Possible continuation of the fault zone in the SW direction was activated by an intensive 2002/2003 aftershock sequence, elongated in the azimuth of N55°E. Beneath the Krakatau earthquake cluster, an aseismic gap exists in the Wadati-Benioff zone of the subducting plate at the depths 100–120 km. We interpret this aseismic gap as a consequence of partial melting inhibiting stress concentration necessary to generate stronger earthquakes, whereas the numerous earthquakes observed in the overlying lithospheric wedge beneath the volcano probably reflect magma ascent in the recent plumbing system of the Krakatau volcano. Focal depth of the deepest events (~100 km) of the Krakatau cluster constrains the location of the primary magma generation to greater depths. The ascending magmatic fluids stress fault segments within the Sunda Strait fault zone and change their friction parameters inducing the observed tectonic earthquakes beneath Krakatau.     

On the aftershock sequence of a 4.6 mb earthquake of the Garhwal Himalaya

Locally recorded data for eighteen aftershocks of a magnitude(m_b) 4.6 earthquake occurring near Ukhimath in the Garhwal Himalaya were analysed. A master event technique was adopted to locate seventeen individual aftershock hypocentres relative to the hypocentre of the eighteenth aftershock chosen as the master event. The aftershock epicentres define an approximately 30 km² rupture zone commensurate with the magnitude of the earthquake. The distribution of epicentres within this zone and the limited amount of first motion data support the view that a group of parallel, sub-vertical, sinistral strike-slip faults oriented N46°, transverse to the regional NW-SE trend of the Garhwal Himalaya, was involved in this seismic episode. Since the estimated focal depth range for aftershocks of this sequence is 3–14 km, we infer that this transverse fault zone extends through the upper crustal layer to a depth of 14 km at least.     

Seismicity pattern in north Sumatra-Great Nicobar region: In search of precursor for the 26 December 2004 earthquake

We analyse the seismicity pattern including b-value in the north Sumatra-Great Nicobar region from 1976 to 2004. The analysis suggests that there were a number of significant, intermediate and short-term precursors before the magnitude 7.6 earthquake of 2 November 2002. However, they were not found to be so prominent prior to the magnitude 9.0 earthquake of 26 December 2004 though downward migration of activity and a 50-day short-term quiescence was observed before the event. The various precursors identified include post-seismic and intermediate-term quiescence of 13 and 10 years respectively, between the 1976 (magnitude 6.3) and 2002 earthquakes with two years (1990–1991) of increase in background seismicity; renewed seismicity, downward migration of seismic activity and foreshocks in 2002, just before the mainshock. Spatial variation in b-value with time indicates precursory changes in the form of high b-value zone near the epicenter preceding the mainshocks of 2004 and 2002 and temporal rise in b-value in the epicentral area before the 2002 earthquake.     

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.     

A procedure to estimate the seismic hazard in an urban area: an application to Acireale (Eastern Sicily)

A “standard procedure” to characterize the seismic hazard of a given area was proposed. It is based on a multidisciplinary approach implying: (1) the knowledge of the seismic history of the area; (2) detailed geological surveys; (3)seismic noise measurements; (4) simulations of earthquake scenarios. The downtown of Acireale, a typical baroque town located on Eastern Sicily, was chosen as the “test area”. A catalog of the local seismogenic faults (able to generate earthquakes in historical times) has been compiled, as well as a seismic catalog for the effects of both local and regional earthquakes. The analysis of both catalogs allowed us to make the following conclusions: (1) the most important seismogenic faults affecting the Acireale municipality do not affect the downtown, while the related local earthquakes attenuate their energy (and intensity) in short (few km) distances; (2) the highest seismic intensity (degree X) experienced in Acireale downtown was caused by the 1693 regional earthquake; (3) over the last 140 years, the downtown has experienced the highest intensity value of VII only once, while six times the intensity was VI. On the whole, this implies a moderate seismic hazard. The estimation of the seismic hazard has been also approached by the experimental method of recording seismic noise. Measurements have been performed at seven different sites, where drills gave detailed information on the shallow subsurface geology to obtain HV (horizontal/vertical) spectral ratios. On the whole, the highest site amplification factor was moderate (about 7). A further investigation based on synthetic seismograms (and spectra) produced by simulating two given earthquake scenarios was also performed. The two scenarios are, respectively, representative of the largest expected earthquake in the area (the 1693 shock) and of a moderate (magnitude ca. 5.5) local earthquake (as the 1818 one). Moderate to strong locally expected accelerations were evidenced.     

Premonitory crustal deformations, strains and seismotectonic features (b-values) preceding Koyna earthquakes

There have been instances of premonitory variations in tilts, displacements, strains, telluric current, seismomagnetic effects, seismic velocities ( V_p, V_s) and their ratio (V_p/V_s), b-values, radon emission, etc. preceding large and moderate earthquakes, especially in areas near epicentres and along faults and other weak zones. Intensity and duration (T) of these premonitory quantities are very much dependent on magnitude (M) of the seismic event. Hence, these quantities may be utilised for prediction of an incoming seismic event well in advance of the actual earthquake. In the recent past, tilts, strain in deep underground rock and crustal displacements have been observed in the Koyna earthquake region over a decade covering pre- and postearthquake periods; and these observations confirm their reliability for qualitative as well as quantitative premonitory indices. Tilt began to change significantly one to two years before the Koyna earthquake of December 10, 1967, of magnitude 7.0. Sudden changes in ground tilt measured in a watertube tiltmeter accompanied an earthquake of magnitude 5.2 on October 17, 1973 and in other smaller earthquakes in the Koyna region, though premonitory changes in tilt preceding smaller earthquakes were not so much in evidence. However, changes in strains in deep underground rock were observed in smaller earthquakes of magnitude 4.0 and above. Furthermore, as a very large number of earthquakes (M = 1–7.0) were recorded in the extensive seismic net in the Koyna earthquake region during 1963–1975, precise b-value variations as computed from the above data, could reveal indirectly the state of crustal (tectonic) strain variations in the earthquake focal region and consequently act as a powerful premonitory index, especially for the significant Koyna earthquakes of December 10, 1967 (M = 7.0) and October 17, 1973 (M = 5.2). The widespread geodetic and magnetic levelling observations covering the pre- and postearthquake periods indicate significant vertical and horizontal crustal displacements, possibly accompanied by large-scale migration of underground magma during the large seismic event of December 10, 1967 in the Koyna region (M = 7.0). Duration (T) of premonitory changes in tilt, strains, etc., is generally governed by the equation of the type logT = A + BM (A and B are statistically determined coefficients). Similar other instances of premonitory evidences are also observed in micro-earthquakes (M = − 1 to 2) due to activation of a fault caused by nearby reservoir water-level fluctuations.     

On the problem of destructive Iranian earthquakes and their causative faults

This article is devoted to evaluating destructive earthquakes (magnitude >6) of Iran and determining properties of their source parameters. First of all, a database of documented earthquakes has been prepared via reliable references and causative faults of each event have been determined. Then, geometric parameters of each fault have been presented completely. Critical parameters such as Maximum Credible Rupture, MCR, and Maximum Credible Earthquake, MCE, have been compiled based on the geometrical parameters of the earthquake faults. The calculated parameters have been compared to the maximum earthquake and the surface rupture which have been recorded for the earthquake faults. Also, the distance between the epicenter of documented earthquake events and their causative faults has been calculated (the distance was less than 20 km for 90% of the data). Then, the distance between destructive earthquakes (with the magnitude more than 6) and the nearest active fault has been calculated. If the estimated distance is less than 20 km and the mechanism of the active fault and the event are reported the same, the active fault will be introduced as a probable causative fault of that earthquake. In the process, all of the available geological, tectonic, seismotectonic maps, aerial geophysical data as well as remote sensing images have been evaluated. Based on the quality and importance of earthquake data, the events have been classified into three categories: (1) the earthquakes which have their causative faults documented, (2) the events with magnitude higher than 7, and (3) the events with the magnitude between 6 and 7. For each category, related maps and tables have been compiled and presented. Some important faults and events have been also described throughout the paper. As mentioned in this paper, these faults are likely to be in high seismic regions with potential for large-magnitude events as they are long, deep and bound sectors of the margins characterized by different deformation and coupling rates on the plate interface.