Geogenic Radon Research Articles

Estimation of the radon production rate in granite rocks and evaluation of the implications for geogenic radon potential maps: A case study in Central Portugal

The goal of this study was to estimate radon gas production rate in granitic rocks and identify the factors responsible for the observed variability. For this purpose, 180 samples were collected from pre-Hercynian and Hercynian rocks in north and central Portugal and analysed for a) 226Ra activity, b) radon (222Rn) per unit mass activity, and c) radon gas emanation coefficient. On a subset of representative samples from the same rock types were also measured d) apparent porosity and e) apparent density. For each of these variables, the values ranged as follows: a) 15 to 587Bqkg-1, b) 2 to 73Bqkg-1, c) 0.01 to 0.80, d) 0.3 to 11.4 % and e) 2530 to 2850kgm-3. Radon production rate varied between 40 to 1386Bqm-3h-1. The variability observed was associated with geologically late processes of low and high temperature which led to the alteration of the granitic rock with mobilization of U and increase in radon 222Rn gas emanation. It is suggested that, when developing geogenic radon potential maps, data on uranium concentration in soils/altered rock should be used, rather than data obtained from unaltered rock.

Quantile regression and Bayesian cluster detection to identify radon prone areas

Albeit the dominant source of radon in indoor environments is the geology of the territory, many studies have demonstrated that indoor radon concentrations also depend on dwelling-specific characteristics. Following a stepwise analysis, in this study we propose a combined approach to delineate radon prone areas. We first investigate the impact of various building covariates on indoor radon concentrations. To achieve a more complete picture of this association, we exploit the flexible formulation of a Bayesian spatial quantile regression, which is also equipped with parameters that controls the spatial dependence across data. The quantitative knowledge of the influence of each significant building-specific factor on the measured radon levels is employed to predict the radon concentrations that would have been found if the sampled buildings had possessed standard characteristics. Those normalised radon measures should reflect the geogenic radon potential of the underlying ground, which is a quantity directly related to the geological environment. The second stage of the analysis is aimed at identifying radon prone areas, and to this end, we adopt a Bayesian model for spatial cluster detection using as reference unit the building with standard characteristics. The case study is based on a data set of more than 2000 indoor radon measures, available for the Abruzzo region (Central Italy) and collected by the Agency of Environmental Protection of Abruzzo, during several indoor radon monitoring surveys.

Correlation analysis of the natural radionuclides in soil and indoor radon in Vojvodina, Province of Serbia

The most dominant source of indoor radon is the underlying soil, so the enhanced levels of radon are usually expected in mountain regions and geology units with high radium and uranium content in surface soils. Laboratory for radioactivity and dose measurement, Faculty of Sciences, University of Novi Sad has rich databases of natural radionuclides concentrations in Vojvodina soil and also of indoor radon concentrations for the region of Vojvodina, Northern Province of Serbia. In this paper we present the results of correlative and multivariate analysis of these results and soil characteristics in order to estimate the geogenic radon potential. The correlative and multivariate analysis were done using Toolkit for Multivariate Analysis software package TMVA package, within ROOT analysis framework, which uses several comparable multivariate methods for our analysis. The evaluation ranking results based on the best signal efficiency and purity, show that the Boosted Decision Trees (BDT) and Multi Layer Preceptor (MLP), based on Artificial Neural Network (ANN), are multivariate methods which give the best results in the analysis. The BDTG multivariate method shows that variables with the highest importance are radio-nuclides activity on 30 cm depth. Moreover, the multivariate regression methods give a good approximation of indoor radon activity using full set of input variables. On several locations in the city of Novi Sad the results of indoor radon concentrations, radon emanation from soil, gamma spectrometry measurements of underlying soil and geology characteristics of soil were analyzed in detail in order to verify previously obtained correlations for Vojvodina soil.

Geogenic and anthropogenic impacts on indoor radon in the Techa River region

