The Impact of Atmospheric and Tectonic Constraints on Radon-222 and Carbon Dioxide Flow in Geological Porous Media - A Dozen-Year Research Summary

Long-term monitoring of Rn-222 and CO 2 at a depth of several dozen meter at the Sde-Eliezer site, located within one of the Dead Sea Fault Zone segments in northern Israel, has led to the discovery of the clear phenomenon that both gases are affected by underground tectonic activity along the Dead Sea Fault Zone. It may relate to pre-seismic processes associated with the accumulation and relaxation of lithospheric stress and strain producing earthquakes. This approach assumes that meteorological in ﬂ uences on physico-chemical parameters are limited at depth since its strength diminishes with the increase of the overlay layer thickness. Hence, the monitoring of natural gases in deep boreholes above the water table enables to reduce the climatic-induced periodic contributions, and thus to identify the speci ﬁ c portion of the radon signals that could be related to regional tectonic pre-seismic activity. The plausible pre-seismic local movement of the two gases at depth - is identi ﬁ ed by the appearance of discrete, random, non-cyclical signals, wider in time duration than 20 h and clearly wider than the sum of the width of the periodic diurnal and semidiurnal signals driven by ambient meteorological parameters. These non-cyclical signals may precede, by one day or more, a forthcoming seismic event. Hence, it is plausible to conclude that monitoring of any other natural gas that is present at depth may show a similar broadening signal and may serve as a precursor too. The necessary technical conditions enabling to distinguish between anomalous signals of gases that may be induced locally by pre-seismic processes at depth, and the relatively low periodic signals that are still established at depth related to external climatic conditions, are presented in detail

The available techniques for the measurement of radon and its decay products (RDPs) in sub-ground media are based on the interaction of either alpha particles or gamma-rays with the sensor material.The most common means for the measurement of radon is through the detection of alpha particle emission from the decay of ,  and , or gamma ray radiation from radon decay products (RDPs) mainly Bi-214 and Pb-214 (e.g., Ivanovich and Harmon, 1982).The detection process includes recording of the alpha or gamma radiation for sequential time intervals (seconds to hours).Electronic automatic radon monitoring station are known to be operated in diverse places such as Japan (e.g., Noguchi and Wakita, 1977), China (e.g., King, 1986), and Europe (e.g., Koch and Heinicke, 1994).Most of these stations measure radon in water or soil using either Lucas cells (based on ZnS for detection of alpha particles) or NaI (for gamma detection).
In continuous long-term radon measurements, the mandatory essential condition is not to violate the natural thermodynamic equilibrium nor the radon and RDPs secular equilibrium, by pumping and grabbing radon samples from the geological porous media below the surface.The measurement technologies that allow working in passive mode are the alpha-particles sensors based on silicon photodiodes or ionization chambers, and gamma detectors based on various scintillation materials.It enables tracking temporal radon fluctuations under natural stable internal conditions, with a high-resolution sampling rate from once every few seconds up to several minutes.
For radon measurement by alpha detectors, the prime concealed problem is the ability of radon gas to move out from the solid grain to the intergranular space in the bedrock (emanation), and then to flow out to the air space (exhalation) in order to be perceived and measured by the detector.By taking into consideration that less than 32 percent of the radon content in rock emanate from the solid grains to the porous media (Hassan et al., 2009) and that the radon exhalation greatly depends on the features of the investigated media, such as its thickness, the internal structure as well as the moisture content of the solid phase, the amount of the radon gas being exhaled to the air and that may temporally be detected is less than 10%.
In addition, when an alpha detector is set inside a borehole, for example, it has to respond to radon particles that first have to enter by diffusion from the air space into the detector's sensing volume (dozens of cm 3 ).Then the alpha radiation which is emitted by the Rn-222 and its RDPs, Po-218 and Po-214, has to reach the internal silicon solid state sensor to produce counting signals.This process greatly reduces the efficiency of the measurement.
On the other side, there is an intrinsic distinction between the alpha and the gamma measurement, related to the sensing capabilities.When a gamma detector is installed, inside a borehole, within a geological formation, it responds to the radon existing within a surrounding spherical body of bedrocks of a radius of R + 35 cm (where R is the borehole radii and the 35 cm is the average gamma absorption length).In addition to the three-scale difference in the size of sensing volume, the gamma detection method enables principally the measurement of radon content directly inside the country rock media.Thus, it is the ideal passive method for measuring radon concentration under natural and undisturbed conditions, suitable for resolving small temporal changes provided that the gamma detectors have high detection sensitivity which is essential to distinguish between minor variations in amplitudes and time.
In light of the above, it is possible to summarize in general the important requirements to be taken into account in choosing radon detection systems for long-time monitoring as follows: • The ability to discern small changes in the intensity of measured signals high resolution capability with high signal to noise ratio • The ability to discern small changes in the variability of the signals at a high level of time separation, at a measurement rate of at least a few seconds high temporal resolution • Reliability in measurements over time-a stable measurement capacity for months to years • Automatic data collection and remote data transmission The rules listed above rule out the use for long-time radon monitoring of the techniques based on alpha particles detection by nuclear track films or active charcoal, or even radon extraction and grabbing from subsurface soil or rocks into ionization or Lucas cells, since the results are not continuous, the measurement times are too long, and the devices may reach saturation.Radon monitoring ionization chambers and silicon photodiodes-based alpha particle detectors that do not disturbed the local environments, and crystal scintillation-based gamma detectors are the most suitable sensors for passive long-term measurements.
If there is a priority to make measurements deeper than the shallow surface and monitor radon and other natural gases at a depth of tens of meters, it is necessary to adapt the detector to this configuration.
Finally, all studies based on radon monitoring are confronted with identifying the mechanisms that control radon temporal changes produced within the geological subsurface media.The models that were developed are based on the conventional assumptions that concentration and pressure gradients within the air-filled pores within soil porous media are the main guidingforces for radon flow by diffusion and advection.For this, a time dependent radon transport equation for the air-filled part of the porous media that combines generation, radioactive decay, diffusion, and advection of radon in porous media, was the best model developed so far (Rogers and Nielson, 1991;Van der Spoel et al., 1998).
This radon transport model for the air-filled part of the porous media as formulated by Van der Spoel (1998) is reproduced below: β is the partition-corrected porosity, Ca is the radon concentrations in the air-phase (Bq/m 3 ), D is the radon diffusion coefficient in air-filled pores (m 2 /s), k is the isotropic air-filled permeability (m 2 ), μ is the air-filled dynamic viscosity (kg/m 1 s 1 ), p is the air pressure (Pa) λ is the decay constant of radon (2.1/10 6 s 1 ), Φ ληρ b C Ra is the radon production rate from the immediate parent Ra-226 (Bq/m 3 s 1 ), ∇ • (D∇Ca) is the diffusive bulk flux density of radon which is proportional to the bulk radon concentration gradient ∇Ca, k μ (∇p • ∇Ca) is the radon advective flux density with velocity u, produced by a local gradient in pressure ∇p (deviation from the aerostatic absolute pressure) calculated from Darcy's law: ∇p − μ k u.But even though the existing model proved to be quite satisfactory for non-isothermal radon transport, it had one major flaw: physical equations do not include temperature, which is one of the most significant parameters affecting gas thermodynamics.
This work builds the experimental basis for extending the above model in order to include explicitly the connection between the temperature gradient and the radon movement.
The data base of the last dozen years will serve to further model radon variability by eventually determining an additional function F(Ca, ΔT, κ, E k , T, t) to effectively describe and quantify the radon thermal flux from the surface heated by the Sun, to the subsurface porous rock and vice versa.It will allow to determine the gas flux, its direction and range, as a function of: Ca-the concentration of radon in the porous media airspaces, ΔT-the temperature gradient, κ-the thermal conductivity of the gas, E k -the transitional kinetic energy of the radon atoms as a function of the temperature T (E k 3/2k B T), and t-the time.

The Goals of This Work
In this article we intend to summarize the research and achievements, carried out in 7 phases over the past dozen years, in connection with the following subjects: The radon and CO 2 temporal periodic variability patterns; Differentiating between the effect of atmospheric temperature and pressure on radon and CO 2 flow through a variety of sub terrains media; Eliminating the climatic-induced periodic contributions in order to identify within the measured time series the portion of the signals that could be associated with the regional geodynamic pre-seismic evolution.
The seven stages of work are described below.

Phase I
Response of subsurface Rn-222 to a 20-ton blast experiment investigated at shallow levels during seismic calibration explosion experiment, simulating a 2.6-M L earthquake, in order to investigate the influence of the explosive blast and the transitory seismic wave fields on the radon transport in the country rock, adjacent to the focus of the explosion (Zafrir et al., 2009).

Phase II
The analysis of the temporal behavior of radon and environmental parameters measured continuously at the Amram tunnel, near Eilat (Gulf of Aqaba).The analysis of the results over a period of several years shows radon variability at multiple time scales, associated with the air temperature outside the tunnel, specifically the temperature gradient between the external environment and the more stable environment inside the tunnel (Barbosa et al., 2010).

Phase III
Study of the suitability of silicon photodiodes alpha particles detectors, gamma scintillation sensors and ionization chamber detectors, for utilization in long-term radon monitoring in sub terrain geological media.A comparison of the efficiency and sensitivity, the capability to resolve signal to noise, background, stability, and reliability of their long-term measurements was presented (Zafrir et al., 2011).

Phase IV
Long-term high-resolution collection of radon time series carried out in the southern part of Israel.It shows that long-term radon monitoring based on simultaneous alpha and gamma measurement enables to differentiate between the impact of ambient temperature and pressure on radon transportation within porous media (Zafrir et al., 2013).

Phase V
Concurrent radon measurements by gamma and alpha detection systems in a deep drilling hole served as a proxy for studying radon movement within the shallow and deep subsurface, as well as for analyzing the effect of various environmental and tectonic parameters on the radon transport pattern.The technique enabled measuring, for the first time, the radon vertical velocity in the subsurface geologic media as 25 m/hr, and the recovery of a preseismic radon anomalous signal apparently associated with a regional geodynamic process (Zafrir et al., 2016).

Phase VI
Open deserted drillings and abandoned wells are an interacting interface between the subsurface hydrosphere, lithosphere, and biosphere to the atmosphere above it.The effect of atmospheric conditions, namely atmospheric pressure and temperature, on air, CO 2 , and Rn-222 transport across the borehole-ambient atmosphere interface, was investigating inside a 110 m deep by 1 m diameter borehole in northern Israel.The air exchange rate of these features was analyzed and, consequently their contribution as sources for greenhouse gas (GHG) emissions to the atmosphere was quantified (Levintal et al., 2020).

Phase VII
Long-term monitoring method of Rn-222 and CO 2 at a depth of several tens of meters at Sde-Eliezer site, located within one of the Dead Sea Fault Zone (DSFZ) segments in northern Israel, has led to the discovery of the phenomenon that both gases are affected by underground tectonic activity along the DSFZ, possible related to the pre-seismic processes associated with the accumulation and relaxation of lithospheric stress and strain producing earthquakes (Zafrir et al., 2019).
All the technological components used in each of the phases described above included the configuration of a system designed for multi-parametric measurement at a depth of tens of meters, which will be presented in the next section.

EXPERIMENTAL SETTING
Each measurement system installed in the above-mentioned sites included at least two radon gamma detectors and, in some sites, also radon alpha detectors.
Several types of gamma detectors have been adopted: standard (2″ diameter) NaI (Tl) gamma scintillation detectors controlled by total count electronic circuit as single channel analyzer (SCA), from two sources-PM-11 system (of Rotem Industries, Israel) and narrow BGO and NaI (Tl) (1.5″ diameter) scintillation detectors (of Scionix, Holland).
These BGO and NaI detectors were chosen after we ascertained, depending on the rules listed above, that they have high sensitivity, stability and suitability for use in a 2inch diameter exploration drilling holes to a depth of several dozen meters.We have also re-defined the energy range of the SCA used in the Scionix's gamma detectors to count gamma rays only between 475 and 3,000 keV, cancelling the measurement of gamma radiation in the field of the Compton scattering, between 50 and 475 keV, enabling to improve the ratio of signal to noise by 70% in comparison to the PM-11 detection systems (Figure 1A).
The basic principle to choose the SCA for total counting of RDPs gamma rays, is the assumption that radon is the only element within the porous media that can vary in time, under the influence of atmospheric variables or tectonic activity.This fluctuation is observed in addition to the natural gamma radiation from other solid radioactive elements as uranium, thorium and potassium that produce a constant local background (Figure 1B).This validates the use of total counting by gamma detection systems equipped with a SCA for monitoring the temporal variation of radon in geological environment.The sensitivity of Scionix's gamma detectors is 0.01-0.02Bq/m 3 per 1 count/1 hr (see Zafrir et al., 2011, based on the comparative test with a calibrated AlphaGUARD detection system, of Saphymo GmbH, Germany).
For radon detection by alpha sensing, the Barasol, BT45 N (Algade Inc., France) is a suitable instrument since it was designed to be used in difficult environments and the calibration of the sensor enables to calculate the volume activity of the Rn-222 in the measured airspace, under undisturbed natural condition.The sensor is an implanted silicon detector with a 400 mm 2 of sensitive area.It enables the counting of Rn-222 and its daughter alpha products by spectrometry of the alpha particles (with energy between 1.5 and 6 MeV) created in the sensing volume.The sensitivity is about 50 Bq/m 3 per 1 pulse/1 hr.
For CO 2 monitoring, the detectors that were lowered down for measurements at depth, were the GMD-20 and GMP-252 infrared gas analyzers of Vaisala, Finland, with measurement ranges of 0-2,000, 0-3,000, and 0-10,000 ppm according to the borehole depth.CO 2 sensor accuracy varied between ±40 ppm with a drift of ±60 ppm/year for the lower range, to ±2% and drift of ±150 ppm/year, for the 0-10,000 ppm range.
The meteorological variables are measured at each site and with the same high time resolution, as well as complementary parameters such as temperature and relative humidity at depth.

THE RESULTS AND THE ACHIEVEMENTS OF THE CONSECUTIVE STUDIES Overview
In Israel, the first radon electronic stations were installed in 1995, in the northwestern sector of the Dead Sea.The measurements, with high temporal resolution of 15 min and with adequate sensitivity, showed the presence of a high and temporally fluctuating radon flux in gravel and water (Vulkan et al., 1995).
A research activity based on cooperation between physicists, geophysics, and geologists has begun in Israel, in 2001, aiming to integrate multi-disciplinary scientific methods in Earth sciences in order to promote the understanding of the processes of earthquake formation and to allow their early prediction (Zafrir et al., 2003).
In the last dozen years, about seven combined and diversified studies have been conducted investigating the impact of temperature, pressure, tectonic and other constrains on radon and CO 2 flow within the subsurface geological media, and they are presented in the following sections.

Phase I: Response of Subsurface Rn-222 to a 20-Ton Seismic Calibration Experiment
An underground 20-ton explosion for the purpose of performing spatial calibration of a seismic systems array in Israel (Gitterman et al., 2005) was utilized to investigate its effect on the radon movement in the shallow geological media, very close to the focal point of the explosion.A network of five radon detectors was setup, with the radon sensors installed a few meters apart along 150 m from the center of the explosion.All the detectors were lowered to a depth of 2 m, protected by a metal pipe and suitable cushioning (Figures 2A,B).
Measurements were initiated 4 days before the explosion and carried on for 9 days after the calibration blast with 15 min integration time (Figure 2C).A higher sampling rate of 10 s was inserted in use for few hours before and after the explosion (Figure 2D).
Although it was possible to assume that in the center of the explosion there might be local changes in the radon concentration due to the crushing of the surrounding bedrock, there was no significant change in the radon content along the measuring line of 150 m for the duration of the expansion of the blast wave along this distance, of about 200-400 milliseconds, under acceleration of 1,700 cm/s 2 (Zafrir et al., 2009).
The calculations showed that in this short period of less than half a second, even under a stronger explosion that would create an earthquake simulation even more powerful than the measured The Gamma spectra of Uranium, Thorium, and Potassium, measured using a 3″ × 3″ NaI scintillator.It displays the Bi-214 energy (RDP) lines of 609, 1,120, 1,764, and 2,204 keV in the gamma-ray spectra from magmatic rocks.The spectrum includes also energy lines of U-238 and Th-232 decay products and the isotope K-40.The Ba-133 line is due to an artificial source added for the stabilization of the electronic amplification of the MCA (Multi Channel Analyzer).(B) A 10-days measurement exhibiting the total daily gamma variations due to the RDPs, while the gamma level of the thorium decay products within the rock is stable.This demonstrates the use of total gamma measurement to track temporal variations of radon in a soil or country rock media.
M 2.6, the radon could not respond and move more than 10 cm.Therefore, no changes were measured in the radon concentration at each of the five detectors.
The original article from 2009 includes a calculation showing that around the gamma detector which was close to 17 m from an explosion of 20 tons of explosive material, developed a huge pressure of about 5 atm that could move radon up to 0.3 m per s.
[In comparison to pressure gradient-controlled variations between 10-100 Pa that were sufficient to constrain advective movement within a sand column-Van der Spoel (1998)].However, the transition time in the explosion that was used for simulating earthquakes with magnitude of 2.6, was 0.25 s.Therefore, the time for the pressure impact on the radon gas, even if the EQ energy would be large in some order of magnitude (M 4.5), would still be too short to move the radon gas more than a few meters away.
Hence it is possible to conclude from this experiment that the effect in the subsurface radon concentration of a sudden pressure increase associated with an occurring earthquake is very reduced, as a result of the very short distance that radon is able to move in such a short time.This does not contradict the possibility of subsurface radon being affected by a preseismic event, as the pressure build-up associated with earthquake preparation processes typically develops over much longer time scales, enabling radon movement over larger distances.

