Pre-seismic variations of atmospheric radon activity as a possible reason for abnormal atmospheric effects ANNALS OF GEOPHYSICS, 58, 5, 2015, A0554; doi:10.4401/ag-6808 A0554 Pre-seismic variations of atmospheric radon activity as a possible reason for abnormal atmospheric effects Vadim V. Surkov National Research Nuclear University MEPhI, Moscow, Russia ABSTRACT Possible atmospheric effects associated with pre-seismic variations of at- mospheric radon concentration occasionally observed prior to earthquake occurrence are subjected to theoretical study in order to assess their po- tential effectiveness. Altitude distribution of atmospheric conductivity is examined as a function of radon and aerosol concentrations. Horizontal components of atmospheric electric fields and currents, which can arise near a boundary of seismo-active region, are estimated. We speculated that changes in radon-induced air ionization can have an impact on fair- weather spectra of ULF electric pulsations. Particular emphasis has been placed on hypothesis for correlation between variations of pre-seismic radon activity and upward thermal radiation measured by satellite over the high seismicity regions. This hypothesis is based on assumption that the radon-induced ionization of atmospheric surface layer can trigger water vapor condensation which in turn results in release of the vapor- ization heat. Our analysis has shown that the upward thermal flux has to be many orders of magnitude lower than that predicted on the basis of this hypothesis. 1. Introduction It is generally believed that the main ionization of the atmospheric surface layer is due to galactic cosmic rays and radiation as well as the naturally occurring ra- dioactive elements naturally found in the air. Shortwave solar radiation and energetic particles of the solar wind predominately affect the upper atmosphere above alti- tudes of 50 km. A mean rate of ion production under the cosmic rays is approximately the same in all points of the Earth surface at all times. At sea level its value is cc = 1.5–1.7 ion pairs per cm 3 per second [e.g., Tverskoy 1962, Chalmers 1967, Israël 1970]. The radioactive ele- ments contained in the surface layer of the soil play in- significant role. For the most the air ionization is due to radioactive atoms which first enter to the atmosphere from the soil and then accumulate in the surface layer. On land, the ionization rate due to radioactive elements that are naturally present in the atmospheric surface layer is estimated to be cr= 8–10 ion pairs/(cm 3 s) al- though this value is slightly lower near the ocean sur- face. According to other estimates both the mechanisms, i.e. the cosmic rays and radioactivity, make approximately equal contribution to the ion production process in the surface layer [e.g., Tverskoy 1962, Zhang et al. 2011]. The ionization rate cc due the cosmic rays increases with increasing altitude and it reaches a peak value at the altitudes 14–15 km [e.g., Tverskoy 1962]. Alterna- tively, the number density of the radioactive atoms de- creases with increasing altitude thereby producing the decrease in the ionization rate cr in such a way that cr becomes practically negligible at the altitudes 4–5 km. A significant factor, which may greatly affect the background radiation, is the heavy gas radon that has the tendency to penetrate permanently along the sys- tems of capillaries and cracks through the rocks into the soil and air. Among many radioactive rocks are ura- nium-containing rock, granite, gneiss, shale, tuff, metal- bearing ore, phosphorite, and potassium-containing rock. The 22286Rn nuclei are efficient sources of the air ioniza- tion since their half-life is about 3.82 days [e.g., Protec- tion against Radon-222 at Home and at Work, 1994]. The half-life of other radon isotopes; that is, 21986Rn and 222 86Rn is much smaller than the above value. The radon con- centration in the air depends on a host of factors such as local species of rocks, the configuration and state of soil surface, the possibility for removal of radioactive ingredients due to rainout, radon flux from depth, the rate of radon dissipation in the atmosphere which in turn depends on the meteorological conditions includ- ing the direction and strength of wind, pressure, air hu- midity and etc. The enhancement of background radiation has been occasionally observed in the regions of tectonic faults of the Earth crust [Toutain and Baubron 1998]. Some researchers have speculated that the abnormal Article history Received June 15, 2015; accepted October 2, 2015. Subject classification: Atmosphere, Aerosols, Radioactivity and isotopes, Earthquake source and dynamics, Heat generation and transport. variations of radon concentration in the atmospheric surface layer can be due to the earthquake preparation processes [Pierce 1976]. A review of anomalous radon precursors have been reported for a very wide range of earthquake magnitudes [Cicerone et al. 2009]. For ex- ample, Virk and Singh [1994] have reported that radon concentration in the soil gases and in the spring water suddenly increased by 2.5 and 1.5 times a week before Uttarkashi earthquake. These measurements have been pursued at the epicentral distance of 300 km. Yasuoka et al. [2009] have detected a gradual increase in the at- mospheric radon emanation two months before Kobe earthquake occurred on January 17, 1995 (Mw = 6.9). The radon volumetric activity was increased up to 20 Bq/m3 just before an earthquake and it was approxi- mately doubled as compared to the background level. It should be noted that the evidence of pre-seismic radon anomalies in the atmosphere are fewer than those detected directly in the soil or in indoor locations [e.g., Cicerone et al. 2009]. However, some researchers have not found any statistically significant changes in radon concentration before earthquakes [Geller 1997, Inan et al. 2008]. Whether the radon concentration variation is a plausible precursor of impending earth- quake or not is still an open question. Also the answer to this question is likely to be dependent from the prop- erties of seismogenetic fault. 2. Local changes of electrical resistance of the at- mosphere The electric conductivity of the air increases ap- proximately exponentially with height if standard con- ditions prevail. In this notation the surface atmospheric layer makes a major contribution to the total electric resistance of the atmosphere. The conductivity of the lower atmosphere is predominately due to small/light ions that have the greatest mobility. The small ions build up as a result of the attachment of a few neutral molecules to the primary ions, which in turn are caused by the cosmic rays and ionizing radiation of the ra- dioactive substances. The small ions many times collide with other particles and take part into enormous num- ber of chemical reactions during their life-time. The composition and number of molecules entering into the composition varies with time. The mean number of these molecules depends on air humidity. On the av- erage, 10–20 neutral molecules are held around each positive ion while a negative small ion contains 5–10 neutral molecules [e.g., Tverskoy 1962, Chalmers 1967, Israël 1970]. To estimate the radon variations impact on the electric resistance of the surface atmospheric layer, we first consider the simplified model of the atmosphere in which there are only two species of opposite-charged small ions. Suppose that the small ions may disappear either by direct recombination with oppositely charged small ions or by attachment to aerosol particles. Tak- ing into account that the concentration ni of both ion species are approximately the same, we come to the fol- lowing equation [e.g., Israël 1970] (1) where a denotes the ion-ion volumetric recombination coefficient, na is aerosol concentration, b is the ion- aerosol attachment coefficient, and J stands for the small ion flux density due to entrainment of small ions with air fluxes, small ion diffusion, small ion drift in elec- tric field and etc. Here we shall examine the quasi-steady small ion distribution in which ^tni = 0 and J is assumed to be constant value. It follows from Equation (1) then [e.g., Harrison et al. 2010] (2) The values of parameters entering Equation (2) can strongly vary with altitude. Whenever the atmos- phere is in thermodynamic equilibrium, the concen- trations of the aerosol particles and radon atoms as well as the rate cr of ion production decrease with altitude according to Boltzmann distribution. By contrast, the rate cc ion production by cosmic rays increases approxi- mately exponentially with altitude z [Wåhlin 1994]. Ac- tually, all these altitude dependences can be approximated by exponential functions as follows where cc0 and cr0 are the ionization rate at the sea level. Here we made use of the following abbreviations lr = kBT/(mr g) and la= kBT/(ma g), where kB is Boltzmann constant, g is acceleration of free fall, mr is a mass of radon atom, ma is an average mass of aerosol particles, and lc is an unknown parameter that will be estimated later. Substituting the temperature T = 293 K into Equation (3), yields lr≈1.1 km. Taking into account that the typical sizes of the aerosols possess value within 0.002–1nm [e.g., Ivlev and Dovgalyuk 1999, Kondratyev et al. 2006], the length la is no greater than several me- ters even for the smallest aerosols. These values can SURKOV 2 (3) 3 serve as rough estimates since the actual altitude dis- tributions of the radon atoms, small ions and aerosols can be subjected to local meteorological conditions, air flow, turbulent and convective gas motion, which can drag these particles, and many other factors. Consider- ing singly charged small ions, the air conductivity is given by v= 2ninie, where e denotes the elementary charge. The small ion mobility ni is inversely propor- tional to the concentration of neutral molecules which in turn exponentially decreases with altitude. Whence it follows that ni =ni0 exp (z/H), where ni0 is the small ion mobility at the sea level and H≈8 km denotes the scale-height of an isothermal atmosphere. Combining Equations (2) and (3) with the above equations, we ob- tain the air conductivity To make an estimate, we use the following nu- merical values [e.g., Pierce 1976, Harrison et al. 2010]: n0≈1.2·10 -4 m2V-1s-1, a= 1.6×10-12 m3s-1, b = 4×10-11 m3s-1 (for the aerosol size of 0.2–0.25 nm). At the ground level the aerosol concentration na0 varies from 200 cm -3 for clean air to 1.5 × 103 cm-3 for dusty air [e.g., Harri- son et al. 2010]. In such a case the condition b2n2a >> 4a(cc+cr) is valid for the atmospheric surface layer. In fact this means that one may neglect the term an2i which describes the small ion recombination in Equation (1). In this extreme case Equation (4) for air conductivity is reduced to (5) With the above parameters one can obtain that for altitudes greater than several tens or hundreds of me- ters, the aerosol concentration na decreases so much that the recombination term becomes greater than bnina. In such a case Equation (4) is simplified to (6) In this altitude range Equation (6) must coincide with the familiar empirical dependence v=v0exp (z/z0) which is valid if z < 3.6 km [Chalmers 1967]. Whence it follows the approximate relation between the parame- ters: H-1+ (2lc) -1 = z0 -1. Taking into account that z0 = 0.82 km we can estimate the unknown parameter lc as fol- lows: lc = 0.5Hz0/(H − z0)≈0.46 km. The variations of radon content may greatly affect the air conductivity and this effect depends on the aerosol concentration. To illustrate this tendency, the altitude dependences of v calculated from Equation (4) are shown in Figure 1 with lines 1 and 2 for the clean air (na0 = 2 ×10 2 cm-3) and with lines 3 and 4 for the dusty air (na0 = 1.5 ×10 3 cm-3). In making the plot of v we have used the above parameters and la =50 m. As is seen from Figure 1, the air conductivity decreases with an increase in the aerosol concentration, which is in turn resulted from the enhancement of the dustiness of the atmos- phere. The implication here is that the majority of the newly formed small ions connect with the aerosol par- ticles which do not contribute to the air conductivity because of their low mobility. In Figure 1, we have plotted functions v(z) for two different rates of the radon-induced ion production in order to illustrate the radon effect on the air conduc- tivity. The lines 1 and 3 correspond to cr0 = 6 cm -3s-1 while the lines 2 and 4 correspond to doubled value of cr0. For example, the similar changes in radon concentra- tion has been observed by Yasuoka et al. [2009] prior to the 1995 Kobe earthquake. Notice that cr varies directly with radon activity concentration which in turn is pro- portional to radon concentration. It is obvious from Figure 1 that the increase of cr0 causes the increase of the air conductivity by 10–20%. The radon effect disappears (that is, of the order of lr) at the altitudes above 1 km in such a way that all the dependences tend to the unified PRE-SEISMIC VARIATIONS OF RADON (4) Figure 1. Model distribution of the air conductivity with altitude. The air conductivity for the aerosol number densities na0 = 2 ×10 2 cm-3 is shown with lines 1 and 2 while the same value for na0= 1.5 ×10 3 cm-3 is shown with lines 3 and 4. The sea level rate cr0 of ion pro- duction for the lines 2 and 4 is twice as great as that for the lines 1 and 3. Calculation from the empirical equation v=v0 exp (z/z0) is shown with line 5. i0 c0 asymptotic dependence given by Equation (6). The local variations of radon concentration in the surface layer may affect the total resistance of the at- mosphere/columnar resistance because the surface layer has the greatest air resistance. To illustrate this ten- dency, we calculate the resistance R per unit of area of the atmospheric surface layer with thickness zs = 3 km Figure 2 shows the dependence of R on ratio k = cr0/cc0. In making this plot we have used the above parameters and Equation (4) for v(z). It is generally be- lieved that the typical value of k can vary from 2 to 5 [e.g., Tverskoy 1962, Zhang et al. 2011]. Setting the background value of k = 4 and using the data by Virk and Singh [1994] and by Yasuoka et al. [2009], we sup- pose that for seismo-active period this value to be twice as large as that for the background value; that is, k ≈ 8. As is seen from Figure 2, this can result in the decrease of the surface layer resistance by 15% for the clean air (line 1) and by 25% for the dusty air (line 2). As would be expected, considering the important role played by the surface layer in the atmospheric resistance, these es- timates hold true for the total columnar resistance of the atmosphere. 3. Radon gas activity and electrical structure of the lower atmosphere Harrison et al. [2010] have suggested that the ab- normal decrease in the atmospheric surface layer re- sistance caused by an enhancement of the radon volumetric activity can result in an increase in the ver- tical fair-weather conduction current ja by several units or tens percent. It follows from Maxwell equation that the density of electric charge t=f0jadv -1/dz, which builds up as a result of the inhomogeneity of the air conductivity, must be proportional to the vertical cur- rent density ja. This implies that the relative changes of t can vary together with ja within the same range. A number of researchers have speculated that the en- hancement of the fair-weather conduction current is followed by an increase in the electric charges of water droplets suspended in the air. This effect is assumed to be favorable for the steam condensation which in turn may explain the observation of the so-called earth- quake clouds [Rycroft et al. 2008, Rycroft et al. 2012, Harrison et al. 2013]. Here we deal with other radon effects which may affect the electrical structure of the area adjacent to ab- normal region. First, we shall estimate a horizontal component of the fair-weather conduction current in the vicinity of the abnormal region with the enhanced radon concentration. Taking into account that the con- duction current density is not a function of altitude, the distribution of electric potential with altitude can be written as where v(z) is given by Equation (4), and {a≈ 250 kV [e.g., Rycroft et al. 2000] is the potential drop between the bottom of the ionosphere z = za and the Earth z = 0. Figure 3 shows the result of numerical integration of Equation (8). As before the lines 1 and 3 correspond to cr0 = 6 cm -3s-1 while the lines 2 and 4 correspond to cr0 = 12 cm -3s-1. It is obvious from Figure 3 that an in- crease of the ionization rate cr0 leads to decrease of po- tential (8) at given altitude. It also can be seen that there are the essential distinctions between background (lines 1 and 3) and abnormal (lines 2 and 4) potential distri- SURKOV 4 (8) (7) Figure 2. Total columnar resistance R of the atmospheric surface layer versus parameter k = cr0/cc0. The lines 1 and 2 correspond to aerosol number densities of 2 × 102 cm-3 and 1.5 × 103 cm-3, respec- tively. The plots are normalized to the value of R taken at k = 1. Figure 3. The same as in Figure 1 but for the atmospheric potential. 5 butions in the altitude range 02–2 km. In order to high- light this peculiarity, the differences between the corre- sponding distributions as a function of altitude are shown in Figure 4. Here the line 1 demonstrates the dif- ference between the first and second distributions shown in Figure 3 whereas the line 2 demonstrates the same but for the third and fourth distributions. This constitutes a mapping of the surface radon activity into space variations of the atmospheric po- tential at the altitudes of several kilometers. The recent measurements [Kabrt et al. 2014] of radon gas activity concentration distributed over a large area can serve as a basis for the searching of the seismic-induced radon anomaly region. As would be expected, the horizontal electric fields and currents occur near the boundary of the abnormal region where the radon concentration changes from the enhanced value to the background one. The reason is that the height-dependence of po- tential in the abnormal area is different from that in the surrounding regions. The horizontal field components can be estimated roughly as follows E||~d{/dl and j||~vE|| where dl is the horizontal space scale at which the radon number density decreases down to the back- ground level while d{ is the corresponding potential difference taken at the fixed height. It follows from Fig- ures 1 and 4 that maximum value d{= 13–25 kV is achieved at such altitudes where the air conductivity v = 8 ×10-14–10-13 S/m. Taking these values and choos- ing dl = 10–50 km, we come to the following estimates: E||≈ 0.26–2.5 V/m and j||≈ 0.02–0.25 pA/m 2. Since the radon anomalies are subjected to temporal variations around seismic event, the above values can only serve as order-of-magnitude estimates of the possible atmos- pheric perturbations resulted from the enhanced radon emanation. A more significant effect can be expected in the vicinity of mesoscale convective systems (MSC) since there is the background atmospheric current well above the fair-weather one. The lateral size and total current of the MSC electrical active region are greater than 150 km and several tens Ampere, respectively, that is much larger than typical thunderstorm parameters [Davydenko et al. 2004]. A portion of this current can close through the abnormal region of the lower at- mosphere with the reduced resistance thereby produc- ing the charge relaxation and changes in electrical structure of MSC. It is pertinent to note that there are a few possible mechanisms which can provide a coupling between ultra low frequency (ULF) electric field and charge den- sity pulsations in the surface atmospheric layer and the variations of radon activity concentration. The fair- weather electric field pulsations have a power-law spec- trum which falls off with frequency as f -m. In the fre- quency range of 0.004–0.06 Hz the spectral index m varies from 2 to 2.4 [Yerg and Johnson 1974] while in the range of 0.01–1 Hz the electrical pulsations and their power spectra can be split into two groups; that is, the so-called “structured spectra” with m = 2.03–3.36 and “unstructured spectra” with m = 1.23–2.89 [Anisi- mov et al. 2002, 2005]. The “structured spectra” has been assumed to be due to aero-electric structures (AESs) in- volved in the neutral-gas turbulence in the surface layer. The enhancement of radon activity concentration, fol- lowed by the production of additional small ions and charged aerosols, can result in the increase of both the space charge and the variance of the charge-density fluctuations v2p. According to the model by [Anisimov et al. 2002] the ULF spectrum of electric pulsations is proportional to v2p and thus the variation of radon con- centration can affect the spectrum amplitude. We can- not come close to exploring these topics in any detail although one may assume that the spectral index is a function of radon concentration as well. The release of radon from the soil to the air de- pends on local rock morphology and genesis including faults, zones of rock fracture and etc. Overall, the radon gas is nonuniformly distributed near the ground surface. This implies that the pre-seismic radon varia- tion in seismo-active region is capable of producing the random spatiotemporal variations of the charge density and electric field in the air. These electric fluctuations appear to have a power-law spectrum in analogous to spectrum of natural ULF magnetic noises produced by random fluctuations of the ionospheric and atmos- pheric currents [Surkov and Hayakawa 2007, 2008, 2014]. Despite this effect seems to be weak, it can influence PRE-SEISMIC VARIATIONS OF RADON Figure 4. Model calculation of the atmospheric potential difference which might arise in horizontal direction. The potentials were cal- culated at the same altitudes but for different ion production rates, cr0, which differ by the factor 2. The lines 1 and 2 correspond to na0 = 2 ×10 2 cm-3 and na0 = 1.5 ×10 3 cm-3, respectively. on the fair-weather spectral index m of natural electric pulsations. So, the changes in the ULF electric spectra amplitude and spectral index can be indicative of pre- seismic tectonic activity. 4. Release of radon from the ground and temperature of the atmospheric surface layer Statistical analysis of satellite thermal imaging data measured over the high seismicity regions and fault sys- tems in the Earth’s crust have shown that the upward infrared (IR) radiation from the atmosphere can exceed the background value as compared to the surrounding regions [Gorny et al. 1993, Tramutoli et al. 2001, Tronin et al. 2002, Ouzounov et al. 2006, Genzano et al. 2009]. This effect is generally interpreted as a by-product of mean temperature increase of the atmospheric surface layer in seismo active regions. The mechanism of the phenomenon was presumed to be due to the release of optically active gases from soil to the air followed by local greenhouse effect [Tronin et al. 2002]. One more a plausible explanation of this phenomenon is the gradual squeezing-out of the groundwater from higher depth towards the ground surface that results from the varia- tions of the tectonic stress in the fault zone. Taking into account that, on the average, the ground temperature increases with depth by 2 K per 100 m, Surkov et al. [2006] have shown that this leads to the convective heat- ing of the ground surface by a few Kelvin since the tem- perature of the deep groundwater is higher than the mean temperature of the ground surface. It was also hypothesized by Pulinets and Ouzounov [2011] and by Pulinets [2011] that the enhancement of the outgoing IR radiation several days before seismic event can be due to the increase in radon concentration in the surface layer. This hypothesis is based on as- sumption that the radon-producing air ionization gives rise to the formation of so-called large cluster ions with typical sizes 1–3 nm; that is, as large as the aerosol par- ticles. These large ions can serve as condensation cen- ters triggering the water vapor condensation which in turn is accompanied with the release of the heat of va- porization. To estimate the maximal value of this effect, we assume that the alpha decay of radon nuclei occurs in a humid air near the ground surface. Let A be the vol- umetric radon activity and Ni be the number of ordi- nary/primary ion pairs resulted from the scattering of single alpha-particle in the air. In a short time the clus- ter ions are formed due to attachment of neutral mol- ecules to the primary ions. The polar molecules of air, including the molecules of H2O, are largely involved in this process. Assuming for the moment that every pri- mary ion gives rise to the formation of cluster ion, the rate of the cluster ion production per unit volume is ANi. The volumetric rate of H2O molecule attachment to the cluster ions is ANiNw, where Nw stands for the mean number of H2O molecules attached to the single cluster ion. The thermal energy released to the atmos- phere per unit time and per unit volume is ANiNwE1, where E1 is the energy release due to attachment of sin- gle H2O molecule to the cluster ion. As a rough esti- mate of E1, one may use the value m/Na, where m≈41 kJ/mol [Gray 1972] is the molar heat of vaporization for water and Na is the Avogadro number. Let hr be the thickness of the near-surface layer which contains the basic part of atmospheric radon. If the heat capacities of the ground and air are both ignored, the maximal value of upward heat flux density can be estimated as follows: qmax = ANiNwhrm/Na. According to the hypothesis by Pulinets and Ouzounov [2011], the so-called big cluster ions (BCI), “which are actually the cluster ions formed after ion- ization and ion’s hydration”, play a crucial role in a water vapor condensation followed by the pre-seismic thermal anomaly. Their numerical estimate is based on assumption that, first, the number of H2O molecules attached to the BCI ranges as high as Nw = 1.3 × 10 11; second, the pre-seismic radon activity reaches the value A = 2 ×103 Bq/m3 and Ni = 4.5×10 5 ion pars [Pulinets and Ouzounov 2011]. Furthermore, it was implicitly supposed that every primary ion is converted to the big ion. Taking hr = 1 m and substituting the above numer- ical values of the parameters into the equation for qmax , gives the estimate: qmax = 16 W/m 2 which coincides with the numerical estimate given by Pulinets and Ouzounov [2011]. Despite this value is believed to be compatible with the observations of abnormal pre-seismic increase of the upward thermal flux [Ouzounov et al. 2006], such an approach has the serious shortcoming that the rate of the BCI production ANi is dramatically overesti- mated. Indeed, the radon activity may affect the BCI rate production through the interaction between the radon-produced small ions and the BCI. Taking into ac- count that, first, about 5–20 neutral molecules are at- tached to the small ion; second, only a part of them are the H2O molecules, it becomes apparent that at least about 1010 small ions needs to create the BCI about which Nw = 1.3 × 10 11 H2O molecules are held. This means that the production rate of the BCI is smaller than that of the small ions by 10 orders of magnitude. In other words, the value ANi = 9 ×10 8 m-3s-1 of the cluster ion production can be applied to the small ions but not to the BCI. Moreover, it seems likely that the parameter A = 2 ×103 Bq/m3 used in the cited papers was also overestimated because this value of radon ac- SURKOV 6 7 tivity concentration is typical for the soil rather than for the atmosphere. It is common knowledge that the world average value of volumetric radon activity, on being measured in the air 1 m above the ground, varies from 7 to 12 Bq/ m3 though this value can reach 50 Bq/ m3 in areas with radon-laden soil [e.g., Protection against Radon-222 at Home and at Work, 1994]. In addition, we recall that the pre-seismic increase in atmospheric radon activity concentration before Kobe earthquake was as high as 20 Bq/m3 [Yasuoka et al. 2009]. If this latter value of A is taken, the above estimate of qmax is further lowered by two orders of magnitude. 5. Discussion and conclusion Our theoretical analysis favours the view [e.g., Harrison et al. 2010, 2013, 2014] that the variations of atmospheric radon concentration may affect the colum- nar resistance of the atmosphere. For example, accord- ing to our estimate the pre-seismic increase in the radon activity concentration by two times can result in the local lowering of the standard atmosphere resist- ance by 10–20% and approximately the same relative increase in the fair-weather/background conduction current. It follows from our model that the greatest changes of the air conductivity can be expected at the altitudes 100–200 m in the case of dusty air of high aerosol content. Although, the radon and aerosols ap- pear to have no effect on the air conductivity above 1 km since the cosmic rays and other mechanisms for air ionization play a major role in this altitudes range. The enhanced radon release from the ground fol- lowed by the variations of the air conductivity may af- fect the electrical structure of the atmosphere in seismo active regions. For example, not only vertical but also the horizontal components of both the background cur- rent density and of the gradient of atmospheric electri- cal potential can arise near the boundaries of seismo active regions and the Earth crust faults. Our estimates have shown that at the fair weather conditions they can reach the values 0.25 pA/m2 and 2.5 V/m, respectively. It appears that the horizontal field variations can be even greater than the above values by several orders of mag- nitude in the presence of strong background current that are typical of mesoscale convective systems and the atmospheric precipitation [Davydenko et al. 2004]. It can be also expected that the increased air ion- ization due to radon decay has an impact on the aero- electrical structures thereby producing the pre-seismic changes in spectral indices and amplitudes of ULF elec- tric pulsations resulted from motion of these struc- tures. A closer look at the electric pulsation spectra before seismic events would help us to understand whether this pre-seismic effect takes place or not. If this is the case, the radon-produced air ionization can serve as a plausible candidate for explanation of this effect. We have analyzed the hypothesis by Pulinets and Ouzounov [2011] that the enhancement of radon re- lease from the ground deeply influences a near-surface process of evaporation and related phenomena such as abnormal variations of upward thermal IR radiation occasionally observed from satellite before earthquake occurrence. The big cluster ions (BCIs) are of great concern in this hypothesis. It should be realized that the properties of the BCI are unusual since they are out of standard classification of ions [e.g., Hirsikko et al. 2011] because of their big size (about 1 nm) and a great num- ber of the attached water molecules (about 1011). Tak- ing into account that the diameter of water molecule is 0.29 nm, it is easy to check that a total volume of these attached molecules are greater than a volume of the in- dividual BCI. In this notation the BCI should be assigned to water droplets rather than to cluster ions. The con- ditions for the generation of BCI are unknown and, perhaps, the BCI seldom if ever occurs in the actual at- mosphere. However the major problem of this hypothesis is that the production rate of the BCI needs to be as high as and even greater than that of small ions in order to explain the observed data. According to our estimates, the air molecule ionization by atmospheric radon decay must lead to the generation of approximately 109 BCI per unit volume and per unit time in order to produce the upward thermal flux about 16 W/m2. This implies that every primary ion originated from alpha particle scattering would end up as the BCI. Besides the atmos- pheric radon activity concentration used by Pulinets and Ouzounov [2011] is two orders of magnitude greater than the actual one. Taken together, it makes this hypothesis very questionable. We can conclude that, in light of our study, the change in the air resistance of surface atmospheric layer by several tens percent is currently the only plausible radon effect which can be predicted with confidence. Although one might expect certain changes in the elec- tric structure of the lower atmosphere. A network of ground-based electric sensors covered seismo-active re- gion is necessary in order to monitor the electrical field perturbations due radon anomaly. Despite these effects seem to be sensitive to meteorological conditions, the electrical perturbations can be detectable at the ground level if it is granted that daily and seasonal variations of fair-weather atmospheric currents and electric fields are eliminated. It appears that further experiments are necessary to provide us with evidence both in support of and against a concept that other pre-seismic radon effect actually exists. PRE-SEISMIC VARIATIONS OF RADON Acknowledgements. This study was supported by the grant no. 13-05-12091 from RFBR. References Anisimov, S.V., E.A. Mareev, N.M. Shikhova and E.M. Dmitriev (2002). Universal spectra of electric field pulsations in the atmosphere, Geophys. Res. Lett., 29 (24), 2217-2220; doi:10.1029/2002GL015765. Anisimov, S.V., E.A. Mareev, N.M. Shikhova, A.E. Sorokin and E.M. Dmitriev (2005). On the electro- dynamical characteristics of the fog, Atmosph. Res., 76 (1-4), 16-28. Chalmers, J.A. (1967). Atmospheric electricity, 2nd ed., Pergamon press, New York. Cicerone, R.D., J.E. Ebel and J. Britton (2009). A sys- tematic compilation of earthquake precursors, Tectonophysics, 476, 371-396. Davydenko, S.S., E.A. Mareev, T.C. Marshall and M. Stolzenburg (2004). Calculation of electric fields and currents of mesoscale convective systems and their influence on the global electric circuit, J. Geophys. Res., D109, 11103-11110; doi:10.1029/2003JD003832. Geller, R.J. (1997). Earthquake prediction: a critical re- view, Geophys. J. Int., 131, 425-450. Genzano, N., C. Aliano, R. Corrado, C. Filizzola, M. Lisi, G. Mazzeo, R. Paciello, N. Pergola and V. Tra- mutoli (2009). RST analysis of MSG-SEVIRI TIR ra- diances at the time of the Abruzzo 6 April 2009 earthquake, Nat. Hazards Earth Syst. Sci., 9, 2073- 2084. Gorny, V.I., B.B. Shilin and G.I. Yasinsky (1993). Ther- mal Aerospace Survey, Nedra, Moscow (in Russian). Gray, D.E. (1972). American Institute of Physics Hand- book, 3rd ed., McGraw-Hill Book Company, 2364 p. Harrison, R.G., K.L. Aplin and M.J. Rycroft (2010). At- mospheric electricity coupling between earthquake regions and the ionosphere, J. Atmos. Sol.-Terr. Phys., 72, 376-381. Harrison, R.G., K.L. Aplin and M.J. Rycroft (2013). Earthquake-cloud coupling through the global at- mospheric electric circuit, Nat. Hazards Earth Syst. Sci. Discuss., 1, 7271-7283; doi:10.5194/nhessd-1- 7271-2013. Harrison, R.G., K.L. Aplin and M.J. Rycroft (2014). Brief Communication: Earthquake-cloud coupling through the global atmospheric electric circuit, Nat. Hazards Earth Syst. Sci., 14, 773-777; doi:10.5194/nhess-14- 773-2014. Hirsikko, A., T. Nieminen, S. Gagné, K. Lehtipalo, H.E. Manninen, M. Ehn, U. Hõrrak, V.-M. Kerminen, L. Laakso, P.H. McMurry, A. Mirme, S. Mirme, T. Petäjä, H. Tammet, V. Vakkari, M. Vana and M. Kul- mala (2011). Atmospheric ions and nucleation: a re- view of observations, Atmos. Chem. Phys., 11, 767- 798; doi:10.5194/acp-11-767-2011. Inan, S., T. Akgül, C. Seyis, R. Saatçılar, S. Baykut, S. Ergintav and M. Bas (2008). Geochemical monitor- ing in the Marmara region (NW Turkey): a search for precursors of seismic activity, J. Geophys. Res., 113, B03401; doi:10.1029/2007JB005206. Israël, H. (1970). Atmospheric Electricity, Israel Pro- gram for Scientific Translations, Jerusalem, 1, 317 p. Ivlev, L.S., Yu.A. Dovgalyuk (1999). Physics of atmos- pheric aerosol systems, St Petersburg, St Petersburg State University, 194 p. (in Russian). Kabrt, F., C. Seidel, A. Baumgartner, H. Friedmann, F. Rechberger, M. Schuff and F.J. Maringer (2014). Radon soil gas measurements in a geological versa- tile region as basis to improve the prediction of areas with a high radon potential, Rad. Prot Dos., 160 (1-3), 217-221; doi:10.1093/rpd/ncu086. Kondratyev, K.Ya., L.S. Ivlev, V.F. Krapivin and C.A. Varotsos (2006). Atmospheric Aerosol Properties: Formation, Processes and Impacts, Springer/PRAXIS, Chichester, UK, 572 p. Ouzounov, D., M. Bryant, T. Logan, S. Pulinets and P. Taylor (2006). Satellite thermal IR phenomena as- sociated with some of the major earthquakes in 1999-2003, Phys. Chem. Earth, 31, 154-163. Pierce, E. (1976). Atmospheric electricity and earth- quake prediction, Geophys. Res. Lett., 3 (3), 185-188. Protection against Radon-222 at Home and at Work (1994). Annals of the ICRP, No. 65, Oxford, Perga- mon. Pulinets, S. (2011). The synergy of earthquake precur- sors, Earthq. Sci., 24, 535-548. Pulinets, S., and D. Ouzounov (2011). Lithosphere-At- mosphere-Ionosphere Coupling (LAIC) model - An unified concept for earthquake precursors valida- tion, J. Asian Earth Sci., 41, 371-382. Rycroft, M.J., S. Israelsson and C. Price (2000), The global atmospheric electric circuit, solar activity and climate change, J. Atmos. Sol.-Terr. Phys., 62, 1563- 1576. Rycroft, M.J., R.G. Harrison, K.A. Nicoll and E.A. Ma- reev (2008). An overview of Earth’s global electric circuit and atmospheric conductivity, Space Sci. Rev., 137 (1-4), 83-105; ISSN 0038-6308, doi: 10.1007/s112 14-008-9368-6. Rycroft, M.J., K.A. Nicoll, K.L. Aplin and R.G. Harrison (2012). Recent advances in global electric circuit cou- pling between the space environment and the tro- posphere, J. Atmos. Sol.-Terr. Phys., 90-91, 198-211. Surkov, V.V., O.A. Pokhotelov, M. Parrot and M. Hayakawa (2006). On the origin of stable IR anom- alies detected by satellite above seismo-active regions, SURKOV 8 9 Phys. Chem. Earth, 31, 164-171. Surkov, V.V., and M. Hayakawa (2007). ULF electromag- netic noise due to random variations of background atmospheric current and conductivity, J. Geophys. Res., 112, D11116; doi:10.1029/2006JD007788. Surkov, V.V., and M. Hayakawa (2008). Natural electro- magnetic ULF noise due to fluctuations of ionos- pheric currents, J. Geophys. Res., 113, A11310; doi:10. 1029/2008JA013196. Surkov, V., and M. Hayakawa (2014). Ultra and Ex- tremely Low Frequency Electromagnetic Fields, Springer Geophysics Series, XVI, Springer, 486 p. Toutain, J.-P., and J.-C. Baubron (1998). Gas geochem- istry and seismotectonics: a review, Tectonophysics, 304, 1-27. Tramutoli, V., D. Bello, G.N. Pergola and S. Piscitelli (2001). Robust satellite technique for remote sens- ing of seismically active areas, Annali di Geofisica, 44, 295-312. Tronin, A.A., M. Hayakawa and O.A. Molchanov (2002). Thermal IR satellite application for earth- quake research in Japan and China, J. Geodynamics, 33, 519-534. Tverskoy, P.N. (1962). Meteorology (Physics of the at- mosphere), Gydrometeoizdat, Leningrad (in Rus- sian). Virk, H.S., and B. Singh (1994). Radon recording of Ut- tarkashi earthquake, Geophys. Res. Lett., 21, 737- 742. Wåhlin, L. (1994). Elements of fair weather electricity, J. Geophys. Res., 99, D10767-D10772. Yasuoka, Y., Y. Kawada, H. Nagahama, Y. Omori, T. Ishikawa, S. Tokonami and M. Shinogi (2009). Pre- seismic changes in atmospheric radon concentration and crustal strain, Phys. Chem. Earth, 34, 431-434. Yerg, D.G., and K.R. Johnson (1974). Short-period fluc- tuations in the fairweather electric field, J. Geophys. Res., 79, 2177-2184. Zhang, K., J. Feichter, J. Kazil, H. Wan, W. Zhuo, A. D. Griffiths, H. Sartorius, W. Zahorowski, M. Ramonet, M. Schmidt, C. Yver, R. E. M. Neubert and E.-G. Brunke (2011). Radon activity in the lower tropo- sphere and its impact on ionization rate: a global es- timate using different radon emissions, Atmos. Chem. Phys., 11, 7817-7838; doi:10.5194/acp-11- 7817-2011. Corresponding author: Vadim Surkov, National Research Nuclear University MEPhI, Moscow, Russia; email: surkovvadim@yandex.ru. © 2015 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. PRE-SEISMIC VARIATIONS OF RADON << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice