SUMMARY Earthquake prediction research has been conducted for over 100 years with no obvious successes. Claims of breakthroughs have failed to withstand scrutiny. Extensive searches have failed to find reliable precursors. Theoretical work suggests that faulting is a non-linear process which is highly sensitive to unmeasurably fine details of the state of the Earth in a large volume, not just in the immediate vicinity of the hypocentre. Any small earthquake thus has some probability of cascading into a large event. Reliable issuing of alarms of imminent large earthquakes appears to be effectively impossible. Key word: Earthquake prediction.

Earthquake prediction: a critical review

The scaling properties of earthquake populations show remarkable similarities to those observed at or near the critical point of other composite systems in statistical physics. This has led to the development of a variety of different physical models of seismogenesis as a critical phenomenon, involving locally nonlinear dynamics, with simplified rheologies exhibiting instability or avalanche-type behavior, in a material composed of a large number of discrete elements. In particular, it has been suggested that earthquakes are an example of a “self-organized critical phenomenon” analogous to a sandpile that spontaneously evolves to a critical angle of repose in response to the steady supply of new grains at the summit. In this stationary state of marginal stability the distribution of avalanche energies is a power law, equivalent to the Gutenberg-Richter frequency-magnitude law, and the behavior is relatively insensitive to the details of the dynamics. Here we review the results of some of the composite physical models that have been developed to simulate seismogenesis on different scales during (1) dynamic slip on a preexisting fault, (2) fault growth, and (3) fault nucleation. The individual physical models share some generic features, such as a dynamic energy flux applied by tectonic loading at a constant strain rate, strong local interactions, and fluctuations generated either dynamically or by fixed material heterogeneity, but they differ significantly in the details of the assumed dynamics and in the methods of numerical solution. However, all exhibit critical or near-critical behavior, with behavior quantitatively consistent with many of the observed fractal or multifractal scaling laws of brittle faulting and earthquakes, including the Gutenberg-Richter law. Some of the results are sensitive to the details of the dynamics and hence are not strict examples of self-organized criticality. Nevertheless, the results of these different physical models share some generic statistical properties similar to the “universal” behavior seen in a wide variety of critical phenomena, with significant implications for practical problems in probabilistic seismic hazard evaluation. In particular, the notion of self-organized criticality (or near-criticality) gives a scientific rationale for the a priori assumption of “stationarity” used as a first step in the prediction of the future level of hazard. The Gutenberg-Richter law (a power law in energy or seismic moment) is found to apply only within a finite scale range, both in model and natural seismicity. Accordingly, the frequency-magnitude distribution can be generalized to a gamma distribution in energy or seismic moment (a power law, with an exponential tail). This allows extrapolations of the frequency-magnitude distribution and the maximum credible magnitude to be constrained by observed seismic or tectonic moment release rates. The answers to other questions raised are less clear, for example, the effect of the a priori assumption of a Poisson process in a system with strong local interactions, and the impact of zoning a potentially multifractal distribution of epicentres with smooth polygons. The results of some models show premonitory patterns of seismicity which could in principle be used as mainshock precursors. However, there remains no consensus, on both theoretical and practical grounds, on the possibility or otherwise of reliable intermediate-term earthquake prediction.

Statistical physics, seismogenesis, and seismic hazard

Pre-earthquake signals: Underlying physical processes

https://ntrs.nasa.gov/api/citations/20010111483/downloads/20010111483.pdf

Charge Generation and Propagation in Igneous Rocks

Despite its extreme importance and years of effort, practical short-term earthquake prediction still remains to be seen. However, research in earthquake-related electromagnetic phenomena have recently shown that such phenomena make for a promising candidate for short-term earthquake prediction. There is a good deal of accumulated evidence of precursory signatures in a wide frequency range (DC-VHF). ULF geomagnetic change is one of the most promising phenomena and it suggests that short-term prediction is realizable. This report first reviews earlier observational facts and presents the latest results in detection of ULF emissions by means of sophisticated signal processing. Finally, we discuss the possible generation and propagation mechanism of earthquake-related ULF signals.

/pdf/ulf-geomagnetic-changes-associated-with-large-earthquakes-3ld4l0pbt7.pdf

ULF Geomagnetic Changes Associated with Large Earthquakes

New observations of magnetic, electric and electromagnetic field variations, possibly related to recent volcanic and seismic events, have been obtained on Mt. Unzen in Japan, Reunion Island in Indian Ocean, the Long Valley volcanic caldera in California, and for faults in China and Russia, California and several other locations. For volcanic events, contributions from different physical processes can be identified during the various eruption stages. Slow processes (weeks to months) include near-surface thermal demagnetization effects, piezomagnetic effects, and effects from rotation/displacement of magnetized material. Rapid processes (seconds to days) include piezomagnetic effects from instantaneous stress redistribution with explosive eruptions and electrokinetic effects from rupture of high pressure fluid compartments commonly encountered in volcanic regions. For seismic events, the observed coseismic offsets are instantaneous, provided care has been taken to ensure sensors are insensitive to seismic shaking and are in regions of low magnetic field gradient. Simple piezomagnetic dislocation models based on geodetically and seismically determined fault parameters generally match the observed signals in size and sign. Electrokinetic effects resulting from rupture of fluid filled compartments at hydrostatic to lithostatic pore pressures can generate transient signals in the frequency band 100 Hz to 0.01 Hz. However, large-scale fluid driven processes are not evident in near-field measurements in the epicentral region minutes to weeks before large earthquakes. The subset of ionospheric disturbances generated by trapped atmospheric pressure waves (also termed gravity waves and/or acoustic waves, traveling ionospheric disturbances or TID's) that are excited by earthquakes and volcanic eruptions are common and propagate to great distances. These are known and expected consequences of earthquakes, volcanic explosions (and other atmospheric disturbances), that must be identified and their effects removed from VLF/ULF electromagnetic field records before associating new observations of ionospheric disturbances with earthquake activity.

M. J. S. Johnston

Papers

Review of electric and magnetic fields accompanying seismic and volcanic activity