Modelling the effect of air exchange on 222Rn and its progeny concentration in a tunnel atmosphere

Radon (222Rn) as a tracer in natural ventilation efficiency assessment in underground workings - an example of "St John Mine" tourist complex in Krobica (the Sudetes, SW Poland)

The authors characterize the use of ²²²Rn as an effective tracer of natural ventilation of an underground site where air circulates within a whole system of workings and ventilation intensity (the number of air exchanges in the space) is determined by atmospheric factors. A radon-related database containing results of measurements conducted at various intervals and at different stages of site accessibility was compiled. During 8 months of the calendar year ²²²Rn activity concentration exceeds the mean annual reference value established by Polish law (300 Bq/m³). These months correspond to periods with low intensity of natural ventilation of the workings and reduced efficiency of air exchange between the site and the atmosphere. They occur in autumn - in the second half of September, in October and November, and in May in spring, and persist for 7 to even 14 days. During these periods, the time spent inside the facility which is considered safe in terms of radiation protection is limited to an average of 6-8 h a day, i.e. from 6 a.m. to 6 p.m. in October, from 11 a.m. to 6 p.m. in November and from 11 a.m. to 5 p.m. in May. The length of a safe stay in the facility is determined by atmospheric factors, mainly the air and ground temperature. The concentrations of other gases in the atmosphere inside the facility comply with Polish mining regulations.

The results of long-term simultaneous measurements of radon exhalation rate, radon concentrations in soil gas and groundwater in the fault zone

The regular monthly radon measurements were carried out in the fault zone on the Western slope of the Beshtau magmatic massif (North Caucasus). The radon exhalation rate from the soil surface, as well as radon concentrations in soil gas at a depth of 0.5 m and in groundwater discharged at a spring located nearby have simultaneously been measured. High seasonal fluctuations in radon exhalation and radon concentration in soil gas, characterized by highs in summer and lows in winter, were registered. In summer, the radon exhalation reached 23.8 Bq m^-2s^-1, and the radon concentration in the soil gas reached 166 kBq m^-3. In winter, both the radon exhalation and the radon concentration in the soil dropped to 0.025 Bq m^-2s^-1 and <3 kBq m^-3, respectively. The concentration of radon in ground water varied over the year in a relatively narrow range (100-210 Bq l^-1), and there were no seasonal fluctuations. A sharp increase in soil radon and radon exhalation in spring and a fall in autumn are timed to the moments when the temperature of the atmospheric air becomes, respectively, higher and lower than the temperature of the rock massif. Both the soil radon concentration and the radon exhalation show a close correlation with the temperature of atmospheric air, but in the first case the relationship is linear, and in the second - exponential. The obtained data confirm the assumption that the seasonal radon variations are caused by atmospheric air circulation in the shallow area of the fault due to the temperature difference between the atmosphere and the rock massif.

The role of atmospheric conditions in CO2 and radon emissions from an abandoned water well

Boreholes and wells are complex boundary features at the earth-atmosphere interface, connecting the subsurface hydrosphere, lithosphere, and biosphere to the atmosphere above it. It is important to understand and quantify the air exchange rate of these features and, consequently their contribution as sources for greenhouse gas (GHG) emissions to the atmosphere. Here, we investigate the effect of atmospheric conditions, namely atmospheric pressure and temperature, on air, CO₂, and radon transport across the borehole-ambient atmosphere interface and inside a 110-m deep by 1-m diameter borehole in northern Israel. Sensors to measure temperature, relative humidity, CO₂, and radon were placed throughout a cased borehole. A standard meteorological station was located above the borehole. Data were logged at a high 0.5-min resolution for 9 months. Results show that climatic driving forces initiated 2 different advective air transport mechanisms. (1) Diurnal and semidiurnal atmospheric pressure cycles controlled daily air transport events (barometric pumping); and (2) There was a correlation between borehole-atmosphere temperature differences and transport on a seasonal scale (thermal-induced convection). Barometric pumping was identified as yielding higher fluxes of vadose zone gases than thermal-induced convection. Air velocities inside the borehole and CO₂ emissions to the atmosphere were quantified, fluctuating from zero up to ~6 m/min and ~5 g-CO₂/min, respectively. This research revealed the mechanisms involved in the process throughout the year and the potential contribution role played by boreholes to GHG emissions.

