Probabilistic assessment of earthquake recurrence in the January 26, 2001 earthquake region of Gujrat, India

Kutch region of Gujrat is one of the most seismic prone regions of India. Recently, it has been rocked by a large earthquake (M _w = 7.7) on January 26, 2001. The probabilities of occurrence of large earthquake (M≥6.0 and M≥5.0) in a specified interval of time for different elapsed times have been estimated on the basis of observed time-intervals between the large earthquakes (M≥6.0 and M≥5.0) using three probabilistic models, namely, Weibull, Gamma and Lognormal. The earthquakes of magnitude ≥5.0 covering about 180 years have been used for this analysis. However, the method of maximum likelihood estimation (MLE) has been applied for computation of earthquake hazard parameters. The mean interval of occurrence of earthquakes and standard deviation are estimated as 20.18 and 8.40 years for M≥5.0 and 36.32 and 12.49 years, for M≥6.0, respectively, for this region. For the earthquakes M≥5.0, the estimated cumulative probability reaches 0.8 after about 27 years for Lognormal and Gamma models and about 28 years for Weibull model while it reaches 0.9 after about 32 years for all the models. However, for the earthquakes M≥6.0, the estimated cumulative probability reaches 0.8 after about 47 years for all the models while it reaches 0.9 after about 53, 54 and 55 years for Weibull, Gamma and Lognormal model, respectively. The conditional probability also reaches about 0.8 to 0.9 for the time period of 28 to 40 years and 50 to 60 years for M≥5.0 and M≥6.0, respectively, for all the models. The probability of occurrence of an earthquake is very high between 28 to 42 years for the magnitudes ≥5.0 and between 47 to 55 years for the magnitudes ≥6.0, respectively, past from the last earthquake (2001).     

Evaluation of seismicity and seismic hazard parameters in Turkey and surrounding area using a new approach to the Gutenberg–Richter relation

In this study, the spatial distributions of seismicity and seismic hazard were assessed for Turkey and its surrounding area. For this purpose, earthquakes that occurred between 1964 and 2004 with magnitudes of M ≥ 4 were used in the region (30–42°N and 20–45°E). For the estimation of seismicity parameters and its mapping, Turkey and surrounding area are divided into 1,275 circular subregions. The b-value from the Gutenberg–Richter frequency–magnitude distributions is calculated by the classic way and the new alternative method both using the least-squares approach. The a-value in the Gutenberg–Richter frequency–magnitude distributions is taken as a constant value in the new alternative method. The b-values calculated by the new method were mapped. These results obtained from both methods are compared. The b-value shows different distributions along Turkey for both techniques. The b-values map prepared with new technique presents a better consistency with regional tectonics, earthquake activities, and epicenter distributions. Finally, the return period and occurrence hazard probability of M ≥ 6.5 earthquakes in 75 years were calculated by using the Poisson model for both techniques. The return period and occurrence hazard probability maps determined from both techniques showed a better consistency with each other. Moreover, maps of the occurrence hazard probability and return period showed better consistency with the b-parameter seismicity maps calculated from the new method. The occurrence hazard probability and return period of M ≥ 6.5 earthquakes were calculated as 90–99% and 5–10 years, respectively, from the Poisson model in the western part of the studying region.     

Fracture characteristics of the 1997 Jiashi, Xinjiang, China, earthquake swarm inferred from source spectra

Broadband P and S waves source spectra of 12 MS5.0 earthquakes of the 1997 Jiashi, Xinjiang, China, earthquake swarm recorded at 13 GDSN stations have been analyzed. Rupture size and static stress drop of these earthquakes have been estimated through measuring the corner frequency of the source spectra. Direction of rupture propagation of the earthquake faulting has also been inferred from the azimuthal variation of the corner frequency. The main results are as follows: ①The rupture size of MS6.0 strong earthquakes is in the range of 10～20 km, while that of MS=5.0～5.5 earthquakes is 6～10 km.② The static stress drop of the swarm earthquakes is rather low, being of the order of 0.1 MPa. This implies that the deformation release rate in the source region may be low. ③ Stress drop of the earthquakes appears to be proportional to their seismic moment, and also to be dependent on their focal mechanism. The stress drop of normal faulting earthquakes is usually lower than that of strike-slip type earthquakes. ④ For each MS6.0 earthquake there exists an apparent azimuthal variation of the corner frequencies. Azimuthally variation pattern of corner frequencies of different earthquakes shows that the source rupture pattern of the Jiashi earthquake swarm is complex and no uniform rupture expanding direction exists.     

On the seismicity of central Kamchatka before the 1997, M = 7.9 Kronotskii earthquake

