Search

Earthquake forecasting is arguably one of the most challenging tasks in Earth sciences owing to the high complexity of the earthquake process. Over the past 40 years, there has been a plethora of work on finding credible, consistent and accurate earthquake precursors. This paper is a cumulative survey on earthquake precursor research, arranged into two broad categories: electromagnetic precursors and radon precursors. In the first category, methods related to measuring electromagnetic radiation in a wide frequency range, i.e., from a few Hz to several MHz, are presented. Precursors based on optical and radar imaging acquired by spaceborne sensors are also considered, in the broad sense, as electromagnetic. In the second category, concentration measurements of radon gas found in soil and air, or even in ground water after being dissolved, form the basis of radon activity precursors. Well-established mathematical techniques for analysing data derived from electromagnetic radiation and radon concentration measurements are also described with an emphasis on fractal methods. Finally, physical models of earthquake generation and propagation aiming at interpreting the foundation of the aforementioned seismic precursors, are investigated.

This paper reports kHz EM observations recorded by the Kardamas station in Ilia, Greece a few days before four near-field earthquakes occurred within a 24 km radius, with epicentral depths below 29 km. This work investigated the fractal features hidden in the EM observations via power-law analysis. All EM signals exhibited characteristic fractal epochs with organisation in space and time. A significant number of accurate fractal segments were delineated in the majority of the EM observations. A significant number of fractal areas corresponded to predictable Class I fBm category with 1 ≤ b ≤ 3 (0 ≤ H ≤ 1). Numerous persistent key-periods are reported with 2 < b ≤ 3 (0.5 ≤ H ≤ 1) which are deemed as signs of impeding earthquakes according to the literature. Numerous segments were found with strong persistent b-values in the range (2.3 ≤ b ≤ 3) (0.65 ≤ H ≤ 1) and b-values corresponding to switching between antipersistency and persistency with (1.7 ≤ b < 2.3) (0.35 ≤ H < 0.65). These are deemed as the most significant precursory signs. Interpretations are given via the asperity model.

The overall seismicity of Pakistan from 1820 to 2020 is analysed in terms of its multifractal behaviour. Seismic events of magnitude ML = 3.0 and above are spatially clustered into four distinct groups, each one corresponding to a different region of high seismic activity. The Multifractal Detrended Fluctuation Analysis (MFDA) method applied on each cluster reveals pronounced inter-cluster heterogeneity in terms of the resulting generalised Hurst exponent and fractality spectrum, possibly due to the particular tectonic characteristics of the regions under investigation. Additional results on the variability of the Gutenberg–Richter b-value across the defined clusters further corroborate the uniqueness of the seismic profile of each region.