Hydrogeological control on carbon dioxide input into the atmosphere of the Chauvet-Pont d'Arc cave

Carbon dioxide (CO₂) concentration (CDC) is an essential parameter of underground atmospheres for safety and cave heritage preservation. In the Chauvet cave (South France), a world heritage site hosting unique paintings dated 36,000 years BP, a high-sensitivity monitoring, ongoing since 1997, revealed: 1) two compartments with a spatially uniform CDC, a large volume (A) (40,000 to 80,000 m³) with a mean value of 2.20 ± 0.01% vol. in 2016, and a smaller remote room (B) (2000 m³), with a higher mean value of 3.42 ± 0.01%; 2) large CDC annual variations with peak-to-peak amplitude of 2% and 1.6% in A and B, respectively; 3) long-term changes, with an increase of CDC and of its annual amplitude since 1997, then faster since 2013, reaching a maximum of 4.4% in B in 2017, decreasing afterwards. While a large effect of seasonal ventilation is ruled out, monitoring of seepage at two dripping points indicated that the main control of CDC seasonal reduction was transient infiltration. During periods of water deficit, calculated from surface temperature and rainfall, CDC systematically increased. The carbon isotopic composition of CO₂, correlated with water excess, is consistent with a time-varying component of CO₂ seeping from above. The CO₂ flux, which is the primary driver of CDC in A and B, inferred using box modelling, was found to confirm the relationship between water excess and reduced CO₂ flux into A, compatible with a more constant flux into B. A buoyancy-driven horizontal CO₂ flow model in the vadose zone, hindered by water infiltration, is proposed. Similarly, pluri-annual and long-term CDC changes can likely be attributed to variations of water excess, but also to increasing vegetation density above the cave. As CDC controls the carbonate geochemistry, an increased variability of CDC raises concern for the preservation of the Chauvet cave paintings.

Variation in natural short-period ionospheric noise, and acoustic and gravity waves revealed by the amplitude analysis of a VLF radio signal on the occasion of the Kraljevo earthquake (Mw = 5.4)

We analyse the lower ionosphere disturbances in the time period around the Mw 5.4 Kraljevo earthquake (EQ), which occurred on 3 November 2010 in Serbia. The results presented herein are based on analysis of the amplitudes of three VLF signals emitted in Italy, UK, and Germany and recorded in Serbia whose variations primarily result from changes in the electrical properties of the lower ionosphere at a distance more than 120 km from the epicentre of the EQ. The primary goals of this study are to reveal specific variations as possible EQ precursors as well as disturbances following the EQ event recorded by the observational equipment, and to investigate whether better time resolution data can affect the analysis of the lower ionosphere disturbances possibly connected to the EQ. We process two sets of data with sampling periods of 1 min and 0.1 s. As the first analysis indicates the absence of long-term disturbances, which can clearly be connected to the Kraljevo EQ, our attention is focused on the study of short-period noise amplitude and the excitation and attenuation of acoustic and gravity waves in the lower ionosphere. Processing of the amplitudes of three VLF signals during the nights of the four EQs with magnitude greater than 4 that occurred in Serbia, as well as EQs of all magnitudes during the three days of 3, 4, and 9 November, indicates that the detected ICV radio signal noise amplitude reduction starting before the Kaljevo EQ is also observed for 13 additional EQ events near the signal propagation path, and occurred during all three days (for all EQs with magnitude greater than 4 and several less intensive events). Excitation and attenuation of acoustic waves are also found for all these EQ events with a magnitude greater than 4.

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 ²²²Rn 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 ²²²Rn 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 ²²²Rn 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/m³. 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 r_s and Kendall's tau. The obtained results showed that the value of the r_s coefficient ranges from 0.38 to even 0.43. The best relationship at the level of r_s = 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 r_s level of 0.37-0.38 between detector no. 5 and 4 and the WT-2 channel. A bit higher than r_s = 0.39 between detector no. 3 and the WT-2 channel. In each case, these were positive correlations. The obtained Spearman's r_s coefficients indicate the correlation between ²²²Rn activity concentration and tectonic activity of the rock mass. The t-statistic, which analyzes the significance of Spearman's coefficient r_s 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: a radioactive therapeutic element

