Radon Signals Research Articles

Radon signals in soil gas associated with earthquake occurrence in Greece: review and perspective

Radon signals in soil gas associated with earthquake occurrence in Greece: review and perspective

A sensitive system based on radon amplification at soil-air interface: Aiming to advance earthquake precursor research

Radon, a natural radioactive gas, serves as a valuable tracer in geophysical research and atmospheric science such as detecting stress induced signal in bedrock. However, the conventional radon monitoring methods often lack the sensitivity required to accurately capture such signals. This limitation, coupled with interference from meteorological effects, poses challenges in distinguishing genuine stress-induced signals. In this study, we propose a novel approach utilizing radon concentration gradients at the soil-air interface to enhance sensitivity and detect stress induced radon signals more effectively. Drawing from pressure diffusion models, we demonstrate how seismic stress accumulation in bedrock alters radon profiles in the sub-soil, providing insights into the mechanisms underlying stress-induced radon variations. Building upon this theoretical framework, we introduce the “Bhabha Radon Observatory for Seismic Application (BhaROSA)," a remote sensing, solar-powered radon observatory designed for widespread deployment and continuous unattended monitoring for big database generation. Field experiments comparing BhaROSA's performance to conventional soil probe techniques validate and confirm the superior sensitivity in line with theoretical predictions. This innovative approach holds promise for improving our understanding of stress dynamics in bedrock and has potential applications in various geophysical and atmospheric science such as earthquake precursory research, geo-genic radon potential and risk assessment. To progress, we propose international alliance and application of deep learning to a big database of precursor signals, which may lead to more informed conclusions on earthquake predictability-an enduring and unsolved challenge for humanity.

The dynamics of Rn-222 cyclic flow within the shallow geological subsurface media as a daily temporal variated source for exhalation into the air

Extensive research on the dynamics of radon gas (Rn-222) originating from the radioactive decay of radium (Ra-226) in geological subsurface media, sheds light on its periodic release into the atmosphere. Radon is a product of the uranium-238 decay chain found within rock and soil grains. While only a fraction of the generated radon escapes (emanates) into porous spaces due to nuclear recoil, it serves as the source for subsurface gas flows and for cyclic exhalation into the soil-atmosphere interface.Ongoing study of radon movement in shallow and deep subsoil, and its emergence at the surface, reveals complete semi-diurnal, diurnal, and seasonal gas flow cycles in the subsoil. Complementary emissions occur nocturnally as radon is released into the atmosphere. Moreover, two natural driving forces govern complex semi-diurnal and diurnal flows below and above the surface. Subsurface gas movement in porous media exhibits nonlinear behavior influenced by surface temperature gradients, resulting in downward flow to depths of up to 100 m. This flow exhibits daily periodicity with depth-dependent time delays, correlating with the diurnal surface temperature cycle. Additionally, pore gas transport into and out of open boreholes responds linearly to semi-diurnal barometric pressure changes, known as barometric pumping.Beyond subsurface phenomena, Europe and Australia increasingly employ nocturnal radon measurements to study atmospheric stability and air quality, assuming that variations in local stationary near-surface radon concentrations reflect atmospheric mixing processes. Recognizing mechanisms governing radon's temporal changes within geological subsurface media highlights the need for continuous underground radon monitoring to validate variations in daily radon exhalation to the surface.On the other hand, monitoring radon at considerable depths minimizes climatic contributions and enhances the ability to discern non-periodic pre-seismic radon signals, independent of atmospheric compulsion. This research offers potential insights into seismic precursors and the complex interplay between subsurface geodynamics and atmospheric conditions.

Time-series analysis of radon monitoring in soil gas in association with earthquakes in Stivos faulting, at Lagadas basin, North Greece

AbstractTime series analysis was applied to the continuous radon level, temperature, pressure, and rainfall to find clear earthquake signals. Radon signals appeared a few days after heavy rains, and radon signals associated with events M = 3.8–4.2 were detected 12 up to 36 days before. The events are complete data recorded from 1983 to 1986, giving discussion and conclusion on M with prediction time and radon anomaly detected in the Stivos faulting near Thessaloniki, N. Greece.

A preliminary study on real-time Rn/Tn discriminative detection using air-flow delay in two ion chambers in series