We report results from a detailed study of seismicity in central Kamchatka for the period from 1960 to 1997 using a modified traditional approach. The basic elements of this approach include (a) segmentation of the seismic region concerned (the Kronotskii and Shipunskii geoblocks, the continental slope and offshore blocks), (b) studying the variation in the rate of M = 4.5–7.0 earthquakes and in the amount of seismic energy release over time, (c) studying the seismicity variations, (d) separate estimates of earthquake recurrence for depths of 0–50 and 50–100 km. As a result, besides corroborating the fact that a quiescence occurred before the December 5, 1997, M = 7.9 Kronotskii earthquake, we also found a relationship between the start of the quiescence and the position of the seismic zone with respect to the rupture initiation. The earliest date of the quiescence (decreasing seismicity rate and seismic energy release) was due to the M = 4.5–7.0 earthquakes at depths of 0–100 km in the Kronotskii geoblock (8–9 years prior to the earthquake). The intermediate start of the quiescence was due to distant seismic zones of the Shipunskii geoblock and the circular zone using the RTL method, combining the Shipunskii and Kronotskii geoblocks (6 years). Based on the low magnitude seismicity (M≥2.6) at depths of 0–70 km in the southwestern part of the epicentral zone (50–100 km from the mainshock epicenter), the quiescence was inferred to have occurred a little over 3 years (40 months) before the mainshock time and a little over 2 years (25 months) in the immediate vicinity of the epicenter (0–50 km). These results enable a more reliable identification of other types of geophysical precursors during seismic quiescences before disastrous earthquakes.     

Time distribution characteristics of regional macroseismic activity in the Sichuan-Yunnan region and its significance to mid-long term prediction

IntroductionThe development and occurrence of macroseismic activity come from the specific structUreenvironment and stress condition. So the space-time distribUtion of strong quakes appears to bevery inhomogeneous. Earthquakes with Ms27.0 in a seismicity period presented often a speeding-up pattern with time in different seismic provinces. That is, the cumulative frequency of earthquakes with Ms27.0 increase with exponent N(t)=ae', in a seismicity period (ZHANG, FU, 1989).It means that t…     

Evolution of cumulative Coulomb failure stress in northeastern Qinghai-Xizang （Tibetan） Plateau and its effect on large earthquake occurrence

We simulate accumulative Coulomb failure stress change in a layered Maxwell viscoelastic media in the north-eastern Qinghai-Xizang(Tibetan)Plateau since 1920.Lithospheric stress/strain evolution is assumed to be drivenby dislocations of large earthquakes(M≥7.0)and secular tectonic loading.The earthquake rupture parameters suchas the fault rupture length,width,and slip are either adopted from field investigations or estimated from their sta-tistic relationships with the earthquake magnitudes and seismic moments.Our study shows that among 20 largeearthquakes(M≥7.0)investigated,17 occurred in areas where the Coulomb failure stress change is positive,with atriggering rate of 85%.This study provides essential data for the intermediate to long-term likelihood estimation oflarge earthquakes in the northeastern Tibetan Plateau.     

Intermediate-term precursors of large (<Emphasis Type="Italic">M</Emphasis> ≥ 6.6) Kamchatka earthquakes for the period from 1987 to 2004: A retrospective assessment of their information content for prediction

Data from the literature were used to systematize intermediate-term (with advance times of 1 month to ∼2.5 years) precursors to the M ≥ 6.6 Kamchatka earthquakes of 1987–2004. The precursors were observed as changes in seismological, geodetic, geophysical, water-level, and hydrochemical parameters. Retrospective assessment of the information content in these intermediate-term precursors for earthquake prediction is in progress. The focus was on estimating the occurrence times of various precursors as functions of earthquake parameters (magnitude M, hypocentral distance R, and epicenter location). In the conditions of the Kamchatka observing network, precursors can be identified by a combination of methods, mostly before M ∼ 7 earthquakes or greater south of the Kronotskii Peninsula, for which M/logR ≥ 3. It is shown that the relative proportion of earthquakes for which precursors have been identified in the observations considered here is 0.43–0.86.     

Periodic spectral characteristics of seismicity before strong earthquakes and their application

IntroductionMaximumentropyspectralmethod(MEM)(Burg,1972)hadbeenamethodusuallyusedinstudyingtheseismicityanditsmainpurposeistofindthedominantspectrainthelong-termseismicityprocessesinthepastyears(Zhu,1985).Inthispaper,themethodisappliedtostudywhethertherearesomespecialspectraofseismicityinsomespecificstagesinearthquake-generatingprocesses.Sowestudyseparatelythenormalandabnormalstageofearthquakeactivity,whoseactiveprocessisregardedasstablestochasticprocess,inordertofindtheirspectracharactersan…     