This paper presents selected issues related to the use of 222Rn in therapeutic treatments. Radon is a radioactive element whose usage in medicine for more than 100 years is based on the radiation hormesis theory. However, owing to the radioactive character of this element and the fact that its alpha-radioactive decay is the source of other radionuclides, its therapeutic application has been raising serious doubts. The author points to potential sources and carriers of radon in the environment that could supply radon for use in a variety of therapies. Except for centuries-long tradition of using radon groundwaters, and later also the air in caves and underground workings, the author would also like to focus on soil air, which is still underestimated as a source of radon. The text presents different methods of obtaining this radioactive gas from groundwaters, the air in caves, mining galleries and soil air, and it presents new possibilities in this field. The author also discusses problems related to the transportation and storage of radon obtained from the environment.

Within radon-prone areas, it is often necessary to de-radon groundwaters that are intended for human consumption and household usage. Also, dry radon wells are used to prevent radon migration from the ground into residential buildings. The author proposes using radon released from radon groundwaters and amassed in dry radon wells for radonotherapy treatments. Thanks to this, it is possible to reduce the cost of radiological protection of people within radon-prone areas while still exploiting the 222Rn obtained for a variety of therapies.

With regard to the ongoing and still unsettled dispute concerning the beneficial or detrimental impact of radon on the human organism, the author puts special emphasis on the necessity of strictly monitoring both the activity concentration of 222Rn in media used for therapeutic treatments and of its radioactive decay products. Monitoring should be also extended to the environments in which such treatments are delivered (inhalatoriums, baths, saunas, showers, pools and other facilities), as well as to the patients – during and after the radonotherapy treatments. It is also essential to monitor the dose of radon and its daughters that is received by persons undergoing radon therapy. This should facilitate the assessment of the effectiveness of these treatments, which may contribute to a fuller understanding of the mechanisms of radon impact, and ionizing radiation in general, on the human organism. This will make it easier to ultimately confirm or reject the radiation hormesis theory. It is also essential to monitor the effective dose that is received by medical and technical staff employed to deliver the radonotherapy treatments.

Predictive analysis of shaft station radon concentrations in underground uranium mine: A case study

This paper presented a method for predicting shaft station radon concentrations in a uranium mine of China through theoretical analysis, mathematical derivation and Monte-Carlo simulation. Based upon the queuing model for tramcars, the average waiting time of tramcars and average number of waiting tramcars were determined, which were further used in developing the predictive model for calculating shaft station radon concentrations. The results exhibit that the extent of variation of shaft station radon concentration in the case study mine is not significantly affected by the queuing process of tramcars, and is always within the allowable limit of 200 Bq m(-3). Thus, the empirical limit of 100,000 T annual ore-hoisting yields has no value in ensuring radiation safety for this mine. Moreover, the developed model has been validated and proved useful in assessing shaft station radon levels for any uranium mine with similar situations.

Investigation of temperature and barometric pressure variation effects on radon concentration in the Sopronbánfalva Geodynamic Observatory, Hungary