Indoor radon concentration was studied in the 14 settlements located near the Techa River, which was contaminated by radioactive wastes in 1950-s. Results of the radon survey were used for analysis of the relationship between the indoor radon and main geologic factors (Pre-Jurassic formations, Quaternary sediments and faults), local geogenic radon potential and anthropogenic factors. Main influencing factors explain 58% of the standard deviation of indoor radon concentration. Association of the air exchange influence over radon concentration with underlying geological media was related to different contributions of geogenic advective and diffusive radon entries. The properties of geological formation to transfer radon gas in interaction with the house can be considered within the radon geogenic potential concept. The study of the radon exposure of the Techa River population can be used to estimate the contribution of natural radon to the overall radiation exposure of the local population during the period of radioactive waste discharges.

Modeling of geogenic radon in Switzerland based on ordered logistic regression

PurposeThe estimation of the radon hazard of a future construction site should ideally be based on the geogenic radon potential (GRP), since this estimate is free of anthropogenic influences and building characteristics. The goal of this study was to evaluate terrestrial gamma dose rate (TGD), geology, fault lines and topsoil permeability as predictors for the creation of a GRP map based on logistic regression. MethodSoil gas radon measurements (SRC) are more suited for the estimation of GRP than indoor radon measurements (IRC) since the former do not depend on ventilation and heating habits or building characteristics. However, SRC have only been measured at a few locations in Switzerland. In former studies a good correlation between spatial aggregates of IRC and SRC has been observed. That’s why we used IRC measurements aggregated on a 10 km × 10 km grid to calibrate an ordered logistic regression model for geogenic radon potential (GRP). As predictors we took into account terrestrial gamma doserate, regrouped geological units, fault line density and the permeability of the soil. ResultsThe classification success rate of the model results to 56% in case of the inclusion of all 4 predictor variables. Our results suggest that terrestrial gamma doserate and regrouped geological units are more suited to model GRP than fault line density and soil permeability. ConclusionOrdered logistic regression is a promising tool for the modeling of GRP maps due to its simplicity and fast computation time. Future studies should account for additional variables to improve the modeling of high radon hazard in the Jura Mountains of Switzerland.

Estimation and mapping of uranium content of geological units in France

In France, natural radiation accounts for most of the population exposure to ionizing radiation. The Institute for Radiological Protection and Nuclear Safety (IRSN) carries out studies to evaluate the variability of natural radioactivity over the French territory. In this framework, the present study consisted in the evaluation of uranium concentrations in bedrocks. The objective was to provide estimate of uranium content of each geological unit defined in the geological map of France (1:1,000,000). The methodology was based on the interpretation of existing geochemical data (results of whole rock sample analysis) and the knowledge of petrology and lithology of the geological units, which allowed obtaining a first estimate of the uranium content of rocks. Then, this first estimate was improved thanks to some additional information. For example, some particular or regional sedimentary rocks which could present uranium contents higher than those generally observed for these lithologies, were identified. Moreover, databases on mining provided information on the location of uranium and coal/lignite mines and thus indicated the location of particular uranium-rich rocks. The geological units, defined from their boundaries extracted from the geological map of France (1:1,000,000), were finally classified into 5 categories based on their mean uranium content. The map obtained provided useful data for establishing the geogenic radon map of France, but also for mapping countrywide exposure to terrestrial radiation and for the evaluation of background levels of natural radioactivity used for impact assessment of anthropogenic activities.

Geographically weighted regression and geostatistical techniques to construct the geogenic radon potential map of the Lazio region: A methodological proposal for the European Atlas of Natural Radiation