Earthquake probabilities and magnitude distribution ( M ≥6.7) along the Haiyuan fault, northwestern China

In recent years, some researchers have studied the paleoearthquake along the Haiyuan fault and revealed a lot of paleoearthquake events. All available information allows more reliable analysis of earthquake recurrence interval and earthquake rupture patterns along the Haiyuan fault. Based on this paleoseismological information, the recurrence probability and magnitude distribution for M≥6.7 earthquakes in future 100 years along the Haiyuan fault can be obtained through weighted computation by using Poisson and Brownian passage time models and considering different rupture patterns. The result shows that the recurrence probability of M _S≥6.7 earthquakes is about 0.035 in future 100 years along the Haiyuan fault. Foundation item: Joint Seismological Science Foundation of China (103034) and Major Research “Research on Assessment of Seismic Safety” from China Earthquake Administration during the tenth Five-year Plan.     

Statistical prediction of the next great earthquake around Tehran,Iran

The analysis of seismic activity variations with space and time is a complex problem. Several statistical methods have been adopted to study these variations. One of the tasks that has attracted the attention of the seismological and statistical community is to explain seismicity patterns by statistical models and apply the results for earthquake prediction. Here the probability distribution of recurrence times as described by Exponential, Gamma, Lognormal, Pareto, Rayleigh and Weibull probability distributions and the idea of conditional probability has been applied to predict the next great (Ms ≥ 6.0 and Ms ≥ 6.5) earthquake around Tehran (r ≤ 200 km). Conditional probability specifies the likelihood that a given earthquake will happen within a specified time. This likelihood is based on the information about past earthquake occurrences in the given region and the basic assumption that future seismic activity will follow the pattern of past activity. The rapid growth of Tehran to approximately 12 million inhabitants has resulted in a much more rapid increase in its vulnerability to natural disasters, especially earthquakes. Several earthquakes affected this region in the past, mostly on the Mosha, Taleqan, Eyvankey and Garmsar faults. The estimated recurrence times for Exponential, Gamma, Lognormal, Pareto, Rayleigh and Weibull distributions has been computed to be 66.64, 14.79, 26.88, 2.37, 67.58 and 80.47, respectively. Accordingly, one may expect that a large damaging earthquake may occur around Tehran approximately every 10 years.     

Present patterns of decelerating–accelerating seismic strain in South Japan

Decelerating generation of preshocks in a narrow (seismogenic) region and accelerating generation of other preshocks in a broader (critical) region, called decelerating–accelerating seismic strain (D-AS) model has been proposed as appropriate for intermediate-term earthquake prediction. An attempt is made in the present work to identify such seismic strain patterns and estimate the corresponding probably ensuing large mainshocks (M ≥ 7.0) in south Japan (30–38° N, 130–138° E). Two such patterns have been identified and the origin time, magnitude, and epicenter coordinates for each of the two corresponding probably ensuing mainshocks have been estimated. Model uncertainties of predicted quantities are also given to allow an objective forward testing of the efficiency of the model for intermediate-term earthquake prediction.     

Possible precursory accelerated Benioff strain in the region of Sistan Suture Zone of Eastern Iran

We investigated whether accelerated seismic strain release precedes large earthquakes occurring in and around the Sistan Suture Zone, Eastern Iran. Online catalogs of teleseismic events occurring post-1960 within the region 27.0°–37.0°N, 55.0°–65.0°E, report five M _w > 7.0 earthquakes, namely, 1968 Dasht-e-Bayaz, 1978 Tabas, 1979 Khuli-Buniabad, 1981 Sirch and 1997 Zirkuh-e-Q’aenat events. We defined four earthquake test episodes, 1968–1978, 1978–1981, 1979–1981, and 1981–1997, with all catalogued intermediate events having magnitudes within 2.0 units that of the final large event. Using the 1968 event as the starting point, we investigated possible increased moderate earthquake activity patterns prior to the large events of 1978, 1981 and 1997 by examining if the cumulative Benioff strain released from such preceding events followed a power law time-to-failure. Our investigation seem to suggest that the 1978, 1981 and 1997 events (i) followed a period of accelerated moderate earthquake activity and (ii) the radius of their optimal critical region, R, scaled with their magnitude, M, according to the scaling law log R ∝ 0.36 M. Our suggestions conform to those proposed by similar investigations in varied seismotectonic regimes.     

Probabilistic Assessment of Earthquake Hazard in Gujarat and Adjoining Region of India