Due to its short half-life, thoron gas has been assumed to have negligible health hazards on humans compared to radon. But, one of the decay products with a long half-life can make it to be transported to a long distance and to cause a severe internal dose through respiration. Since most commercial radon detectors can not discriminate thoron signals from radon signals, it is very common to overestimate radon doses which in turn result in biased estimation of lung cancer risk in epidemiological studies. Though some methods had been suggested to measure thoron and radon separately, they could not be used for real-time measurement because of CR-39 or LR-115. In this study, an effective method was suggested to measure radon and thoron separately from the free air. It was observed that the activity of thoron decreases exponentially due to delay time caused by a long pipe between two chambers. Therefore from two ion chambers apart in time, it was demonstrated that thoron and radon could be measured separately and simultaneously. We also developed a collimated alpha source and with this source and an SBD, we could convert the ion chamber reading to count rate in cps.

Analysis of Radon Measurements in Relation to Daily Seismic Activity Rates in the Vrancea Region, Romania.

Many previous research studies have shown how local and even regional earthquakes can significantly affect the release of radon in the soil. The aim of this work is to investigate the relationship between radon measurements and the daily seismic activity rate and develop a methodology that allows estimating the seismic activity rate using only radon measurements. To carry out this study, the earthquake catalogue of the Vrancea region (Romania) has been used to estimate the daily seismic activity rate during a given time period, in which radon measurements were also recorded, from January 2016 to September 2020. The Vrancea zone represents the most active seismic zone in Europe and is located on the eastern edge of the strongly bent Carpathian arc. In the case of the radon measurements, seasonal behaviours and linear trends due to non-seismic factors have been identified and subsequently removed. The discrete wavelet transform has been used to analyse the radon signal at two different scales: long and short periods. From the analysis carried out on a long-period scale, an approximate linear relationship has been obtained between the radon series and the daily seismic activity rate, which provides insights into the behaviour of the seismic activity in the study region with only the radon information. In addition, the study reveals certain characteristics that could be used as precursors of earthquakes at different scales: weeks in the case of the estimated daily seismic activity rate, and days in the case of the short-period signal obtained by the wavelet analysis. The results obtained for this region allow us to hope that the analysis of the radon time series can become an effective complement to the conventional seismic analysis used in operational earthquake forecasting.

Study of Surface Emissions of 220Rn (Thoron) at Two Sites in the Campi Flegrei Caldera (Italy) during Volcanic Unrest in the Period 2011–2017

The study concerns the analysis of 220Rn (thoron) recorded in the surface soil in two sites of the Campi Flegrei caldera (Naples, Southern Italy) characterized by phases of volcanic unrest in the seven-year period 1 July 2011–31 December 2017. Thoron comes only from the most surface layer, so the characteristics of its time series are strictly connected to the shallow phenomena, which can also act at a distance from the measuring point in these particular areas. Since we measured 220Rn in parallel with 222Rn (radon), we found that by using the same analysis applied to radon, we obtained interesting information. While knowing the limits of this radioisotope well, we highlight only the particular characteristics of the emissions of thoron in the surface soil. Here, we show that it also shows some clear features found in the radon signal, such as anomalies and signal trends. Consequently, we provide good evidence that, in spite of the very short life of 220Rn compared to 222Rn, both are related to the carrier effect of CO2, which has significantly increased in the last few years within the caldera. The hydrothermal alterations, induced by the increase in temperature and pressure of the caldera system, occur in the surface soils and significantly influence thoron’s power of exhalation from the surface layer. The effects on the surface thoron are reflected in both sites, but with less intensity, the same behavior of 222Rn following the increasing movements and fluctuations of the geophysical and geochemical parameters (CO2 flux, fumarolic tremor, background seismicity, soil deformation). An overall linear correlation was found between the 222−220Rn signals, indicating the effect of the CO2 vector. The overall results represent a significant step forward in the use and interpretation of the thoron signal.

Singular spectrum analysis on soil radon time series (222Rn) in Kachchh, Gujarat, India: detection of periodic oscillations and earthquake precursors