Phase II: The Analysis of the Long-Term Temporal Variability of Radon in Research Tunnel
Continuous radon monitoring using co-located alpha and gamma detectors was performed over a period of more than 10 years at the 170 m research tunnel in Mt.Amram.Radon was measured under stable conditions more than 100 m from the tunnel's entrance and below a rock burden of about 100 m (see Figure 4, "Phase III: A Comparative Assessment of Radon Sensing Systems for Long-Term Monitoring in Sub Terrain Geological Media" section, Figure 5, and "Phase IV: Impact of Atmospheric Aspects on Radon Transportation Within Porous Media" section).This research tunnel in the Arava valley, which is part of the Negev desert, is particularly suitable for the investigation of the factors inducing radon variability, as some of the effects known to influence radon mobility, such as rainfall and soil moisture, are negligible in such a dry (<20 mm/year) environment.The natural environmental conditions in the tunnel are constant, specifically indoor humidity and temperature are kept stable throughout the year.In addition to the radon and environmental parameters which were measured inside the tunnel, air temperature and barometric pressure were also measured outside the tunnel.
The time series resulting from the continuous radon monitoring by both alpha and gamma systems display corresponding temporal variations (Figure 3A).Thus, the simultaneous use of alpha and gamma monitoring systems demonstrated that the detailed signals visible in both time series are genuine signals reflecting radon variability, as they are obtained from different detection systems based on distinct measurement principles.
The temporal variability is characterized by a clear seasonal pattern, with maximum values in summer.As expected, the temperature and atmospheric pressure measured outside the tunnel also exhibit a clear seasonal pattern, in contrast the temperature measured inside the tunnel, at a distance of 140 m from the tunnel's entrance, which is very stable along the year with a value of 28 ± 0.3 °C.The humidity is also fairly constant and close to 90%.
The temporal variability displayed by the time series of radon and environmental parameters covers multiple time scales, from less than 1 h to years.The decomposition of the time series using the maximal overlap discrete wavelet transform (Percival, 2008) allows to extract the dominant modes of variability and to isolate seasonal, intra-seasonal, and diurnal time scales.The results show that at all the time scales temperature is the main factor responsible for the observed temporal patterns (Barbosa et al., 2010).
At the seasonal scale the variability of radon and temperature is very similar, with both variables displaying a maximum in summer; however, a detailed analysis shows that the temperature outside the tunnel increases since mid-January, while radon concentrations remain low, only increasing slightly in mid-March and then faster by mid-May.However, the similarity between the two seasonal patterns becomes evident when analyzing in detail the difference between the temperature outside the tunnel and the constant temperature inside the tunnel (hereafter noted as temperature gradient), as displayed in Figure 3B, left).
The temporal variability of the temperature gradient also drives intra-seasonal variability.radon anomalies spanning multiple days are only observed when the temperature gradient is positive.This dependence of radon on the temperature gradient is non-linear, as the amplitude of the radon anomaly is not directly related to the amplitude of the temperature gradient, but dependent on the sign of the temperature gradient (Figure 3B,right).
As for the case of seasonal and intra-seasonal variability, the daily variations are closely related to the temperature gradient between the internal part of the tunnel and the exterior.Daily variations are absent in the winter, corresponding to a negative temperature gradient.When the external temperature rises and becomes higher than the temperature in the inner part of the tunnel, the radon time series display daily cycles with maximum at night.
Temperature influences natural ventilation as well as radon transport within the porous media.Thus, the mechanisms through which temperature can influence radon variability include temperature-driven air ventilation, through the air medium, and temperature-driven radon migration by thermal gradient dependent function-in the porous media (Nield and Bejam, 2006;Zafrir et al., 2007;Zafrir et al., 2008).

Phase III: A Comparative Assessment of Radon Sensing Systems for Long-Term Monitoring in Sub Terrain Geological Media
A comparative study between alpha particles radon detectors based on silicon photodiode or passive ionization cells, and gamma-rays radon detectors based on scintillation materials, were performed in order to determine their suitability for long-term monitoring in terms of sensitivity, resolution in the parameters of signals intensity and time, as well as long-term operating reliability (Zafrir et al., 2011).
The test system was installed in the research tunnel at Mt. Amram near Eilat, within a closed room at a depth of 140 m (Rn1 and Rn2 in Figure 4A).The constant internal temperature of 28 ± 0.3 °C in the tunnel's rooms have enabled the comparison of the various radon detectors under the same environmental conditions.The setup included: Five gamma-ray systems with different dimensions (noted A-E, Rotem, Israel, and Scionix, Holland), one silicon photodiode alpha detector (noted F, Algade, France) and one calibrated ionization chamber (noted G, Saphymo, Germany).All the detectors were exposed over several days to the same natural radioactivity levels under undisturbed thermodynamic equilibrium between the air space and the tunnel walls, recording count rates with high time resolution of every 15 or 30 min (Figure 4B).
A conversion scale that enable to evaluate the counting rate (CR) per 15 or 30 min of the various detectors sensors (A-F) at peaks in Bq/m 3 , is included in Figure 4 (Figure 4C).
The main result was that gamma detectors have a high sensitivity between 2 and 4 times in relation to the alpha detectors, enabling the detection of low sub-diurnal radon fluctuations.
In relation to the alpha detector, it was found that the statistical error during the measurements was higher than 20%, because of the low sensitivity (see Figure 4B), even in sampling intervals of 30 min [0.022 counts per 1 h for 1 Bq/m 3 by the Barasol (F)].Therefore, in order to identify variations with a required statistical accuracy of at least 5% (which can be obtained easily with each one of the gamma detectors), it is necessary to take measurements with this alpha detector in time periods of at least 8 h or even more.
There is another advantage of using a gamma detector for the measurement of radon temporal fluctuations within the geologic environment: gamma detectors do not violate environmental thermodynamic conditions or the radon secular equilibrium, as happens when grabbing or pumping radon samples.Therefore, the gamma detectors have the highest priority regarding the alternative detectors that were examined in this work.Furthermore, even mechanically it is also harder to handle and lower ionization chambers to typical depths of narrow exploration boreholes.

Phase IV: Impact of Atmospheric Aspects on Radon Transportation Within Porous Media
Long-term high-resolution Rn-222 time series, were collected in the subsurface of the southern area of Israel, at the Mt.Amram research tunnel near Eilat and in shallow and deep drillings at Gevanim site in Makhtesh Ramon (Steinitz and Piatibratova, 2009;Zafrir et al., 2013).
In general, if the radon measurements are performed at a tightly closed site, as in the Amram tunnel, the radon within the air space will be in undisturbed environmental conditions and will reach temporal thermodynamic equilibrium with the radon in the bedrocks as well as secular equilibrium with its RDPs, as encountered in the previous "Phase III: A Comparative Assessment of Radon Sensing Systems for Long-Term Monitoring in Sub Terrain Geological Media" section, Figure 4.
Then, the radon time series in the tunnel, exhibit daily and seasonal radon levels varying periodically and depending on the climatic temperature gradient outside the rock surface, 100 m above the horizontal tunnel (Figures 3B, 5B, 6A,B).The Fourier amplitude spectrum of temperature and of the "summer Radon" is greatly similar (Figures 6C,E), and the radon diurnal periodical variations are canceled during the winter (Figure 6F).
Hence, the tunnel measurement system reveals a radon signal that pushed down via the country rock porous media, by the daily heating of the surface at a height of 100 m above the tunnel.The heating engine is the gradient temperature ΔT a (down) pointing downward (Figure 5B).In a delay of about 10 h after the maximum exterior heating, the radon signal reaches the tunnel (Figure 6B).Ten hours is the time it takes for the gas to flow from the surface to the tunnel.When the surface heating is fading with the external temperature, the size of the temperature gradient ΔT a is reset and its direction turned upward ΔT a (up).Then also the radon daily signal fades over time.In winter, when the temperature inside the tunnel is stable and higher than the surface temperature, the interior temperature gradient ΔT i is always directed toward the surrounding rock, and inhibits the radon daily flow into the tunnel space (Figures 3B,  5B, 6A).
The temperature gradient is the only geophysical driver available to force the radon via the air-filled space of the surrounding bedrock, as it will be demonstrated also in the next phase ("Phase V: Differentiating Between Radon Anomalies Produced by Climatic Constituent and Pre Seismic Tectonic Activity" section).The alternative mechanism, pressure gradient on the surface, is too weak to compel radon to move along the pore's media within the bedrock to 100 m depth.
In contrast, the radon concentration in open boreholes in Makhtesh Ramon, as measured in the boreholes air space by the alpha particle's detectors at 1.2 and 85 m depth, suggests an association with the diurnal (24 h) and semidiurnal (12 h) periodic components and even with multiday barometric pressure variations (Figure 7B).The FFT spectra of the radon signals exhibits a similar component which is different from the temperature periodicity (Figures 7C, 6E).
The appearance of double daily signals like the one circled in drawing 7b is analogous to the pressure's half-day cycle S2 and S1 (Figure 7C) as presented by Dai and Wang (1999).Those peaks appear locally at 8-10 a.m. and 12-2 night (Steinitz and Piatibratova, 2009), and they demonstrate the connection within the open drilling between the barometric pressure and the radon movement in the inner space (see also "Phase VI: CO 2 and Rn-222 Emissions From an Abandoned Water Well Under the Influence of Climatic Pressure and Temperature" section and "Phase VII: Radon and CO 2 Monitoring Technique in Deep Subsurface, as a Proxy for Investigating Tectonic Pre-seismic Activity" section).
Accordingly, when the semi-daily pressure reaches the minimum (negative pressure in Figure 7B, bottom), the radon can flow into the drilling air space in order to achieve equilibrium.On the other side, when the semi-daily increase in pressure occurs (positive pressure gradient) it prevents the migration of radon from the country rock into the borehole air space.This mechanism was demonstrated by Van der Spoel (1998), showing that even pressure gradient controlled variations between 10 and 100 Pa were sufficient to constrain advective movement within a sand column.
At the multiday scale, the barometric pressure variations produce a long multi-day radon signals at a borehole depth of 85 m, with the same time span as in the shallow depth of 1.2 m (Figure 7B, top).But the radon daily and half a day periodicity that was observed at 1.2 m air space is concealed at the depth of 85 m most of the time except for few days when the pressure varies around 950 mbar level.
The understandings and knowledge that came out of this study, about utilizing in parallel two different types of radon measurement methods in closed and open-air spaces, led to the idea to operate these combined methods of radon measurement within deep abandoned boreholes as a proxy for underground geodynamic research.

