Testing earthquake prediction algorithms: Statistically significant advance prediction of the largest earthquakes in the Circum-Pacific, 1992-1997

V.G. Kossobokov; L.L. Romashkova; V. I. Keilis-Borok; J. H. Healy

doi:10.1016/s0031-9201(98)00159-9

Home/Archive/USGS

USGS · 70021777

Testing earthquake prediction algorithms: Statistically significant advance prediction of the largest earthquakes in the Circum-Pacific, 1992-1997

V.G. Kossobokov·L.L. Romashkova·V. I. Keilis-Borok·J. H. Healy

View original source ↗Export metadata ↓

Abstract

Algorithms M8 and MSc (i.e., the Mendocino Scenario) were used in a real-time intermediate-term research prediction of the strongest earthquakes in the Circum-Pacific seismic belt. Predictions are made by M8 first. Then, the areas of alarm are reduced by MSc at the cost that some earthquakes are missed in the second approximation of prediction. In 1992-1997, five earthquakes of magnitude 8 and above occurred in the test area: all of them were predicted by M8 and MSc identified correctly the locations of four of them. The space-time volume of the alarms is 36% and 18%, correspondingly, when estimated with a normalized product measure of empirical distribution of epicenters and uniform time. The statistical significance of the achieved results is beyond 99% both for M8 and MSc. For magnitude 7.5 + , 10 out of 19 earthquakes were predicted by M8 in 40% and five were predicted by M8-MSc in 13% of the total volume considered. This implies a significance level of 81% for M8 and 92% for M8-MSc. The lower significance levels might result from a global change in seismic regime in 1993-1996, when the rate of the largest events has doubled and all of them become exclusively normal or reversed faults. The predictions are fully reproducible; the algorithms M8 and MSc in complete formal definitions were published before we started our experiment [Keilis-Borok, V.I., Kossobokov, V.G., 1990. Premonitory activation of seismic flow: Algorithm M8, Phys. Earth and Planet. Inter. 61, 73-83; Kossobokov, V.G., Keilis-Borok, V.I., Smith, S.W., 1990. Localization of intermediate-term earthquake prediction, J. Geophys. Res., 95, 19763-19772; Healy, J.H., Kossobokov, V.G., Dewey, J.W., 1992. A test to evaluate the earthquake prediction algorithm, M8. U.S. Geol. Surv. OFR 92-401]. M8 is available from the IASPEI Software Library [Healy, J.H., Keilis-Borok, V.I., Lee, W.H.K. (Eds.), 1997. Algorithms for Earthquake Statistics and Prediction, Vol. 6. IASPEI Software Library]. ?? 1999 Elsevier Science B.V. All rights reserved.

Explore related subjects

Physics of the Earth and Planetary Interiors Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

V.G. Kossobokov, L.L. Romashkova, V. I. Keilis-Borok, J. H. Healy. 1999. Testing earthquake prediction algorithms: Statistically significant advance prediction of the largest earthquakes in the Circum-Pacific, 1992-1997. https://doi.org/10.1016/s0031-9201(98)00159-9

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2021

Electrical conductivity of the lithosphere-asthenosphere system

Electromagnetic geophysical methods image the electrical conductivity of the subsurface. Electrical conductivity is an intrinsic material property that is sensitive to temperature, composition, porosity, volatile and/or melt content, and other physical properties relevant to the solid Earth. Therefore, imaging the electrical structure of the crust and mantle yields valuable information on the physical and chemical state of the lithosphere-asthenosphere system. Here we explore the viability of the passive magnetotelluric (MT) method for constraining upper mantle properties. We approach this problem in four successive steps: 1) review the electrical conductivity behavior of relevant materials; 2) predict the bulk electrical conductivity structure of oceanic and continental lithosphere for a suite of representative physical states; 3) generate synthetic MT data from the conductivity predictions; 4) compare and discuss the conductivity predictions and the synthetic data with select case studies from oceanic and continental settings. Our aim is to clarify the uncertainties associated with drawing inferences from electrical conductivity observations and ultimately to provide a basis for assigning confidence levels to interpretations.