The influence of various meteorological parameters on the observed time series of soil radon (222Rn) in the form of periodic oscillations is profound and well established. An attempt is made in the present study to identify and remove the periodic oscillations using singular spectrum analysis (SSA) to identify earthquake precursors in the time series of soil radon (222Rn) recorded at Badargadh, Kachchh region of Gujarat, India, for the period from February 21 to June 10, 2011. The temperature has shown diurnal oscillations, whereas radon and pressure show the mixture of diurnal and semi-diurnal oscillations. A Butterworth digital filter has been applied to the observed data of all parameters (222Rn, pressure, temperature, and humidity) to remove the long periodic trends (more than 30 days). The periodic and aperiodic components in the time series of soil radon are identified and separated using SSA. While correlating the radon with other parameters, we found that the aperiodic radon has negative correlation with temperature and the radon has positive correlation with humidity and pressure. Two earthquakes of magnitudes M3.7 and M4.2 which occurred on March 26, 2011, and May 17, 2011, respectively, are being investigated to identify anomalies before these earthquakes. These earthquakes are located within a distance of 100 km from Badargadh. Based on the statistical analysis of mean+2*sigma, we observed the anomalous radon signal before the occurrence of these earthquakes with relative percentage of 94–96% at Badargadh.

Long‐Term Monitoring and Characterization of Soil Radon Emission in a Seismically Active Area

AbstractRadon (222Rn) is a radioactive gas, which originates everywhere in rocks of Earth's crust through a process of natural decay of other radioactive elements. The monitoring of soil radon can be useful to address a multitude of geological and environmental issues. However, some difficulties arise when trying to use radon as earthquake precursor because the earthquake‐related anomaly cannot be easily and univocally discriminated from other anomalies of different origin. To explore the relationship between soil radon emission and seismicity, a local network is operating in southeastern Sicily (Italy). Its peculiarities are the uniform geological condition of the monitoring sites and the long data record with simultaneous measures of radon and the main climate variables. In this paper, we applied continuous wavelet transformation to a ≈ 3.5 years long soil radon time series to detect periodic variations. Results indicate the occurrence of cycles at annual, semiannual, and monthly periodicity, which are ascribable to the effects of the climatic‐environmental parameters. The periodic components have been modeled and the signal conveniently filtered. We show as the characterization of the long‐term behavior of radon signals is essential to recognize anomalies in radon emission, which can be related to geological/environmental phenomena. The methodology proposed in this work provides a reliable characterization of the radon time series and can be applied at various spatial/temporal scales, depending on the scientific objectives to be achieved. This approach can also represent the base for further analysis like the investigation of the modulation between the periodic components and short‐term forecasts.

Continuous radon monitoring during seven years of volcanic unrest at Campi Flegrei caldera (Italy)

This is a seven-year study (1/7/2011-31/12/2017) of radon monitoring at two sites of Campi Flegrei caldera (Neaples, Southern Italy) that in the last 70 years experienced repeated phases of volcanic unrest. The sites are equipped with devices for radon detection, based on the spectrometry analysis of the α-particles of radon daughters. A hybrid method, as combination of three known methods, is applied for the identification of residuals (anomalies) and trends of the time series of Radon. The results are compared with the following indicators of current caldera unrest: the tremor caused by the major fumarolic vent registered by a seismic station; the cumulative of background seismicity; the maximum vertical deformation acquired by GPS networks during the current phase of uplift; the temperature-pressure of the hydrothermal system estimated based on gas geo-indicators. The comparisons show strong correlation among independent signals and suggest that the extension of the area affected by current Campi Flegrei crisis is larger than the area of seismicity and of intense hydrothermal activity from which the radon stations are 1–4 km away. These results represent an absolute novelty in the study of a such calderic area and mark a significant step forward in the use and interpretation of the radon signal.

Detecting time series anomalies using hybrid methods applied to Radon signals recorded in caves for possible correlation with earthquakes

Anomalies in the Radon activity concentration time series recorded in five European caves (Czech Republic, Slovakia, Slovenia) are detected using three hybrid methods: (1) multiple linear regression and autoregressive integrated moving average statistical methods, (2) Empirical Mode Decomposition with Support Vector Regression techniques and (3) the Singular Spectrum Analysis composed with a predicting methodology. Results coming from the three methods are compared and the best hybrid method is selected based on statistical evaluation criteria of the uncertainty. Radon anomalies occur ± 30 days from earthquake occurrence, selected according to the Dobrovolsky’s earthquake preparation zone formula and to seismic events (with magnitude ≥ 4) occurred in the neighboring European Countries to the monitoring caves. The anomalies detection furnishes results consistent across the used methodologies, as proven by the calculation of a statistical parameter that search the presence of anomalies coming from the hybrid methods within ± 30 days from earthquake event.

Testing of 222Rn application for recognizing tectonic events observed on water-tube tiltmeters in underground Geodynamic Laboratory of Space Research Centre at Książ (the Sudetes, SW Poland)

