Radon in indoor air of primary schools: a systematic survey to evaluate factors affecting radon concentration levels and their variability.

UNLABELLED: In order to optimize the design of a national survey aimed to evaluate radon exposure of children in schools in Serbia, a pilot study was carried out in all the 334 primary schools of 13 municipalities of Southern Serbia. Based on data from passive measurements, rooms with annual radon concentration >300 Bq/m(3) were found in 5% of schools. The mean annual radon concentration weighted with the number of pupils is 73 Bq/m(3), 39% lower than the unweighted 119 Bq/m(3) average concentration. The actual average concentration when children are in classrooms could be substantially lower. Variability between schools (CV = 65%), between floors (CV = 24%) and between rooms at the same floor (CV = 21%) was analyzed. The impact of school location, floor, and room usage on radon concentration was also assessed (with similar results) by univariate and multivariate analyses. On average, radon concentration in schools within towns is a factor of 0.60 lower than in villages and at higher floors is a factor of 0.68 lower than ground floor. Results can be useful for other countries with similar soil and building characteristics. PRACTICAL IMPLICATIONS: On average, radon concentrations are substantially higher in schools in villages than in schools located in towns (double,on average). Annual radon concentrations exceeding 300 Bq/m3 were found in 5% of primary schools (generally on ground floors of schools in villages). The considerable variability of radon concentration observed between and within floors indicates a need to monitor concentrations in several rooms for each floor. A single radon detector for each room can be used provided that the measurement error is considerable lower than variability of radon concentration between rooms.

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.

An episode of Ru-106 in air over Europe, September-October 2017 - Geographical distribution of inhalation dose over Europe.

Between the end of September and early October 2017, 106Ru was recorded by air monitoring stations across parts of Europe. In the environment, this purely anthropogenic radionuclide can be detected very rarely only. As far as known, 106Ru is only used in radiotherapy and possibly in radiothermal generators. Therefore, the episode drew considerable interest in the monitoring community, although the activity concentrations and resulting exposure were far below radiological concern. Health consequences can be practically excluded except possibly near the source. 106Ru in aerosols could be detected for several weeks and in some regions of Central and Eastern Europe tens, up to over 100 mBq/m³ were measured as one-day means. Discussions about a possible source continue until today (early 2019). Atmospheric back-modelling led to trajectories likely originating in the Southern to Northern Ural region of Russia and possibly Northern Kazakhstan. Suspiciously, no other anthropogenic radionuclides have been observed alongside, except minute concentrations of comparatively short-lived 103Ru (half life 39 d vs. 376 d for 106Ru). Due to the absence of other anthropogenic radionuclides, a reactor accident can be excluded, although both Ru isotopes are fission products generated in nuclear reactors. The exposure resulting from 106Ru activity concentration in air exceeded 200 mBq × d/m³ in some parts of Central and Eastern Europe. This leads to inhalation doses of up to about 0.3 µSv regionally, assuming the radiologically most efficient speciation, lacking better information, and inhalation dose conversion factors from ICRP 119. We show an interpolated map of the dose distribution over parts of Europe where sufficient measurements are available to us. Overlaying population density, we give an estimate of collective dose. The opportunity is also used to give a short review of origin, properties and use of 106Ru, as well as of accidents which involved release of this radionuclide.

Machine learning methods as a tool to analyse incomplete or irregularly sampled radon time series data.

