Fractal Patterns in Groundwater Radon Disturbances Prior to the Great 7.9 Mw Wenchuan Earthquake, China
Abstract
:1. Introduction
2. Materials and Methods
2.1. Experimental Aspects
2.1.1. Earthquake Activity
2.1.2. Measurement Setup
3. Mathematical Aspects
3.1. Fractal and Long Memory
3.2. Hurst Exponent
- (i)
- If , then the series has a positive long-range autocorrelation. A series’ high value is followed by a series’ low value, and vice versa. High Hurst exponents suggest persistent interactions that are predicted to occur in the series’ far future;
- (ii)
- If , then the low values follow high values in the time series, and vice versa. There is an ongoing exchange between low and high values for low H values in the time series’ future (this is known as anti-persistency);
- (iii)
- If associated processes are random, then the time series are totally uncorrelated.
3.3. Detrended Fluctuation Analysis (DFA)
3.3.1. Application of DFA
- (i)
- The initial time series is, first, integrated:The entire average value of the time series is denoted in Equation (1) by the symbol <…>, and k stands for the various time scales.
- (ii)
- The integrated time series, , is then separated into equal, non-overlapping bins of length n.
- (iii)
- The function that represents the trend in the bin is then fitted. Simple linear trends or polynomials of second order or higher order may be used. Here, the linear function is used. This linear function’s y coordinate is denoted by the notation in each box n.
- (iv)
- The local linear trend, , is then subtracted from the integrated time series, , which is detrended in each box of length n. The detrended time series, , is determined in this manner and for each bin as follows:
- (v)
- The integrated and detrended time series’ fluctuations’ root mean square (rms) is then computed for each bin of size n as follows:
- (vi)
- For various sizes of the scale boxes, the method steps (i)–(v) are repeated. This reveals the specific sort of connection between and n. If there are long-term relationships in the time series, then and n have an exponential relationship.The DFA scaling exponent of Equation (4) assesses the strength of the time series’ long-term relationships.
- (vii)
- A linear association between and is found via the logarithmic transformation of Equation (4). A strong linear connection suggests that the accompanying variations are long-lasting and, consequently, have a long memory. The square of Spearman’s () is used in this paper to measure the accuracy of the linear fit. Good linear fits are defined as having 0.95 or above [1,15,39,57,84].
3.3.2. Sliding Window DFA
- (a)
- The time series was segmented into windows of 64 samples. This segmentation approximately yields a two-month series’ part for the PZHS, SPS, GS, and MSS LSR stations, which record one measurement per day. The 64-sample window was also employed for fractal analysis of the data from three monitoring stations of urban air pollution with precisely the same measurement recording rate, namely, one measurement per day [39]. In a recent paper for the PZHS, a 256 segmentation DFA was employed [25], whereas, for radon in soil measurements, an approach of 128-sample window was utilised [85]. Nevertheless, since the windows are shifted 1 sample forward (sliding window technique), the whole signal is analysed, except from a 64-sample window, which was the final one. On the other hand, the 64-sample windowing yields a 64 h window for the HSR station of KDS, i.e., an analysis of about 2.5 days. Despite this differentiation, it is noteworthy that for a radon station in Pakistan, with the same recording rate as the one for KDS, a 64-window analysis was also utilised [16]. DFA from the data of KDS was analysed with 64 sample windows for consistency.
- (b)
- Every window was fitted using the least-squares fit of vs. in accordance with Equation (4). The data were fitted to a straight line without seeking cross-overs, as in the related literature [1,25,39], with the restriction that the slope of the fit displays a square of Spearman’s correlation coefficient above or equal to 0.95.
- (c)
- The window was advanced by one sample, and the steps (a) and (b) were repeated until the signal’s end.
- (d)
- DFA slopes, , were plotted against time, and the associated plot data were exported to ASCII output files for further use.
3.4. Fractal Dimension Analysis
3.4.1. Katz’s Method
3.4.2. Higuchi’s Method
3.4.3. Sevcik’s Method
3.4.4. Computational Methodology of Fractal Dimension
- (i)
- As in Section 3.3.2, the time series was segmented into windows of 64 samples. As mentioned, this segmentation approximately corresponds to, for the LSR stations (PZHS, SPS, GS, and MSS), a two-month signal. The 64-sample windowing was also employed in the fractal dimension calculation (with the same methods) from the data of the three LSR urban pollution stations with identical rates of measurements, i.e., one measurement per day [39]. As also mentioned in Section 3.3.2, for the HSR station KDS, the 64-sample segmentation corresponds to approximately 2.5 days. In a previous fractal dimension analysis (with the same methods), a 256-window approach was implemented for the PZHS [25]; however, in a very recent fractal dimension analysis with an identical methodology for an HSR radon station in Pakistan with the same rate of measurements as the one for KDS (one measurement per hour), a 64-window approach was utilised [16]. Finally, as in Section 3.3.2 and for consistency with the windowing of the other stations, a 64-sample window was chosen here as well for the KDS station.
- (ii)
- The fractal dimensions of each method were calculated as follows:
- Higuchi’s method: Equal to the slope, D, of the first-order least-squares fit of the logarithmic transformation of Equation (13), namely, the relation of versus , for . In the aforementioned analysis for the urban air pollution stations [39], , whereas in the analysis of radon in Pakistan [16] and of the electromagnetic disturbances of the Ileia station, Greece, the approach was used. Based on the two latter papers, was also selected here.
- Sevcik’s method: Equal to the Hausdorff dimension of Equation (16) () for N = 64, namely, equal to the number of samples in each window, which constitutes parameter L.
- (iii)
- Each window was forwarded one sample (sliding window technique), and procedures (i)–(ii) were iterated until the end of the time series.
- (iv)
- Time evolution plots of the fractal dimensions in accordance with the Katz’s, Higuchi’s, and Sevcik’s methods were generated, and the partial data were extracted to ASCII files for further use.
3.5. Power-Law Analysis
Computational Methodology of Power-Law Analysis
- (a)
- The time series was separated into 128-sample windows. This is a double window size in comparison to the other two methods. This is performed because power-law analysis does not work well with small-sized windows. For the LSR stations (PZHS, SPS, GS, and MSS), this segmentation corresponds to a 4-month signal and, roughly, a 5-day signal for the KDS. In previous publications, a 128-sample window was employed in parameter estimations [85] for recordings of similar recording rates, whilst in others, a 512-sample window [99] with a recording rate of one measurement every 10 min was employed.
- (b)
- The power spectrum, , based on the Morlet wavelet, as well as the central Morlet frequency, f, were calculated in each segment.
- (c)
- The parameters and were fitted via least squares. Exponents, , and power amplification, , were computed for every window under the constraint that Spearman’s .
- (d)
- Steps (a) through (c) were iterated to the end of the time series. At each iteration, the window was shifted one sample forwards. As with the other techniques, the whole time series was covered except the last window.
- (e)
- The and data were tabulated and saved in ASCII format for further use.
3.6. Further Issues
3.6.1. Formation of Analysis Classes
- (a)
- Class I: This class comprises windows that are associated with DFA least-squares log–log fits with a Spearman’s coefficient of and, simultaneously, a scaling exponent between , i.e., modelled by the fBm class [84].Class I segments:
- (b)
- Class II: This class includes time series segments with DFA’s (i.e., they do not adhere to the prominent fBm class) or (i.e., they adhere to the fractional Gaussian noise (fGn) class).Significantly, Class II segments:
- They are complements of Class I segments.
3.6.2. Comparisons of the Fractal Results
- (1)
- From DFA’s exponent as:
- (2)
- From fractal dimension (D) as:(Berry’s equation)
- (3)
- From power-law as:
3.7. Meta-Analysis
- (a)
- (Step 1): According to user-defined thresholds, each ASCII output results file is computationally scanned for out-of-threshold values. The ASCII files carrying the fractal dimension values are searched for under threshold values, whilst the ASCII files containing DFA’s exponents and the power-law values are searched for over threshold values. New ASCII step 1 files are generated that contain the out-of-threshold values.
- (b)
- (Step 2): Under the restriction that each segment’s first sample date is arbitrarily considered as the date of the whole segment, the step 1 ASCII files are computationally filtered to find areas with common dates. The above computational process results in the full coverage of all dates, except the one of the last window. The whole procedure is iterated over the results of all possible combinations of the following:
- DFA versus fractal analysis or versus at least two fractal dimension calculation techniques (six combinations);
- Fractal analysis versus at least two fractal dimension calculation techniques (four combinations);
- One fractal dimension calculation technique versus the other two (three combinations).
4. Results and Discussion
- If , then the associated time series is a temporal fractal and follows the Class I category:
- If , then the time series follows anti-persistent paths;
- If , then the time series follows persistent paths.
- If , then the time series is of low predictability and follows the Class II category:
- If , then the fluctuations in the related processes are not growing and, hence, a stationary system describes the series;
- If , then the underlying dynamics are random and the related system has no memory.
- (a)
- All over- or under-threshold results of all fractal methods (step 1, meta-analysis) for the KDS, MSS, and GS. The threshold results of each station are combined per 2, 3, 4, and 5 methods (step 2, meta-analysis; a total of 13 combinations) versus all 19 earthquakes of Table 2 and Figure 17.As mentioned in the previous paragraph, it is important not only to identify footprints using one or more techniques (already conducted here) but more importantly to link the different techniques focusing on similar aspects of the problem at hand. To achieve this:
- (1)
- The exact over- or under-threshold dates were located computationally from the fractal outputs of each station (step 1, meta-analysis). These dates are year, month, day, and hour for the HSR and KDS, and year, month, and day for the MSS and GS. This is conducted through a serial search.
- (2)
- The common threshold dates from all different techniques were found through an incremental computational search. The outputs used are from all methods and with up to 13 different combinations of these. All outputs were generated through special software and were stored in a computer for use.
- (3)
- The earthquake data from USGS [116] were transformed into an adequate ASCII file for the generation of the final plot.
- (b)
- Wherever the symbols of different methods coincide in time, this means that the signs of seismicity are provided by more than one method. If all 13 methods coincide, this means that the evidence is maximised. The more the techniques point to similar findings, the more rigid the evidence is. It should be emphasised to the reader that this coinciding is conducted on the step 1 results, that is, on the fractal outputs.
5. Conclusions
Author Contributions
Funding
Data Availability Statement
Conflicts of Interest
References
- Nikolopoulos, D.; Petraki, E.; Yannakopoulos, P.H.; Priniotakis, G.; Voyiatzis, I.; Cantzos, D. Long-lasting patterns in 3 kHz electromagnetic time series after the ML = 6.6 earthquake of 2018-10-25 near Zakynthos, Greece. Geosciences 2020, 10, 235. [Google Scholar] [CrossRef]
- Cicerone, R.; Ebel, J.; Britton, J. A systematic compilation of earthquake precursors. Tectonophysics 2009, 476, 371–396. [Google Scholar] [CrossRef]
- Hough, S. The Great Quake Debate: The Crusader, the Skeptic, and the Rise of Modern Seismology; University of Washington Press: Washigton, DC, USA, 2020. [Google Scholar]
- Hayakawa, M.; Hobara, Y. Current status of seismo-electromagnetics for short-term earthquake prediction. Geomat. Nat. Hazards Risk 2010, 1, 115–155. [Google Scholar] [CrossRef]
- Molchanov, O.A.; Hayakawa, M. Seismo-Electromagnetics and Related Phenomena: History and Latest Results; Number A8; Terrapub: Tokyo, Japan, 2008; p. 189. [Google Scholar]
- Ouzounov, D.; Pulinets, S.; Hattori, K.; Taylor, P. Pre-Earthquake Processes: A Multidisciplinary Approach to Earthquake Prediction Studies Posters; Wiley: Hoboken, NJ, USA, 2018; p. 384. [Google Scholar]
- Petraki, E.; Nikolopoulos, D.; Nomicos, C.; Stonham, J.; Cantzos, D.; Yannakopoulos, P.; Kottou, S. Electromagnetic Pre-earthquake Precursors: Mechanisms, Data and Models-A Review. J. Earth Sci. Clim. Chang. 2015, 6, 250. [Google Scholar] [CrossRef]
- Petraki, E.; Nikolopoulos, D.; Panagiotaras, D.; Cantzos, D.; Yannakopoulos, P.; Nomicos, C.; Stonham, J. Radon-222: A Potential Short-Term Earthquake Precursor. J. Earth Sci. Clim. Chang. 2015, 6, 282. [Google Scholar] [CrossRef]
- Uyeda, S.; Nagao, T.; Kamogawa, M. Short-term earthquake prediction: Current status of seismo-electromagnetics. Tectonophysics 2009, 470, 205–213. [Google Scholar] [CrossRef]
- Conti, L.; Picozza, P.; Sotgiu, A. A Critical Review of Ground Based Observations of Earthquake Precursors. Front. Earth Sci. 2021, 9, 676766. [Google Scholar] [CrossRef]
- Ghosh, D.; Deb, A.; Sengupta, R. Anomalous radon emission as precursor of earthquake. J. Appl. Geophys. 2009, 187, 245–258. [Google Scholar] [CrossRef]
- Liu, C.; Dong, P.; Zhu, B.; Shi, Y. Stress Shadow on the Southwest Portion of the Longmen Shan Fault Impacted the 2008 Wenchuan Earthquake Rupture. J. Geophys. Res. Solid Earth 2018, 123, 9963–9981. [Google Scholar] [CrossRef]
- Yin, Y.; Wang, F.; Sun, P. Landslide hazards triggered by the 2008 Wenchuan earthquake, Sichuan, China. Landslides 2009, 6, 139–152. [Google Scholar] [CrossRef]
- Xu, Y.; Koper, K.D.; Sufri, O.; Zhu, L.; Hutko, A.R. Rupture imaging of the Mw 7.9 12 May 2008 Wenchuan earthquake from back projection of teleseismic P waves. Geochem. Geophys. Geosystems 2009, 10, Q04006. [Google Scholar] [CrossRef]
- Nikolopoulos, D.; Matsoukas, C.; Yannakopoulos, P.H.; Petraki, E.; Cantzos, D.; Nomicos, C. Long-Memory and Fractal Trends in Variations of Environmental Radon in Soil: Results from Measurements in Lesvos Island in Greece. J. Earth Sci. Clim. Chang. 2018, 9, 1–11. [Google Scholar] [CrossRef]
- Rafique, M.; Iqbal, J.; Shah, S.A.A.; Alam, A.; Lone, K.J.; Barkat, A.; Shah, M.A.; Qureshi, S.A.; Nikolopoulos, D. On fractal dimensions of soil radon gas time series. J. Atmos. Sol.-Terr. Phys. 2022, 227, 105775. [Google Scholar] [CrossRef]
- Shi, Z.; Wang, G.; Wang, C.y.; Manga, M.; Liu, C. Comparison of hydrological responses to the Wenchuan and Lushan earthquakes. Earth Planet. Sci. Lett. 2014, 391, 193–200. [Google Scholar] [CrossRef]
- Yin, X.Z.; Chen, J.H.; Peng, Z.; Meng, X.; Liu, Q.Y.; Guo, B.; Li, S.C. Evolution and Distribution of the Early Aftershocks Following the 2008 Mw 7.9 Wenchuan Earthquake in Sichuan, China. J. Geophys. Res. Solid Earth 2018, 123, 7775–7790. [Google Scholar] [CrossRef]
- Audemard, F.; Azuma, T.; Baiocco, F.; Baize, S.; Blumetti, A.M.; Brustia, E.; Clague, J.; Comerci, V.; Esposito, E.; Guerrieri, L.; et al. Earthquake Environmental Effect for Seismic Hazard Assessment: The ESI Intensity Scale and the EEE Catalogue; ISPRA—Servizio Geologico d’Italia: Roma, Italy, 2015; Volume XCVII, pp. 5–8. [Google Scholar]
- Keller, E.A.; Pinter, N. Active Tectonics, Earthquakes, Uplift and Landscape, 2nd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2002; 362p. [Google Scholar]
- Zhou, X.; Du, J.; Chen, Z.; Cheng, J.; Tang, Y.; Yang, L.; Xie, C.; Cui, Y.; Liu, L.; Yi, L.; et al. Geochemistry of soil gas in the seismic fault zone produced by the Wenchuan Ms 8.0 earthquake, southwestern China. Geochem. Trans. 2010, 11, 5. [Google Scholar] [CrossRef] [PubMed]
- Liu, J.Y.; Chen, Y.I.; Chen, C.H.; Liu, C.Y.; Chen, C.Y.; Nishihashi, M.; Li, J.Z.; Xia, Y.Q.; Oyama, K.I.; Hattori, K.; et al. Seismoionospheric GPS total electron content anomalies observed before the 12 May 2008 Mw7.9 Wenchuan earthquake. J. Geophys. Res. Space Phys. 2009, 114. [Google Scholar] [CrossRef]
- Huayong, N.; Hua, G.; Yanchao, G.; Blumetti, A.M.; Comerci, V.; Di Manna, P.; Guerrieri, L.; Vittori, E. Comparison of Earthquake Environmental Effects and ESI intensities for recent seismic events in different tectonic settings: Sichuan (SW China) and Central Apennines (Italy). Eng. Geol. 2019, 258, 105149. [Google Scholar] [CrossRef]
- Ren, H.; Liu, Y.; Yang, D. A preliminary study of post-seismic effects of radon following the Ms 8.0 Wenchuan earthquake. Radiat. Meas. 2012, 47, 82–88. [Google Scholar] [CrossRef]
- Alam, A.; Wang, N.; Zhao, G.; Mehmood, T.; Nikolopoulos, D. Long-lasting patterns of radon in groundwater at Panzhihua, China: Results from DFA, fractal dimensions and residual radon concentration. Geochem. J. 2019, 53, 341–358. [Google Scholar] [CrossRef]
- Alam, A.; Wang, N.; Zhao, G.; Barkat, A. Implication of radon monitoring for earthquake surveillance using statistical techniques: A case study of Wenchuan earthquake. Geofluids 2020, 2020, 2429165. [Google Scholar] [CrossRef]
- Alam, A.; Wang, N.; Petraki, E.; Barkat, A.; Huang, F.; Shah, M.A.; Cantzos, D.; Priniotakis, G.; Yannakopoulos, P.H.; Papoutsidakis, M.; et al. Fluctuation Dynamics of Radon in Groundwater Prior to the Gansu Earthquake, China (22 July 2013: M s = 6.6): Investigation with DFA and MFDFA Methods. Pure Appl. Geophys. 2021, 178, 3375–3395. [Google Scholar] [CrossRef]
- Ma, T.; Wu, Z. Precursor-Like Anomalies prior to the 2008 Wenchuan Earthquake: A Critical-but-Constructive Review. Int. J. Geophys. 2012, 2012, 583097. [Google Scholar] [CrossRef]
- Mandelbrot, B.B.; Ness, J.W.V. Fractional Brownian motions, fractional noises and applications. J. Soc. Ind. Appl. Math 1968, 10, 422–437. [Google Scholar] [CrossRef]
- Morales, I.O.; Landa, O.; Fossion, R.; Frank, A. Scale invariance, self-similarity and critical behaviour in classical and quantum system. J. Phys. Conf. Ser. 2012, 380, 012020. [Google Scholar] [CrossRef]
- Musa, M.; Ibrahim, K. Existence of long memory in ozone time series. Sains Malays. 2012, 41, 1367–1376. [Google Scholar]
- Vadrevu, K.P. Fractal analysis revealed persistent correlations in long-term vegetation fire data in most South and Southeast Asian countries. Environ. Res. Commun. 2023, 5, 011001. [Google Scholar] [CrossRef]
- May, R.M. Simple mathematical models with very complicated dynamics. Nature 1976, 261, 459–467. [Google Scholar] [CrossRef]
- Sugihara, G.; May, R. Nonlinear forecasting as a way of distinguishing chaos from measurement error in time series. Nature 1990, 344, 734–741. [Google Scholar] [CrossRef]
- Liu, C.; Liang, J.; Li, Y.; Shi, K. Fractal analysis of impact of PM2.5 on surface O3 sensitivity regime based on field observations. Sci. Total Environ. 2023, 858, 160136. [Google Scholar] [CrossRef]
- Pastén, D.; Pavez-Orrego, C. Multifractal time evolution for intraplate earthquakes recorded in southern Norway during 1980–2021. Chaos Solitons Fractals 2023, 167, 113000. [Google Scholar] [CrossRef]
- Nikolopoulos, D.; Moustris, K.; Petraki, E.; Cantzos, D. Long-memory traces in PM _ 10 time series in Athens, Greece: Investigation through DFA and R/S analysis. Meteorol. Atmos. Phys. 2021, 133, 261–279. [Google Scholar] [CrossRef]
- Chelidze, T.; Matcharashvili, T.; Mepharidze, E.; Dovgal, N. Complexity in Geophysical Time Series of Strain/Fracture at Laboratory and Large Dam Scales: Review. Entropy 2023, 25, 467. [Google Scholar] [CrossRef]
- Nikolopoulos, D.; Moustris, K.; Petraki, E.; Koulougliotis, D.; Cantzos, D. Fractal and long-memory traces in PM10 time series in Athens, Greece. Environmnets 2019, 6, 29. [Google Scholar] [CrossRef]
- Hurst, H. Long term storage capacity of reservoirs. Trans. Am. Soc. Civ. Eng. 1951, 116, 770–808. [Google Scholar] [CrossRef]
- Hurst, H.; Black, R.; Simaiki, Y. Long-Term Storage: An Experimental Study; Constable: London, UK, 1965. [Google Scholar]
- Lopez, T.; Martınez-Gonzalez, C.; Manjarrez, J.; Plascencia, N.; Balankin, A. Fractal Analysis of EEG Signals in the Brain of Epileptic Rats, with and without Biocompatible Implanted Neuroreservoirs. AMM 2009, 15, 127–136. [Google Scholar] [CrossRef]
- Eftaxias, K.; Balasis, G.; Contoyiannis, Y.; Papadimitriou, C.; Kalimeri, M. Unfolding the procedure of characterizing recorded ultra low frequency, kHZ and MHz electromagnetic anomalies prior to the L’Aquila earthquake as pre-seismic ones—Part 2. Nat. Hazards Earth Syst. Sci. 2010, 10, 275–294. [Google Scholar] [CrossRef]
- Kilcik, A.; Anderson, C.; Rozelot, J.; Ye, H.; Sugihara, G.; Ozguc, A. Nonlinear Prediction of Solar Cycle 24. Astrophys. J. 2009, 693, 1173–1177. [Google Scholar] [CrossRef]
- Chattopadhyay, A.; Khondekar, M. An investigation of the relationship between the CME and the Geomagnetic Storm. Astron. Comput. 2023, 43, 100695. [Google Scholar] [CrossRef]
- Granero, M.S.; Segovia, J.T.; Perez, J.G. Some comments on Hurst exponent and the long memory processes on capital Markets. Phys. A Stat. Mech. Its Appl. 2008, 387, 5543–5551. [Google Scholar] [CrossRef]
- Musaev, A.; Makshanov, A.; Grigoriev, D. The Genesis of Uncertainty: Structural Analysis of Stochastic Chaos in Finance Markets. Complexity 2023, 2023, 1302220. [Google Scholar] [CrossRef]
- Pérez-Sienes, L.; Grande, M.; Losada, J.C.; Borondo, J. The Hurst Exponent as an Indicator to Anticipate Agricultural Commodity Prices. Entropy 2023, 25, 579. [Google Scholar] [CrossRef] [PubMed]
- Vogl, M. Hurst exponent dynamics of S&P 500 returns: Implications for market efficiency, long memory, multifractality and financial crises predictability by application of a nonlinear dynamics analysis framework. Chaos Solitons Fractals 2023, 166, 112884. [Google Scholar] [CrossRef]
- Dattatreya, G. Hurst Parameter Estimation from Noisy Observations of Data Traffic Traces. In Proceedings of the 4th WSEAS International Conference on Electronics, Control and Signal Processing, Miami, FL, USA, 17–19 November 2005; pp. 193–198. [Google Scholar]
- Wang, F.; Wang, H.; Zhou, X.; Fu, R. A Driving Fatigue Feature Detection Method Based on Multifractal Theory. IEEE Sens. J. 2022, 22, 19046–19059. [Google Scholar] [CrossRef]
- Zhou, H.; Chang, F. The long-memory temporal dependence of traffic crash fatality for different types of road users. Phys. A Stat. Mech. Appl. 2022, 607, 128210. [Google Scholar] [CrossRef]
- Li, X.; Polygiannakis, J.; Kapiris, P.; Peratzakis, A.; Eftaxias, K.; Yao, X. Fractal spectral analysis of pre-epileptic seizures in terms of criticality. J. Neural Eng. 2005, 2, 11–16. [Google Scholar] [CrossRef] [PubMed]
- Escobar-Ipuz, F.; Torres, A.; García-Jiménez, M.; Basar, C.; Cascón, J.; Mateo, J. Prediction of patients with idiopathic generalized epilepsy from healthy controls using machine learning from scalp EEG recordings. Brain Res. 2023, 1798, 148131. [Google Scholar] [CrossRef] [PubMed]
- Wijayanto, I.; Humairani, A.; Hadiyoso, S.; Rizal, A.; Prasanna, D.L.; Tripathi, S.L. Epileptic seizure detection on a compressed EEG signal using energy measurement. Biomed. Signal Process. Control 2023, 85, 104872. [Google Scholar] [CrossRef]
- Rehman, S.; Siddiqi, A. Wavelet based Hurst exponent and fractal dimensional analysis of Saudi climatic dynamics. Chaos Solitons Fractals 2009, 39, 1081–1090. [Google Scholar] [CrossRef]
- Petraki, E. Electromagnetic Radiation and Radon-222 Gas Emissions as Precursors of Seismic Activity. Ph.D. Thesis, Department of Electronic and Computer Engineering, Brunel University London, London, UK, 2016. [Google Scholar]
- Fujinawa, Y.; Takahashi, K. Electromagnetic radiations associated with major earthquakes. Phys. Earth Planet. Inter. 1998, 105, 249–259. [Google Scholar] [CrossRef]
- Hayakawa, M.; Ida, Y.; Gotoh, K. Multifractal analysis for the ULF geomagnetic data during the Guam earthquake. In Proceedings of the IEEE 6th International Symposium on Electromagnetic Compatibility and Electromagnetic Ecology, Saint Petersburg, Russia, 21–24 June 2005; pp. 21–24, 239–243. [Google Scholar]
- Hayakawa, M. VLF/LF radio sounding of ionospheric perturbations associated with earthquakes. Sensors 2007, 7, 1141–1158. [Google Scholar] [CrossRef]
- Nikolopoulos, D.; Alam, A.; Petraki, E.; Papoutsidakis, M.; Yannakopoulos, P.; Moustris, K.P. Stochastic and self-organisation patterns in a 17-year PM10 time series in Athens, Greece. Entropy 2021, 23, 307. [Google Scholar] [CrossRef] [PubMed]
- Skordas, E.S. On the increase of the “non-uniform” scaling of the magnetic field variations before the M(w)9.0 earthquake in Japan in 2011. Chaos 2014, 24, 023131. [Google Scholar] [CrossRef] [PubMed]
- Stanley, H.E. Powerlaws and universality. Nature 1995, 378, 597–600. [Google Scholar] [CrossRef]
- Sarlis, N.; Skordas, E.; Varotsos, P.; Nagao, T.; Kamogawa, M.; Tanaka, H.; Uyeda, S. Minimum of the order parameter fluctuations of seismicity before major earthquakes in Japan. Proc. Natl. Acad. Sci. USA 2013, 110, 13734–13738. [Google Scholar] [CrossRef]
- Becker, M.; Karpytchev, M.; Hu, A. Increased exposure of coastal cities to sea-level rise due to internal climate variability. Nat. Clim. Chang. 2023, 13, 367–374. [Google Scholar] [CrossRef]
- Ivanova, K.; Ausloos, M. Application of the detrended fluctuation analysis (DFA) method for describing cloud breaking. Phys. A Stat. Mech. Appl. 1999, 274, 349–354. [Google Scholar] [CrossRef]
- Koscielny-Bunde, E.; Bunde, A.; Havlin, S.; Roman, H.E.; Goldreich, Y.; Schellnhuberet, H. Indication of a Universal Persistence Law Governing Atmospheric Variability. Phys. Rev. Lett. 1998, 81, 729–732. [Google Scholar] [CrossRef]
- Rahmani, F.; Fattahi, M.H. Climate change-induced influences on the nonlinear dynamic patterns of precipitation and temperatures (case study: Central England). Theor. Appl. Climatol. 2023, 152, 1147–1158. [Google Scholar] [CrossRef]
- Linhares, R.R. Fractional poisson process: Long-range dependence in DNA sequences. Model Assist. Stat. Appl. 2023, 18, 33–43. [Google Scholar] [CrossRef]
- Peng, C.; Mietus, J.; Havlin, S.; Stanley, H.; Goldberger, A. Long-range anti-correlations and non-Gaussian behavior of the heartbeat. Phys. Rev. Lett. 1993, 70, 1343–1346. [Google Scholar] [CrossRef] [PubMed]
- Buldyrev, S.V.; Goldberger, A.L.; Havlin, S.; Mantegna, R.N.; Matsa, M.E.; Peng, C.K.; Simons, M.; Stanley, H.E. Long-range correlation properties of coding and noncoding DNA sequences: GenBank analysis. Phys. Rev. E Stat. Phys. Plasmas Fluids Relat. Interdiscip. Top. 1995, 51, 5084–5091. [Google Scholar] [CrossRef] [PubMed]
- Ivanov, P.C.; Rosenblum, M.G.; Peng, C.K.; Mietus, J.E.; Havlin, S.; Stanley, H.E.; Goldberger, A.L. Multifractality in human heartbeat dynamics. Nature 1999, 399, 461–465. [Google Scholar] [CrossRef]
- Mateo-March, M.; Moya-Ramón, M.; Javaloyes, A.; Sánchez-Muñoz, C.; Clemente-Suárez, V.J. Validity of detrended fluctuation analysis of heart rate variability to determine intensity thresholds in elite cyclists. Eur. J. Sport Sci. 2023, 23, 580–587. [Google Scholar] [CrossRef] [PubMed]
- Rogers, B.; Schaffarczyk, M.; Gronwald, T. Improved Estimation of Exercise Intensity Thresholds by Combining Dual Non-Invasive Biomarker Concepts: Correlation Properties of Heart Rate Variability and Respiratory Frequency. Sensors 2023, 23, 1973. [Google Scholar] [CrossRef] [PubMed]
- Eftaxias, K.; Balasis, G.; Contoyiannis, Y.; Papadimitriou, C.; Kalimeri, M.; Athanasopoulou, L.; Nikolopoulos, S.; Kopanas, J.; Antonopoulos, G.; Nomicos, C. Unfolding the procedure of characterizing recorded ultra low frequency, kHZ and MHz electromagnetic anomalies prior to the L’Aquila earthquake as pre-seismic ones-Part 1. Nat. Hazards Earth Syst. Sci. 2009, 9, 1953–1971. [Google Scholar] [CrossRef]
- Gotoh, K.; Hayakawa, M.; Smirnova, N.; Hattori, K. Fractal analysis of seismogenic ULF emissions. Phys. Chem. Earth 2004, 29, 419–424. [Google Scholar] [CrossRef]
- Hayakawa, M.; Ida, Y.; Gotoh, K. Fractal (mono- and multi-) analysis for the ULF data during the 1993 Guam earthquake for the study of prefracture criticality. Curr. Dev. Theory Appl. Wavelets 2008, 2, 159–174. [Google Scholar]
- Varotsos, P.; Sarlis, N.; Skordas, E. Natural Time Analysis: The New View of Time. Precursory Seismic Electric Signals, Earthquakes and Other Complex Time-Series; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
- Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Order Parameter and Entropy of Seismicity in Natural Time before Major Earthquakes: Recent Results. Geosciences 2022, 12, 225. [Google Scholar] [CrossRef]
- Hu, K.; Ivanov, P.; Chen, Z.; Carpena, P.; Stanley, H. Effect of trends on Detrended Fluctuation Analysis. Phys. Rev. E 2001, 64, 1–19. [Google Scholar] [CrossRef]
- Peng, C.; Buldyrev, S.; Simons, M.; Havlin, S.; Stanley, H.; Goldberger, A. On the mosaic organization of DNA sequences. Phys. Rev. E 1994, 49, 1685–1689. [Google Scholar] [CrossRef] [PubMed]
- Peng, C.; Havlin, S.; Stanley, H.; Goldberger, A. Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series. Chaos 1995, 5, 82–87. [Google Scholar] [CrossRef]
- Peng, C.; Hausdor, J.; Havlin, S.; Mietus, J.; Stanley, H.; Goldberger, A. Multiple-time scales analysis of physiological time series under neural control. Phys. A Stat. Mech. Appl. 1998, 249, 491–500. [Google Scholar] [CrossRef] [PubMed]
- Nikolopoulos, D.; Yannakopoulos, P.H.; Petraki, E.; Cantzos, D.; Nomicos, C. Long-Memory and Fractal Traces in kHz-MHz Electromagnetic Time Series Prior to the ML = 6.1, 12/6/2007 Lesvos, Greece Earthquake: Investigation through DFA and Time-Evolving Spectral Fractals. J. Earth Sci. Clim. Chang. 2018, 9, 1–15. [Google Scholar]
- Nikolopoulos, D.; Petraki, E.; Nomicos, C.; Koulouras, G.; Kottou, S.; Yannakopoulos, P.H. Long-Memory Trends in Disturbances of Radon in Soil Prior ML=5.1 Earthquakes of 17 November 2014 Greece. J. Earth Sci. Clim. Chang. 2015, 6, 244. [Google Scholar]
- Katz, M. Fractals and the analysis of waveforms. Comput. Biol. Med. 1988, 18, 145–156. [Google Scholar] [CrossRef] [PubMed]
- Raghavendra, B.; Dutt, D.N. Computing Fractal Dimension of Signals using Multiresolution Box-counting Method. Int. J. Electr. Comput. Energetic Electron. Commun. Eng. 2010, 4, 183–198. [Google Scholar]
- Higuchi, T. Approach to an irregular time series on basis of the fractal theory. Phys. D Nonlinear Phenom. 1988, 31, 277–283. [Google Scholar] [CrossRef]
- de la Torre, F.C.; Ramirez-Rojas, A.; Pavia-Miller, C.; Angulo-Brown, F.; Yepez, E.; Peralta, J. A comparison between spectral and fractal methods in electrotelluric time series. Rev. Mex. Fis. 1999, 45, 298–302. [Google Scholar]
- de la Torre, F.C.; Gonzaalez-Trejo, J.; Real-Ramírez, C.; Hoyos-Reyes, L. Fractal dimension algorithms and their application to time series associated with natural phenomena. J. Phys. Conf. Ser. 2013, 475, 1–10. [Google Scholar] [CrossRef]
- Sevcik, C. On fractal dimension of waveforms. Chaos Solit. Fract. 2006, 27, 579–580. [Google Scholar] [CrossRef]
- Cantzos, D.; Nikolopoulos, D.; Petraki, E.; Yannakopoulos, P.H.; Nomicos, C. Earthquake precursory signatures in electromagnetic radiation measurements in terms of day-to-day fractal spectral exponent variation: Analysis of the eastern Aegean 13/04/2017–20/07/2017 seismic activity. J. Seismol. 2018, 22, 1499–1513. [Google Scholar] [CrossRef]
- Ida, Y.; Hayakawa, M. Fractal analysis for the ULF data during the 1993 Guam earthquake to study prefracture criticality. Nonlin. Process. Geophys. 2012, 13, 409–412. [Google Scholar] [CrossRef]
- Ida, Y.; Yang, D.; Li, Q.; Sun, H.; Hayakawa, M. Fractal analysis of ULF electromagnetic emissions in possible association with earthquakes in China. Nonlin. Process. Geophys. 2012, 19, 577–583. [Google Scholar] [CrossRef]
- Smirnova, N.; Hayakawa, M. Fractal characteristics of the ground-observed ULF emissions in relation to geomagnetic and seismic activities. J. Atmos. Sol. Ter. Phy. 2007, 69, 1833–1841. [Google Scholar] [CrossRef]
- Yonaiguchi, N.; Ida, Y.; Hayakawa, M.; Masuda, S. Fractal analysis for VHF electromagnetic noises and the identification of preseismic signature of an earthquake. J. Atmos. Sol. Ter. Phy. 2007, 69, 1825–1832. [Google Scholar] [CrossRef]
- Eftaxias, K. Footprints of non-extensive Tsallis statistics, self-affinity and universality in the preparation of the L’Aquila earthquake hidden in a pre-seismic EM emission. Phys. A Stat. Mech. Appl. 2010, 389, 133–140. [Google Scholar] [CrossRef]
- Kapiris, P.; Peratzakis, J.P.A.; Nomikos, K.; Eftaxias, K. VHF-electromagnetic evidence of the underlying pre-seismic critical stage. Earth Planets Space 2002, 54, 1237–1246. [Google Scholar] [CrossRef]
- Nikolopoulos, D.; Petraki, E.; Marousaki, A.; Potirakis, S.; Koulouras, G.; Nomicos, C.; Panagiotaras, D.; Stonhamb, J.; Louizi, A. Environmental monitoring of radon in soil during a very seismically active period occurred in South West Greece. J. Environ. Monit. 2012, 14, 564–578. [Google Scholar] [CrossRef]
- Telesca, L.; Lasaponara, R. Vegetational patterns in burned and unburned areas investigated by using the detrended fluctuation analysis. Phys. A Stat. Mech. Appl. 2006, 368, 531–535. [Google Scholar] [CrossRef]
- Eftaxias, K.; Contoyiannis, Y.; Balasis, G.; Karamanos, K.; Kopanas, J.; Antonopoulos, G.; Koulouras, G.; Nomicos, C. Evidence of fractional-Brownian-motion-type asperity model for earthquake generation in candidate pre-seismic electromagnetic emissions. Nat. Hazard Earth Syst. 2008, 8, 657–669. [Google Scholar] [CrossRef]
- Varotsos, P.; Alexopoulos, K. Physical properties of the variations of the electric field of the earth preceding earthquakes, I. Tectonophysics 1984, 110, 73–98. [Google Scholar] [CrossRef]
- Varotsos, P.; Alexopoulos, K. Physical properties of the variations of the electric field of the earth preceding earthquakes, II. Tectonophysics 1984, 110, 99–125. [Google Scholar] [CrossRef]
- Varotsos, P.; Sarlis, N.; Lazaridou, M.B.N. Statistical evaluation of earthquake prediction results. Comments Success Rate Alarm. Rate Acta Geophys. Pol. 1996, 44, 329–347. [Google Scholar]
- Varotsos, P.; Sarlis, N.; Skordas, E. Magnetic field variations associated with SES. The instrumentation used for investigating their detectability. Proc. Jpn. Acad. Ser. B 2001, 77, 87–92. [Google Scholar] [CrossRef]
- Petraki, E.; Nikolopoulos, D.; Chaldeos, Y.; Koulouras, G.; Nomicos, C.; Yannakopoulos, P.H.; Kottou, S.; Stonham, J. Fractal evolution of MHz electromagnetic signals prior to earthquakes: Results collected in Greece during 2009. Geomat. Nat. Hazards Risk 2016, 7, 550–564. [Google Scholar] [CrossRef]
- Pinault, J.; Baubron, J. Signal processing of soil gas radon, atmospheric pressure, moisture, and soil temperature data: A new approach for radon concentration modeling. J. Geoph. Res. Sol. EA 1996, 101, 3157–3171. [Google Scholar] [CrossRef]
- Eftaxias, K.; Panin, V.; Deryugin, Y. Evolution-EM signals before earthquakes in terms of mesomechanics and complexity. Phys. Chem. Earth 2007, 29, 445–451. [Google Scholar] [CrossRef]
- Eftaxias, K.; Sgrigna, V.; Chelidze, T. Mechanical and electromagnetic phenomena accompanying preseismic deformation: From laboratory to geophysical scale. Tectonophysics 2007, 341, 1–5. [Google Scholar] [CrossRef]
- Gotoh, K.; Hayakawa, M.; Smirnova, N. Fractal analysis of the ULF geomagnetic data obtained at Izu Peninsula, Japan in relation to the nearby earthquake swarm of June–August 2000. Nat. Haz. Earth Syst. 2003, 3, 229–234. [Google Scholar] [CrossRef]
- Hayakawa, M.; Kawate, R.; Molchanov, O.; Yumoto, K. Results of ultra-low-frequency magnetic field measurements during the Guam earthquake of 8 August 1993. Geophys. Res. Lett. 1996, 23, 241–244. [Google Scholar] [CrossRef]
- Hayakawa, M.; Itoh, T.; Hattori, K.; Yumoto, K. ULF electromagnetic precursors for an earthquake at Biak, Indonesia on 17 February 1996. Geophys. Res. Lett. 2000, 27, 1531–1534. [Google Scholar] [CrossRef]
- Kapiris, P.; Eftaxias, K.; Nomikos, K.; Polygiannakis, J.; Dologlou, E.; Balasis, G.; Bogris, N.; Peratzakis, A.; Hadjicontis, V. Evolving towards a critical point: A possible electromagnetic way in which the critical regime is reached as the rupture approaches. Nonlinear Proc. Geoph. 2003, 10, 1–14. [Google Scholar] [CrossRef]
- Smirnova, N.; Hayakawa, M.; Gotoh, K. Precursory behavior of fractal characteristics of the ULF electromagnetic fields in seismic active zones before strong earthquakes. Phys. Chem. Earth 2004, 29, 445–451. [Google Scholar] [CrossRef]
- Smirnova, N.A.; Kiyashchenko, D.A.; Troyan, V.N.; Hayakawa, M. Multifractal Approach to Study the Earthquake Precursory Signatures Using the Ground-Based Observations. Rev. Appl. Phys. 2013, 2, 3. [Google Scholar]
- USGS. Available online: https://earthquake.usgs.gov/earthquakes/map/?extent=22.41103,91.51611&extent=35.56798,113.57666&range=search&baseLayer=terrain&timeZone=utc&search=%7B%22name%22:%22Search%20Results%22,%22params%22:%7B%22starttime%22:%222008-01-01%2000:00:00%22,%22endtime%22:%222009-01-01%2000:00:00%22,%22maxlatitude%22:36.844,%22minlatitude%22:22.999,%22maxlongitude%22:116.323,%22minlongitude%22:96.504,%22minmagnitude%22:5.5,%22orderby%22:%22time%22%7D%7D (accessed on 26 May 2023).
Station Code | Station Name | Latitude | Longitude | Distance (km) |
---|---|---|---|---|
KDS | Kangding station | 30.12 | 102.17 | 152.2 |
GS | Ganzi station | 31.61 | 100.01 | 325.5 |
MSS | Mingshan station | 30.1 | 103.1 | 105.6 |
PZHS | Panzhihua station | 26.51 | 101.74 | 526.0 |
SPS | Sonpan station | 32.65 | 103.6 | 182.5 |
i/i | Year | Month | Day | Hour | Minute | Second | Latitude | Longitude | Depth (km) | |
---|---|---|---|---|---|---|---|---|---|---|
1 | 2008 | 8 | 31 | 8 | 31 | 10 | 5.6 | 26.232 | 101.97 | 10 |
2 | 2008 | 8 | 30 | 8 | 30 | 53 | 6.0 | 26.241 | 101.889 | 11 |
3 | 2008 | 8 | 21 | 12 | 24 | 30 | 6.0 | 25.039 | 97.697 | 10 |
4 | 2008 | 8 | 5 | 9 | 49 | 17 | 6.0 | 32.756 | 105.494 | 6 |
5 | 2008 | 8 | 1 | 8 | 32 | 43 | 5.7 | 32.033 | 104.722 | 11 |
6 | 2008 | 7 | 24 | 9 | 30 | 9 | 5.7 | 32.747 | 105.542 | 10 |
7 | 2008 | 7 | 23 | 19 | 54 | 44 | 5.5 | 32.752 | 105.498 | 4 |
8 | 2008 | 5 | 27 | 8 | 37 | 51 | 5.7 | 32.71 | 105.54 | 10 |
9 | 2008 | 5 | 25 | 8 | 21 | 49 | 6.1 | 32.56 | 105.423 | 18 |
10 | 2008 | 5 | 17 | 8 | 25 | 48 | 5.8 | 32.24 | 104.982 | 9 |
11 | 2008 | 5 | 16 | 5 | 25 | 47 | 5.6 | 31.355 | 103.351 | 3 |
12 | 2008 | 5 | 13 | 7 | 7 | 8 | 5.8 | 30.89 | 103.194 | 9 |
13 | 2008 | 5 | 12 | 20 | 8 | 50 | 5.6 | 31.413 | 103.889 | 21.7 |
14 | 2008 | 5 | 12 | 11 | 11 | 2 | 6.1 | 31.214 | 103.618 | 10 |
15 | 2008 | 5 | 12 | 9 | 42 | 24 | 5.5 | 31.527 | 104.092 | 10 |
16 | 2008 | 5 | 12 | 6 | 43 | 14 | 5.8 | 31.211 | 103.715 | 10 |
17 | 2008 | 5 | 12 | 6 | 42 | 8 | 5.7 | 31.342 | 104.682 | 10 |
18 | 2008 | 5 | 12 | 6 | 61 | 56 | 5.7 | 31.586 | 104.032 | 10 |
19 | 2008 | 5 | 12 | 6 | 28 | 1 | 7.9 | 31.002 | 103.322 | 19 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Alam, A.; Nikolopoulos, D.; Wang, N. Fractal Patterns in Groundwater Radon Disturbances Prior to the Great 7.9 Mw Wenchuan Earthquake, China. Geosciences 2023, 13, 268. https://doi.org/10.3390/geosciences13090268
Alam A, Nikolopoulos D, Wang N. Fractal Patterns in Groundwater Radon Disturbances Prior to the Great 7.9 Mw Wenchuan Earthquake, China. Geosciences. 2023; 13(9):268. https://doi.org/10.3390/geosciences13090268
Chicago/Turabian StyleAlam, Aftab, Dimitrios Nikolopoulos, and Nanping Wang. 2023. "Fractal Patterns in Groundwater Radon Disturbances Prior to the Great 7.9 Mw Wenchuan Earthquake, China" Geosciences 13, no. 9: 268. https://doi.org/10.3390/geosciences13090268
APA StyleAlam, A., Nikolopoulos, D., & Wang, N. (2023). Fractal Patterns in Groundwater Radon Disturbances Prior to the Great 7.9 Mw Wenchuan Earthquake, China. Geosciences, 13(9), 268. https://doi.org/10.3390/geosciences13090268