Research on relationships between variation in 222Rn activity concentration and tectonic events recorded using the instruments of the Geodynamic Laboratory of SRC PAS at Książ (the Sudetes, SW Poland) had been conducted since 2014. The performed analyses of variation have demonstrated the spatial character of changes in 222Rn activity concentration. Their time-course is comparable in all parts of the underground laboratory. This means that gas exchange between the lithosphere and the atmosphere occurs not only through fault zones but also through all surfaces of the underground workings: the floors, the sidewalls and the roofs. Further, some relationships between 222Rn activity concentration and tectonic activity of the orogen have been demonstrated with the use of Pearson's linear correlation coefficient. The comparison between temporal distribution (times series) of radon activity concentration and water-tube tiltmeters (WTs) demonstrated that radon data have regular oscillations which can be approximated using the sine function with a 12 month cycle (seasonal changes) and amplitude in the range of 1000–1500 Bq/m3. To compare the collected radon signal data and tectonic activity, we used linear function as the simplest method of trend assessment. Pearson's correlation coefficient r cannot be accepted as appropriate for assessing the interdependencies between variables because they do not have a normal distribution, and the relationship between them is not linear.It was noted that each series of data, namely radon activity concentration and tectonic activity determine the series of deviations above and below the trend function. Because of the non-fulfillment of the above assumptions, we used nonparametric equivalents such as Spearman's rank correlation coefficient rs and Kendall's tau. The obtained results showed that the value of the rs coefficient ranges from 0.38 to even 0.43. The best relationship at the level of rs = 0.43 was determined between the radon activity concentration recorded by detector no. 3 and the tectonic activity of the rock mass registered on the WT-2 channel. Similar at the rs level of 0.37–0.38 between detector no. 5 and 4 and the WT-2 channel. A bit higher than rs = 0.39 between detector no. 3 and the WT-2 channel. In each case, these were positive correlations. The obtained Spearman's rs coefficients indicate the correlation between 222Rn activity concentration and tectonic activity of the rock mass. The t-statistic, which analyzes the significance of Spearman's coefficient rs is a descriptive measure of the accuracy of regression matching to empirical data. It takes values in the range of percentage and provides informations about which part of the total variability of the radon activity concentration (Y) observed in the sample has been explained (determined) by regression in relation to tectonic activity of the rock mass (X). In our case, approximately f 40% to more than 50% of the radon activity concentration (Y) was explained by regression in relation to the tectonic activity of the rock mass. We obtained similar results with the use of Kendall's tau coefficient.Precise description of the character of this relationship requires further, more detailed analyses, such as comparing characteristics of the distributions based on trend variation like Monte Carlo simulation, Multivariate Adaptive Regression Splines or neural networks.

Radon-based atmospheric stability classification in contrasting sub-Alpine and sub-Mediterranean environments

A recently-developed radon-based technique is used to investigate relative changes in summertime atmospheric stability at two sites in Slovenia with contrasting geographical settings. Although atmospheric stability for both sites (50 km apart) was shown to be governed by similar synoptic conditions, their contrasting settings caused differences in mixing conditions for each stability category. At the urban sub-Alpine site Ljubljana, situated within a topographic basin, wind speeds associated with the most stable conditions were 0.2–0.3 m s−1. By comparison, corresponding wind speeds for the near-coastal sub-Mediterranean site Ajdovščina, located at the foothills of the Trnovski gozd barrier, were 0–0.2 m s−1. The wind direction at Ljubljana under stable conditions (∼80°) was consistent with drainage flow into the basin along the Sava River valley. The corresponding wind direction at Ajdovščina was 20–40°, consistent with gentle katabatic drainage from the flanks of the Trnovski gozd barrier. After removing fetch effects on radon variability at each site, a large contrast in local contributions to the radon signal was noted: the diurnal amplitude of the local radon signal increased from ∼24 Bq m−3 at Ljubljana to ∼47 Bq m−3 at Ajdovščina. This difference was attributed to a greater nocturnal radon accumulation rate at Ajdovščina (3.5 Bq m−3 h−1 vs 2.1 Bq m−3 h−1) due to higher radon fluxes from flysch and carbonate rocks compared to the sea and lake sediments in the Ljubljana Basin. The ability of radon to consistently distinguish subtle changes in atmospheric mixing at sites with contrasting topographic settings indicates that it will be a powerful tool for characterising air quality in these complex environments. Specifically, diurnal radon cycles indicate that the capability of the atmosphere to dilute primary pollutants is considerably less in the basin environment.