Machine learning is a class of statistical techniques which has proven to be a powerful tool for modelling the behaviour of complex systems, in which response quantities depend on assumed controls or predictors in a complicated way. In this paper, as our first purpose, we propose the application of machine learning to reconstruct incomplete or irregularly sampled data of time series indoor radon (222Rn). The physical assumption underlying the modelling is that Rn concentration in the air is controlled by environmental variables such as air temperature and pressure. The algorithms "learn" from complete sections of multivariate series, derive a dependence model and apply it to sections where the controls are available, but not the response (Rn), and in this way complete the Rn series. Three machine learning techniques are applied in this study, namely random forest, its extension called the gradient boosting machine and deep learning. For a comparison, we apply the classical multiple regression in a generalized linear model version. Performance of the models is evaluated through different metrics. The performance of the gradient boosting machine is found to be superior to that of the other techniques. By applying learning machines, we show, as our second purpose, that missing data or periods of Rn series data can be reconstructed and resampled on a regular grid reasonably, if data of appropriate physical controls are available. The techniques also identify to which degree the assumed controls contribute to imputing missing Rn values. Our third purpose, though no less important from the viewpoint of physics, is identifying to which degree physical, in this case environmental variables, are relevant as Rn predictors, or in other words, which predictors explain most of the temporal variability of Rn. We show that variables which contribute most to the Rn series reconstruction, are temperature, relative humidity and day of the year. The first two are physical predictors, while "day of the year" is a statistical proxy or surrogate for missing or unknown predictors.

Cryoconites from Alpine glaciers: Radionuclide accumulation and age estimation with Pu and Cs isotopes and 210Pb.

Cryoconites ("cold dust", derived from the Greek) are aeolian sediments accumulated on glacier surfaces. In cryoconites from the surface of the Stubacher Sonnblickkees, a temperate Austrian glacier, extremely high activity concentrations of artificial and natural radionuclides were found. Artificial radionuclides stem from two clearly distinguishable sources, global fallout from the nuclear weapons testing era deposited over a period of years until roughly 1966 and the fallout from Chernobyl in 1986, which was essentially deposited as a single input during one week. Anthropogenic radionuclides identified were 137Cs, 134Cs, 238Pu, 239+240Pu, 90Sr, 241Am, 60Co, 125Sb, 154Eu, and 207Bi. The naturally occurring radionuclides detected were the long-lived radon decay product 210Pb, the primordial radionuclide 4 K and the cosmogenic 7Be. Isotopic ratios of 134Cs/137Cs and 239+240Pu/238Pu were used to separate the nuclide inventory into the contributions of the two aforementioned sources, which show varying degrees of mixing and provide information on the mixing age of the cryoconites. Since isotopic ratios of Pu often have high uncertainties due to low absolute concentrations, age estimation based on this method can be quite inaccurate. Additional information about the age of cryoconites was obtained through analysis of 210Pb, which is constantly deposited over time.

A geostatistical autopsy of the Austrian indoor radon survey (1992-2002).

Indoor radon (Rn) concentrations have been measured intensively in many countries to assess the burden of diseases associated with exposure to this radioactive gas. So-called radon-risk maps have consequently been produced to delineate areas with high levels. Geostatistical techniques are commonly used nowadays to map a range of environmental variables, in particular to generate probability maps of exceeding a given threshold. However, very few case studies in which indoor radon measurements have been investigated using geostatistical techniques have been published so far. By analyzing around 12,000 Rn measurements made in Austrian ground floors during a 10-year survey, we aim here to review and discuss the potential of geostatistics for mapping an environmental variable that shows very strong local variability. In particular, we show how kriging of the scale components can shed new light on various factors that affect the very high spatial variability of the variable.

137Cs and 90Sr transfer to milk in Austrian alpine agriculture.

The alpine regions of Austria were among the most contaminated territories outside of the former USSR after the Chernobyl accident. In the investigated province of Salzburg the median (137)Cs surface deposition was 31.4 kBq m(-2) with maximum values exceeding 90 kBq m(-2) (May 1986). To quantify the transfer of (137)Cs and (90)Sr from vegetation to milk in these seminatural conditions nine seasonally grazed alps were identified and vegetation and milk sampled during summer 2002 and summer 2003. Mean+/-SD milk transfer coefficients (fm) for (137)Cs and (90)Sr were 0.0071+/-0.0009 d l(-1), and 0.0011+/-0.0004 d l(-1), respectively; which for (137)Cs is markedly higher than those fm values found in intensive agricultural systems. Transfer kinetics for (137)Cs into cow milk were approximated using a 2-compartment model with a short and a long-term component. Fitting the model to empirical data results in reliable estimates of the time constant of the short-term component, biological half-life 1.06+/-0.28 d, whereas the estimates of the long-term component are subject to high uncertainties.