Phase V: Differentiating Between Radon Anomalies Produced by Climatic Constituent and Pre Seismic Tectonic Activity
The combination of long-term, radon measurement by gamma and alpha methods within deep abandoned boreholes enabled to distinguish between radon anomalous signals at depth induced by the periodical climatic parameters (temperature and barometric pressure), and those that may be driven by deep tectonic processes along active seismic faults.The concurrent identification of radon signals in various depths with 50 m between them, allowed for the first time to calculate the radon vertical velocity in the bedrock under the influence of the temperature gradient on the surface.It also revealed a radon anomaly that apparently occurred ahead of a tectonic event, the Nuweiba earthquake, which occurred within the Dead Sea Fault Zone on June 27, 2015, with magnitude of M 5.5 (Zafrir et al., 2016).
The first deep dedicated combined radon monitoring system was installed in January 2015 at an abandoned water well (110 m depth and about 1 m diameter) at the Sde Eliezer site, which is located within the Hula Valley western border fault (Figure 8A).
Lowering simultaneously gamma and alpha radon detectors to different depths within the deep borehole (Figure 8B) and monitoring their temporal variations, enabled to differentiate in the radon time series signals, between the daily periodical influence of the ambient temperature (stars in Figure 9A) and the semi-daily effects of the barometric pressure (arrows in Figure 9A).
The gamma detectors at 10 and 60 m depth at Sde-Eliezer site (Figures 8B, 9A,B) are narrow BGO detectors, presented in "Experimental Setting" section and "Phase III: A Comparative Assessment of Radon Sensing Systems for Long-Term Monitoring in Sub Terrain Geological Media" section.These detectors measure contributions from a source volume of 1.5 m 3 of surrounding bed rock (Figure 8B).
According to Figure 9A the maximum signal in the raw data for the BGO gamma sensor at the depth of 60 m (for example, in blue in 27/4/15), is about 53,000 counts per 15 min minus the background of 4,000 counts per 15 min.It is thus equal to a net rate of 49,000 counts per 15 min, which is equivalent to 16,333 Bq/m 3 according to the above calibration scale (Figure 4c).This peak is reached, with a delay of about 2 h from the peak of the radon in the borehole airspace, as measured by the alpha sensor at 40 m (see Figure 9C).This radon peak of 125 counts per 15 min, with a background of five counts, corresponds, according to the stated calibration scale (Barasol sensor) to 120/0.0186,666 Bq/m 3 .This fact proves that the concentration of the radon within the sealed iron casing airspace (see Figure 8B) corresponds to less than 40% of the radon content that is measured by the gamma detector.Thus more than half of the gamma counts result from the contribution of the radon in the surrounding bedrock in the form of a background from the radon that remains stuck within the rock grains, and a variable part from the radon able to move within the rock porous media, mainly driven by the influence of the Sun heating on the surface creating a temperature gradient ("Phase II: The Analysis of the Long-Term Temporal variability of Radon in Research Tunnel" section, and "Phase IV: Impact of Atmospheric Aspects on Radon Transportation Within Porous Media" section).The background portion is a further proof that this contribution originates from outside the iron pipe, as it does not decrease to near zero as the radon signal within the casing airspace does under the effect of the positive barometric pressure gradient (Figures 9B,C).
The semi-daily cycle measured by the gamma and alpha detectors (Figure 9A) results from the radon which is pumped from the water into the well airspace by the barometric extraction (minima) at early morning, as measured by the alpha detector at 40 m.The secondary peak of the three radon signals appear at around the same hour, 4 a.m.(Figure 9C).The contribution from the mobile radon within the surrounding bedrock to the radiation measured by the gamma detectors is very low at that time.
In addition, the pronounced time delays between the radon peeks measured by the three sensors (Figure 9C), prove that if all the radon variations would happen in the same borehole airspace then we would have to look for fast streaming from the middle (the first peak at 40 m) upward to 10 m and then back downwards to reach the 60 m level 2 h later.It seems very unlikely as there is a single direction of CO 2 streaming, from 70 m level to 40, 20, and 5 m, and vice versa, as described in "Phase VI: CO 2 and Rn-222 Emissions From an Abandoned Water Well Under the Influence of Climatic Pressure and Temperature" section and in Figures 10B,D, with a faster air velocity between 2 and 6 m/min (120-360 m/h).Now, considering the radon temporal variation in the bedrock around the pipe.From 10 m downward, we would find after about 2 h a similar concentration peak of about several thousand Bq/m 3 at a depth of 60 m, high above the radon concentration within the borehole airspace that reduced at that moment to at least half level (Figure 9C).The temperature-driven downward movement of radon outside the sealed casing, brings radon into the saturation zone at 80 m.Since the borehole casing is not perforated down to the water table, radon can only enter the casing from the groundwater.Considering that Ra-226 was not found in the Sde Eliezer water, radon streamed downward into the groundwater via the porous rock, driven down by the surface temperature gradient ΔT.It means that the groundwater is the source for supplying radon to the open borehole air space.
In a shallow, open, abandoned wells such as the Sde-Eliezer well, long-term equilibrium between the radon within the porous media and inside the open borehole air space is not necessarily established.Therefore, the radon measured by the gamma and the alpha detectors is different.The negative barometric pressure gradient is the driving force that pumps the radon upward via the open borehole (as named Barometric Pumping in "Phase VI: CO 2 and Rn-222 Emissions From an Abandoned Water Well Under the Influence of Climatic Pressure and Temperature" section).(Zafrir et al., 2013).In this figure, three counts per 15 min of the gamma sensors peaks at 10, and 60, are equivalent to concentration of 1 Bq/m 3 within the surrounding country rock, and 0.018 counts per 15 min of the alpha particle detector peak is equivalent also to 1 Bq/m 3 .
The assumption that radon can flow downward via the porous media is also suggested from other measurements at diverse sites, including various time delay measurements (see "Discrimination Between the Influence of Climatic Temperature and Pressure on Radon Flow in Geological Media" section).
Analyzing the time series in detail, (Figure 9B) an interesting relationship between the radon signals and the gradients of both temperature and pressure parameters appears: the disappearance of the radon signals at some time intervals, as occurred on 22 March, 30 March, 10-14 April, or 23-24 Aapril.This phenomenon is related to the inversion of the temperature gradient ΔT T − T 10m , from positive to negative (when T 10m is the stable temperature in the well at a 10 m depth).This situation happens when the outside temperature becomes lower than the constant temperature of 21.29 °C at the 10 m depth (as described in "Phase II: The Analysis of the Long-Term Temporal variability of Radon in Research Tunnel" section and "Phase IV: Impact of Atmospheric Aspects on Radon Transportation Within Porous Media" section).It is also happened when the inversion from negative to positive occurred of the pressure gradient ΔP dP/dt 9h , from negative to positive occurred, leading to the same effect.The disappearance of radon daily signals because of the inversion in the temperature gradient was already observed (Barbosa et al., 2010;Choubey et al., 2011).
The whole process is very clear: the higher the external temperature, the more gas will flow down the bedrock.If more radon gas reaches the water, then more radon is sucked into the well airspace.At noon, because parallel to the rising of the FIGURE 10 | Analysis of the relationship between air velocity (V air ) and atmospheric pressure changes (dP atm /dt).(A) Experimental setting.(B) Demonstrates the method using a single inflow event.(C) CO 2 and radon responses to changes in barometric pressure (P atm ).Labels 1, 2, and 3 refer to the time when P atm begins to decrease, becomes stable (i.e., between points 2 and 3), and increases, respectively.(D) Plot displays V air as a function of dP atm /dt using 59 analyzed inflow events over three months (March-May 2018).Black arrow in (A) denotes the time interval between the minimum decrease value of CO 2 (dCO 2 /dt) at −5 and −70 m, from which V air was quantified.Red area represents the 95% confidence range according to the calculated linear regression (black line) (Levintal et al., 2020).
external temperature during the day, the barometric pressure is decreasing accordingly, and the barometric pumping is getting stronger.Then, more radon is measured by the alpha detectors at a depth of 40 m (Figure 9A).
The time lag of about 2 h between the maximum values of the two gamma detectors separated vertically by 50 m, defines the radon vertical velocity as 25 m/h in the local subsurface porous media (Figure 9C).This is the first determination of the downward radon transportation velocity within porous rock media, by direct field measurement.It also proves that radon in the Amram tunnel moves downwards to the tunnel (100 m depth) in 10 h (Figure 6B).The most important observation in the study of the radon time series obtained during the first 10 months of 2015 at the Sde Eliezer site, was the discovery of a broadened pre-seismic radon signal (wider than the sum of the expected width of the diurnal and semidiurnal radon signals combined) that have been ahead by one day (Figure 12 at "Phase VII: Radon and CO 2 Monitoring Technique in Deep Subsurface, as a Proxy for Investigating Tectonic Pre-seismic Activity" section) the Nuweiba M 5.5, 4 km deep earthquake, which happened on June 27, 2015, near Eilat (GE) about 480 km from Sde-Eliezer site at Hula Valley.It was the first evidence that a pre-seismic earth strain, stress or deformation along an active fault like the DSFZ, may produce a faraway local pressure gradient (Aharonov and Scholtz 2018) to induce, non-periodic, radon anomalies.
Since no experimental method proved or contradicted, the assumption that a long-lasting stress or strain produced during pre-seismic rupture processes can exist and operate along active faults, far away to hundreds of kilometers, we can raise the hypothesis that a pressure gradient in the magnitude of a few pascal (10 -5 atm) can be produced by pre-seismic tectonic driving forces in the subsurface, and induce radon and CO 2 anomalous flow within the porous media, even for strains below the tidal magnitude of ∼10 −8 .Models to assess the total strain that can develop in geological media before seismic events, to a distance of a few hundred kilometers from the center of the earthquake, as a function of its magnitude, were presented in the literature (Cicerone et al., 2009;Wang and Manga, 2010;Skelton et al., 2015;Woith et al., 2018).
Phase VI: CO 2 and Rn-222 Emissions From an Abandoned Water Well Under the Influence of Climatic Pressure and Temperature Boreholes and wells are a complex interplaying interface between the subsurface hydrosphere, lithosphere, and biosphere to the atmosphere above it.Quantitative computing the air exchange rate between these features is therefore very relevant, as these are potential sources for greenhouse gas (GHG) emissions to the atmosphere.The influence of atmospheric conditions, including barometric pressure and ambient temperature, on the air, CO 2 , and radon transport across the borehole-ambient atmosphere interface and groundwater-borehole bottom interface is studied based on measurements inside a 110 m deep by 1 m diameter borehole in Sde-Eliezer site in northern Israel (Levintal et al., 2020).
In parallel the setup included thermocouples (T type) and Relative Humidity probes (Humitter 50 U, also of Vaisala, Finland) at almost the same depths, and radon alpha-particle detector (Barasol of Algade Inc., France) as described previously ("Experimental setting" section, "Phase III: A Comparative Assessment of Radon Sensing Systems for Long-Term Monitoring in Sub Terrain Geological Media" section, and "Phase IV: Impact of Atmospheric Aspects on Radon Transportation Within Porous Media" section).
• The high resolution of CO 2 measurements enables to determine in-situ the air velocities within the borehole.
During barometric pumping (BP), when there is inflow of atmospheric air from the top to the bottom direction, the CO 2 measured by the upper sensor (−5 m) is expected to decrease before the CO 2 variation will be measured at the bottom (−70 m).Then the distance of 65 m between the two sensors can be divided by the time interval between sensor responses, in order to get the vertical component of the CO 2 velocity within the borehole (Vair), assuming that CO 2 transport is advective-carried as part of the air mass movement within the borehole.The basic assumption was based on the photographed proof that the borehole casing above the 80 m depth is a sealed iron pipe, and the only CO 2 inlet and outlet are the bottom or the top entries.Therefore, the minimum value of dCO 2 /dt provides an easily identifiable signal.The time difference between sensors' minimum value used to calculate the velocity of the atmospheric air pulse between sensors.This analysis is illustrated in Figures 10B-D.Preliminary calculations indicated that the average vertical CO 2 velocity inside the borehole equal ∼1.3 m/s.
Calculating the air exchange rate at the well from the measured air/CO 2 velocities, using two independent methods as detailed in Levintal et al., 2020, (fluctuating up to ∼6 m/min) gave a mass flux rate to the atmosphere of ∼5 g-CO 2 /min CO 2 emissions, becoming higher during the summer with two daily peaks as presented for example in Figure 15 at the next "Phase VII: Radon and CO 2 Monitoring Technique in Deep Subsurface, as a Proxy for Investigating Tectonic Pre-seismic Activity" section.This reveals the time-dependent behavior of boreholes as an additional source of CO 2 emissions.On a yearly scale, the calculated CO 2 emissions from Sde-Eliezer site were similar in magnitude to a wheat field in the range of 100 to 6,000 m 2 .
The most notable results are summarized as follow: • Barometric pressure changes were the governing mechanism for air movement along the well.The thermal internal instability played a minor role.• CO 2 , water vapor, and radon, moved along the air space in a similar manner.• The transport of gases from the vadose zone to the atmosphere via the borehole airspace depends on two processes: the transport of the pore gas into the borehole lower entry, and the movement of the gas out of the borehole and into the atmosphere.
Phase VII: Radon and CO 2 Monitoring Technique in Deep Subsurface, as a Proxy for Investigating Tectonic Pre-seismic Activity

