Seasonal Variation of Radon Concentrations in Russian Residential High-Rise Buildings

: Assessment of the annual radon concentration is often required in indoor radon surveys of territories and individual dwellings for comparison with reference levels, studying factors affecting radon accumulation in dwellings, assessment of exposure in epidemiological studies, etc. The indoor radon surveys were carried out in multistorey buildings in eight Russian cities using solid state nuclear track detectors with an exposure period of three months. For these surveys, the estimation of annual indoor radon concentration was required to compare radon levels in buildings of high-and low-energy-efﬁciency classes located in different cities. To develop approaches to seasonal normalization in high-rise buildings, long-term one-hour radon concentration series obtained applying radon-monitors in 20 ﬂats were analyzed. The dependency of indoor radon concentration on the indoor–outdoor temperature difference was studied taking into account the known natural, technogenic and anthropogenic factors affecting radon levels. The developed model of seasonal variations in multistorey buildings includes winter, summer, and demi-season periods, which differ both in ventilation intensity and dependency of radon concentration on the temperature difference. The developed model allows to estimate annual radon concentration taking into account the actual distribution of outdoor temperatures during the exposure of the track detectors. Radon entry from building materials a Outdoor radon concentration


Introduction
The relevance of studying radon ( 222 Rn) in indoor atmosphere is connected with findings on association of radon progeny inhalation with human health effects. Based on the review of uranium miners' study, BEIR VI Committee [1] has proved that exposure to radon gas and its decay products causes lung cancer. A strong causal association between residential radon exposure and lung cancer mortality is supported by pooled analysis of case control studies performed in Europe and North America [2,3]. A combined evaluation of estimates of residential case control studies meta-analysis and geographically aggregated data on lung cancer mortality and average indoor radon concentration yielded an excess risk of lung cancer of about 14% (90% CI: 0.10-0.18) per radon concentration 100 Bq/m 3 linked to an exposure period of 25-30 years [4]. The health risk assessment has demonstrated that radon is the second cause of cancer after smoking [1,5]. According to a special study of the Lancet journal, radon is ranked as the fifth environmental risk factor for human health in 2019 [6].
With regard to the ICRP [7,8] and WHO [9] recommendations, the IAEA has developed a strategy for protecting humans against indoor radon [10]. In particular, ICRP recommended to establish national indoor radon concentration reference levels in the range from 100 to 300 Bq/m 3 [8]. Protection of the population in the situation of exposure to radon in buildings requires reliable information on levels of indoor radon concentration. National, regional, and single home surveys with direct measurements of indoor radon concentration provide the basis for future decisions on control over the public exposure due to radon [8,10]. Radon surveys estimate the population radon exposure, identify radon-prone areas, and assess compliance with the reference level.
Seasonal variations of indoor radon concentration are mainly affected by relationship between diffusion and advection radon entry mechanisms, and dependency of ventilation rate and ventilation practice on climatic conditions. As shown by Arvela [30], in low-rise residential houses with natural ventilation, where the pressure difference-driven radon flow from soil is the dominant radon source, winter concentrations are normally higher than summer concentrations. Bossew et al. [15] consider the tendency to open windows in summer and keep them closed in winter, which leads to higher indoor accumulation of geogenic radon in winter. The normal seasonality is also observed in temperate and continental climate countries when radon diffusion from building materials is the dominant source of radon entry and the "normal" ventilation is plasticized [31,32]. Model of radon entry and accumulation in multi-flat buildings developed by Yarmoshenko et al. [33] for Russian climatic conditions predicts deviation from the normal seasonality in flats where low ventilation frequency is maintained. A wide range of magnitudes of seasonal variation between homes caused by variability of influencing factors is emphasized in many studies [18][19][20].
Seasonality of indoor radon concentration has to be considered for planning and analyzing the radon survey results. To know the annual average radon concentrations, measurements ideally should be made over a total period of twelve consecutive months. However, the measurements taken over the course of an entire year are impractical because of a long delay for the results and other reasons. Shorter exposure periods (typically three months) are widely used and included to indoor radon survey protocols in many countries [26,27,[34][35][36][37]. In the case of the shorter duration of measurements some efforts have to be made to estimate the annual radon concentrations.
The most common approach for season normalization is application of seasonal correction factors estimated by results of special study. The seasonal correction factors allow estimates of the annual mean radon concentrations in homes surveyed over several months by multiplying the measured value with the appropriate empirically determined factor. In the UK, seasonal correction factors were first derived by Wrixon et al. [24] and lately updated with estimation depending on the starting month and the measurement period [19,34]. A more simplified approach consists in estimating the winter/summer ratios, which is appropriate if integrated measurements are conducted during six months [38][39][40] or one season; usually winter or the heating period is recommended for radon surveys [12,22]. Tsapalov et al. [22] proposed to assess the compliance with reference level taking to account the uncertainty of annual average radon concentration estimation depending on the measurement duration. In some studies, only a moderate variation of seasonal radon averages from the annual average was found for measurements with duration of three or more months and season normalization was not applied [37,[41][42][43]. However annual average radon concentration is required for various purposes, in particular, studying the influence of geogenic and anthropogenic factors on indoor radon concentration.
Further investigation of radon seasonal variations is relevant for various combinations of radon sources and entry driving forces (diffusion and advection, soil, and building materials etc.,), climate, and living habits. In particular, there are only a few studies of the radon seasonal variations where multistorey houses (mostly below 5 storeys) were included in the consideration [15,30,38].
Multi-flat middle-and high-rise buildings are common in Russian cities and modern trends are such that more than 2/3 of the houses under construction are multistorey.
From the radiation protection point of view the attention to the new high-rise houses is especially attracted due to energy-saving measures implemented in modern building construction [44]. In buildings of high-energy-efficiency classes, the conditions are provided for lower air exchange rate and higher indoor radon accumulation [45]. In previous studies in Ekaterinburg, large city in the Urals region of Russia, application of radon-monitors with obtaining long-term, up to six months, radon series demonstrated great advantages for researching indoor radon accumulation in high-rise buildings. In particular the role of building materials as a source of indoor radon at high storeys was revealed [46]. Important information on ventilation patterns was obtained as well [31,47].
In this study, the long-term radon series together with data of integrated measurements were analyzed to study the seasonal variations of radon concentration and develop approaches to seasonal normalization in multistorey buildings.

Available Data on Radon Temporal Variation
The following data were used as experimental material to study the temporal and seasonal variations of radon concentration in multistorey buildings:

1.
Radon series with duration of more than two months, obtained in flats of multistorey buildings under normal operation (applying radon monitors); 2.
Results of a survey in a sample of flats in multistorey buildings in Chelyabinsk carried out in two seasons (applying solid-state nuclear track detectors); 3.
Results of a survey of flats in a multistorey building in Ekaterinburg, carried out in two seasons (applying solid-state nuclear track detectors).
Radon series measurements were performed in flats of staff of Institute of Industrial Ecology and their relatives, and colleagues from other cities. Two-season radon measurements were conducted in the flats of some participants of the large radon survey, earlier performed in Chelyabinsk [44]. The detailed indoor radon study in the multistorey building in Ekaterinburg was previously described by Vasilyev et al. [48]; the flats, where both winter and summer three-months measurements were performed, are included in this study.
Long-term serial one-hour measurements of radon concentration were performed using four AlphaGuard PQ2000Pro radon monitors (SAPHYMO GmbH, Germany), which were calibrated regularly using a NIST calibration source (SRM 4973-16). The characteristics of the flats in which the measurements were carried out are presented in the Table 1. Most of the measurements were carried out in Ekaterinburg, while single series were obtained in St. Petersburg and Chelyabinsk. The Figure 1 shows the duration of the series using the atmospheric temperature scale during the measurement period. The series included in this analysis cover temperature ranges between highest and lowest outdoor temperatures from 20 to 50 • C. Most of these series were previously used in studies of factors influencing the radon entry and accumulation in multistorey buildings [31]. This analysis also includes series that were not previously considered in the mentioned studies due to the low radon concentration and the inability to determine the intervals of controlled and uncontrolled air exchange. Some of the measurements were intentionally carried out in flats of high-energyefficiency buildings built over the past 10 years. High energy efficiency of a multistorey building was determined taking into account such characteristics as the type of building structure, low permeability of the building envelope, double glazing windows, balcony glazing, and others [44]. The results of serial measurements of radon concentration were used to assess the dependence of radon concentration on atmospheric temperature. For these purposes average radon concentrations were estimated in 3 • C outdoor temperatures intervals for each radon series.
Integrated measurements in Chelyabinsk were carried out in a sample of ten multistorey flat buildings in the periods of August-November 2019 and November 2020-January 2021. The houses in which the measurements were made belong to different types of building structure. The survey results are presented in Table 2. As can be seen from the table, the winter radon concentrations, on average, are more than two times higher than those obtained in the warmer period of the year. The coefficient of seasonal variation in 10 flats varied significantly from 1.44 to 3.51. The integrated measurements in flats in the multistorey building in Ekaterinburg were carried out in 2017. Twenty three flats were examined in two seasons, which was 23% of all flats in the building. The object was built in 2012 using monolithic concrete, with a low permeability of the shell and other solutions that provide a high class of energy efficiency. The data obtained at this object were previously reported in [47,48]. In particular, there was no significant tendency to radon concentration falling off with a height for flats above ground floor [48]. The average results of radon concentration measurements in the building are presented in Table 3. On average, radon concentrations in cold and warm seasons differ more than 1.5-1.6 times. In 26% of flats, the seasonal coefficient was in the range 0.96-1.1, that can be considered as low seasonal variations. As shown in the study [47,48], the variability of indoor radon concentrations and seasonality in this sample of flats are associated with different ventilation habits of dwellers.  The results of radon concentration-integrated measurements in two seasons in Chelyabinsk and Ekaterinburg were used to verify the model dependences of radon concentration on outdoor temperature.
Application of the developed model of temperature seasonal normalization is demonstrated in the example of a radon survey in Nizhny Novgorod. The survey is a continuation of the study Yarmoshenko et al. [44] with application of the same approach to sample formation and indoor radon measurements. The characteristics of the sample of flats and the results obtained are presented in the Table 4. Radon concentration measurements in Ekaterinburg, Chelyabinsk, and Nizhny Novgorod were carried out using solid-state nuclear track detectors RSKS type provided by Radosys [49]. To support the quality of measurements, RSKS detectors were calibrated both by the producer and within the framework of internal laboratory control of the Institute of Industrial Ecology using a system similar to that described by Mostafa et al. [50] placing RSKS detectors in a radon chamber with AlphaGuard PQ2000Pro. The detection limit for RSKS detectors for three-months exposure is defined as 5 Bq/m 3 .

Climate
In this study, data obtained in four cities of Russia located in different terrestrial and climatic zones were used. The main climatic characteristics are presented in the Table 5 and in Figure 2. The cities of Ekaterinburg and Chelyabinsk are located in a zone of continental climate, which is characterized by a significant difference between summer and winter temperatures. St. Petersburg is located on the shores of the Gulf of Finland of the Baltic Sea; seasonal variations of temperature and other climate parameters are more moderate in this city. Nizhny Novgorod is located in the center of the European part of Russia, which is the most populated Russian region. Therefore, Nizhny Novgorod can be considered the most typical Russian city in terms of climate.  Atmosphere temperature distributions in Ekaterinburg, Chelyabinsk, and Nizhny Novgorod are presented in Figure 3 according to the data of meteorological stations located in these cities (www.rp5.ru, accessed on 10 June 2021). The distributions of outdoor temperature are estimated for the whole year and the period of measurements performed in respective cities. The temperature distributions were used for temperature normalization of three-months indoor radon concentration measurements.

Seasonal Variation Model
When constructing a model of the dependence of radon concentration on the outdoor temperature, a basic physical model of radon entry into a flat in a multistorey building was used, whose separate parts are described in the previous studies [31,33,47]. The model is based on a theoretical concept of the mechanisms of indoor radon entry appeared in the 1980s and 1990s when the scientific community developed a general strategy on protection against indoor radon [51][52][53]. The key points of the model are as follows.

1.
The main mechanism of the radon entry into flats in multistorey buildings is diffusion from building materials. The main characteristic of the building material, which determines the rate of diffusion of radon entry, is the Ra-226 concentration. The variability of this characteristic within a large country like Russia can be quite significant. At the same time, the variability of other characteristics of the building materials that determine the diffusion is less significant.

2.
The main mechanism limiting the accumulation of radon is the ventilation of the building. There are two main modes of ventilation in a building: controlled and uncontrolled. In uncontrolled mode, all ventilation devices are closed and inactivated, while the air exchange rate (AER) depends on the permeability of the building envelope and the temperature difference. The controlled mode involves opening windows or turning on mechanical ventilation systems in order to achieve the AER desired by residents. The physical model describing ventilation takes into account the leakage area, temperature differences, and the frequency and duration of controlled ventilation [54,55].

3.
In continental climate during cold winter, below certain negative temperatures, controlled ventilation stops. As the outdoor temperature rises, the frequency and ventilation duration gradually increases. At summer temperatures, the ventilation frequency is maximum. At the same time, the temperature difference in the summer period is close to zero, so the AER is mainly determined by the wind pressure. 4.
In large multistorey buildings, there may be radon flow between the rooms of the building, caused by pressure differences. As a result, advective transport of radon from the non-living spaces with low air exchange to living rooms can be observed.

5.
From the point of view of radon concentration modeling, modern buildings of a high-energy-efficiency class have a lower AER in an uncontrolled mode.
Taking into account the above theoretical considerations, a model of the dependence of radon concentration in a multistorey building on the outdoor air temperature is suggested ( Figure 4). The temperature dependence of radon concentration differs for three temperature ranges. At low temperatures (T < T1), radon concentration is determined by the AER in an uncontrolled mode, which depends on the temperature difference [33,54,55]: where k leak is the temperature-independent leakage coefficient, T ind and T are the air temperatures indoor and outside the building, respectively. At T = T1, the maximum radon concentration is achieved due to the lowest AER for the entire period when air exchange occurs in an uncontrolled mode only. At high temperatures T > T3, the minimum radon concentration, Rn min , is achieved due to the maximum frequency and duration of ventilation. One of the factors determining the value of Rn min is the radon concentration in the atmosphere.
In the intermediate temperature range (spring/autumn), there is a gradual change in the frequency and duration of ventilation in the controlled mode, which increase with outdoor temperature. To describe the behavior of radon concentration in this range, the following equation is used: where 4.394 is the normalizing factor at which the conditions Rn (T2) = 0.9 (Rn max − Rn min ) + Rn min and Rn (T3) = 0.1 (Rn max − Rn min ) + Rn min are met. Taking into account the previous studies, some values were accepted for the model parameters included in Equation (1) ( Table 6). To move from a descriptive model to a physical one, it is necessary to determine the relationship between Rn min , radon in the atmosphere, and the AER. As shown in Table 6, taking into account that in the summer period the air temperature can be quite high, in this study it is assumed that AER (T3) = 1 h −1 (taking into account the air exchange is partly associated with the wind pressure). Table 6. Accepted parameters or the model.

Parameter
Value or Equation The radon concentration, which is expected to be obtained during the measurement period with known temperature frequency distribution, is equal to: where Rn mod (T i , {m}) is radon concentration at temperature T i , calculated by the model, {m} is a set of model parameters.
The set of parameters {M} appropriate for describing the dependence of radon concentration on outdoor temperature in a certain room was obtained using the condition: where RnOBS is the measured radon concentration. Additional boundary conditions, assumptions and approximations can also be taken into account to choose some of the model parameters.
After the function Rn mod (T i , {M}) is determined, the annual average radon concentration is calculated using the equation: Rn an = ∑ i Rn mod (T i , {M})·P(T i ), (5) where Rn an is the annual average radon concentration.
Seasonal correction factors can be estimated for regions with certain function P (T). For example, for the three-months integrated radon concentration measurements, 12 values can be calculated for each months of the measurement starting, by a manner similar to that implemented in the UK [19,34]. If the temperature distribution over the measurement period starting in j-th month is described by Pm j (T i ), then the expected radon concentration for a measurement period of three months RnEXP j is defined by equation: (6) and the seasonal correction factor k j is equal to the ratio: kj = Rn an RnEXP j (7)

Estimating the Parameters of The Descriptive Model
The results of the fitting of the parameters of the model dependence of radon concentration on temperature according to the analysis of the radon concentration time series in flats of multistorey buildings are presented in Table 7. For 20 series, the correlation coefficient between the observed radon concentration and the fitted model values is in the range from 0.75 to 0.99. Two series were too short in temperature range and were excluded from further analysis. The coverage of the temperature range during the monitoring measurements made it possible to determine the full set of parameters (T1, T2, T3, Rn min , Rn max ) describing seasonal variations in 10 flats. In 14 flats, sufficiently long series were obtained to determine all the parameters except the temperature T1. In all 20 series included in further analysis, the parameters T2 and Rn max were determined.
The analysis of the obtained time series and the fitted parameters of the model made it possible to distinguish two main patterns and several subtypes of seasonal variations. The identified patterns and subtypes are schematically shown in Figure 5. The first type of seasonal variation includes flats in which the radon concentration is higher in the cold period than in the warm one (pattern A, including subtypes A1-A3). Subtypes A1-A3 differ in either temperature T2 or Rn min parameters. In flats assigned to subgroup A1, the radon concentration Rn min was found corresponding to the radon entry from building materials aD, taking into account the radon concentration in the atmosphere. In subgroup A2, Rn min exceeds the value that can be expected under diffusion entry domination. In subgroup A3, there is a significant shift in temperature T2 to higher values. Pattern B includes the cases when radon concentration in summer is higher than in winter, which, within the framework of a descriptive model, can be described by assigning a value above 5 • C (5 • C and 10 • C in our examples) to parameter T1. Pattern B-2 flats (10%). For subgroups A1 and A2, which together make up 70% of flats, the overall average is T2 = −1.1 C.
Thus, for the flats of the most common type, the following rounded temperature parameters of the seasonal variation model can be accepted: T1 = −5 • C, T2 = −1 • C, T3 = 15 • C.

Verification of the Model
To verify the obtained model of seasonal variations, independent radon survey data in two samples of flats, carried out in the summer and winter periods of the year were used. For each sample of seasonal data, the model parameters Rn max and Rn min were estimated, at which the average annual radon concentration, estimated based on summer and winter measurements, are equal ( Table 8). The table also shows the values of the radon diffusion entry rate from building materials aD. As it can be seen in Table 8, the values of Rn max and Rn min obtained for the flat building in Ekaterinburg generally correspond to those observed in other buildings for which there is a series of measurements (Table 7). For Chelyabinsk, a negative value of radon concentration was obtained, that has no objective meaning. It can be assumed that the negative value is the result of an error in the radon measurement, high sampling heterogeneity, and other reasons. Taking into account the known factors affecting the radon entry and accumulation, the minimum radon concentration is derived by the entry from building materials (aD) and the dilution with atmospheric air with a certain radon content at the AER taken equal to 1 h −1 . If more justified Rn min is accepted, then different average annual radon concentrations are estimated applying summer and winter measurements: 34 and 44 Bq/m 3 , respectively. The relative deviations of the values calculated in this way from that obtained in the descriptive model are 10% and 16%.
For the building in Ekaterinburg, a relatively high Rn min value is observed. With the corresponding parameters aD, AER = 1 h −1 and radon concentration in the atmosphere of 5 Bq/m 3 , in order to achieve Rn min = 108 Bq/m 3 , it is necessary to take an additional entry equal to about 80 Bq/ (m 3 h). As shown in the study of indoor radon temporal variations in multistorey buildings [31], this entry may be associated with advective flows inside the building from service premises with low air exchange to living rooms. High advective entry due to flows between different parts of a house can be observed in buildings with low air exchange in uncontrolled mode and high Ra−226 content in building materials [46]. A high Rn max value is an indicator of the presence of a significant advective redistribution of air with a high radon content between the parts of a building. In Chelyabinsk buildings, there are no conditions for the formation of such entry. Table 9 shows the calculation of the average annual radon concentration in a sample of buildings in Nizhny Novgorod, taking into account the results of measuring radon concentration during three months starting in November. The table shows the observed temperature distributions during the year, calculated from the archived observational data for two years from 7 April 2019 to 6 April 2021, and the distributions of temperatures during the exposure period (from 24 October 2020 to 4 February 2021). The parameters of the model of radon concentration seasonal variations for Nizhny Novgorod were selected taking into account the requirement (4) such as: Rn min = 8.6 Bq/m 3 , Rn max = 40 Bq/m 3 , a D = 3.6 Bq/(m 3 h). After seasonal normalization, the average annual radon concentration for a sample of flats in multistorey buildings in Nizhny Novgorod was estimated as 23 Bq/m 3 . Table 9. Calculation of the average annual radon concentration in a sample of buildings in Nizhny Novgorod. For the climatic conditions of the city of Nizhny Novgorod, seasonal correction factors were estimated using Equation (7) for calculating the average annual radon concentration based on three-months measurements using integrated methods ( Figure 6). For comparison, the figure also shows the data obtained for the UK [34]. As can be seen in the figure, the range of seasonal correction factors in the Russian city is higher, which is apparently associated with greater temperature variation, including higher temperatures in summer and lower in winter. On the whole, the character of the seasonal variations obtained in Russia and Great Britain is consistent.

Discussion
This paper presents the model of seasonal variations of radon concentration in flats of multistorey buildings based on the results of analysis of radon concentration long-term one-hour series obtained using radon monitors. The developed model makes it possible to use an approach for assessing the average annual radon concentration based on the temperature normalization. Temperature normalization has certain advantages over other methods, as shown in [19,34], which are easily realizable with the availability of electronic weather archives.
The analysis of radon series gained clear and physically understandable dependence of radon concentration on the atmosphere temperature. When air exchange is maintained at a level corresponding to comfortable microclimate for most people, radon concentration in flats is higher in winter than in summer. This character of seasonal variations in radon concentration in multistorey buildings corresponds to the predominant radon entry from building materials. The parameters of the temperature dependence of radon concentration in general do not depend on the climatic conditions within the zones of continental climates, in which cold winters and warm summers occur. The T1 parameter determines the temperature below which the vents in dwellings either do not open at all, or open slightly for a few minutes. Under these conditions, the use of a simple model is justified, in which radon concentration is proportional to indoor outdoor temperature difference. As shown in this paper, on average, this temperature T1= −5 • C. The middle of the demiseason interval T2-T3 is approximately 7 • C, which corresponds to the end of the heating season. The average daily temperature of 8 • C is the value at which local government regulatory decision to turn off the municipal central heating system is issued in Russia. At a temperature above T3, which, according to this study, is approximately 15 • C, in flats of multistorey buildings the frequency and duration of ventilation approach maximum values. The adopted parameters of seasonal variations depend on the preferences of residents about a comfortable temperature in dwellings and the norms for regulating the heating season that are the same for the whole of Russia, therefore, the presented model of seasonal variations can be used throughout Russia.
A life style in which ventilation remains at a low level in the summer and radon concentration in the warm period may be higher than in the cold one is not typical (pattern B). Therefore, a model that takes into account the more common preferences on ventilation during a warm period (pattern A) can be accepted to describe the average seasonal variations in large randomized samples of buildings.
Less significant deviations from the typical character of the dependence of radon concentration on temperature may be associated with other features of the life style or living habits. In the available sample of flats, in which measurements of long-term series of radon concentrations were carried out, there was a relative increase of the parameter T2 in every fifth flat. This parameter shows the ambient air temperature at which the transition between winter and demi-season ventilation modes occurs. It can be assumed that at temperatures above −1 • C, residents of these flats prefer to control the indoor air temperature by adjusting a heating system. Currently, only small parts of buildings are equipped with the technical ability to regulate the temperature of the heating batteries in Russia. Therefore, in the future, with a greater prevalence of energy-saving systems in buildings, the typical value of T2 parameter in the seasonal model may increase.
The applicability of the proposed model of radon seasonal variations is shown in flats of multistorey buildings with the results of radon concentration measurements by integrated methods. The discrepancies that were obtained with the results of the survey of a sample of ten flats in Chelyabinsk are insignificant. Using the example of energy-efficient multistorey building in Ekaterinburg, it was shown that taking into account the features of the advective transport of radon inside the building, full consistency of the model and observed results can be obtained.
The developed model of seasonal variations and the corresponding approach to the temperature normalization of radon concentration measurements in flats of multistorey buildings in Russia are generally consistent with the approaches developed in other studies [12,15,19,25]. The nature of seasonal variations is common for most countries, in which the seasons are defined as warm and cold (with dedicated heating period).
The proposed approach to the seasonal normalization of radon concentrations is planned to be used in the future to analyze the results of radon study in buildings of various types in large Russian cities [44]. This study aims to determine the effect of a building's energy efficiency on radon accumulation and compare radon levels in buildings located in different regions. Given that more than 80% of the Russian population of large urbanized region lives in multistorey buildings, this model can be applied to randomized samples surveyed in such cities.

Conclusions
The main conclusions of the study are as follow: 1.
Indoor radon seasonal variation can be described taking into account the natural, technogenic, and anthropogenic factors affecting the radon levels. 2.
The model indoor radon seasonal variation developed for the situation of radon entry from building materials in multistorey buildings associates the seasonal radon concentration to frequency and duration of controlled ventilation. For continental climate, the radon temperature dependence can be typically modeled using three temperatures which determine periods of different ventilation pattern as follow: summer ventilation pattern at outdoor temperature above 15 • C, demi-season ventilation at outdoor temperature between 15 and −1 • C, total insulation with temperature below −5 • C.

3.
Estimation of the annual average indoor radon concentration using the results of the integrated radon measurements during period shorter than twelve months using the developed model is performed by comparing the frequency distributions of atmospheric temperature during the period and a whole year. 4.
Developed model of indoor radon seasonal variations allows the temperature normalization to estimate average annual radon concentration obtained in the radon surveys in Russian cities in continental climate zones. Institutional Review Board Statement: Not applicable.