Activity ratios of 137Cs, 90Sr and 239+240Pu in environmental samples.

Both global and Chernobyl fallout have resulted in environmental contamination with radionuclides such as 137Cs, 90Sr and 239+240Pu. In environmental samples, 137Cs and 239+240Pu can be divided into the contributions of either source, if also the isotopes 134Cs and 238Pu are measurable, based on the known isotopic ratios in global and Chernobyl fallout. No analogous method is available for 90Sr. The activity ratios of Sr to Cs and Pu, respectively, are known for the actual fallout mainly from air filter measurements; but due to the high mobility of Sr in the environment, compared to Cs and Pu, these ratios generally do not hold for the inventory many years after deposition. In this paper we suggest a method to identify the mean contributions of global and Chernobyl fallout to total Sr in soil, sediment and cryoconite samples from Alpine and pre-Alpine regions of Austria, based on a statistical evaluation of Sr/Cs/Pu radionuclide activity ratios. Results are given for Sr:Cs, Sr:Pu and Cs:Pu ratios. Comparison with fallout data shows a strong depletion of Sr against Cs and Pu.

INDOOR RADON, THORON AND THEIR PROGENY CONCENTRATIONS IN HIGH THORON RURAL SERBIA ENVIRONMENTS.

This article deals with the variation of radon (Rn), thoron (Tn) and their progeny concentrations expressed in terms of equilibrium equivalent concentrations (EERC and EETC), in 40 houses, in four villages of Sokobanja municipality, Southern Serbia. Two types of passive detectors were used: (1) discriminative radon-thoron detector for simultaneous Rn and Tn gases measurements and (2) direct Tn and Rn progeny sensors (DRPS/DTPS) for measuring Rn and Tn progeny concentrations. Detectors were exposed simultaneously for a single period of 12 months. Variations of Tn and EETC appear higher than those of Rn and EERC. Analysis of the spatial variation of the measured concentrations is also reported. This work is part of a wider survey of Rn, Tn and their progeny concentrations in indoor environments throughout the Balkan region started in 2011 year.

The relation between radon in schools and in dwellings: A case study in a rural region of Southern Serbia.

Recognized as a significant health hazard, radon (Rn) has been given increasing attention for years. Surveys of different kinds have been performed in many countries to assess the intensity and the geographical extent of possible Rn problems. Common surveys cover mainly dwellings, the indoor place with highest occupancy, and schools, where people spend a large fraction of their lifetime and which can also be considered exemplary for Rn exposure at workplaces; it has however been observed that relating them is difficult. It was unclear whether residential Rn at a location, or in a region, can be predicted by Rn at a school of that location, or vice versa. To current knowledge, no general rule seems applicable, as few models to describe the relationship between Rn in dwellings and in schools have been developed. In Southern Serbia, a Rn survey in a predominantly rural region was based on measurements in primary schools. The question arose whether or to which degree the results can be considered as indicative or even representative for residential Rn concentrations. To answer the question an additional survey of indoor Rn concentrations in dwellings was initiated, designed and performed in Sokobanja district in 2010-2012 in a manner to be able to detect a relationship if it exists. In the study region, 108 dwellings in 12 villages and towns were selected, with one primary school each. In this paper, we investigate how a relation between Rn in schools and dwellings could be identified and quantified, by developing a model and using experimental data from both the above main and additional surveys. The key criterion is the hypothesis that the relation dwellings - schools, if it exists, is stronger for dwellings closer to a school than for those dwellings further away. We propose methods to test the hypothesis. As result, the hypothesis is corroborated at 95% significance level. More specifically, on town level (typical size about 1 km), the Rn concentration ratio dwelling/school is about 0.8 (geometrical mean), with geometrical standard deviation (GSD) about 1.9. For dwelling and school hypothetically in the same location, the ratio is estimated about 0.7 with GSD about 1.5. We think that the methodology can be applied to structurally similar problems. The results could be used to create "conditional maps" of Rn concentration in dwellings, i.e., for example a map of probabilities that indoor Rn concentrations in dwellings exceed 100 Bq/m3, as function of Rn concentration in the local school.

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.