The Gujarat and adjoining region falls under all four seismic zones V, IV, III and II of the seismic zoning map of India, and is one of the most seismically prone intracontinental regions of the world. It has experienced two large earthquakes of magnitude M _w 7.8 and 7.7 in 1819 and 2001, respectively and several moderate earthquakes during the past two centuries. In the present study, the probability of occurrence of earthquakes of M ≥ 5.0 has been estimated during a specified time interval for different elapsed times on the basis of observed time intervals between earthquakes using three stochastic models namely, Weibull, Gamma and Lognormal. A complete earthquake catalogue has been used covering the time interval of 1819 to 2006. The whole region has been divided into three major seismic regions (Saurashtra, Mainland Gujarat and Kachchh) on the basis of seismotectonics and geomorphology of the region. The earthquake hazard parameters have been estimated using the method of maximum likelihood. The logarithmic of likelihood function (ln L) is estimated and used to test the suitability of models in three different regions. It was found that the Weibull model fits well with the actual data in Saurashtra and Kachchh regions, whereas Lognormal model fits well in Mainland Gujarat. The mean intervals of occurrence of earthquakes are estimated as 40.455, 20.249 and 13.338 years in the Saurashtra, Mainland Gujarat and Kachchh region, respectively. The estimated cumulative probability (probability that the next earthquake will occur at a time later than some specific time from the last earthquake) for the earthquakes of M ≥ 5.0 reaches 0.9 after about 64 years from the last earthquake (1993) in Saurashtra, about 49 years from the last earthquake (1969) in Mainland Gujarat and about 29 years from the last earthquake (2006) in the Kachchh region. The conditional probability (probability that the next earthquake will occur during some specific time interval after a certain elapsed time from last earthquake) is also estimated and it reaches about 0.8 to 0.9 during the time interval of about 57 to 66 years from the last earthquake (1993) in Saurashtra region, 31 to 51 years from the last earthquake (1969) in Mainland Gujarat and about 21 to 28 years from the last earthquake (2006) in Kachchh region.     

Evidence for a Tang-Song Dynasty great earthquake along the Longmen Shan Thrust Belt prior to the 2008 <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 7.9 Wenchuan earthquake,China

The magnitude (M _w) 7.9 Wenchuan earthquake occurred on 12 May 2008 in the Longmen Shan region of China, the transition zone between the Tibetan Plateau and the Sichuan Basin, resulting in widespread damage throughout central and western China. The steep, high-relief eastern margin of the Tibetan Plateau has undergone rapid Cenozoic uplift and denudation accompanied by folding and thrusting, yet no large thrust earthquakes are known prior to the 2008 M _w 7.9 Wenchuan earthquake. Field and excavation investigations reveal that a great historical earthquake occurred in the Sichuan region that ruptured a >200-km-long thrust fault within the Longmen Shan Thrust Belt, China, which also triggered the 2008 M _w 7.9 Wenchuan earthquake. The average co-seismic slip amount produced by this historical earthquake is estimated to be 2–3 m, comparable with that caused by the 2008 Wenchuan earthquake. Paleoseismic and archaeological evidence and radiocarbon dating results show that the penultimate great earthquake occurred in the Sichuan region during the late Tang-Song Dynasty, between AD 800 and 1000, suggesting a recurrence interval of ~1,000–1,200 years for Wenchuan-magnitude (M = ~8) earthquakes in the late Holocene within the Longmen Shan Thrust Belt. This finding is in contrast with previous estimates of 2,000–10,000 years for the recurrence interval of large earthquakes within the Longmen Shan Thrust Belt, as obtained from long-term slip rates based on the Global Positioning System and geological data, thereby necessitating substantial modifications to existing seismic-hazard models for the densely populated region at the eastern marginal zone of the Tibetan Plateau.     

Temporal distribution of earthquakes using renewal process in the Dasht-e-Bayaz region

Temporal distribution of earthquakes with M _w > 6 in the Dasht-e-Bayaz region, eastern Iran has been investigated using time-dependent models. Based on these types of models, it is assumed that the times between consecutive large earthquakes follow a certain statistical distribution. For this purpose, four time-dependent inter-event distributions including the Weibull, Gamma, Lognormal, and the Brownian Passage Time (BPT) are used in this study and the associated parameters are estimated using the method of maximum likelihood estimation. The suitable distribution is selected based on logarithm likelihood function and Bayesian Information Criterion. The probability of the occurrence of the next large earthquake during a specified interval of time was calculated for each model. Then, the concept of conditional probability has been applied to forecast the next major (M _w > 6) earthquake in the site of our interest. The emphasis is on statistical methods which attempt to quantify the probability of an earthquake occurring within a specified time, space, and magnitude windows. According to obtained results, the probability of occurrence of an earthquake with M _w > 6 in the near future is significantly high.     