In many countries, assessment programmes are carried out to identify areas where people may be exposed to high radon levels. These programmes often involve detailed mapping, followed by spatial interpolation and extrapolation of the results based on the correlation of indoor radon values with other parameters (e.g., lithology, permeability and airborne total gamma radiation) to optimise the radon hazard maps at the municipal and/or regional scale. In the present work, Geographical Weighted Regression and geostatistics are used to estimate the Geogenic Radon Potential (GRP) of the Lazio Region, assuming that the radon risk only depends on the geological and environmental characteristics of the study area. A wide geodatabase has been organised including about 8000 samples of soil-gas radon, as well as other proxy variables, such as radium and uranium content of homogeneous geological units, rock permeability, and faults and topography often associated with radon production/migration in the shallow environment. All these data have been processed in a Geographic Information System (GIS) using geospatial analysis and geostatistics to produce base thematic maps in a 1000 m × 1000 m grid format. Global Ordinary Least Squared (OLS) regression and local Geographical Weighted Regression (GWR) have been applied and compared assuming that the relationships between radon activities and the environmental variables are not spatially stationary, but vary locally according to the GRP. The spatial regression model has been elaborated considering soil-gas radon concentrations as the response variable and developing proxy variables as predictors through the use of a training dataset. Then a validation procedure was used to predict soil-gas radon values using a test dataset. Finally, the predicted values were interpolated using the kriging algorithm to obtain the GRP map of the Lazio region. The map shows some high GRP areas corresponding to the volcanic terrains (central-northern sector of Lazio region) and to faulted and fractured carbonate rocks (central-southern and eastern sectors of the Lazio region). This typical local variability of autocorrelated phenomena can only be taken into account by using local methods for spatial data analysis. The constructed GRP map can be a useful tool to implement radon policies at both the national and local levels, providing critical data for land use and planning purposes.

Indoor radon, geogenic radon surrogates and geology – Investigations on their correlation

The indoor radon concentration was measured in most houses in a couple of municipalities in Austria. At the same time the activity concentration of radium in soil, the soil gas radon concentration, the permeability of the ground and the ambient dose equivalent rate were also measured and the geological situations (geological units) were recorded too. From the indoor radon concentration and different house and living parameters a radon potential (Austrian radon potential) was derived which should represent the radon concentration in a standard room. Another radon potential (Neznal radon potential) was calculated from the soil gas radon concentration and the permeability. The aim of the investigation was to correlate all the different variables and to test if the use of surrogate data (e.g. geological information, ambient dose equivalent rate, etc.) can be used to judge the radon risk for an area without performing numerous indoor measurements.

Estimating the terrestrial gamma dose rate by decomposition of the ambient dose equivalent rate

An extensive network of dose rate monitoring stations continuously measures ambient dose rate across Europe, as part of the EURDEP system. Its purpose is early warning in radiological emergencies and documenting its temporal and spatial evolution. In normal conditions, when there is no contribution to the dose rate signal coming from fresh anthropogenic contamination, the data represent the radiation “background”, i.e. the combined natural radiation and existing anthropogenic contamination (by global and Chernobyl fallout). These data are being stored, but have so far not been evaluated in depth, or used for any purpose. In the framework of the EU project ‘European Atlas of Natural Radiation’ the idea has emerged to exploit these data for generating a map of natural terrestrial gamma radiation. This component contributes to the total radiation exposure and knowing its geographical distribution can help establishing local ‘radiation budgets’. A further use could be found in terrestrial dose rate as a proxy of the geogenic radon potential, as both quantities are related by partly the same source, namely uranium content of the ground. In this paper, we describe in detail the composition of the ambient dose equivalent rate as measured by the EURDEP monitors with respect to its physical nature and to its sources in the environment. We propose and compare methods to recover the terrestrial component from the gross signal. This requires detailed knowledge of detector response. We consider the probes used in the Austrian, Belgian and German dose rate networks, which are the respective national networks supplying data to EURDEP. It will be shown that although considerable progress has been made in understanding the dose rate signals, there is still space for improvement in terms of modelling and model parameters. An indispensable condition for success of the endeavour to establish a Europe-wide map of terrestrial dose rate background is progress in harmonising the European dose rate monitoring network.

Mapping geogenic radon potential by regression kriging