Radon concentration variation has been monitored since 2009 in the artificial gallery of the Sopronbánfalva Geodynamic Observatory, Hungary. In the observatory, the radon concentration is extremely high, 100-600 kBq m(-3) in summer and some kBq m(-3) in winter. The relationships between radon concentration, temperature and barometric pressure were separately investigated in the summer and winter months by Fast Fourier Transform, Principal Component Analysis, Multivariable Regression and Partial Least Square analyses in different frequency bands. It was revealed that the long-period radon concentration variation is mainly governed by the temperature (20 kBq m(-1) °C(-1)) both in summer and winter. The regression coefficients between long-period radon concentration and barometric pressure are -1.5 kBq m(-3) hPa(-1) in the summer and 5 kBq m(-3) hPa(-1) in the winter months. In the 0.072-0.48 cpd (cycles per day) frequency band the effect of the temperature is about -1 kBq m(-3) °C(-1) and that of the barometric pressure is -5 kBq m(-3) hPa(-1) in summer and -0.5 kBq m(-3) hPa(-1) in winter. In the high frequency range (>0.48 cpd) all regression coefficients are one order of magnitude smaller than in the range of 0.072-0.48 cpd. Fast Fourier Transform of the radon concentration, temperature and barometric pressure time series revealed S1, K1, P1, S2, K2, M2 tidal constituents in the data and weak O1 components in the radon concentration and barometric pressure series. A detailed tidal analysis, however, showed that the radon tidal components are not directly driven by the gravitational force but rather by solar radiation and barometric tide. Principal Component Analysis of the raw data was performed to investigate the yearly, summer and winter variability of the radon concentration, temperature and barometric pressure. In the summer and winter periods the variability does not change. The higher variability of the radon concentration compared to the variability of the temperature and the barometric pressure shows that besides the temperature and barometric pressure variations other agents, e.g. natural ventilation of the observatory, wind, etc. also play an important role in the radon concentration variation.

Short-term radon activity concentration changes along the Underground Educational Tourist Route in the Old Uranium Mine in Kletno (Sudety Mts., SW Poland)

Short-term (222)Rn activity concentration changes along the Underground Educational Tourist Route in the Old Uranium Mine in Kletno were studied, based on continuous measurements conducted between 16 May 2008 and 15 May 2010. The results were analysed in the context of numbers of visitors arriving at the facility in particular seasons and the time per day spent inside by staff and visitors. This choice was based on partially published earlier findings (Fijałkowska-Lichwa and Przylibski, 2011). Results for the year 2009 were analysed in depth, because it is the only period of observation covering a full calendar year. The year 2009 was also chosen for detailed analysis of short-term radon concentration changes, because in each period of this year (hour, month, season) fluctuations of noted values were the most visible. Attention has been paid to three crucial issues linked to the occurrence and behaviour of radon and to the radiological protection of workers and visitors at the tourist route in Kletno. The object of study is a complex of workings in a former uranium mine situated within a metamorphic rock complex in the most radon-prone area in Poland. The facility has been equipped with a mechanical ventilation system, which is turned on after the closing time and at the end of the working day for the visitor service staff, i.e. after 6 p.m. Short-term radon activity concentration changes along the Underground Educational Tourist Route in the Old Uranium Mine in Kletno are related to the activity of the facility's mechanical ventilation. Its inactivity in the daytime results in the fact that the highest values of (222)Rn activity concentration are observed at the time when the facility is open to visitors, i.e. between 10 a.m. and 6 p.m. The improper usage of the mechanical ventilation system is responsible for the extremely unfavourable working conditions, which persist in the facility for practically all year. The absence of appropriate radiological protection (i.e. preventive measures like shortening working day, dosimetric measurements in the workplace) is a serious problem in the Kletno adit.

Estimating the importance of factors influencing the radon-222 flux from building walls