Consistent CMT solutions from Harvard University before great earthquakes in Kurile Islands and its significance for earthquake prediction

Introduction In the development of seismic science,the improvement of observational instruments is prone to produce new parameters and promote rapid growth of seismology.For example,after the global digital seismograph network is established,some studies that could not be made or easily made in the past are now being carried out successively.As a result,analytical studies based on digital data have been made one after another and digital seismology has formed.Up to now,Harvard University has …     

Statistical testing of the relation between burst of aftershocks and strong earthquakes

According to geological tectonics and seismic activites this paper devided North China (30°–45°N, 105°–130°E) into four areas. We analyzed the North China earthquake catalogue from 1970 to 1986 (from 1965 to 1986 for Huabei, the North China, plain region) and identified forty-two bursts of aftershock. Seven of them occurred in aftershock regions of strong earthquakes and seventeen of them in the seismic swarm regions. The relation between strong earthquakes with the remaining eighteen bursts of aftershocks has been studied and tested statistically in this paper. The result of statistical testing show that the random probabilityp of coincidence of bursts of aftershock with subsequent strong earthquakes is less than six percent. By Xu’sR scoring method the efficacy of predicting strong earthquake from bursts of aftershock is estimated greater than 39 percent. Following the method proposed in the paper we analyzed the earthquake catalogue of China from 1987 to June, 1988. The results show that there was only one burst of aftershock occurred on Jan. 6, 1988 withM=3.6 in Xiuyan of Northeast China. It implicates that a potential earthquake withM _S⩽5 might occur in one year afterwards in the region of Northeast China. Actually on Feb. 25, 1988 an earthquake withM _S=5.3 occurred in Zhangwu of Northeast China. Another example is Datong-Yanggao shock on October 18, 1989 which is a burst of aftershock. Three hours after an expected shock withM =6.1 took place in the same area. Two examples above have been tested in practical prediction and this shows that bursts of aftershocks are significant in predicting strong earthquakes. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 273–280, 1991. Part of earthquake catalogue is from Jinbiao Chen, Peiyan Chen and Quanlin Li.     

On the relation of moderate-short term anomaly of earth resistivity to earthquake

IntroductionMany anomalies due to earthquake have been recorded in observation of earth-resistivity for30 years and over, which showed that there objectively existed the anomalies of each-resistivity.The crustal strUcture and medium conditions are quite complex, so the complexity of the temporal,spatial and intensive development of the anomalies is inevitable. Both of time and amplitUde ofanomalies among some stations near an epicenter are different (even among different observational directi…     

Diagnosis of Time of Increased Probability (TIP) for Volcanic Earthquakes at Mt. Vesuvius

The possibility of intermediate-term earthquake prediction at Mt. Vesuvius by means of the CN algorithm is explored. CN was originally designed to identify the Times of Increased Probability (TIPs) for the occurrence of strong tectonic earthquakes, with magnitude M ≥ M₀, within a region a priori delimited. Here the CN algorithm is applied, for the first time, to the analysis of volcanic seismicity. The earthquakes recorded at Mt. Vesuvius during the period from February 1972 to June 2004 are considered, and the magnitude threshold M₀ selecting the events to be predicted is varied within the range: 3.0–3.3. Satisfactory prediction results are obtained, by retrospective analysis, when a time scaling is introduced. In particular, when the length of the time windows is reduced by a factor 2.5–3, with respect to the standard version of CN algorithm, more than 90% of the events with M ≥ M₀ occur within the TIP intervals, with TIPs occupying about 30% of the total time considered. The control experiment ``Seismic History' demonstrates the stability of the obtained results and indicates that the CN algorithm can be applied to monitor the preparation of impending earthquakes with M ≥ 3.0 at Mt. Vesuvius.     

Occurrence probability of moderate to large earthquakes in Italy based on new geophysical methods

We develop new approaches to calculating 30-year probabilities for occurrence of moderate-to-large earthquakes in Italy. Geodetic techniques and finite-element modelling, aimed to reproduce a large amount of neotectonic data using thin-shell finite element, are used to separately calculate the expected seismicity rates inside seismogenic areas (polygons containing mapped faults and/or suspected or modelled faults). Thirty-year earthquake probabilities obtained from the two approaches show similarities in most of Italy: the largest probabilities are found in the southern Apennines, where they reach values between 10% and 20% for earthquakes of M _W ≥ 6.0, and lower than 10% for events with an M _W ≥ 6.5.