Overview
Long-term monitoring method of Rn-222 and CO 2 at a depth of several tens of meters at Sde-Eliezer site, located within one of the DSFZ segments in northern Israel (Phase V: Differentiating Between Radon Anomalies Produced by Climatic Constituent and Pre Seismic Tectonic Activity), has led to the discovery of the phenomenon that both gases are affected by underground tectonic activity along the DSFZ, possible related to the preseismic processes associated with the accumulation and relaxation of lithospheric stress and strain producing earthquakes (Zafrir et al., 2019).
The Results of the First Period From February 2015 to May 2017 The first period of the long-term radon monitoring carried out at the Sde-Eliezer site lasted for two and a half years (844 days) from February 2015 to May 2017.In parallel, about 260 earthquakes with magnitude M ≥ 2.0, that have been recorded during the 844 days along the DSFZ area, was added to the acquired data (Figure 11A).The DSFZ area in Figure 11B is defined by the Geophysical Survey of Israel as the region from 27 °to 36 N in latitude and longitude from 32 °to 38 E (see https://earthquake.co.il/en/earthquake/searchEQS.php).The largest variations in the radon temporal flux are mostly controlled by the temperature gradient on the surface as can be observed in Figure 11A at different depths.
The role of the temperature gradient on radon temporal behavior is probably not linear since the ratio between the radon level in winter to summer varies by a factor of 3-10 (Figure 11A), while the temperature varies within 10% span (28 °change vs. 285K).In parallel, the pressure effects, according to the existing physical model, arise linearly through the negative pressure gradient term.Therefore, it is impossible to predict so far, the total impact of the temperature and pressure at particular time, by the current contemporary model.
Other climatic variables, such as rain and snow, were found unimportant according to the results between winter and spring (December to end of March) and the long Israeli summer (April to November) during which there are no visible signs of rain or snow.
In addition, three, very pronounced radon signals that were different in shape and preceded earthquakes with magnitude >4.5, appeared during disparate events which occurred during the above-mentioned two-and-a-half-year, and marked by a, b and c in Figure 11A as presented in the following section.

The Data Analysis Methods to Identify the Nature of the Gases Anomalies
The procedure to exclude periodical radon signals.
Since the existing radon model and the physical equations for radon transport in the subsurface geologic media do not include temperature, which is one of the most significant parameters affecting gas thermodynamics, it is difficult to predict the radon temporal behavior induced by this parameter as was done for the barometric pressure (Perrier and Girault, 2013).On the other hand, in order to identify the nature of the residual anomalies, their specific contribution must be excluded from the measured time series.
One of the possible methods for isolating the components belonging to the periodic signals produced by the climatic parameters in the measured time series is to apply the Fourier domain filtering and component isolation procedure (AutoSignal ™ , Systat Software Inc.).
Utilizing these existing algorithms that work in the frequency domain for low pass truncation, frequency components above a chosen ideal cutoff frequency be removed from the original radon time series, leaving the frequencies below it unaltered as well as any non-periodic discrete events.
Choosing a frequency cutoff less the diurnal one (1 cycle per day) excludes all the unwanted diurnal, semi diurnal and higher periodical signals (maybe harmonics) that are induced by the atmospheric parameters.
Another extraction method that relies on the Modulation Spectral Analysis is the spectrogram phase space representation in which one may see simultaneously the temporal as well as the spectral information of a given signal.
The main advantage of using a spectrogram is its capability of being a better visualization tool showing in a more emphasized manner the relevant temporal anomaly.
In any case, we have proved by the two methods, based on FFT low-pass truncation filtering and the time domain reconstruction analysis, that indeed a radon anomaly is identified and that it lasts for more than 20 h.
The Recovery of a Pre-seismic Radon Anomalous Signal Apparently Associated With a Regional Geodynamic Process These are the three disparate events that occurred during the above-mentioned two-and-a-half-year, in which very pronounced radon signals seem to preceded earthquakes with magnitude >4.5:   ,F.
Utilizing an existing algorithm that work in the frequency domain for low-pass truncation, frequency components above a chosen ideal cutoff frequency can be removed from the original radon time series, leaving the frequencies below it unaltered as well as any non-periodic discrete events.Choosing a frequency cutoff less the diurnal one (1 cycle per day) excludes all the unwanted diurnal, semidiurnal, and higher periodical signals (maybe harmonics) that are induced by the atmospheric parameters (e.g., Figures 12G,H).
The Second Period From May 2017 to End of 2018-Radon and CO 2 in Depth Overview In the middle of 2017, the setup at the Sde-Eliezer site was expanded to study the impact of atmospheric variables, (ambient pressure and temperature) on air flow within the deep dry borehole at the site-110 m deep and 1 m wide (Levintal et al., 2018).Temperature, relative humidity, CO 2 (at 5, 20, 40, and 70 m), and radon detectors (measured by gamma detectors at 10, 60, and 88 m and by alpha at 40 m) were placed along the cased boreholes (Figures 8B and 10A) and barometric pressure probe was installed at the surface.Since then to present, all borehole data has been logged at a high 30 s temporal resolution.
Beyond that initially anticipated, the results show that CO 2 in the open borehole air space, as measured by four CO 2 detectors at different depths, follows the radon (measured by alpha detector at 40 m) as well as the radon temporal variations at the surrounding bedrocks (measured by gamma detectors at 10, 60, and 88 m), and both are driven by the climatic parameters as analyzed for the radon in our previous sections.Therefore, the CO 2 reveals the same daily and semi daily temperature and pressure periodicities in the frequency time domain (Zafrir et al., 2019).Improving the sampling rate by 30 times (from 15 min to high 30 s resolution) and implementing four CO 2 detectors at different depths enabled measuring the CO 2 velocity in the Sde-Eliezer borehole under the influence of barometric pressure that seems to be more effective than temperature within the air space inside the open well (see "Phase VI: CO 2 and Rn-222 Emissions From an Abandoned Water Well Under the Influence of Climatic Pressure and Temperature" section).
Radon, CO 2 , and Pre-Seismic Processes-The 2018 Lake Kinneret Scenario In order to understand whether there is a link between earthquakes in the DSFZ and the gas signal measurements at the depth in the abandoned well, located within the northern part of the fault belonging to the DSFZ (WMF in Figure 13), we have chosen to focus on the emergence of 103 seismic events in Lake Kinneret (Sea of Galilee) that occurred in the summer of 2018.
The temporal variations of the radon measurements at Sde-Eliezer site (located on the WMF about 35 km north of Lake Kinneret) and the rest of the environmental parameters, acquired simultaneously with all the earthquakes that occurred in the Lake Kinneret, are presented in Figure 14.
Despite the situation that the number of tremors were over 100, only two of them had magnitude of M 4.2 and M 4.5, and all the rest had lower magnitude between M 2 and M 4. Therefore, there was no early expectation that tectonic activity in the subsurface would be followed by a pre-seismic radon anomalous signals, as occurred before (The Recovery of a Preseismic Radon Anomalous Signal Apparently Associated With a Regional Geodynamic Process section).Nonetheless, about nine non-periodic radon signals, broader than 20 h each, have been exposed during 75 consecutive day from 10 June to 20 August (Figure 14).
The encircled segments designated as a, b, c, d, e, f, g, h, i in Figure 14, occurred within the time series in which the signals of radon are exceeding the daily and semi daily periodic effect of atmospheric temperature and pressure as analyzed by the FFT time domain reconstruction with low-pass truncation filter and short time Fourier transform (STFT) methods (for example see Figure 15).The results of the FFT filtering and the time domain reconstruction of the time section between 10 June and 12 July pinpointed the two anomalous peaks (a)) and (b) in Figure 14, which started before the first rupture.It shows that two radon, clear, discreet, non-periodic and broadened with time duration that exceeds 20 h, on 12 and 22 June, preceding by 24 and 12 days the more powerful tremors of July 4, 2018, with magnitude M 4.5.
A selected example in Figure 15, shows the analysis of the signals of radon and CO 2 during 24 consecutive days at the Sde-Eliezer site, from 26 July to August 18, 2018 including (f), (g), (h)) and (i).The tremors in this time interval, fall below the magnitude of M 4, and the relations between them and the gases temporal variations become non-diagnostic.There is a mixture of earthquakes on one hand and radon and CO 2 broadening signals, on the other.The broadened signals could probably be related to geodynamic activity as we have shown in the cases of earthquakes with magnitude above 4.5 (Figures 12A-F) but it is not possible to associated them to a specific event.
At least, based on our previous investigation on "Response of radon in a seismic calibration explosion" ("Overview" section), we suggest that there will be no an apparent response in the concentration of radon below the surface as a direct effect of short-time transient excess pressure produced by any one of the swarm's earthquake, including the highest one of M 4.5 on 4 July (Figure 14).It is also possible to assume that the phenomenon of dilation in the time of the CO 2 signals that occurs simultaneously with the extension of the radon signals (Figure 15), and of course independent with it, is genuine.
In general, any geodynamic process that causes stress or strain release affects subsurface fluids, irrespectively if the process will finally produce an earthquake or not.