Radon (222Rn) gas is produced in the radioactive decay chain of uranium (238U) which is an element that is naturally present in soils. Radon is transported mainly by diffusion and convection mechanisms through the soil depending mainly on the physical and meteorological parameters of the soil and can enter and accumulate in buildings. Health risks originating from indoor radon concentration can be attributed to natural factors and is characterized by geogenic radon potential (GRP). Identification of areas with high health risks require spatial modeling, that is, mapping of radon risk. In addition to geology and meteorology, physical soil properties play a significant role in the determination of GRP. In order to compile a reliable GRP map for a model area in Central-Hungary, spatial auxiliary information representing GRP forming environmental factors were taken into account to support the spatial inference of the locally measured GRP values. Since the number of measured sites was limited, efficient spatial prediction methodologies were searched for to construct a reliable map for a larger area. Regression kriging (RK) was applied for the interpolation using spatially exhaustive auxiliary data on soil, geology, topography, land use and climate. RK divides the spatial inference into two parts. Firstly, the deterministic component of the target variable is determined by a regression model. The residuals of the multiple linear regression analysis represent the spatially varying but dependent stochastic component, which are interpolated by kriging. The final map is the sum of the two component predictions. Overall accuracy of the map was tested by Leave-One-Out Cross-Validation. Furthermore the spatial reliability of the resultant map is also estimated by the calculation of the 90% prediction interval of the local prediction values. The applicability of the applied method as well as that of the map is discussed briefly.

From the European indoor radon map towards an atlas of natural radiation.

In 2006, the Joint Research Centre of the European Commission launched a project to map radon at the European level, as part of a planned European Atlas of Natural Radiation. It started with a map of indoor radon concentrations. As of May 2014, this map includes data from 24 countries, covering a fair part of Europe. Next, a European map of geogenic radon, intended to show ‘what earth delivers’ in terms of radon potential (RP), was started in 2008. A first trial map has been created, and a database was established to collect all available data relevant to the RP. The Atlas should eventually display the geographical distribution of physical quantities related to natural radiation. In addition to radon, it will comprise maps of quantities such as cosmic rays and terrestrial gamma radiation. In this paper, the authors present the current state of the radon maps and the Atlas.

Towards a Brazilian radon map: consortium radon Brazil.

Recently, the idea of generating radon map of Brazil has emerged. First attempts of coordinating radon surveys--carried out by different groups across the country--and initial discussions on how to proceed on a larger scale were made at the First Brazilian Radon Seminary, Natal, September 2012. Conventionally, it is believed that indoor radon is no major problem in Brazil, because the overall benign climate usually allows high ventilation rates. Nevertheless, scattered measurements have shown that moderately high indoor radon concentrations (up to a few hundred Bq m⁻³) do occur regionally. Brazilian geology is very diverse and there are regions where an elevated geogenic radon potential exists or is expected to exist. Therefore, a Brazilian Radon Survey is expected to be a challenge, although it appears an important issue, given the rising concern of the public about the quality of its environment.

Determination of radon prone areas by optimized binary classification

Geogenic radon prone areas are regions in which for natural reasons elevated indoor radon concentrations must be expected. Their identification is part of radon mitigation policies in many countries, as radon is acknowledged a major indoor air pollutant, being the second cause of lung cancer after smoking. Defining and estimating radon prone areas is therefore of high practical interest.In this paper a method is presented which uses the geogenic radon potential as predictor and thresholds of indoor radon concentration for defining radon prone areas, from which thresholds for the geogenic radon potential are deduced which decide whether a location is flagged radon prone or not, in the absence of actual indoor observations.The overall results are different maps of radon prone areas, derived from the geogenic radon map, and depending (1) on the criterion which defines what a radon prone area is; and (2) on the choice of score whose maximization defines the optimal classifier. Such map is not the result of a transfer model (geogenic to indoor radon), but of the optimization of a classification rule. The method is computationally simple but has its caveats and statistical traps, some of which are also addressed.

Mapping the geogenic radon potential: methodology and spatial analysis for central Hungary