Radiation hazard in dwellings is dominated by the contribution of radon-222 released from soil and bedrock, but the contribution of building materials can also be important. Using a simple air mixing model in a 2-story house with an attic and a basement, it is estimated that a significant risk arises when the Wall Radon exhalation Flux (WRF) exceeds 10×10(-3) Bq·m(-2)·s(-1). WRF is studied using a multiphase advection-diffusion 3-layer analytical model with advective flow, possibly induced by a pressure deficit inside the house compared with the outside atmosphere. To first order, in most circumstances, the WRF is proportional to the wall thickness and to the radon source term, the effective radium concentration EC(Ra), which is the product of the radium-226 concentration by the emanation coefficient E. The WRF decreases with increasing material porosity and exhibits a maximum for water saturation of about 50%. For EC(Ra)=10 Bq·kg(-1), in many instances, WRF is larger than 10×10(-3) Bq·m(-2)·s(-1) and, therefore, EC(Ra)=10 Bq·kg(-1) can be considered as the typical limit not to be exceeded by building materials. An upper limit of the WRF is obtained in the purely advective regime, independent of porosity or moisture content, which can thus be used as a robust safety guideline. The sensitivity of WRF to temperature, due to the temperature sensitivity of EC(Ra) or the temperature sensitivity of radon Henry constant can be larger than 5% for the seasonal variation in the presence of slight pressure deficit. The temperature sensitivity of EC(Ra) is the dominant effect, except for moist walls. Temperature and moisture variation effects on the WRF potentially can account for most observed seasonal variations of radon concentration in houses, in addition to seasonal changes of air exchange, suggesting that the contribution of walls should be considered when designing remediation strategies and studied with dedicated experiments.

The distribution of radon in tunnels with different geological characteristics in China

In China, as the economy is developing and the population is expanding, some underground buildings have been used as supermarkets, restaurants and entertainment places. Tunnels in mountains are one type of underground building, and the radon ((222)Rn) level in tunnels is an important issue. Radon levels in different type tunnels appear to differ, and relatively higher levels of (222)Rn are associated with particular types of bedrock. The (222)Rn levels in tunnels in five different geological characteristics were analyzed. Those built in granite had the highest (222)Rn levels with a geometric mean (GM) of 280Bqm(-3), while those built in limestone (GM: 100Bqm(-3)) and andesitic porphyry (GM: 96Bqm(-3)) were lower. The sequence of (222)Rn concentrations was: granite>tuff>quartz sandstone>limestone>andesitic porphyry, and the (222)Rn in granite was statistically significantly higher than in limestone and andesitic porphyry. Tunnels built in granite, tuff, quartz sandstone, limestone tended to have higher (222)Rn concentrations in summer than in winter, while the reverse tendency was true in andesitic porphyry tunnels. Only the difference in limestone was statistically significant.

Spatiotemporal variation of radon and carbon dioxide concentrations in an underground quarry: coupled processes of natural ventilation, barometric pumping and internal mixing

Radon-222 and carbon dioxide concentrations have been measured during several years at several points in the atmosphere of an underground limestone quarry located at a depth of 18 m in Vincennes, near Paris, France. Both concentrations showed a seasonal cycle. Radon concentration varied from 1200 to 2000 Bq m(-3) in summer to about 800-1400 Bq m(-3) in winter, indicating winter ventilation rates varying from 0.6 to 2.5 x 10(-6) s(-1). Carbon dioxide concentration varied from 0.9 to 1.0% in summer, to about 0.1-0.3% in winter. Radon concentration can be corrected for natural ventilation using temperature measurements. The obtained model also accounts for the measured seasonal variation of carbon dioxide. After correction, radon concentrations still exhibit significant temporal variation, mostly associated with the variation of atmospheric pressure, with coupling coefficients varying from -7 to -26 Bq m(-3) hPa(-1). This variation can be accounted for using a barometric pumping model, coupled with natural ventilation in winter, and including internal mixing as well. After correction, radon concentrations exhibit residual temporal variation, poorly correlated between different points, with standard deviations varying from 3 to 6%. This study shows that temporal variation of radon concentrations in underground cavities can be understood to a satisfactory level of detail using non-linear and time-dependent modelling. It is important to understand the temporal variation of radon concentrations and the limitations in their modelling to monitor the properties of natural or artificial underground settings, and to be able to assess the existence of new processes, for example associated with the preparatory phases of volcanic eruptions or earthquakes.