In situ Monitoring of Subsurface Radon-Experimental Point of View
The long-term geophysical investigation of radon in shallow and deep boreholes, based on passive measuring systems (no pumping and gas circulation that disturb the local environmental equilibrium) was able to detect radon temporal changes, with very high sensitivity and high time resolution.It enabled the researchers to eliminate climatic-induced periodic contributions from the radon time series, and thus extract the part of the signals that could be correlated with regional geodynamic processes.The various monitoring systems that were found to be suitable for use in a passive way included various types of gamma and alpha detectors.
The basic principle to choose the gamma detection by SCA technique for total counting of RDPs' (radon decay products) gamma rays, is the assumption that radon is the only radioactive element that can vary in terms of its amount in the porous media under the impact of upper climatic or deeper geodynamic activities.These fluctuations are measurable as an additional increment to the measured gamma radiation from natural elements as uranium, thorium and potassium that form a constant background in time.This background also includes the contributions of all the RDPs atoms that remain confined inside the solid material and are not emanated to the porous media neither free to move under the influence of various driving forces.
Therefore, the variable contribution in the measurement of radon by the gamma detector does not actually exceed 10-25% of the total radiation that is measured.Thus, the first test for a system's performance is its capability to recover the daily changes in the temporal radon transport variation within a soil media, while the detector is lowered into ground depth of at least 2-5 m in an open area.It also needs to take into consideration that too lower radon contents do not allow to measure significant temporal variations in the concentration of radon, under the daily periodic behavior.Experimentally, it should be higher than 1 kBq/m 3 , few times above 200 Bq/m 3 , the average outdoor radon levels.
BGO and NaI gamma detectors were selected after being ascertained that they have high sensitivity, stability and suitability for use in a 2-inch diameter drilling holes to a depth of several dozen meters.We have also re-defined the energy range the SCA of the gamma detectors to count gamma rays between 475 and 3,000 keV, cancelling the measurement of gamma radiation in the field of the Compton scattering, between 50 and 450 keV, enabling to improve the ratio of signal to noise by 70% in comparison to the common detection systems (Figure 1A).
The main result was that gamma detectors have a higher sensitivity between 2 and 4 times in relation to the alpha detector, which enables them reveal to detect low sub-diurnal fluctuations.The construction of hermetic and water-protected housing enables the use of these radon gamma radiation detectors in a variety of environments such as soil, country rock, ground and seawater, and any subterranean spaces including drillings.