A note on 207Bi in environmental samples.

Traces of the radionuclide (207)Bi were identified in soil and cryoconite (glacier sediment) samples from Alpine regions of Austria. This nuclide has been produced in thermonuclear explosions mainly in the early 1960s and subsequently dispersed in the atmosphere. Activity concentrations up to 22 Bq/kg d.m. have been found. The ratio (207)Bi:(137)Cs(global fallout) equals (1.70+/-0.12)10(-3), which is in accordance with literature data. When low levels of (207)Bi are assessed by gamma spectrometry, corrections must be made for a gamma line produced in the lead shield by neutron activation due to cosmic neutrons.

Variance of indoor radon concentration: Major influencing factors.

Variance of radon concentration in dwelling atmosphere is analysed with regard to geogenic and anthropogenic influencing factors. Analysis includes review of 81 national and regional indoor radon surveys with varying sampling pattern, sample size and duration of measurements and detailed consideration of two regional surveys (Sverdlovsk oblast, Russia and Niska Banja, Serbia). The analysis of the geometric standard deviation revealed that main factors influencing the dispersion of indoor radon concentration over the territory are as follows: area of territory, sample size, characteristics of measurements technique, the radon geogenic potential, building construction characteristics and living habits. As shown for Sverdlovsk oblast and Niska Banja town the dispersion as quantified by GSD is reduced by restricting to certain levels of control factors. Application of the developed approach to characterization of the world population radon exposure is discussed.

A very long-term HPGe-background gamma spectrum.

Background gamma spectra, taken with one HPGe detector between 1987 and 2003, are summed up to represent a total collect time of 1.03x10(8) s. Standard background gamma lines belonging to the natural decay series as well as more uncommon ones are identified. A high temporal variability of the background lines related to radon decay products is found in comparison with lines related to other radionuclides contained in the environment and in the structural material of the detector.

Estimation of radon prone areas through binary classification, part 2: radon prone geologies.

A radon prone geology is one for which the probability is increased that in a house built on it, elevated indoor Rn concentration will be encountered, or that its Rn potential will be increased. Labelling geological units as Rn prone or not can be an important support in deciding whether a geographical or administrative region in which that geological unit occurs, should be called Rn prone area, possibly in absence of other predictors. In this article a method is proposed which, given a set of geological classes, sorts the classes into Rn prone and non-Rn prone classes depending on a classification criterion which one can choose according the purpose. The method is computationally simple and is demonstrated on the example of Germany.

Status of the European Atlas of Natural Radiation.

According to the EURATOM (European Atomic Energy Community) Treaty, one of the missions of the Joint Research Centre (JRC) of the European Commission (EC) is to collect, process, evaluate and present data on environmental radioactivity. In 2006, the JRC started the 'European Atlas of Natural Radiation' project, in order to give an overview of the geographic distribution of sources of, and exposures to, natural radiation. As a first task, a map of indoor radon concentration was created, because in most cases this is the most important contribution to exposure, and since it could be expected that data collection would take quite some time, because radon (Rn) surveys are very differently advanced between European countries. The authors show the latest status of this map. A technically more ambitious map proved the one of the geogenic Rn potential (RP), due to heterogeneity of data sources across Europe and the need to develop models to estimate a harmonised quantity which adequately measures or classifies the RP. Further maps currently in the making include those of secondary cosmic radiation, of terrestrial gamma radiation and of the concentrations of the elements U, Th and K that are its source. In this article, the authors show the progress of some of these maps.

Distribution of uranium, thorium and some stable trace and toxic elements in human hair and nails in Niska Banja Town, a high natural background radiation area of Serbia (Balkan Region, South-East Europe).

