The overall seismicity of Pakistan from 24 B.C 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.

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 space borne 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