Progressive failure on the North Anatolian fault since 1939 by earthquake stress triggering

10 M ≥ 6.7 earthquakes ruptured 1000 km of the North Anatolian fault (Turkey) during 1939–1992, providing an unsurpassed opportunity to study how one large shock sets up the next. We use the mapped surface slip and fault geometry to infer the transfer of stress throughout the sequence. Calculations of the change in Coulomb failure stress reveal that nine out of 10 ruptures were brought closer to failure by the preceding shocks, typically by 1–10 bar, equivalent to 3–30 years of secular stressing. We translate the calculated stress changes into earthquake probability gains using an earthquake-nucleation constitutive relation, which includes both permanent and transient effects of the sudden stress changes. The transient effects of the stress changes dominate during the mean 10 yr period between triggering and subsequent rupturing shocks in the Anatolia sequence. The stress changes result in an average three-fold gain in the net earthquake probability during the decade after each event. Stress is calculated to be high today at several isolated sites along the fault. During the next 30 years, we estimate a 15 per cent probability of a M ≥ 6.7 earthquake east of the major eastern centre of Ercinzan, and a 12 per cent probability for a large event south of the major western port city of Izmit. Such stress-based probability calculations may thus be useful to assess and update earthquake hazards elsewhere.     

An improved region–time–length algorithm applied to the 1999 Chi-Chi, Taiwan earthquake

By means of the region–time–length (RTL) algorithm, which is widely used for investigating the precursory seismicity changes in China, Italy, Japan, Russia and Turkey, we examine the precursory seismic activity occurred prior to the 1999, M _w = 7.6, Chi-Chi earthquake around its epicentre. Based on our calculation of the RTL values, the epicentral area has been found to strongly exhibit the signature of anomalous activity, associated with the seismic quiescence and activation, before the main shock. Also proposed in this study is a helpful method for determining two important parameters used in the RTL analysis, the characteristic time and distance. Such method will largely reduce the ambiguity in the original RTL algorithm.     

Statistical analysis of strong-motion accelerograms and its application to earthquake early-warning systems

In view of increasing damage due to earthquakes, and the current problems of earthquake prediction, real-time warning of strong ground motion is attracting more interest. In principle, it allows short-term warning of earthquakes while they are occurring. With warning times of up to tens of seconds it is possible to send alerts to potential areas of strong shaking before the arrival of the seismic waves and to mitigate the damage, but only if the seismic source parameters are determined rapidly. The major problem of an early-warning system is the real-time estimation of the earthquake's size.
We investigated digitized strong-motion accelerograms from 244 earthquakes that occurred in North and Central America between 1940 and 1986 to find out whether their initial portions reflected the size of the ongoing earthquake. Applying conventional methods of time-series analyses we calculate appropriate signal parameters and describe their uncertainties in relation to the magnitude and epicentral distance. The study reveals that the magnitude of an earthquake can be predicted from the first second of a single accelerogram within ±1.36 magnitude units. The uncertainty can be reduced to about ±0.5 magnitude units if a larger number (≥8) of accelerograms are available, which requires a dense network of seismic stations in areas of high seismic risk.     

Seismic quiescence before the M 7, 1988, Spitak earthquake, Armenia

A detailed analysis of the 35  yr of seismicity between 1962 and 1997 using a gridding technique shows that the M 7, Spitak earthquake of 1988 December 7 was preceded by a quiescence anomaly that started at approximately 1984±0.5, and lasted about 5±0.5  yr, up to the main shock. This quiescence anomaly had a radius of about 20±3  km, estimated from circular areas with 75 per cent rate decrease, centred at the point of maximum significance of the anomaly. The quiescence was clearly present in the aftershock volume during the 5  yr before the 1988 main shock, but its statistically strongest expression was located 30  km NW of the epicentre. This anomaly fulfills the association rules between precursory quiescence anomalies and main shocks, even for a tight definition, and is therefore proposed as a case of precursory quiescence. The largest value of the standard deviate Z , found by random selection of samples by gridding, was Z =14 for a time window of T _w=3  yr, using a sample size of N =300 events. This makes this anomaly the strongest observed so far, and it is the first documented in an environment of continental collision. There are no false alarms exceeding in significance the precursor. The Armenian earthquake catalogue used for this study had 4600 earthquakes with M ≥ M _min=2.2 in the area bounded by 39.5° to 42°N/42.5° to 47°E. From the point of view of homogeneous reporting this is the best catalogue we have analysed so far. The limits of the data used and the density of the grid are dictated by the data, and have no influence on the results. The choice of free parameters does not influence the results significantly within the following limits: 100≤ N ≤500, 2≤ T _w≤7, 2.2≤ M _min≤2.8.