Physics of the Earth and Planetary Interiors↗

USGS2021

Slip distribution and rupture history of the August 11, 2012, double earthquakes in Ahar – Varzaghan, Iran, using joint inversion of teleseismic broadband and local strong motion data

We use combined teleseismic and strong motion data sets to investigate finite-fault slip models for a double of earthquakes that occurred on August 11, 2012, in northwestern Iran near the cities of Ahar and Varzaghan. The data include teleseismic P-waveforms retrieved from broadband seismic stations located between 30°–94° from the earthquakes and local strong motion data recorded by the Iran Strong Motion Network, installed and operated by the Building and Housing Research Centre. We first invert teleseismic P-waveforms and local strong motion data separately. For the first event (12:23 UTC), the teleseismic broadband inversion yields a somewhat deeper and simpler distribution of slip than the local strong motion inversion. The strong motion inversion results in a more complex distribution because of higher frequency content but can also be influenced by complexities in the propagation path. For the second event (12:34 UTC), the slip distribution from strong motion data is more similar to the teleseismic result and shows a simple slip area with a small relative movement to the west. To resolve the differences between the results of these two data sets and obtain a better constrained slip model, we perform a joint inversion of teleseismic broadband and local strong motion data. The joint inversion for the first event shows two asperities with a maximum slip of 3.9 m up- dip from the hypocenter and extending to the west between depths of 1 and 5 km. A second narrower high-slip area is seen just above the hypocenter from 6 to 10 km depth. The total moment for this earthquake is calculated to be M o = 3.8 × 10 25 dyn-cm (3.8 × 10 18 N.m) (M w 6.4). For the second event, the results of the joint inversion show a simple slip distribution that is mainly confined in a single patch around the hypocenter with a depth range from about 10 to 13 km and maximum slip of 1.9 m. We compute a total seismic moment of M o = 1.6 × 10 25 dyn-cm (1.6 × 10 18 N.m) (M w 6.1) for the second event. The largest stress drops for the first event occur above the hypocenter with an average stress drop over the rupture area of 120 bar (12 Mpa). For the second event, the maximum stress drop occurs at the reported focal depth with an average stress drop over the rupture area of 80 bar (8 Mpa).

East Anatolian Fault↗

USGS2018

The 1952 Kern County, California earthquake: A case study of issues in the analysis of historical intensity data for estimation of source parameters

Seismic intensity data based on first-hand accounts of shaking give valuable insight into historical and early instrumental earthquakes. Comparing an observed intensity distribution to intensity-prediction models based on modern calibration events allows the magnitude to be estimated for many historic earthquakes. Magnitude estimates can also potentially be refined for earthquakes for which limited instrumental data are available. However, the complicated nature of macroseismic data and the methods used to collect and interpret the data introduce significant uncertainties. In this paper, we illustrate these challenges and possible solutions using the 1952 Kern County, California, earthquake as a case study. Published estimates of its magnitude vary from M W 7.2–7.5, making it possibly the second largest in California during the 20th century. We considered over 1100 first-hand reports of shaking, supplemented with other data, and inferred the magnitude in several ways using intensity prediction equations, yielding a preferred intensity magnitude M I 7.2 ± 0.2, where the uncertainty reflects our judgement. The revised intensity distribution reveals stronger shaking on the hanging wall, south of the surface expression of the White Wolf fault, than on the footwall. Characterizing the magnitude and shaking distribution of this early instrumental earthquake can help improve estimation of the seismic hazard of the region. Such reinterpreted intensities for historic earthquakes, combined with U.S. Geological Survey (USGS) Did You Feel It? data for more recent events, can be used to produce a uniform shaking dataset with which earthquake hazard map performance can be assessed.

California↗