Discrimination Between the Influence of Climatic Temperature and Pressure on Radon Flow in Geological Media
Combined long-term radon measurement by gamma and alpha methods within deep abandoned boreholes empowered the capability to distinguish between radon anomalous signals at depth induced by temperature and barometric pressure.The main evidences of our research are: 1)The radon present in country rock formations that is measured by gamma radiation detectors is propelled by the surface temperature gradient to at least a proven depth of 100 m.The radon reveals a daily periodicity similar to the one per day cycle of the surface temperature.The gamma detector at each level, presents very sharp, clear and accurate peaks as a result of a high counting rate and low error, with a specific time lag between each other.It was found that the time lag depends on the downward radon velocity within the bedrock type.
The list of results that includes measurement from different monitoring sites is as follows: • At Sde Eliezer site, Israel, a time lag of 2 h was measured between 10 and 60 m depth, implying a downward radon velocity of 25 m per h (Zafrir et al., 2016).• At Nachal Mor site, Israel, a time lag of 2 h was measured between 10 and 40 m depth implying a downward radon velocity of 15 m per h (Zafrir et al., 2019).• At Amram tunnel site, Israel, a time lag of 10 h was measured, between 0 and 100 m depth implying a downward radon velocity of 10 m per h (Zafrir et al., 2013).• At Bentonite mine, in Makhtesh Ramon, Israel, a time lag of 5 h was measured between 0 and 40 m depth, implying a downward radon velocity of 8 m per h (Zafrir and Malik, 2010).• At Ghuttu site in Garhwal Himalaya, India, a time lag of 3 h was measured between 0 and 10 m depth, implying a downward radon velocity of 3.3 m per h (Choubey et al., 2011).
In general, from a geophysical point of view, there was no need to find any other gases that would bear radon up or down (e.g.Etiope and Martinelli, 2002;Neri et al., 2016).Once the Sun surface heating provides energy to the radon atoms their ballistic moment (mv) downward is at least 4 times larger than any other natural gas atoms.The same process can lead to the emergence of radon and other gases, close to the surface, under the effects of underground heating associated with volcanic events.
2) The radon concentration within any open underground space, such as an abandoned drilling, as measured by a detector of its alpha particles, is controlled by the negative barometric pressure gradient (in anti-correlation).
3) Most of the radon daily signals that are recognized as induced by climatic parameters, appear once or twice a day at specific times: early morning and afternoon (e.g.,, 4-6 AM and 5-7 PM, in 60 m depth at Sde-Eliezer site, Figure 9C).
4) The amplitudes of the radon periodical signals are controlled by the intensity of the climatic driving force, in linear dependency with the pressure gradient according to the existing physical model, and with largest non-linear variations induced by the ambient temperature gradient, that according to the ratio between the radon level in winter to summer varies by a factor of 3-10 while the temperature varies within 10% span (28 change vs. 285 K).
These phenomenological indications cause difficulties in accepting the results of any statistical analysis that is based on daily, multi-day, seasonal and multi-year averages calculations, of the radon gas levels in any media.Since the half-life of radon is short and constant (3.8 days), and naturally it always tends to be in secular equilibrium, only external driving forces will cause the radon to move within the geological media.Hence, such calculations actually reveal and reflect the temporal variations of the external climatic parameters, the ambient temperature and pressure, which undulate a combined change of the amplitude of the radon signal during daily scale, multi-day variation and significant changes between seasons too.This assumption is undoubtedly correct when the radon gas detector is installed at a shallow depth of a few dozen centimeters.
In addition to the difference in the intensity of the radon signals between the summer and winter, it is also possible to discern in all the above works, the multi day heat waves, mainly in the spring and summer, that are simultaneously producing waves of radon that are prolonging within the same time (e.g.,, Figures 6A,B).Over these heat waves there are exposed signals in a daily cycle with a peak that follows the maximum of the daily temperature, with a delay in time depending on the depth of the detector.The same phenomena can be observed with barometric pressure waves lasting few days that also affect the continuous level of the radon under the surface in addition to the semidaily and diurnal signals (Figure 7B).It is important to emphasize that the two phenomena with daily and multi-day variations, are not always dependent on each other, as for example, the semi-daily cycle that exists only in the periodicity of barometric pressure, prompted by an atmospheric tide (Zafrir et al., 2013).
Therefore, it is impossible to separate the influence of the climatic driving force by additive arithmetic operations.For example, the subtraction of the moving average of 7 days from the daily average as part of the statistical manipulation to calculate anomalies higher than twice the standard deviation is uncertain and problematic (e.g., Ghosh et al., 2009;Cannelli et al., 2016;Singh et al., 2017;Barkat et al., 2018;Alam et al., 2020).Even the attempt to adapt regression equations to radon time series utilizing every environmental variable and generating a best linear regression model to correlate between radon variations and meteorological parameters, does not enable to resolve and defined the subsurface contribution to the anomalous radon temporal variation (Neri et al., 2016).
5) There can be a number of ways for analyzing, identifying the nature of the anomalies in the measured time series, dismount to various components, and isolating the parameters operating as driving forces for the radon movements beneath the surface.The convolution is a mathematical procedure of combining different signals to form an another signal, or isolate components and detect signals with filtering and reconstruction techniques.It is the single and most important technique in Digital Signal Processing (e.g.,, Smith, 2002).
Although, there is no yet explicit mathematical expression for the thermal gradient rule in the radon transportation process, to be insert within the convolution integral, one can use the Fourier analysis to convert signals from their original time domain to a presentation in the frequency domain and vice versa.
We apply the Fourier domain filtering and component isolation procedure (by AutoSignal ™ of Systat Software Inc.).
Utilization of an existing algorithm that work in the frequency domain for low-pass truncation of frequency components above a chosen ideal cutoff frequency, enabled to remove them from the original radon time series, leaving the frequencies below it unaltered as well as any non-periodic discrete events.Choosing a frequency cutoff less the diurnal (one cycle per day), excludes all the unwanted diurnal, semidiurnal, and higher periodical signals (maybe harmonics) that are induced by the atmospheric parameters and proved that indeed we measured radon anomaly that its emergence lasts for more than 20 h and that these non-cyclical signals may precede by several hours or more a forthcoming seismic events even if they are weak.
It is important to note that according to our ongoing experience, the broad signals associated with tectonic activity almost cannot be revealed even by the measurements of radon temporal variation at 10 m depth.It was easy to observe within the long-time series acquired during the last 4 years, that the entire radon measured signals at the shallow depth of 10 m, used as reference for the deep data acquisition, are mostly radon periodical signals produced under the direct effect of the climatic driving forces on the ground surface.
As long as investigations continue to rely on detectors installed at a depth that does not exceed 10 m, the likelihood of discovering non-cyclical signals associated with tectonic activity in depth remains slim.
Currently, it is clear that there are gaps in our knowledge, including depths, time scale, range and magnitude of the energy release by pre-seismic geophysical processes and the relationship between radon temporal behavior and these various tectonic processes.Moreover, it seems that there is no proof for the assumption that the radon signal induced by seismo-tectonic processes has to appear with an intensity higher than the climaticinduced periodic signal, except in the case of a high-magnitude earthquake and a measuring system that is located near the event, such as occurred in Kobe, Japan (M 7.2 at 30 km distance, Igarashi et al., 1995).
Radon, CO 2 and Other Natural Gas in Depth as a Proxy for Investigating Tectonic Pre-seismic Activity As observed in the Sde Eliezer results, radon signals that are induced by the ambient temperature and pressure have diurnal and semidiurnal periodicity and amplitudes which are variable and not uniform.In general, they preserve the multiday profiles of the atmospheric temperature and barometric pressure as well.Monitoring radon at a depth of several dozen meters, substantially attenuates the climatic contribution and increases the possibility of resolving from the radon temporal spectrum the preseismic radon signals that are not periodic and are independent from the atmospheric driving forces.
In parallel it was observed that CO 2 , within the internal airspace of the borehole at Sde Eliezer, within various depth (5, 20, 40, and 70 m) follows the radon (measured by alpha detector at 40 m) as well as the radon temporal variations at the surrounding bedrock (measured by gamma detectors at 10, 60, and 88 m), and both are driven by the climatic parameters as analyzed for radon (e.g., Figure 15).Thus, CO 2 reveals the same daily and semi daily temperature and pressure periodicities in the frequency time domain and clearly responds, simultaneously with radon, to the apparently tectonic pre-seismic driving force, with the same rise, fall, and broadening times.
Improving the sampling rate by 30 times (from 15 min to high 30 s resolution) and implementing four CO 2 detectors at different depths enabled measuring the CO 2 velocity in the Sde-Eliezer borehole under the influence of barometric pressure that seems to be more effective than temperature within the air space inside the open well (Levintal et al., 2020).
Similar to radon, the high CO 2 concentration above the atmospheric content (about 400 ppm) and up to 2,800 ppm at Sde-Eliezer site, entered the sealed borehole's iron casing through the perforated lower sections below 80 m, and then degassing into the sealed casing from the water table.The source of the high CO 2 content could be the microbial activity within the vadose zone, or the presence of an organic substance (peat/lignite) trapped in the Hula Valley's bedrock.Hence, it is reasonable to assume the same thermodynamic mechanism supplies CO 2 as well as radon to the groundwater.
The above conclusion means that the two gases, radon and CO 2 are driven down within the geological media by the combined temperature ΔT and pressure ΔP gradients.Therefore, they reveal the same daily and semi-daily periodicity (Figure 15).Now, it is expected that recovering of non-periodic broadening of radon and CO 2 signals, enables eliminating the climaticinduced periodic contributions and extract the residual portion of the radon and CO 2 signals related to the regional tectonic preseismic processes, and could be used as earthquake precursors before the main rupture occurs.It is not essential that each one of these geodynamic processes will end as an earthquake.
In order to use the CO 2 temporal variations as a proxy for preseismic activity its concentration within the geological media, should be at least twice as high as the atmospheric contents of about 400 ppm.Otherwise, any change will only be a fluctuation of the atmospheric concentration by few tens of ppm under the influence of the barometric pressure, as measured at Nachal Mor site (Zafrir et al., 2019) or even reaching lower levels as 200 ppm in isolated tunnels as in Amram site (Zafrir et al., 2013).The results at Nachal Mor for example, prove that there are no organic materials in rocks and soil around the Dead Sea and they are not the source of radon there either (Kronfeld et al., 1991).
Another important result that emerged is our deep gas monitoring technology.The technology may become a useful tool for the investigation of seismic precursors since similarly to radon and CO 2 , the existence of any natural gas such as nitrogen, oxygen, methane, hydrogen sulfide, carbon monoxide and helium within deep subsurface media can serve as a proxy for pre-seismic precursory phenomena.
In general, any geodynamic process that causes stress or strain release affects subsurface fluids, irrespectively if the process will finally produce an earthquake or not.