A detailed geogenic radon potential (GRP) mapping based on field soil gas radon and soil gas permeability measurements was carried out in this study. A conventional continuous variable approach was used in this study for GRP determination and to test its applicability to the selected area of Hungary. Spatial pattern of soil gas radon concentration, soil permeability and GRP and the relationship between geological formations and these parameters were studied by performing detailed spatial analysis.Exploratory data analysis revealed that higher soil gas radon activity concentration and GRP characterizes the mountains and hills than the plains. The highest values were found in the proluvial–deluvial sediments, rock debris on the downhill slopes eroded from hills. Among the Quaternary sediments, which characterize the study area, the fluvial sediment has the highest values, which are also located in the hilly areas. The lowest values were found in the plain areas covered by drift sand, fluvioeolic sand, fluvial sand and loess. As a conclusion, radon is related to the sediment cycle in the study area.A geogenic radon risk map was created, which assists human health risk assessment and risk reduction since it indicates the potential of the source of indoor radon. The map shows that low and medium geogenic radon potential characterizes the study area in central Hungary. High risk occurs only locally. The results reveal that Quaternary sediments are inhomogeneous from a radon point of view, fluvial sediment has medium GRP, whereas the other rock formations such as drift sand, fluioeolic sand, fluvial sand and loess, found in the study area, have low GRP.

Geographical distribution of the annual mean radon concentrations in primary schools of Southern Serbia – application of geostatistical methods

Between 2008 and 2011 a survey of radon (222Rn) was performed in schools of several districts of Southern Serbia. Some results have been published previously (Žunić et al., 2010; Carpentieri et al., 2011; Žunić et al., 2013). This article concentrates on the geographical distribution of the measured Rn concentrations. Applying geostatistical methods we generate “school radon maps” of expected concentrations and of estimated probabilities that a concentration threshold is exceeded. The resulting maps show a clearly structured spatial pattern which appears related to the geological background. In particular in areas with vulcanite and granitoid rocks, elevated radon (Rn) concentrations can be expected. The “school radon map” can therefore be considered as proxy to a map of the geogenic radon potential, and allows identification of radon-prone areas, i.e. areas in which higher Rn radon concentrations can be expected for natural reasons.It must be stressed that the “radon hazard”, or potential risk, estimated this way, has to be distinguished from the actual radon risk, which is a function of exposure. This in turn may require (depending on the target variable which is supposed to measure risk) considering demographic and sociological reality, i.e. population density, distribution of building styles and living habits.

A statistical evaluation of the influence of housing characteristics and geogenic radon potential on indoor radon concentrations in France

Radon-222 is a radioactive natural gas produced by the decay of radium-226, known to be the main contributor to natural background radiation exposure. Effective risk management needs to determine the areas in which the density of buildings with high radon levels is likely to be highest. Predicting radon exposure from the location and characteristics of a dwelling could also contribute to epidemiological studies. Beginning in the nineteen-eighties, a national radon survey consisting in more than 10,000 measurements of indoor radon concentrations was conducted in French dwellings by the Institute for Radiological Protection and Nuclear Safety (IRSN). Housing characteristics, which may influence radon accumulation in dwellings, were also collected. More recently, the IRSN generated a French geogenic radon potential map based on the interpretation of geological features. The present study analyzed the two datasets to investigate the factors influencing indoor radon concentrations using statistical modeling and to determine the optimum use of the information on geogenic radon potential that showed the best statistical association with indoor radon concentration. The results showed that the variables associated with indoor radon concentrations were geogenic radon potential, building material, year of construction, foundation type, building type and floor level. The model, which included the surrounding geogenic radon potential (i.e. the average geogenic radon potential within a disc of radius 20 km centered on the indoor radon measurement point) and variables describing house-specific factors and lifestyle explained about 20% of the overall variability of the logarithm of radon concentration. The surrounding geogenic radon potential was fairly closely associated with the local average indoor radon concentration. The prevalence of exposure to radon above specific thresholds and the average exposures to radon clearly increased with increasing classes of geogenic radon potential. Combining the two datasets enabled improved assessment of radon exposure in a given area in France.

Dynamics of soil gas radon concentration in a highly permeable soil based on a long-term high temporal resolution observation series

