Abstract
Time series analysis was applied to the continuous radon level, temperature, pressure, and rainfall to find clear earthquake signals. Radon signals appeared a few days after heavy rains, and radon signals associated with events M = 3.8–4.2 were detected 12 up to 36 days before. The events are complete data recorded from 1983 to 1986, giving discussion and conclusion on M with prediction time and radon anomaly detected in the Stivos faulting near Thessaloniki, N. Greece.
Similar content being viewed by others
Introduction
Greek seismicity is among the highest worldwide due to the positions of the tectonic plates crashing at the east Mediterranean, giving options for radon studies in soil and groundwaters as an earthquake indicator (Fig. 1a). These experiments have been performed in Greece since the 80s and continue using the advanced monitoring technology of the time. The results contribute to the comprehensive worldwide database obtained over a half-century derived from experiments. Active faults have been investigated in Japan, Himalayan-Altai and Caucasus regions, Italy, Spain, and France [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. In addition, the seismic fault of St. Andreas in California gave the first results after long-term radon monitoring and the models between earthquake magnitude and the time passed from the radon signal, followed by models based on Himalayan studies [15,16,17,18,19,20,21,22,23,24,25].
The Greek peninsula has mostly been concerned with earthquake prediction since 2000, and integrated and continuous radon monitoring has detected abnormal changes in subsurface radon concentrations before earthquakes. The radon studies began after the catastrophic earthquake M = 6.5 occurred on June 20, 1978, in the Thessaloniki area and M = 6.6 on February 24, 1981, in the Aigion area. The Stivos fault at the Lagadas basin, near the Thessaloniki metropolitan area, has been periodically monitored from 1983 to 2002 [26, 27]. The Alkyonides and Aigion Faults at the Gulf of Corinth near the urban area of Athens and Patras were also occasionally studied from 1995 until today [28,29,30]. Another study in the Pyrgos area of West Greece near the Gulf of Corinth [31], on Lesvos Island and Crete [32], and an experiment in the Ioannina area in Northwest Greece [33] were also performed.
This work reports on data of radon variation in soil gas at the Stivos fault in the Lagadas basin associated with meteorological parameters and with data of seismic events to find radon signals predicting earthquakes. A time series analysis was performed at the continuous radon level data, temperature, and pressure. Radon signals after heavy rains were recorded, while radon anomalies to earthquakes were also detected. The results complete the events recorded at this fault since 1983, giving conclusions for modelling M with prediction time.
Instrumentation
The study was performed in the Thessaloniki area, North Greece, especially at the Lagadas Basin, monitoring the Stivos faulting in the Serbomacedonian zone in North Greece (Fig. 1b). Three stations for radon monitoring were placed at Gerakarou, Stivos and Scholari creating a triangle. Soil gas radon has been monitored since 1983–1985 and 1996–2000 using solid-state nuclear track detectors in specially made plastic devices 30 cm in height (see Fig. 2a). The device was set inside an outer plastic tube 1 m long with a 7 cm diameter [27]. From 1999 to 2002, silicon diode detectors in the BARASOL device 50 cm in height (ALGADE, France, see Fig. 2b) replaced the SSNTD, having a functional area of 450 mm2 of Si with 100 μm depth, a radon detection limit of 50 Bq m−3 and a saturation activity of 3 MBq m−3. The space 50 cm above the device was filled with Styrofoam material in granule form to isolate the devices thermally. Atmospheric parameters, temperature (°C), barometric pressure (mbar) and rain precipitation (mm) were recorded simultaneously by the sensors also included in probes [4, 7]. The devices were calibrated using a radon chamber of 110 l with a leakage rate 10–6 h−1, two orders of magnitude lower than the 7.5 10–4 h−1 radon decay constant, and a Ra-Rn source of 20 kBq [34].
Results and discussion
The results of radon concentrations in soil gas carried out with BARASOL devices for August 1999 to December 2002 at the Gerakarou, Stivos, and Sholari stations are shown in Fig. 3. The gap in the data collected from all radon stations (December 9, 1999–February 14, 2000) was due to software updates due to the new millennium. After 2000, technical problems occurred in the device installed at the Gerakarou station, causing the gaps observed. In the area of the M = 6.5 event in 1978, the radon levels at Gerakarou and Sholari were 25 ± 5 kBq.m−3 and at the Stivos 50 ± 15 kBq m−3 higher enough than the uncertainties involved in the measurements.
From August to December 1999, integrated measurements using SSNTD and continuous by Baracol device were in the field, presenting a good coherence between them (Fig. 4). The impact of rainfall and earthquakes is also clear, especially in the Stivos station, while no significant peaks and drops appeared at Scholari and Gerakarou (Fig. 3). A clear radon peak appeared at Stivos station due to rainfall in November 1999 and reached its maximum of 80 kBq.m−3 during two days after rains. The radon increases in the device after heavy rainfall due to increased barometric pressure inside the porous ground around the device that causes radon to flow through the device, which is not affected by the rain. After the radon peak, following the most interesting result, a radon signal reached its minimum of < 10 kBq.m−3 on November 26. Then it returned to normal levels on December 5, 1999—an inverse behaviour followed by the pressure giving a significant increase when radon was decreased dramatically. One week after the signal's end, an M = 4.5 earthquake occurred on December 12, 1999, predicting 36 days from the first sign and 7 days after the end of the radon anomaly.
The radon peak due to rainfall was 2σ higher, while the radon drop detected before a seismic event was -3σ lower than the radon level standard deviation of 15 kBq.m−3. The effect of rainfall is shown in all stations, with radon peaks appearing during a few days after heavy rains. A rainfall radon peak + 2σ higher than the radon deviation of 5 kBq.m−3 appeared in Scholari station in April 2001 due to continuous weekly heavy rainfall. A peak was also detected in Stivos, but this time, the radon level was inside the radon deviation of 15 kBq.m−3. The above observation referred to any radon anomaly detected at Stivos station after 2000 instead of the clear anomalies detected during 1999. In the case of the Gerakarou station, no significant radon anomalies were also detected during 1999–2002 that could be associated with rainfall or earthquakes. At both stations, radon levels detected have variation lower than the radon standard deviation σ, in contrast to the Scholari data (Fig. 3).
Following the results after 2000, an inverse behaviour of radon level with temperature was observed at Gerakarou, while no influence appears with pressure, indicating that diffusion through soil porosity could be the dominant mechanism of soil gas transfer in a common fracture underground [16,17,18]. On the contrary, Stivos and Scholari's results show an inverse behaviour between radon and pressure, while temperature does not affect the radon level in soil gas, suggesting that advection through fractures and cracks in the underground is the primary process [15, 21]. Interestingly, at Sholari station, after the earthquakes occurred in the Autumn of 2001, the pressure and radon followed inverse behaviour, with pressure appearing similar to the Stivos station, where it was always detected. A possible origin of the lower variability of soil Rn concentration observed at Stivos station could be that the station is located at the junction of the two faulting directions presented in Fig. 1. In contrast, Gerakarou and Scholari stations are located at each one of the faulting directions, with the Scholari to be the closet station to the epicentre of the 1978 event.
However, it is more important and valuable to compare the results obtained from all three stations to evaluate the effects of earthquakes on soil radon concentration, but the lower than the standard deviations of radon variability in soil gas detected at Gerakarou and Stivos stations gives Scholari measurements the advantage for further time series analysis. The 15-min raw data collected are entered into a Python program and averaged daily to exclude the daily radon variation. After that, the daily data were averaged every 30 days (monthly). The program performed ANOVA analysis, and the difference between the daily and monthly results produced by ANOVA analysis was input into Mathematica to analyze the time series of radon, temperature and pressure curves, see Fig. 5. A fitting process was applied, defining the cosine curve as the standard, and the results were for \({\text{Radon}} = - 71.14 \cdot \cos \left( {1.06 \cdot t - 9.01} \right)\), for \({\text{Temperture}} = 0.015 \cdot \cos \left( { - 4.64 \cdot t + 1.01} \right)\) and for \({\text{Pressure}} = - 0.196 \cdot \cos \left( {0.99 \cdot t + 1.21} \right)\). The x-axis is expressed in days after the analysis started (14 February 2000). Significant variation in radon levels was observed during the Autumn of 2001 and 2002 when the pressure and temperature variations were the lowest (Fig. 5). Because radon levels at Sholari station seem to be influenced by the pressure, as illustrated in Fig. 3, the fit radon curve and pressure combined in a system with unknowns the time (t) and the radon level, which were estimated considering only the pressure known. The difference between the measured in the field and radon estimated values obtained by solving the system is presented in Fig. 5.
The figures show the signals of radon drops caused by precipitation and seismic events at the peaks. The signals related to earthquakes, since the variations were 2–3 σ, the radon standard deviation of 5 kBq −3 in the Scholari station. A long and complex radon signal appeared with a drop observed after rains in July 2001, followed by an increase of 10 kBq.m−3 before an M = 3.9 occurred on 10 August 2001 (543 days from 14 Feb 2000), predicting 17 days before. Then, another drop of 10 kBq m−3 due to rainfall in late July, followed by a significant boost to 15 kBq.m−3, appears in a period of low rainfall. The variation was detected 23 days before an event swan M = 4.0 and 4.2 happened on 8 October 2001 (602 days from 14 Feb 2000). Another short radon variation has a clear peak of 15 kBq.m−3 detected 12 days back to ML = 3.8 occurring on 8 Mars 2002 (753 days from 14 Feb 2000) (Fig. 6).
The association of the seismic events with the radon signals is based on the relation proposed between the magnitude M and the epicentral distance R (km). The Dobrovolsky (1979) equation, log R = 0.43 M, transformed by Hauksson and Goddard [35] as M = 2.4 log R – 0.43 and by Morelli [36] as M ≥ (2.3 ± 0.2) log (R) − (0.4 ± 0.3) gives the maximum possible distance, R, where a radon station should be located from the epicentral to detect precursory radon signals. Earthquakes with magnitudes M3–4 can be detected by stations between 20 and 70 km away, as in the Lagadas basin during the study, while for M4.5 -5.5, the distance increased to 90–300 km.
The results of continuous radon monitoring support the data obtained using integrated measurements with SSNTD in the previous years (Table 1). Then, radon peaked 30 and 45 days after two earthquakes, M = 4.7 (August 26, 1983) and M = 5.2 (February 19, 1984), with epicentral distances of 15 and 150 km, respectively [26]. Also, radon signal registration 30 days before an earthquake ML = 4.5 occurred on December 12, 1999 [27], the signal remained high, while a secondary radon peak appeared 15 days before an earthquake M = 3.4 occurred on January 31, 2000, and after returning to the normal level. The epicentral distance was around 20–30 km. The predicting time estimated for events lower than M = 4 agrees with data detected in the Ioannina area, where predicting time ranged from 4 to 15 days before earthquakes with M = 4–3, and the radon variations were 3 times higher than the normal [33].
The results collected over two decades from the Stivos area allow comparison with equations proposed between precursor time (tp, days) and earthquake magnitude (M). The proposed equations suggest log tp, = A + B*M, where A and B are statistically determined parameters. Based on that, Rikitake 1976 calculated log tp, = − 1.83 + 0.7 * M [22], although [25] proposed log tp, = 0.32 * M [25]. The data obtained from continuous radon monitoring show a higher predicting time than the integrated measurements for the same magnitude (see Table 1). If only the continuous radon data were fitted, then the Rikitake 1976 equation seems more applicable. However, the Virk equation is superior when it fits all the data obtained from this area (Fig. 7) since an equation log tp, = 0.36 * M was fitted to the experimental data with R2 = 89%, almost identical to the Virk 1995 proposal. It has to be considered that the equations proposed above have been based on earthquakes with M > 5 but are still applicable to data associated with events with M < 5.
The specific equations could be referred to as every seismic fault. Still, each fault has a particular personality structure. Many explanations have been proposed regarding the earthquake occurrence, as the crack-avalanche model [1] proposes an underground weathering formation due to accelerating cracking in the region of a future earthquake focal area, so stress corrosion in the rock matrix may be associated with abnormal behaviour of radon level. A similar proposition was published based on the dilatancy diffusion model; the earthquake forces crack the rocks, increasing the total porosity and radon entries into the porous, giving a signal under steady conditions. At the same time, pressure decreases in the entire preparation zone due to increased cracking and radon driven to the surface [9].
Based on this theory, another proposed equation includes, among the others, the relative shift detected in radon signal level ΔR = (Rp − Rb)/Rb, where Rp is the radon level of the signal, and Rb is the average radon level. Then, considering the radon lifetime, the following equation gives the prediction time before earthquake magnitude ML = 2 * log (ΔR/α*T) − 15.3, where α is statistically estimated, including the radon decay constant [23]. The seismic events at the Stivos faulting, North Greece, associated with radon signals detected by integrated and continuous radon monitoring, as presented in Table 1, with radon shift from 2 up to 6 times, demonstrate that the value α is around 9.1 × 10–11 d−1. The equation results are illustrated in Fig. 8. This value approaching 1 × 10–10 d−1 is one order of magnitude higher than that estimated using data from the Himalayas and reported as 1.8 × 10–11 d−1 [23], maybe due to different geomorphological structure and seismological behaviour between the two areas, one sited in a strong energy density (J m−3) region while the other in a lower one.
Founded on the Ramola et al. [23] proposed equation (Fig. 8), low radon signals were detected long before the event corresponds to low magnitude (< 4). At the same time, high radon signals occur shortly before a catastrophic earthquake. The Stivos data shows 2 to 6 times the average radon registered 12- 45 days before M = 3.8–5.2 events occurred. If a radon signal > 12 times the average detected, catastrophic earthquake M > 8 shall appear in < 1–3 days or severe M = 6–8 in 5–1 weeks, respectively. To resolve that issue and handle extremely high radon signals, long-term radon monitoring for a decade or more regarding improving this method. Moreover, more stations should be employed to improve statistics since radon varies through space and time significantly.
Conclusions
A time series analysis was conducted at the continuous radon level data, temperature, and pressure to find clear radon signals from rainfall and earthquakes. The background radon level was 25 ± 5 kBq m−3 at the Gerakarou and Sholari stations. A higher level of 50 ± 15 kBq m−3 appeared at Stivos, at the junction of faulting directions in the area. Clear radon peaks appeared due to 50–100 mm rainfall during a few days. A radon anomaly associated with events M 3.8–4.2 was detected 12 up to 36 days before. The data was complemented with previous results obtained from the same area. The seismic events were associated with the radon signal following that the epicentral distance from the radon station was in the frame of the well-adopted equations previously published. These and previous data of the prediction time (tp day) associated with earthquake M following the equation log T = 0.36*M, almost identical to the [25] proposal. Since the radon signal level ΔR considered, the [23] equation M = 2 * log (ΔR/α*T) − 15.3 fits the data well with α = 9.1 × 10–11. The Stivos data shows 2 to 6 times more than average radon levels registered 12–45 days before ML = 3.8–5.2 events occurred. However, more experimental data must be associated with earthquakes to reach safe conclusions.
References
Dobrovolsky IP, Zubkov SI, Miachkin VI (1979) Estimation of the size of earthquake preparation zones. Pure Appl Geophys 117(5):1025–1044
Iida T, Ikebe W (1986) A low-energy beta-particle imaging system for measuring tritium distributions. Nucl Inst Methods Phys Res A253:119–127
Khan HA, Tufail M, Qureshi AA (1990) Radon signals for earthquake prediction and geological prospection. Med J Islamic World Acad Sci 3:229–231
Monnin MM, Seidel JL (1992) Radon in soil and groundwater related to major geophysical events: a survey. Nucl Inst Methods Phys Res A 314:316–330
Igarashi G, Saeki S, Takahata N, Sumikawa K, Tasaka S, Sasaki Y, Takahashi M, Sano Y (1995) Groundwater radon anomaly before the Kobe earthquake in Japan. Science 269:60–61
Wyss M (1997) Cannot earthquakes be predicted? Science 278:487–488
Trique M, Richon P, Perrier F, Avouac JP, Sabroux JC (1999) Radon emanation and electric potential variations associated with transient deformation near reservoir lakes. Nature 399:137–141
Yang T, Walia V, Chyi L, Fu C, Chen C, Liu T, Song S, Lee C, Lee M (2005) Variations of soil radon and thoron concentrations in a fault zone and prospective earthquakes in SW Taiwan Rad. Measurem 40:96–502
Ghosh D, Deb A, Sengupta R, Bera S, Patra KK (2007) Pronounced soil radon anomaly—precursor of recent earthquakes in India. Radiat Meas 42:466–471
Ll F, Baixeras C, Moreno V, Bach J (2008) Soil radon levels across the Amer fault. Rad Meas 43:S319–S323
Tsunomori F, Shimodate T, Ide T, Tanaka H (2017) Radon concentration distributions in shallow and deep groundwater at the Tachikawa fault zone. J Environ Rad 172:106–112
De Luca G, Di Carlo G, Tallini MA (2018) Record changes in the Gran Sasso groundwater before, during and after the 2016 Amatrice earthquake in central Italy. Sci Rep 8:1–16
Barberio M, Gori F, Barbieri M et al (2020) New Observations in Central Italy of groundwater responses to the worldwide seismicity. Sci Rep 10:17850
Feng X, Zhong J, Yan R, Zhou Z, Tian L, Zhao J, Yuan Z (2022) Groundwater radon precursor anomalies identification by EMD-LSTM model. Water 14:69
Clements WE, Wilkening MH (1974) Atmospheric pressure effects on radon-222 transport across the earth–air interface. J Geophys Res 79:5025–5029
King CY (1978) Radon emanation on San Andreas fault. Nature 271:516–519
Mogro-Campero A, Fleischer R, Likes R (1980) Changes in subsurface radon concentration associated with earthquakes. J Geophys Res 85(B6):3053–3057
King CY (1980) Episodic radon changes in subsurface soil gas along active faults and possible relation to earthquakes. J Geophys Res 85:3065–3078
Fleischer RL (1981) Dislocation model for radon response to distant earthquakes. Geophys Res Lett 8(5):477–480
Fleischer RL (1983) Theory of alpha recoil effects on radon release and isotopic disequilibrium. Geochim Cosmochim Acta 47:778–784
Shapiro MH, Rice A, Mendenhall MH, Melvin D, Tombrello TA (1985) Recognition of environ-mentally caused variations in radon time series. Pure Appl Geophys 122:309–326
Rikitake T (1987) Earthquake prediction. Elsevier, Amsterdam
Ramola RC, Singh S, Virk HS (1988) A model for the correlation between radon anomalies and magnitude of earthquakes. Int J Radiat Appl Instrum Part D Nuclear Tracks Radiat Meas 15(1–4):689–769
Singh B, Singh S, Virk HS (1993) Earthquake prediction studies in Kangra Valley using plastic track recorders. Nucl Tracks Radiat Meas 22(1–4):459
Virk HS (1996) A critique of empirical scaling relationship between earthquake magnitude, epicentral distance and precursor time for interpretation of radon data. J Earthq Predict Res 5:574–583
Papastefanou C, Manolopoulou M, Savvides E, Charalambous S (1989) Radon monitoring at the Stivos Fault following the ML=6.5 earthquake at Thessaloniki, Greece, on June 20, 1978. Nucl Geophys 3:49–56
Papastefanou C, Manolopoulou M, Stoulos S, Ioannidou A, Gerasopoulos E (2001) Radon measurements along active faults in the Langadas Basin, northern Greece. Nat Hazards Earth Syst Sci 1:159–164
Richon P, Bernard P, Labed V, Sabroux LC, Beneïto A, Lucius D, Abbad S, Robe MC (2007) Results of monitoring 222Rn in soil gas of the Gulf of Corinth region, Greece. Radiat Meas 42:87–93
Katsanou K, Stratikopoulos K, Zagana E, Lambrakis N (2010) Radon changes along main faults in the broader Aigion region, NW Peloponnese. Bull Geol Soc Greece 43(4):1726–1736
Karastathis VK, Eleftheriou G, Kafatos M, Tsinganos K, Tselentis GA, Mouzakiotis E, Ouzounov D (2022) Observations on the stress-related variations of soil radon concentration in the Gulf of Corinth. Greece Sci Rep 12(1):5442
Nikolopoulos D, Petraki E, Marousaki A, Potirakis S, Koulouras G et al (2012) Environmental monitoring of radon in the soil during a very seismically active period occurred in Southwest Greece. J Environ Monit 14:564–578
Nevinskya I, Tsvetkova T, Dogrub M, Aksoyc E, Inceozc M et al (2018) Results of the simultaneous radon measurements around the Black Sea for seismological applications. J Environ Radiat 192:48–66
Papachristodoulou C, Stamoulis K, Ioannides K (2020) Temporal variation of soil gas radon associated with seismic activity: a case study in NW Greece. Pure Appl Geophys 177:821–836
Stoulos S, Manolopoulou M, Papastefanou C (2003) Assessment of natural radiation exposure and radon exhalation from building materials in Greece. J Environ Radioact 69:225–240
Hauksson E, Goddard JG (1981) Radon earthquake precursor studies in Iceland. J Geophys Res Solid Earth 86:7037–7054
Morelli D (2012) Radon as earthquake precursor. Earthquake research and analysis—statistical studies observations and planning. InTech, Rijeka
Acknowledgements
The authors would like to express their gratitude to Prof. Constantin Papastefanou for his continuous support during his duty of experiments involving radon signals in soil gas associated with earthquakes in the specific region.
Funding
Open access funding provided by HEAL-Link Greece.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Stylianos, S., Alexandra, I. Time-series analysis of radon monitoring in soil gas in association with earthquakes in Stivos faulting, at Lagadas basin, North Greece. J Radioanal Nucl Chem 332, 4581–4590 (2023). https://doi.org/10.1007/s10967-023-09170-2
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10967-023-09170-2