Radon signals in geological (natural) geogas and in a simultaneous enhanced confined mode simulation experiment.

An enhanced confined mode (ECM) radon simulation experiment, tested in the laboratory in Jerusalem, was relocated to a subsurface geophysical observatory located 400 km apart, at a depth of 150 m and with a stable temperature. Five gamma sensors are placed around the ECM canister and lead shielding minimizes the influence of natural local gamma radiation. Simultaneous measurement of the geological radon and from radon in the ECM system indicates that the temporal variation of gamma radiation from radon in the ECM system contains annual, multi-day and daily signals, that correspond to signals in the local geological radon. This implies that a common external driver influences the radiation pattern of the geological radon and from radon inside the ECM canister. Once activated at BGO the typical variation pattern of the experimental system occurring in the laboratory changed to that occurring at the observatory. This is interpreted to indicate that the overall style of the temporal patterns of radiation from radon is site dependent. The outcome of this investigation conforms and further substantiates the recent suggestion that a component in solar radiation is driving the annual and daily periodic components in the variation of radon. New geophysical research potential is indicated.

The imprint of thermally induced devolatilization phenomena on radon signal: implications for the geochemical survey in volcanic areas

The imprint of thermally induced devolatilization phenomena on radon signal: implications for the geochemical survey in volcanic areas

Transient radon signals driven by fluid pressure pulse, micro-crack closure, and failure during granite deformation experiments

In seismically active fault zones, various crustal fluids including gases are released at the surface. Radon-222, a radioactive gas naturally produced in rocks, is used in volcanic and tectonic contexts to illuminate crustal deformation or earthquake mechanisms. At some locations, intriguing radon signals have been recorded before, during, or after tectonic events, but such observations remain controversial, mainly because physical characterization of potential radon anomalies from the upper crust is lacking.Here we conducted several month-long deformation experiments under controlled dry upper crustal conditions with a triaxial cell to continuously monitor radon emission from crustal rocks affected by three main effects: a fluid pressure pulse, micro-crack closure, and differential stress increase to macroscopic failure. We found that these effects are systematically associated with a variety of radon signals that can be explained using a first-order advective model of radon transport. First, connection to a source of deep fluid pressure (a fluid pressure pulse) is associated with a large transient radon emission increase (factor of 3–7) compared with the background level. We reason that peak amplitude is governed by the accumulation time and the radon source term, and that peak duration is controlled by radioactive decay, permeability, and advective losses of radon. Second, increasing isostatic compression is first accompanied by an increase in radon emission followed by a decrease beyond a critical pressure representing the depth below which crack closure hampers radon emission (150–250 MPa, ca. 5.5–9.5 km depth in our experiments). Third, the increase of differential stress, and associated shear and volumetric deformation, systematically triggers significant radon peaks (ca. 25–350% above background level) before macroscopic failure, by connecting isolated cracks, which dramatically enhances permeability.The detection of transient radon signals before rupture indicates that connection of initially isolated cracks in crustal rocks may occur before rupture and potentially lead to radon transients measurable at the surface in tectonically active regions. This study offers thus an experimental and physical basis for understanding predicted or reported radon anomalies.

Radon surveys and monitoring at active volcanoes: learning from Vesuvius, Stromboli, La Soufrière and Villarrica

Abstract Understanding the behaviour of fluids in hydrothermal systems is a key factor in volcano monitoring. Measuring gas emissions in volcanic areas is strategic for detecting and interpreting precursory signals of variations in volcanic activity. The role of radon as a potential precursor of earthquakes has been extensively debated. However, radon anomalies appear to be better suited to forecast eruptive episodes as we know the loci of volcanic eruptions and we can follow the evolution of volcanic activity. Radon mapping is an effective tool in assessing diffuse and concentrated degassing at the surface. We hereby summarize the in-soil radon emissions collected worldwide and further discuss a collection of data on our key targets. These are closed-conduit and open-conduit volcanoes: Vesuvius (Italy) and La Soufrière (Guadeloupe, Lesser Antilles), Stromboli (Italy) and Villarrica (Chile), respectively. In all the above volcanoes, faults and fracture systems control radon degassing. Automatic and real-time measurements help us to detect major changes in volcanic activity. We present and discuss the radon time series associated with the last effusive eruption at Stromboli. Spectral analyses reveal diurnal and semi-diurnal cycles being probably modulated by atmospheric variations. Multiple linear regression (MLR) analyses have been performed by filtering the radon signals from the effects of local environmental parameters. The residuals do not show particular variations or precursory peaks as the gases have been released from this open-conduit volcano before the onset of the effusive phase (7 August 2014). It is finally emphasized that radon is not the sole precursor, and we should also rely on other geochemical and geophysical parameters. In this perspective, we propose a methodological procedure that can contribute to improving volcano surveillance in an attempt to mitigate volcanic risk.