Human hair and nails can be considered as bio-indicators of the public exposure to certain natural radionuclides and other toxic metals over a long period of months or even years. The level of elements in hair and nails usually reflect their levels in other tissues of body. Niska Banja, a spa town located in southern Serbia, with locally high natural background radiation was selected for the study. To assess public exposure to the trace elements, hair and nail samples were collected and analyzed. The concentrations of uranium, thorium and some trace and toxic elements (Mn, Ni, Cu, Sr, Cd, and Cs) were determined using inductively coupled plasma mass spectrometry (ICP-MS). U and Th concentrations in hair varied from 0.0002 to 0.0771 µg/g and from 0.0002 to 0.0276 µg/g, respectively. The concentrations in nails varied from 0.0025 to 0.0447 µg/g and from 0.0023 to 0.0564 µg/g for U and Th, respectively. We found significant correlations between some elements in hair and nails. Also indications of spatial clustering of high values could be found. However, this phenomenon as well as the large variations in concentrations of heavy metals in hair and nail could not be explained. As hypotheses, we propose possible exposure pathways which may explain the findings, but the current data does not allow testing them.

Long-term measurements of radon, thoron and their airborne progeny in 25 schools in Republic of Srpska.

This article reports results of the first investigations on indoor radon, thoron and their decay products concentration in 25 primary schools of Banja Luka, capital city of Republic Srpska. The measurements have been carried out in the period from May 2011 to April 2012 using 3 types of commercially available nuclear track detectors, named: long-term radon monitor (GAMMA 1)- for radon concentration measurements (C(Rn)); radon-thoron discriminative monitor (RADUET) for thoron concentration measurements (C(Tn)); while equilibrium equivalent radon concentration (EERC) and equilibrium equivalent thoron concentrations (EETC) measured by Direct Radon Progeny Sensors/Direct Thoron Progeny Sensors (DRPS/DTPS) were exposed in the period November 2011 to April 2012. In each school the detectors were deployed at 10 cm distance from the wall. The obtained geometric mean concentrations were C(Rn) = 99 Bq m(-3) and C(Tn) = 51 Bq m(-3) for radon and thoron gases respectively. Those for equilibrium equivalent radon concentration (EERC) and equilibrium equivalent thoron concentrations (EETC) were 11.2 Bq m(-3) and 0.4 Bq m(-3), respectively. The correlation analyses showed weak relation only between C(Rn) and C(Tn) as well as between C(Tn) and EETC. The influence of the school geographical locations and factors linked to buildings characteristic in relation to measured concentrations were tested. The geographical location and floor level significantly influence C(Rn) while C(Tn) depend only from building materials (ANOVA, p ≤ 0.05). The obtained geometric mean values of the equilibrium factors were 0.123 for radon and 0.008 for thoron.

Distribution of radiocaesium in an Austrian forest stand.

Within an Austrian spruce stand, vertical distribution of radiocaesium in soil as well as 137Cs concentration in different forest ecosystem compartments including spruce and surface water were investigated 10 years after the Chernobyl accident. The total 137Cs inventory in the forest was estimated to be 46 kBq m-2 (ref. date: 86-05-01). From the collected data annual input rates via litterfall of 0.48% per year and output rates through waterflows of only 0.02% per year were derived. The results identify the high importance of forest soils as a sink for radiocaesium. The estimated ecological residence half-times turned out to be highest in the organic soil horizons (1-3 years per cm), whereas in mineral horizons the values decrease significantly. As a consequence, soil inventory represents more than 95% of the total, whereas only approximately 3.3% of the 137Cs inventory is stored in the living biomass of spruce trees and a further 0.5% in the phytomass of understorey vegetation.

Multifractal analysis of the 137Cs fallout pattern in Austria resulting from the Chernobyl accident.

The cumulative deposition of the 137Cs fallout in Austria resulting from the passage of the Chernobyl cloud has been investigated by applying correlation dimension and hyperbolic frequency distribution methods. For the analysis, a total of 1,881 deposition values were used, which were collected by the Federal Environmental Agency of Austria and the Federal Ministry of Health, representing all available measurements of 137Cs in soil made in Austria after the Chernobyl accident. From these data a hyperbolic exponent for the frequency distribution of 4.0 and a set of fractal correlation dimensions, which decrease from 1.426 +/- 0.022 (for the whole network) to 0.706 +/- 0.047 (for 137Cs values > or = 100 kBq m(-2)), were derived, thus confirming that the fallout pattern can be described as a multifractal.