PERFORMANCE OF HANDHELD NAI(TL) SPECTROMETERS AS DOSIMETERS BY LABORATORY AND FIELD DOSE RATE MEASUREMENTS.

In the framework of the IAEA Coordinated Research Project (CRP) J02012 on 'Advancing Radiation Detection Equipment for Detecting Nuclear and Other Radioactive Material Out of Regulatory Control', the properties of two commercial instruments (1) InSpector 1000 analyzer (Canberra), with a 2″ × 2″ NaI(Tl) scintillator and (2) RIIDEYE M-G3 analyzer (Thermo Scientific), with a 3″ × 3″ NaI(Tl) scintillator, were evaluated as dosimeters by laboratory and field measurements. In the Ionizing Radiation Calibration Laboratory (IRCL) of the Greek Atomic Energy Commission, the NaI(Tl) spectrometers were tested in order to measure Ambient gamma Dose Equivalent Rate (ADER). The NaI(Tl) scintillators were irradiated in a homogeneous field with 662 keV photons with different ADER values from 0.17 to 100 µSv h-1 at 0° incidence (irradiation field perpendicular to the detector's front window) and at 90° incidence. For each irradiation, the measured ADER by the spectrometers and the 'true' ADER values (provided by the IRCL) were compared. In addition, the angular dependence (0-359°) of the ADER response of the spectrometers was studied with a 152Eu source placed at 1, 2 and 3 m from the spectrometers. The ADER dependence as function of the distance from the 152Eu source (at 0° incidence) measured by the two detectors was compared with the theoretical one. In the field studies, ADER was measured by the spectrometers at seven locations belonging to the Greek Early Warning System Network (which is based on Reuter-Stokes ionization chambers). These locations have different ADER values ranging from 20 to 120 nSv h-1. In these locations, gamma ADER were also deduced (1) by in situ gamma spectrometry measurements with portable Germanium HPGe detectors and (2) by the Reuter-Stokes ionization chambers (by subtraction of the cosmic radiation). Gamma dose measurements were also performed with the InSpector 1000 and RIIDEYE M-G3 detectors in 25 locations (beaches) of Northern Greece and compared with the ADER values deduced by sand sample analysis with gamma spectroscopy. Beaches with sand are good candidates for such type of measurements since they are commonly flat and in principle the natural radionuclides are homogenously distributed.

RESPONSE OF THE GREEK EARLY WARNING SYSTEM REUTER-STOKES IONIZATION CHAMBERS TO TERRESTRIAL AND COSMIC RADIATION EVALUATED IN COMPARISON WITH SPECTROSCOPIC DATA AND TIME SERIES ANALYSIS.

The Telemetric Early Warning System Network of the Greek Atomic Energy Commission consists mainly of a network of 24 Reuter-Stokes high-pressure ionization chambers (HPIC) for gamma dose rate measurements and covers all Greece. In the present work, the response of the Reuter-Stokes HPIC to terrestrial and cosmic radiation was evaluated in comparison with spectroscopic data obtained by in situ gamma spectrometry measurements with portable hyper pure Germanium detectors (HPGe), near the Reuter-Stokes detectors and time series analysis. For the HPIC detectors, a conversion factor for the measured absorbed dose rate in air (in nGy h-1) to the total ambient dose equivalent rate H*(10), due to terrestrial and cosmic component, was deduced by the field measurements. Time series analysis of the mean monthly dose rate (measured by the Reuter-Stokes detector in Thessaloniki, northern Greece, from 2001 to 2016) was performed with advanced statistical methods (Fast Fourier Analysis and Zhao Atlas Marks Transform). Fourier analysis reveals several periodicities (periodogram). The periodogram of the absorbed dose rate in air values was compared with the periodogram of the values measured for the same period (2001-16) and in the same location with a NaI (Tl) detector which in principle is not sensitive to cosmic radiation. The obtained results are presented and discussed.

RADON MIGRATION IN SOIL AND ITS RELATION TO TERRESTRIAL GAMMA RADIATION IN DIFFERENT LOCATIONS OF THE GREEK EARLY WARNING SYSTEM NETWORK.

Radon concentration as function of the soil depth was measured during the years (2011-2015), in a location of the Aristotle University campus. Radium distribution in soil was found constant. On the contrary, as expected, radon concentration increases with soil depth. The experimental distribution was reproduced by solving the general transport equation (diffusion and advection). From the general radon migration (diffusion and advection) equation it was indirectly deduced, from the measured radon profile in the soil, the radon exhalation rate from the soil (26.7 ± 4.5 Bq m-2 h-1). In the same location, during 2010-2015, 113 direct radon exhalation measurements were performed and give an average value for the 6 years of 21.1 ± 3 Bq m-2 h-1. The comparison between the radon exhalation rates deduced by the indirect and direct method indicates the validity of the diffusion advection model predictions concerning the radon exhalation rate from soil. The relation between radon migration in soil and terrestrial gamma radiation was studied. In particular, in the present study was investigated the correlation between gamma radiation 1 m above soil and radon exhalation in six locations of the Greek early warning system network. A positive correlation between gamma dose rate in air and radon exhalation rate from soil was found.

Environmental impact of phosphogypsum stockpile in remediated Schistos waste site (Piraeus, Greece) using a combination of γ-ray spectrometry with geographic information systems.

From 1979 to 1989, ten million tons of phosphogypsum, a waste by-product of the Greek phosphate fertilizer industry, was disposed into an abandoned limestone quarry in Schistos former waste site, Piraeus (Greece). The quarry has been recently closed and remediated using geomembranes and thick soil cover with vegetation. A part of the deposited phosphogypsum has been exposed due to intense rainfall episodes leading to concerns about how could potentially released radioactivity affect the surrounding environment. This study seeks to assess the environmental impact of the phosphogypsum deposited in the Schistos quarry, using laboratory-based γ-ray spectrometry measurements and geographical information systems. Radioactivity concentrations were mapped onto spatial-data to yield a spatial-distribution of radioactivity in the area. The data indicate elevated (226)Ra concentrations in a specific area on the steep south-eastern cliff of the remediated waste site that comprises uncovered phosphogypsum and is known to be affected by local weather conditions. (226)Ra concentrations range from 162 to 629 Bq/kg, with an average activity being on the low side, compared to the global averages for phosphogypsum. Nevertheless, the low environmental risk may be minimized by remediating this area with geomembranes and thick soil cover with vegetation, a technique, which has worked successfully over the remainder of the remediated quarry.

PERIODICITY ANALYSIS OF GAMMA RADIATION MEASUREMENTS IN THESSALONIKI, NORTHERN GREECE, IN THE PERIOD 1988-2015.

Gamma radiation measurements were performed during the last 27 y, starting from 1988, with a NaI(Tl)-based Xetex 501A radiation monitor located outside the Nuclear Technology Laboratory of the Aristotle University of Thessaloniki in Northern Greece, and a time series was created. Measurements were also performed in the same place during 1995-98 and 2013-15 with portable high purity germanium (HPGe) detector. The total absorbed dose rate in air decreases exponentially with time. The total absorbed dose rate in air is the sum of the gamma dose rates due to (1) uranium series, (2) thorium series, (3) 40K and (4) 137Cs (due to the Chernobyl accident). In addition, a small contribution due to cosmic radiation is measured by the radiation monitor. From the time-dependence measurements with the HPGe detector, it was found that the time dependence of the absorbed gamma dose rate in air due to (1) uranium series, (2) thorium series and (3) 40K is quite constant. On the contrary, gamma dose rate due to 137Cs decreases exponentially with an effective half-life (t½) of ∼13.5 y, stronger than expected due to the natural decay of 137Cs. Time series analysis of the mean monthly total absorbed dose rate in air was performed. Fourier analysis reveals several periodicities, and applying Zhao-Atlas-Marks transform unravels the time distribution of those periodicities. There are three main discernible periodicities: 12 ± 0.2, 42.3 ± 2.9 and 53.2 ± 3.2 months. One of them is of a seasonal character (annual cycle) and can be linked to seasonal atmospheric variations and is strongly visible from 1988 to 2002 and 2008 to 2014. The other two (42.3 ± 2.9 and 53.2 ± 3.2 months) were found to be also related to meteorological parameters (air temperature), and they were very intense during the years 2002-4 when the annual periodicity was weak. Apart from the three main periodicities, there are also four others (14.7, 18.6, 21.3 and 27.3 months) with lower magnitudes; of which, three agree well with literature data periodicities in solar activity. Different possible mechanisms that can influence the gamma radiation measurements, due to solar activity, were discussed.

Measurements and modelling of 137Cs distribution on ground due to the Chernobyl accident: a 27-y follow-up study in Northern Greece.

Following the Chernobyl accident, an area of ∼1000 m(2) in the University farm of the Aristotle University of Thessaloniki was considered as a test ground for radioecological measurements. The radiocesium deposition in this area, due to the Chernobyl accident, was 20 kBq m(-2). The profile of (137)Cs in the soil of this area was measured systematically from 1987 to 2012. The form of the profile has changed over the years. During the 1987-2000 period the (137)Cs distribution was reproducible by a sum of two exponentials. However, at least since 2005 the (137)Cs distribution can be successfully fitted by a single exponential function. The long-time (∼27 y) evolution study of the (137)Cs distribution in soil permit one to extract with the use of a simple compartment model, the mean vertical migration velocity of (137)Cs. Vertical migration of (137)Cs in soil is a very slow process. The mean vertical migration velocity is estimated to be 0.14 cm y(-1).The relative good comparison between the time dependence of the (137)Cs distribution in soil and the model predictions indicate that the simple model used is realistic.

Implementation of CFD modeling in the performance assessment and optimization of secondary clarifiers: the PVSC case study.

The water industry and especially the wastewater treatment sector has come under steadily increasing pressure to optimize their existing and new facilities to meet their discharge limits and reduce overall cost. Gravity separation of solids, producing clarified overflow and thickened solids underflow has long been one of the principal separation processes used in treating secondary effluent. Final settling tanks (FSTs) are a central link in the treatment process and often times act as the limiting step to the maximum solids handling capacity when high throughput requirements need to be met. The Passaic Valley Sewerage Commission (PVSC) is interested in using a computational fluid dynamics (CFD) modeling approach to explore any further FST retrofit alternatives to sustain significantly higher plant influent flows, especially under wet weather conditions. In detail there is an interest in modifying and/or upgrading/optimizing the existing FSTs to handle flows in the range of 280-720 million gallons per day (MGD) (12.25-31.55 m(3)/s) in compliance with the plant's effluent discharge limits for total suspended solids (TSS). The CFD model development for this specific plant will be discussed, 2D and 3D simulation results will be presented and initial results of a sensitivity study between two FST effluent weir structure designs will be reviewed at a flow of 550 MGD (∼24 m(3)/s) and 1,800 mg/L MLSS (mixed liquor suspended solids). The latter will provide useful information in determining whether the existing retrofit of one of the FSTs would enable compliance under wet weather conditions and warrants further consideration for implementing it in the remaining FSTs.

Follow-up study of indoor radon in Greek buildings.

The Nuclear Technology Laboratory of the Aristotle University of Thessaloniki (NTL-AUTh) and the Greek Atomic Energy Commission (GAEC) have a continuous collaboration on indoor radon measurements in Greek buildings since 1999. In the present work, the existing database was enriched with 590 indoor radon measurements in 295 houses and 76 indoor radon measurements in 38 workplaces. In total in the present work, 1948 indoor radon measurements in 974 buildings performed by the NTL-AUTh and GAEC from 1999 to 2012 in 8 of the 13 administrative regions of Greece are presented and discussed.

Environmental measurements and inspections on imported foods and feedstuffs in Greece after the Fukushima accident.

The radionuclides released during the accident at the Fukushima Daichii nuclear power plant following the Tohoku earthquake and tsunami on 11 March 2011 were dispersed in the whole north hemisphere. Traces of (131)I, (134)Cs and (137)Cs reached Greece and were detected in air, grass, sheep milk, ground deposition, rainwater and drainage water. Members of Six Greek laboratories of the national network for environmental radioactivity monitoring have collaborated with the Greek Atomic Energy Commission (GAEC) and carried out measurements during the time period between 11 March 2011 and 10 May 2011 and reported their results to GAEC. These laboratories are sited in three Greek cities, Athens, Thessaloniki and Ioannina, covering a large part of the Greek territory. The concentrations of the radionuclides were studied as a function of time. The first indication for the arrival of the radionuclides in Greece originating from Fukushima accident took place on 24 March 2011. After 28 April 2011', concentrations of all the radionuclides were below the minimum detectable activities (<10 µBq m(-3) for (131)I). The range of concentration values in aerosol particles was 10-520 µBq m(-3) for (131)I, 10-200 µBq m(-3) for (134)Cs and 10-200 µBq m(-3) for (137)Cs and was 10-2200 µBq m(-3) for (131)I in gaseous phase. The ratios of (131)I/(137)Cs and (134)Cs/(137)Cs concentrations are also presented. For (131)I, the maximum concentration detected in grass was 2.2 Bq kg(-1). In the case of sheep milk, the maximum concentration detected for (131)I was 2 Bq l(-1). Furthermore, more than 200 samples of imported foodstuff have been measured in Greece, following the EC directives on the inspection of food and feeding stuffs.

Gamma radiation measurements and Monte Carlo computations following myocardial perfusion imaging with 201Tl.

In the current study, the time-dependent retention of (201)Tl-thallous chloride (111 MBq) was measured in a 56-y-old man undergoing myocardial perfusion imaging. For 23 d following the (201)Tl injection, total-body retained activity was measured by (i) in situ gamma spectrometry using a portable high-purity germanium (HPGe) detector and (ii) ex vivo urine radioassay using a shielded HPGe detector. The time-dependent decrease in total-body activity followed a monoexponential function, exp(-0.011 t), with an excellent correlation (R(2)=0.9988) between the experimental data and the fitted values. The effective half-life, Teff, of (201)Tl (physical half-life, Tph: 72.9 h) was therefore 63 h and the biological half-life, Tb, 463 h=19.3 d, identical to those measured in the same patient in 1997 (i.e. 14 y ago). The time-dependent decrease in the urine activity concentration, which followed a monoexponential function, exp(-0.0115 t), corroborated the foregoing results. The correlation (R(2)=0.9939) between the experimental data and the fitted values was again excellent. The effective half-life, Teff, was 60.26 h and the biological half-life, Tb, 348 h=14.5 d. Monte Carlo simulation using a simple model of the patient as a unit-density cylinder filled with water and containing a uniform distribution of (201)Tl yielded photon flux results in reasonable agreement with the measured data.

Critical modeling parameters identified for 3D CFD modeling of rectangular final settling tanks for New York City wastewater treatment plants.

New York City Environmental Protection is in the process of incorporating biological nitrogen removal (BNR) in its wastewater treatment plants (WWTPs) which entails operating the aeration tanks with higher levels of mixed liquor suspended solids (MLSS) than a conventional activated sludge process. The objective of this paper is to discuss two of the important parameters introduced in the 3D CFD model that has been developed by the City College of New York (CCNY) group: (a) the development of the 'discrete particle' measurement technique to carry out the fractionation of the solids in the final settling tank (FST) which has critical implications in the prediction of the effluent quality; and (b) the modification of the floc aggregation (K(A)) and floc break-up (K(B)) coefficients that are found in Parker's flocculation equation (Parker et al. 1970, 1971) used in the CFD model. The dependence of these parameters on the predictions of the CFD model will be illustrated with simulation results on one of the FSTs at the 26th Ward WWTP in Brooklyn, NY.

Environmental radioactivity measurements in Greece following the Fukushima Daichi nuclear accident.

Since the double disaster of the 9.0 magnitude earthquake and tsunami that affected hundreds of thousands of people and seriously damaged the Fukushima Daichi power plant in Japan on 11 March 2011, traces of radioactive emissions from Fukushima have spread across the entire northern hemisphere. The radioactive isotope of iodine (131)I that was generated by the nuclear accident in Fukushima arrived in Greece on 24 March 2011. Radioactive iodine is present in the air either as gas or bound to particles (aerosols). The maximum (131)I concentrations were measured between 3 and 5 April 2011. In aerosols the maximum (131)I values measured in Southern Greece (Athens) and Northern Greece (Thessaloniki) were 585±70 and 408±61 µΒq m(-3), respectively. (131)I concentrations in gas were about 3.5 times higher than in aerosols. Since 29 April 2011, the (131)I concentration has been below detection limits. Traces of (137)Cs and (134)Cs were also measured in the air filters with an activity ratio of (137)Cs/(134)Cs equal to 1 and (131)I/(137)Cs activity ratio of about 3. Since 16 May 2011, the (137)Cs concentration in air has been determined to be about the same as before the Fukushima accident. Traces of (131)I were also measured in grass and milk. The maximum measured activity of (131)I in sheep milk was about 2 Bq l(-1) which is 5000 times less than that measured in Greece immediately after the Chernobyl accident. The measured activity concentrations of artificial radionuclides in Greece due to the Fukushima release, have been very low, with no impact on human health.

A comparison study between radon concentration in schools and other workplaces.

The Nuclear Technology Laboratory of the Aristotle University of Thessaloniki has since 1999 an open research project of indoor radon measurements in Greek workplaces. Since now 1380 measurements in 690 workplaces have been performed. Most (75 %) of the workplaces were offices in schools. The remaining 25 % were offices, mainly in public buildings. In the present study, a possible correlation between radon concentration in schools and other workplaces is investigated and discussed.

Indoor radon levels in Greek schools.

Radon and gamma dose rate measurements were performed in 512 schools in 8 of the 13 regions of Greece. The distribution of radon concentration was well described by a lognormal distribution. Most (86%) of the radon concentrations were between 60 and 250 Bq m(-3) with a most probable value of 135 Bq m(-3). The arithmetic and geometric means of the radon concentration are 149 Bq m(-3) and 126 Bq m(-3) respectively. The maximum measured radon gas concentration was 958 Bq m(-3). As expected, no correlation between radon gas concentration and indoor gamma dose rate was observed. However, if only mean values for each region are considered, a linear correlation between radon gas concentration and gamma dose rate is apparent. Despite the fact that the results of radon concentration in schools cannot be applied directly for the estimation of radon concentration in homes, the results of the present survey indicate that it is desirable to perform an extended survey of indoor radon in homes for at least one region in Northern Greece.

Development of a flocculation sub-model for a 3-D CFD model based on rectangular settling tanks.

To assess performance and evaluate alternatives to improve the efficiency of rectangular Gould II type final settling tanks (FSTs), New York City Department of Environmental Protection and City College of NY developed a 3D computer model depicting the actual structural configuration of the tanks and the current and proposed hydraulic and solids loading rates. Fluent 6.3.26™ was the base platform for the computational fluid dynamics (CFD) model, for which sub-models of the SS settling characteristics, turbulence, flocculation and rheology were incorporated. This was supplemented by field and bench scale experiments to quantify the coefficients integral to the sub-models. The 3D model developed can be used to consider different baffle arrangements, sludge withdrawal mechanisms and loading alternatives to the FSTs. Flocculation in the front half of the rectangular tank especially in the region before and after the inlet baffle is one of the vital parameters that influences the capture efficiency of SS. Flocculation could be further improved by capturing medium and small size particles by creating an additional zone with an in-tank baffle. This was one of the methods that was adopted in optimizing the performance of the tank where the CCNY 3D CFD model was used to locate the in-tank baffle position. This paper describes the development of the flocculation sub-model and the relationship of the flocculation coefficients in the known Parker equation to the initial mixed liquor suspended solids (MLSS) concentration X0. A new modified equation is proposed removing the dependency of the breakup coefficient to the initial value of X0 based on preliminary data using normal and low concentration mixed liquor suspended solids values in flocculation experiments performed.

Experimental and theoretical study of radon distribution in soil.

Radon concentration as a function of the soil depth (0-2.6 m) was measured during the years 2002-2003 and 2003-2004 on the Aristotle University campus. Radium distribution in soil was found constant. On the contrary, as expected, radon concentration increased with soil depth. However, the radon concentration did not follow the well known monotonous increase, which levels off to a saturation value. In both radon distributions, radon concentration increased up to a soil depth of about 80 cm, seemed to remain constant at depths of 80-130 cm, and then increased again. The experimental distribution was reproduced by solving the general transport equation (diffusion and advection). The main finding of the numerical investigation is that the aforementioned, experimentally observed, profile of radon concentration can be explained theoretically by the existence of two soil layers with different diffusion-advection characteristics. Soil sample analysis verified the existence of two different soil layers. Different boundary conditions of the radon concentration at the soil surface were used for the solution of the diffusion-advection equation. It was found that the calculated radon concentration in the soil is, away from the soil surface, the same for the two boundary conditions used. However, from the (frequently used) boundary condition of zero radon concentration at the soil surface, the experimental profile of the radon concentration at the soil surface cannot be deduced. On the contrary, with more appropriate boundary conditions the radon concentration at the soil surface could be deduced from the experimental profile. The equivalent diffusion coefficient could be uncovered from the experimental profile, which can then be used to estimate the radon current, which is important, for example, for the estimation of radon entrance to dwellings.

Indoor radon measurements in areas of northern Greece with relatively high indoor radon concentrations.

Indoor radon concentrations were measured in 77 schools of the prefecture of Xanthi in northern Greece. The arithmetic mean radon concentration is 231 Bq m(-3) with a range between 45 and 958 Bq m(-3). Thirty five schools have mean radon concentration above 200 Bq m(-3) and nine schools have mean radon concentration above 400 Bq m(-3). From continuous radon gas measurements in the school with a relative higher radon concentration (958 Bq m(-3)) was deduced the 'true' radon concentration, defined as the radon concentration in the school during the hours of the presence of teachers and scholars. The mean 'true' radon concentration for a time period of about 2 weeks was 104 Bq m(-3). The mean radon concentration for the same 2 weeks was seven times higher (700 Bq m(-3)). Greek and International regulations for radon in workplaces refer to only the mean annual radon concentration. It would be preferable for schools to replace the mean annual radon concentration with the 'true' radon concentration.

Measurements and Monte Carlo calculations of photon energy distributions in MAYAK PA workplaces.

Photon energy distributions were measured in different workplaces of the Mayak Production Association (MPA), which was the first plutonium production plant in the former Soviet Union. In situ gamma spectrometry measurements were performed with a portable germanium detector. The spectral stripping method is used for the conversion of the in situ gamma-ray spectra to photon fluence rate energy distribution. This method requires the simulation of the portable germanium detector, which has been performed based on the MCNP code of Los Alamos. Measured photon fluence rate energy distributions were compared with calculated photon energy distributions (with the MCNP code) in two different workplaces: in the first workplace the geometry exposure was known. On the contrary, in the second workplace, as in most workplaces of MPA, the exposure geometry was unknown. The results obtained from the comparison between the experimental and calculated photon fluence rate energy distributions are presented and discussed.

Follow up study of radiocesium contamination in a Greek forest ecosystem.

Radiocesiun distribution in the different parts of a Quercus conferta Kit ecosystem in Northern Greece was measured in 2005-2006, twenty years after the Chernobyl accident. The comparison between the results of this study and those previously measured (1993-1995) in the same ecosystem gives information about the long-term behavior of 137Cs in forest ecosystems. The major part of the 137Cs inventory is still in the upper layers of the soil. The radiocesium distribution in soil is fixed and has been in equilibrium at least since 1993, when the first measurements were performed. The major contamination mechanism of leaves and wood is root uptake.

Twenty-year follow-up study of radiocesium migration in soil.

The profile of (137)Cs present in undisturbed soil due to the Chernobyl accident was measured repeatedly for approximately 20 y. The vertical migration of (137)Cs in soil is a very slow process. The mean vertical migration velocity is estimated at approximately 0.1-0.2 cm y(-1). A method based on in situ gamma spectrometry measurements and Monte Carlo computations, aimed at estimating the profile of (137)Cs without performing any soil sampling, is investigated.