FIGURE 2 |
FIGURE 2 | (A) Map (Israel grid) showing relation of radon measurement sites relative to 20-ton explosion zone.(B) Radon monitoring was performed at a depth of about 2 m and a "Sensing Volume" of about 0.65 m 3 .(C) Overview of 15-min resolution time series of five gamma detector systems during the 13-days experiment.It presents the variations in radon concentration in the porous media in Bq/m 3 around sensors (By total gamma counting rate per 15 min minus the background counting rate).The time of the explosion is marked with a vertical line.All sensors show build-up of Radon in the porous media within the "Sensing Volume", and typical daily variations.(D) Radon concentration measured 15 h before and after the explosion, including the high resolution 10-s data averaged to 15 min (in squares).Explosion time is indicated.

FIGURE 3 |
FIGURE 3 | Time-series of measurements at the Amram tunnel (15-min temporal resolution), from top to bottom panels: (A) Radon by gamma, and radon by alpha) during 6 years.(B) (left) The similarity between the two seasonal patterns becomes evident when considering the difference between the external temperature and the constant temperature inside the tunnel (temperature gradient).(B) (right) The intra-seasonal variability is also associated with changes in the temperature gradient between the inner part of the tunnel and the outside environment.

FIGURE 4 |
FIGURE 4 | (A) Radon setup in the Amram tunnel.(B) Variations of radon from 20 February to March 18, 2009, as measured simultaneously by seven different radon sensors.The AlphaGUARD system(G), which stopped measuring on March 12, presents the radon concentration in Bq/m 3 per 60 min measurement cycle time, while the rest present the CR (Counting Rate) per 30 min (A), or per 15 min (B-F).Two radon anomalies occurred: a small (in amplitude) and short (in time) one for almost two days from 26 to 28 February and a predominant one from 4 to March 18, 2009.(C) a conversion scale that enable to evaluate the counting rate, CR, per 15 or 30 min of the various detectors sensors (A-F) at peaks in Bq/m 3 .

FIGURE 7 |
FIGURE 7 | (A) The radon setup in Makhtesh Ramon at Gevanim valley including alpha detectors at 1.2 and 85 m. (B) Radon time series as measured by alpha detectors in the boreholes air space near the surface at 1.2 m depth and at depth of 85 m. (C) The FFT amplitude spectrum of the ambient pressure and radon temporal concentration changes at the Gevanim Valley, as measured by alpha counting at 1.2 m depth.

FIGURE 8 |
FIGURE 8 | (A) Sde Eliezer site in the Hula basin western border fault-HWBF, and map of the main segments of the Dead Sea Fault Zone (DSFZ) in northern Israel (Weinberger et al., 2009).(B) The radon monitoring setup.(C) The geological column section.

FIGURE 9 |
FIGURE 9 | (A) Example of time interval of 4 days in which the radon at different depths consistently follows the daily temperature profile with a variable lag between each other.The dotted lines in the figure mark daily noon and the full lines the midnight.Radon early morning semi-diurnal signals are exhibited by the three radon detectors with anti-correlations to the pressure (see arrows) in addition to their response to the daily temperature periodicity after noon (see stars).The parameters are: ambient barometric pressure, radon by gamma detectors (at 10 and 60 m), radon by alpha detector (at 40 m), ambient temperature and internal temperature (at 10 m).(B) The radon time series at 60 m depth in relation to the gradients of the climatic parameters: pressure (ΔP, dark purple) and temperature (ΔT, red).(C) The daily normalized values of each parameter's (maximum for radon at 10 and 60 m, ambient temperature, radon at 40 m and minimum for barometric pressure) summarized for a period of 150 days: 16 January to June 14, 2015.(D) The basic principles for the behavior of radon movement under the climatic influence.Radon in a rock media (as measured by gamma detector) is driven by the surface temperature gradient ΔT.Radon in the open borehole air space (as measured by alpha detector) is driven by the pressure variation ΔP(Zafrir et al., 2013).In this figure, three counts per 15 min of the gamma sensors peaks at 10, and 60, are equivalent to concentration of 1 Bq/m 3 within the surrounding country rock, and 0.018 counts per 15 min of the alpha particle detector peak is equivalent also to 1 Bq/m 3 .

FIGURE 11 |
FIGURE 11 | (A) Continuous time series of the measured parameters at the Sde-Eliezer site (with 15-min temporal resolution) collected during the first two and a half years (844 days).The parameters measured including: barometric pressure (dark red), ambient temperature (green), borehole temperature at 10 m depth (pink), radon within the surrounding rocks at 10 m depth as measured by gamma rays (red), radon within the borehole air space at 40 m depth as measured by alpha particles (dark yellow), radon within the surrounding rock at 60 m depth as measured by gamma rays (blue), radon within the groundwater at 88 m depth (since November 2015, cyan).(B) More than 260 earthquakes with magnitude greater than M 2) recorded by the Israel Seismic Network of The Geophysical Institute of Israel [http://seis.gii.co.il].For three time periods marked in letters a, b, c, there were three earthquakes preceded by an unusual radon anomalous signal.

FIGURE 12 |
FIGURE 12 | Case a: (A)The first-time interval, June 2015, that reveals anomalous broadened radon signals while most of the radon signals are recognized as induced by climatic parameters by appearing once or twice a day at specific times, 3-4 A.M. and 4-8 P.M. (B) The 26 June radon anomalous signal is different from the radon periodic signals appearing once or twice a day.Case b: (c) This time interval reveals anomalous broadened radon signals on 25 april, 1 and 15 May 2016, preceding by days and hours the Nuweiba, M 5.1 earthquake (NU) of May 16, 2016.(D) It happens while most of the radon signals are recognized as induced by climatic parameters, by appearing once or twice a day at specific times.Case c: (E) and (F) These time intervals between 30 June to July 26, 2016, reveals two anomalous broadened radon signals of at least 20 h each, precede the Palmyra Syria M 4.5 earthquake (Pa) on 25 July.It occurred 300 km northeast of the Hula Valley on one of the branches of the main fracture into Syria called Palmyrides, (Al-Zoubi et al., 2006).The signal from the deepest detector at 88 m depth became the (Continued )

FIGURE 12 |
FIGURE 12 | most pronounced since the water table within the borehole of the abundant well, dropped down below 100 m, by unique event on June 25, 2016 (see Figure11) as a result of massive pumping in the Hula Valley area, since February 2016, because of a drought.The EQs locations are Palmyra (Pa), Gulf of Eilat (GE) and East Sinai (ES).(G) and (H) are the FFT low-pass filtering and the reconstruction spectrogram of the radon's 60 m time series in the time interval of 30 June to July 27, 2016, shows the discrete radon broadened signals that are marked in the with black arrows.The results of the FFT filtering and the time domain reconstruction of the above time intervals (E) show that although the procedure does not disable the entire cyclical signature in the original data, it keeps the seasonal fluctuated background under 30,000 counts per 15 min at the 60 m radon gamma detector.(H) The reconstruction STFT spectrogram of the radon 60 m time series in the same time interval.The procedure pinpoints the two broadened signals 6 and 2 days before the earthquake on July 25, 2016.The FFT analysis have done to each one of the above cases.In this figure, three counts per 15 min of the gamma sensors at 10, 60 and 88 m, are equivalent to concentration of 1 Bq/m 3 within the surrounding country rock, and 0.018 counts per 15 min of the alpha particle detector is equivalent also to 1 Bq/m 3 .

FIGURE 13
FIGURE 13 | (A) Accumulation of 103 earthquakes in the Dead Sea seismogenic active faults at the bottom of Lake Kinneret, which were recorded during the 2018 summer between 4 July and 18 August.(B) The tectonic map showing the internal structure of the lake(Hurwitz et al., 2002).EMF and WMF are the eastern and western marginal faults, KNF is the Kfar Nahum fault, and JF and AF are the Jordan and Almagor faults.

FIGURE 14 |
FIGURE 14 | Time series interval of 75 consecutive days at the Sde-Eliezer site.The Sde-Eliezer site is located about 35 km north of Lake Kinneret that was under weak seismic activity during July and August 2018.The non-periodic, broader than 20 h signals, on 12 (a) and 22 (b) June, 4 (c), 13 (d), 22 (e), and 28 & 29 (f) July, and 2 (g), 7& 8 (h), and 12&13 (i) August, may reveal the pre-seismic progression of several weak geodynamic events in depth (details in Figures24, and 25).In this figure, 0.1 counts per 30 s of the gamma sensors at 10, 60 and 88 m, are equivalent to concentration of 1 Bq/m 3 within the surrounding country rock, and 0.0015 counts per 30 s of the alpha particle detector is equivalent also to 1 Bq/m 3 .