This paper studies the temporal variation of soil gas radon activity concentration in a highly permeable (k = 2.0E-11 m2) sandy-gravelly soil in order to understand if temporal variation of soil gas radon activity concentration can affect geogenic radon potential determination. Geogenic radon potential provides information about the potential risk from radon. Its calculation takes into account the equilibrium, saturated at infinite depth, soil gas radon activity concentration (c∞). This concentration may vary at annual time scale due to the environmental conditions. A long-term (yearly) and high temporal resolution (15 min) observation, applied in this study, reveal various temporal features such as long-term trend, seasonality, daily periodicity and sudden events in soil gas radon time series.Results show seasonal and daily periodical variation of the measured soil gas radon activity concentration (csoilRn) in a highly permeable sandy-gravelly soil with definite seasons without obvious long transitional periods. The winter (from October 2010 to April 2011) is characterized by 2.5 times higher average soil gas radon activity concentration (median is 7.0 kBq m−3) than the summer (August, September 2010 and May, June, July 2011) (median is 2.8 kBq m−3). Daily periodicity, which is much less than the seasonal one, controls the soil gas radon activity concentration mainly in the summer season. Average (AM) value of csoilRn is higher at night than in the daytime with about 18% and 3.8% in summer and in winter, respectively. As a conclusion, in case of single csoilRn measurement on a highly permeable (k ≥ 2.0E-11 m2) soil, similar to our test site, csoilRn should be corrected according to the seasons for calculating the equilibrium activity concentration c∞ value.

The European map of the geogenic radon potential

As part of its projected European Atlas of Natural Radiation (EANR), the Joint Research Centre (JRC) of the European Commission, in cooperation with research institutions and radioprotection authorities all over Europe, is currently developing a map of the geogenic radon potential. In an accompanying report the state of knowledge, mapping approaches and problems are discussed. We explain the rationale and the legal situation in Europe and present an overview on the main problems stemming from the heterogeneity of input datasets between participating countries and from the definition of input variables and their differently implemented sampling procedures or protocols. Further topics are definition of the target variable which quantifies the geogenic radon potential and its estimation from heterogeneous input and proxy variables, as well as problems specific to mapping, such as choice of mapping support and resolution. The geogenic map was preceded by a European map of indoor radon concentrations, which is still growing as ever more countries decide to participate, and which served as training for harmonisation problems occurring in the European data realm. We shall also briefly discuss its main results and implications for the geogenic map.

Indoor radon probability calculated from the Czech soil gas radon data in a grid net for the European Geogenic Radon Map construction: test of feasibility

The construction of the European Geogenic Radon Map in a proposed grid system 10 × 10 km requires the data test to derive the probability of exceeding the indoor action level 200 Bq m−3 from the geologically based data. The Czech Republic disposes both indoor and soil gas data sets to test the real probability to exceed 200 Bq m−3 from indoor radon measurements and to compare it with the probability calculated from soil gas radon concentrations. Comparison of real and calculated probability enables to delineate the areas, where under- or overestimation can be expected. The results of data processing show minor differences between processing the raw data in generalised polygons of geological units and in a grid net, when using the generalised geological characteristics of grid cells.

Towards a european atlas of natural radiation: Goal, status and future perspectives

One of the tasks of the European Commission (EC) under the Euratom Treaty is to collect, validate and provide information about the levels of radioactivity in the environment. In order to offer to the public a more balanced view on the annual dose that it may receive from environmental radioactivity, a few years ago we decided to explore the feasibility of preparing a European Atlas of Natural Radiation. By now 21 European countries have already provided indoor radon data, and efforts continue to extend this information to more countries. In addition, we started to investigate the feasibility of a “geogenic radon map”, which would show “what earth delivers” in terms of potential radon hazard. In this paper we present the current state of the art of preparing the Atlas and provide detailed statistics for the results already obtained. The current efforts are still focussed on radon, but also progress on other components, such as cosmic rays and terrestrial gamma radiation, will be presented and discussed.