Seasonal variations of natural ventilation and radon-222 exhalation in a slightly rising dead-end tunnel

The concentration activity of radon-222 has been monitored, with some interruptions, from 1997 to 2005 in the end section of a slightly rising, dead-end, 38-m long tunnel located in the Phulchoki hill, near Kathmandu, Nepal. While a high concentration varying from 6 x 10(3) Bq m(-3) to 10 x 10(3) Bq m(-3) is observed from May to September (rainy summer season), the concentration remains at a low level of about 200 Bq m(-3) from October to March (dry winter season). This reduction of radon concentration is associated with natural ventilation of the tunnel, which, contrary to expectations for a rising tunnel, takes place mainly from October to March when the outside air temperature drops below the average tunnel temperature. This interpretation is supported by temperature measurements in the atmosphere of the tunnel, a few meters away from the entrance. The temporal variations of the diurnal amplitude of this temperature indeed follow the ventilation rate deduced from the radon measurements. In the absence of significant ventilation (summer season), the radon exhalation flux at the rock surface into the tunnel atmosphere can be inferred; it exhibits a yearly variation with additional transient reductions associated with heavy rainfall, likely to be due to water infiltration. No effect of atmospheric pressure variations on the radon concentration is observed in this tunnel. This experiment illustrates how small differences in the location and geometry of a tunnel can lead to vastly different behaviours of the radon concentration versus time. This observation has consequences for the estimation of the dose rate and the practicability of radon monitoring for tectonic purposes in underground environments.

Anthropogenic CO2-flux into cave atmosphere and its environmental impact: A case study in the Císarská Cave (Moravian Karst, Czech Republic)

The evolution of CO2 levels was studied in the ventilated and unventilated Nagel Dome chamber (the Císarská Cave) with- and without human presence. Based on a simplified dynamic model and CO2/Rn data (222Rn considered as a conservative tracer), two types of CO2-fluxes into the chamber were distinguished: (1) the natural input of (2-4) x 10(-6) m3 s(-1), corresponding to a flux of (8.5-17) x 10(-10) m3 m(-2) s(-1) and (2) an anthropogenic input of (0.6-2.5) x 10(-4) m3 s(-1), corresponding to an average partial flux of (4.8-7.7) x 10(-6) m3 s(-1) person(-1). The chamber ventilation rates were calculated in the range from 0.033 to 0.155 h(-1). Comparison of the chamber CO2-levels with chamber dripwater chemistry indicates that the peak CO2-concentrations during stay of persons (log p(CO2) approximately -2.97, -2.89, and -2.83) do not reach the theoretical values at which dripwater carbonate species and air CO2 are at equilibrium (log p(CO2[DW]) approximately -2.76 to -2.79). This means that CO2-degassing of the dripwaters will continue, increasing supersaturation with respect to calcite (dripwater saturation index defined as SI(calcite) = a(Ca2+)a(CO3(2-))/10(-8.4) varied in the range from 0.76 to 0.86). The p(CO2[DW]) values, however, would easily be exceeded if the period of person stay in the chamber had been slightly extended (from 2.85 to 4 h under given conditions). In such case, the dripwater CO2-degassing would be inverted into CO2-dissolution and dripwater supersaturation would decrease. Achieving the threshold values at which water become aggressive to calcite (log p(CO2[EK]) approximately -1.99, -2.02, and -1.84) would require extreme conditions, e.g., simultaneous presence of 100 persons in the cave chamber for 14 h. The study should contribute to a better preservation of cave environment.