Indications for influence of artificial (man-made) activity on radon signals, in simulation experiments

Radon (Rn-222; a radioactive noble gas) is characterized by large temporal variations that differ significantly from variations of (i) other trace elements in geogas (noble gases); (ii) variation patterns of other dynamic geophysical systems (atmospheric, tidal). Consensus exists that there is no simple and straightforward understanding of the phenomena and its behaviour. This lacuna in the understanding of the underlying principles hampers the development of applications—such as radon as a proxy of processes in the seismogenic context. Using results from field investigations and simulation experiments the GSI suggested that an unidentified extraterrestrial component, probably in solar radiation, drives periodic radon signals in the diurnal and annual frequency bands. Recent findings from experimental investigations shed additional perspectives allowing a new evaluation of the issue. Particular transient signals, measured with alpha and gamma detectors, are interpreted to reflect the influence of artificial activity. Criteria are (i) signals lasting several hours that occur around midday on workdays (Sunday–Thursday); (ii) signals composed of a train of around 10 strong pulses, each lasting less than 15 min, occurring within several hours once a week, from Wednesday afternoon/evening to Thursday morning. A first interpretation is that an unidentified artificial activity of some sort (industrial?) generates and emits an unidentified agent that reaches enhanced confined mode experiments at the GSI laboratory, which respond to the incoming agent in the form of radon signals. Developing the capability of identification of such an earth-bound source generating an influencing agent is a key step towards understanding of external influence on radioactivity of radon.

Establishing the range of background for radon variations in groundwater along the Serghaya fault in southwestern Syria

Long-term groundwater radon measurements were carried out in two selected karstic springs emerging from the Serghaya fault zone in southwestern Syria. The work is aimed at determining the range of radon background along the concerned fault, which represents a prominent branch segment of the Dead Sea Fault System in the region. The obtained data was statistically analyzed and the mean radon values have been standardized in terms of probability of magnitude in order to enable the separation between normal radon variations from other anomalous or geotectonic related values. The results revealed a general background range of radon concentration varying between 5000 and 13000 Bqm-3, and thus all radon values lying outside this range were considered anomalous. However, remarkable clusters of high radon concentrations with peak values more than twice times the background level, were observed in both monitoring springs through the period (1992-1994). These abnormal radon values were positively correlated with simultaneous time of high precipitations, while no significant earthquake activities were recorded in the region during that period. Therefore, such radon signals do not seem to be a reflection of regional changes in crustal strain, but rather they mostly indicate evidences of radon response to the groundwater table fluctuations due to high precipitations. These consequences may confirm the usefulness of radon application as an important natural tracer in hydrogeological investigations.

Delayed signatures of underground nuclear explosions.

Radionuclide signals from underground nuclear explosions (UNEs) are strongly influenced by the surrounding hydrogeologic regime. One effect of containment is delay of detonation-produced radioxenon reaching the surface as well as lengthening of its period of detectability compared to uncontained explosions. Using a field-scale tracer experiment, we evaluate important transport properties of a former UNE site. We observe the character of signals at the surface due to the migration of gases from the post-detonation chimney under realistic transport conditions. Background radon signals are found to be highly responsive to cavity pressurization suggesting that large local radon anomalies may be an indicator of a clandestine UNE. Computer simulations, using transport properties obtained from the experiment, track radioxenon isotopes in the chimney and their migration to the surface. They show that the chimney surrounded by a fractured containment regime behaves as a leaky chemical reactor regarding its effect on isotopic evolution introducing a dependence on nuclear yield not previously considered. This evolutionary model for radioxenon isotopes is validated by atmospheric observations of radioxenon from a 2013 UNE in the Democratic People’s Republic of Korea (DPRK). Our model produces results similar to isotopic observations with nuclear yields being comparable to seismic estimates.