Blog

Novel environmental monitoring detector for discriminating fallout and airborne radioactivity | Scientific Reports

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  13, Article number: 22550 (2023 ) Cite this article Labr3 Scintillator

Novel environmental monitoring detector for discriminating fallout and airborne radioactivity | Scientific Reports

Early warning networks are used for detecting abnormal radioactivity levels in the environment. State-of-the-art networks are equipped with both dose rate detectors and spectrometric stations. Current networks don’t automatically discriminate between radioactivity on the ground and in the air. A novel directional sensing gamma radiation detector utilizing a collimated phoswich scintillator was developed. The signals from the two scintillator materials are separated using a pulse shape discrimination. The separated signals are employed to determine the radioactivity concentrations on the ground and in the air assuming specific concentration distributions. Limitations related to imperfect directional sensing and dead time are discussed.

Early warning networks are an important part of nuclear accident preparedness. One of their tasks is timely detection of abnormal radioactivity levels in the environment. Networks were established1 in many European countries after the Chernobyl accident2, but e.g. in Finland the history of the external radiation monitoring network goes back to the 1960s3. Traditionally, dose rate detectors are used, but current state-of-the-art early warning networks are equipped with spectrometric stations to identify radionuclides. An effort to harmonize European early warning dosimetry networks by presenting recommendations and requirements for both dose rate detectors and spectrometers is presented in1. Requirements identified include sensitivity to changes in environmental dose rate, dose rate measurement range and readout frequency.

Present in-situ measurement stations have the inherent restriction of not being able to distinguish between airborne radioactivity and ground-deposited fallout. Operational intervention levels (OIL) are often expressed in terms of external dose rates, fallout radioactivity concentration or airborne radioactivity concentration4,5. E.g. in the report Protective Measures in Early and Intermediate Phases of a Nuclear or Radiological Emergency4, jointly published by the Nordic radiation protection and nuclear safety authorities, it is suggested that the intervention levels for adults taking iodine tablets are an external dose rate exceeding 100 µSv/h or an iodine concentration in inhaled air exceeding 10,000 Bq/m3 for two days. With present early warning networks using dose rate detectors or spectrometers, it is only possible to directly determine whether the external dose rate limit is exceeded.

The amount of radioactive nuclides in the air and on the ground can be determined by several methods. The most accurate approach is to collect soil and air samples that are then measured in a laboratory. Since this approach is both labor-intensive and slow, its usage is limited, especially in the early phase of a nuclear accident when results are required rapidly. Another approach is to perform in-situ measurements with collimated gamma-ray detectors that mainly record either the gamma-rays originating from the ground or the gamma-rays originating from the air. In earlier research this has been done by manually repeating the measurement with and without a gamma-ray shield6,7. Alternatively, the measurement station has two different gamma-ray spectrometers with different collimators8,9. A disadvantage of these solutions for airborne and fallout measurements is that they either require some manual work or double the number of gamma-ray detectors per measurement station.

This paper presents a prototype of a novel directional sensing gamma radiation detector utilizing a collimated phoswich scintillator. A phoswich scintillator is a combination of several scintillators. Different applications of phoswiches have been studied, such as atmospheric radioxenon monitoring10, hard X-ray astronomy11 and low-energy gamma-ray detection12. The in-situ method presented here does not require actively changing the measurement geometry nor performing multiple measurements. Because of this it could automatically perform some of the tasks of emergency responders and thus lessen their workload and accumulated radiation dose.

Figure 1 shows the prototype detector. The detector contains a phoswich scintillator package custom-made by Scionix Holland B.V. The package consists of vertically stacked thallium-doped sodium iodide NaI(Tl) and sodium-doped cesium iodide CsI(Na) scintillation crystals, each of 38 mm diameter and 25 mm length, optically coupled to a Hamamatsu R6231 photomultiplier tube (PMT). The diameter of the PMT is 51 mm and it is surrounded by a solid mu-metal magnetic shielding. The phoswich package is hermetically sealed with a light-tight 0.5 mm thick aluminum housing.

Cross-sectional view of the detector. 1—NaI(Tl) crystal, 2—Lead collimator, 3—CsI(Na) crystal, 4—Ø51 mm Hamamatsu R6231 PMT, 5—Magnetic shield, 6—Bridgeport Instruments usbBase, 7—USB port.

The choice of a NaI(Tl)/CsI(Na) phoswich was made after extensive laboratory testing of different combinations of scintillator crystals. This combination was selected for three reasons. Firstly, the light outputs of these scintillator materials feature distinct decay times so that the signals are easily discriminated. Secondly, these scintillator materials are transparent to each other’s scintillation light, enabling good light collection efficiency which is essential for good energy resolution. Thirdly, in NaI(Tl) and CsI(Na), the energy resolution is not highly degraded by the small losses in light collection efficiency, which are unavoidable in a phoswich design 13,14. The NaI(Tl)/CsI(Na) phoswich is also relatively inexpensive.

The detector is instrumented with a usbBase multichannel analyzer (MCA) made by Bridgeport Instruments that is plugged onto the 14-pin PMT socket. In addition to the MCA, the usbBase incorporates a high-voltage supply including the voltage divider. The MCA has a 12-bit ADC that operates with a speed of 40 megasamples per second. The usbBase is powered and controlled over USB.

The phoswich crystals are surrounded by a lead collimator that shields the top scintillation crystal from gamma radiation coming from the ground and vice versa shields the bottom scintillation crystal from gamma radiation coming from above the detector. The collimator has a cross-sectional shape of an isosceles trapezoid with a shorter base closer to the scintillators. The short and long base are 10 and 50 mm respectively. The inner and outer diameters are 43 and 100 mm respectively. The collimator design was optimized by Monte Carlo simulations made with the GEANT4 toolkit15.

The signals from the two scintillator materials can be separated by using pulse shape discrimination (PSD). Typically, the NaI(Tl) scintillator light pulse has a primary decay time of 0.23 µs16 at room temperature and CsI(Na) has a primary decay time of 0.63 µs17. Due to the very different decay times, the signals from these two scintillator materials are easy to separate.

To enable PSD, a custom firmware was loaded onto the usbBase MCA. With this firmware, the MCA records the scintillator events in list mode. Each entry in the list contains the time stamp, pulse integration over a long time window and pulse integration over a short time window for one scintillator event.

The long integration time was used to determine the energy of the gamma ray, and therefore, it was set long enough to capture most of the scintillation light emitted. After careful optimisation, an integration time of 4 μs was selected. Increasing the integration time further only slightly improved the CsI(Na) energy resolution and had no impact on the NaI(Tl) energy resolution. The slight improvement was considered not to justify the increased detector dead time that the longer integration time causes.

The short integration time was optimised for PSD. PSD performs best when the short integration time maximises the difference in integrated signals between the two scintillator materials. After testing various parameter values, the short integration time was set to 0.35 μs.

The parameter used for pulse shape discrimination was defined as

where \(I_{short}\) and \(I_{long}\) are the pulse integration over the short and long time windows, respectively. The offset parameter \(C\) does not have a clear physical meaning, but it was empirically noticed that if a constant value of 13 is subtracted from the pulse integration, the PSD parameter obtained does not depend on the pulse energy. Commonly used PSD models PSD = Ishort/Ilong and PSD = 1 – Ishort/Ilong were also tested but eventually not used, since the PSD parameters were energy depended18,19.

Figure 2 presents an example of a PSD plot for a measurement with a Cs-137 source that emits 662 keV gamma rays, placed on the centre axis above the NaI(Tl) scintillator (angle 0°; see Fig. 4 for angle). The two clusters at channels 580 and 660 are caused by full 662 keV deposition in the CsI and NaI scintillator, respectively (the detector was calibrated so that channel 662 corresponded to 662 keV for the NaI scintillator). The events left from the clusters are Compton scattered gamma rays with partial energy deposition in one detector. The diagonal line connecting the clusters is caused by complete 662 keV absorption shared between the two scintillator materials. Since the detector recorded only 1420 counts per second, random coincidences are hardly visible.

2D histogram of PSD versus energy of the gamma ray events recorded from a Cs-137 source at angle 0° (see Fig. 4 for angle). The horizontal axis represents the parameter used for energy reconstruction (signal amplitude integrated over 4 μs) and the vertical axis the parameter used for pulse-shape discrimination (calculated with Eq. 1). The blue horizontal line shows the minimum PSD value for NaI(Tl) events and the orange horizontal line the maximum PSD value for CsI(Na) events.

Events with the PSD parameter larger than 0.5 were considered to originate from the NaI(Tl) scintillator whereas events with the PSD parameter smaller than 0.3 were considered to originate from the CsI(Na) scintillator. The energy spectra obtained by applying pulse-shape cuts to the same data are presented in Fig. 3.

Undiscriminated (tot) and discriminated energy spectra of a Cs-137 source measured with the phoswich detector. The spectra are derived from the data presented in Fig. 2. The full absorption of 662 keV gamma rays from the Cs-137 source creates a peak at channel 660 in the NaI(Tl) scintillator spectrum and at channel 580 in the CsI(Na) scintillator. The different position of the peaks is mainly due to the difference in the light output of the two materials.

By PSD and gamma spectrum analysis, peak count rates can be extracted from the measurement data. Calibration factors that convert peak count rates to radioactivity concentrations need to be determined for specific fallout or airborne radioactivity geometries. Although the phoswich scintillators have a primary measurement direction and object (air and ground respectively), the collimator and the scintillators do not completely shield either scintillator from the opposing direction. The photopeak count rate S of the upper and lower scintillator in the phoswich can thus be described as

where εV is a calibration factor for the radioactivity in the air for a given cloud geometry and εA is the calibration factor for the radioactivity deposited on the ground. The indexes u and l represent the upper and lower scintillator. Ac represents the radioactivity concentration in the air (Bq/m3), and Ad represents the deposition density on the ground (Bq/m2). In the simplest case the concentrations are constant over the specified volume and area. In a more general case, the radioactivity concentrations Ac and Ad are linear coefficients of concentration distributions for which the corresponding calibration factors have been determined.

From Eqs. (2) and (3) we get

The calibration factors can be calculated at a gamma energy E for specific airborne and fallout geometries. Next, we calculate the calibration factors for an airborne activity in the shape of an infinite half sphere and an infinite fallout surface.

We assume that the detector efficiency does not depend on the azimuth angle due to symmetry. Let’s further assume that we know the absolute peak efficiency εp(E, r0, θ) for a point source at a source-detector distance r0 as a function of the polar angle θ. This absolute peak efficiency can be determined experimentally as

where A is the activity of the point source, I(E) is the yield of photons with energy E and S(E, r0, θ) is the peak count rate of the scintillator.

For a detector at height h from the ground, the calibration factor can be calculated as a sum of two sets of integrals: one for the volume above the center of the detector and one for the volume below the center of the detector. This gives us

where µ(E) is the attenuation constant in air and θ is the azimuthal angle.

The integral over θ can be estimated by measuring the absolute peak efficiency at distance r0 at m different polar angles θi from 0 to π (i.e. straight down to straight up) and interpolating εp(E, r0, θ) linearly between the measured values at angles θi. With θ1 = 0, θn+1 = π/2 and θm = π we get

The remaining integral can be solved numerically.

The fallout calibration factor is calculated similarly. We assume that the detector is at height h from a flat circular surface

and k is the distance from the center of the circle. By integrating by parts and letting the surface be infinite, we get

where E1 refers to the exponential integral:

The angular response of the detector was measured with a 185 MBq Cs-137 source at a distance of 4 m. Cs-137 emits 662 keV photons with a yield of 0.85120 and has a significant role in releases during nuclear power plant accidents21,22. Measurements were performed at 15° intervals around the centre of rotation, i.e. the intersection point between the NaI(Tl) and the CsI(Na) crystals. Figure 4 presents the measured angular response of the 662 keV peak area. For operational use, a peak efficiency calibration should be done over the full energy region of interest.

Polar plot of the phoswich response, measured with a 185 MBq Cs-137 source at 4 m distance and at different polar angles (θ). A cross section view of the detector is depicted in the background. The peak area values represent the number of counts in the 662 keV full energy peak obtained in 1-min measurements.

The two scintillators have a clearly different angular response, especially in the downward direction. The smaller peak areas at angles 150°–210° are explained by the electronics below the scintillators (see Fig. 1).

Using the measured photopeak areas, Eqs. (6), (8) and (11), and a detector position of 1.5 m above ground level, we get the calibration factors presented in Table 1.

For example, if we measure Su = 100 s−1 and Sl = 200 s−1, using the values in Table 1 and Eqs. (4) and (5) we can calculate that the airborne radioactivity concentration is 12 kBq/m3 and the concentration deposited on the ground is 500 kBq/m2 (this example assumes a 100% photon yield).

According to IAEA safety standards, operational intervention levels (OILs) for initiating different parts of emergency plans should be established5. The Nordic radiation protection and nuclear safety authorities have published joint guidelines and recommendations for protective measures in emergencies4. In STUK et al.4, the OIL for sheltering the population indoors due to airborne intense gamma and beta emitters such as Cs-137 is 10 kBq/m3. Similarly, STUK et al.4 determines the OIL for continuing sheltering indoors because of deposited intense gamma and beta emitters as 10 MBq/m2. This level of contamination is defined as “extreme” in4. According to4 with Cs-137, this fallout activity concentration corresponds to an external dose rate of 25 µSv/h.

A detector should be optimized for its planned use. There are several factors limiting the use of any spectrometer in a given situation, such as deadtime and energy resolution. These factors have to be taken into account also in the design of, and operational plans for, the phoswich detector.

The deadtime of the present detector without the collimator was measured at an angle of 90° (see Fig. 4) with different dose rates caused by a Co-60 source. With an ambient dose equivalent rate of 22 µSv/h, the dead time was around 20%. The collimator and energy spectrum have an effect on the deadtime, but the measurement result hints at a detector limitation because of unreasonable deadtime close to the OIL fallout activity defined in4. There is a trade-off between detection efficiency and deadtime at high dose rates.

The limitations of discriminating fallout from airborne radioactivity with the phoswich detector can be studied. Figure 5 presents the effect of increasing fallout concentrations on the uncertainty of the calculated airborne activity concentration. With a 10 MBq/m2 Cs-137 fallout and 10 kBq/m3 Cs-137 air concentration, the relative uncertainty of the estimated air concentration is around 18%. Because of stronger attenuation, the discrimination capabilities should be better for lower gamma energies.

Combined standard uncertainty of the calculated airborne Cs-137 activity vs. airborne activity, for different fallout Cs-137 concentrations. The uncertainties are calculated according to 23, assuming a 10-min measurement, using the calibration factors and respective uncertainties in Table 1, and using the square-root of the peak area as an estimate for the peak area uncertainty.

One factor limiting the use of a stationary in-situ gamma spectrometer in a fallout situation is contamination of the detector surface. This is an issue also for the detector presented here and limits its discriminating capabilities. The effect of contamination can be reduced by covering the in-situ detector with a plastic cover that can be changed. The possible contamination could also be measured and discriminated from the fallout and airborne radioactivity by attaching a third scintillator to the phoswich package24.

It is worth noting that the uncertainty related to the radioactivity distribution introduces significant uncertainties in the calibration factors. This is also a general problem of in-situ measurements of radioactivity in the air or on the ground. In reference25, methods for assessing the representativeness of a measurement site are presented. In reference7, calibration factors correcting for the differences between a semi-infinite cloud and slabs of finite thickness are calculated. Reference7 states that the semi-infinite approximation introduces an error of 20% compared to a 180 m thick slab source. Similar geometry calculations could be done to study the effect on the discriminating capabilities of the phoswich detector. The calibration factors of the detector can be calculated for different geometries to improve the accuracy and precision of the activity estimates and thus improve the response in emergency management.

The phoswich detector presented here demonstrates that it is possible to automatically discriminate fallout and airborne radioactivity using only one detector and one in-situ measurement geometry. If needed, the design and performance can be further optimized to specific operational tasks of the early warning network. Using the phoswich detector, the network can automatically provide information on the physical quantities used in the operational intervention levels.

The datasets used and/or analyzed during the current study can be made available under a collaboration agreement with The Radiation and Nuclear Safety Authority of Finland. Please contact the corresponding author for more information.

Dombrowski, H. et al. Recommendations to harmonize European early warning dosimetry networks. J. Instrum. 12, P12024 (2017).

UNSCEAR, Report to the General Assembly: Sources and effects of ionizing radiation. Volume II: Effects, Annex J. (United Nations, 2001).

Lahtinen, J. & Koivukoski, J. External dose rate monitoring in Finland: History, experiences and a glimpse at the future. Radiat. Prot. Dosim. 187(2), 249–261 (2019).

STUK et al. Protective measures in early and intermediate phases of a nuclear or radiological emergency (2014).

International Atomic Energy Agency, Preparedness and Response for a Nuclear or Radiological Emergency. IAEA Safety Standards Series No. GSR Part 7 (2015).

Honkamaa, T., Toivonen, H. & Nikkinen, M. Monitoring of airborne contamination using mobile equipment. STUK-A130, Finnish centre for radiation and nuclear safety (1996).

Sowa, W. Direct measurement of homogeneously distributed radioactive air contamination with germanium detectors. Radiat. Prot. Dosim. 32(3), 171–176 (1990).

Article  MathSciNet  CAS  Google Scholar 

Hatakka, J., Paatero, J., Viisanen, Y. & Mattson, R. Variations of external radiation due to meteorological and hydrological factors in Central Finland. Radiochemistry 40(6), 534–538 (1998).

Casanovas, R., Morant, J. J. & Salvado, M. Development and calibration of a real-time airborne radioactivity monitor using direct gamma-ray spectrometry with two scintillation detectors. Appl. Radiat. Isot. 89, 102–108 (2014).

Article  CAS  PubMed  Google Scholar 

Luo, F. et al. A well-type NaI(Tl)/BC404 phoswich detector for atmospheric radioxenon monitoring. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 986, 164715 (2021).

Li, H., Ji, J., Feng, H., Zhang, Z. & Han, D. A large area LaBr 3/NaI phoswich for hard X-ray astronomy. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 671, 24–28 (2012).

Article  CAS  ADS  Google Scholar 

Yamasoto, K. et al. CsI(Tl)/plastic phoswich detector enhanced in low-energy gamma-ray detection. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 550(3), 609–615 (2005).

Article  CAS  ADS  Google Scholar 

Mengesha, W., Taulbee, T. D., Rooney, B. D. & Valentine, J. D. Light yield nonproportionality of CsI(Tl), CsI(Na), and YAP. IEEE Trans. Nucl. Sci. 45(3), 456–461 (1998).

Article  CAS  ADS  Google Scholar 

Valentine, J. D., Rooney, B. D. & Li, J. The light yield nonproportionality component of scintillator energy resolution. IEEE Trans. Nucl. Sci. 45(3), 512–517. https://doi.org/10.1109/23.682438 (1998).

Article  CAS  ADS  Google Scholar 

Agostinelli, S. et al. Geant4: A simulation toolkit. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 506(3), 250–303 (2003).

Article  CAS  ADS  Google Scholar 

Moszyński, M. et al. Temperature dependences of LaBr 3(Ce), LaCl3(Ce) and NaI(Tl) scintillators. Nucl. Instrum. Methods Phys. Res. A Accel. Spectrom. Detect. Assoc. Equip. 568, 739–751 (2006).

Gamage, K. A. A., Joyce, M. J. & Hawkes, N. P. A comparison of four different digital algorithms for pulse-shape discrimination in fast scintillators. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 642(1), 78–83 (2011).

Article  CAS  ADS  Google Scholar 

Min, S. et al. Phoswich detectors in sensing applications. Sensors 21, 4047 (2021).

Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

CEA/LNE-LNHB. Nucléide - Lara, Library for gamma and alpha emissions. http://www.lnhb.fr/Laraweb/ (2023).

Voutilainen, M. et al. Potential consequences of hypothetical nuclear power plant accidents in Finland. STUK (2022).

Evangeliou, N., Balkanski, Y., Cozic, A. & Møller, A. P. Global and local cancer risks after the Fukushima Nuclear Power Plant accident as seen from Chernobyl: A modeling study for radiocaesium (134Cs & 137Cs). Environ. Int. 64, 17–27 (2014).

Article  CAS  PubMed  Google Scholar 

Joint Committee for Guides in Metrology (JCGM) Evaluation of measurement data—Guide to the expression of uncertainty in measurement. https://www.bipm.org/documents/20126/2071204/JCGM_100_2008_E.pdf/cb0ef43f-baa5-11cf-3f85-4dcd86f77bd6 (2008).

Bogdanoff, V. Detector contamination monitoring system for Finnish radiation early warning network. Retrieved from http://urn.fi/URN:NBN:fi:jyu-202011106585 (2020).

Zähringer, M. & Sempau, J. The assessment of the representativeness of data from dose rate monitoring stations. Radiat. Prot. Dosim. 64(4), 275–282 (1996).

This work has been partly supported by the Spanish Funding Agency for Research (AEI) through Projects No. PID2019-104390GB-I00. The authors are also grateful for STUK and HIP for allocating resources for this project.

Present address: NeutronGate Oy, 11710, Riihimäki, Finland

Present address: Department of Physics, University of Jyvaskyla, P.O. Box 35, 40014, Jyvaskyla, Finland

Radiation and Nuclear Safety Authority, 01370, Vantaa, Finland

Philip Holm, Sakari Ihantola, Ville Bogdanoff, Kari Peräjärvi & Maarit Muikku

Helsinki Institute of Physics, University of Helsinki, P.O. Box 64, 00014, Helsinki, Finland

Institute of Structure of Matter, CSIC, Serrano 113-Bis, 28006, Madrid, Spain

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

P.H. mainly contributed to the development of detailed data analysis process and writing the manuscript. S.I. and V.B. (V.B. under S.I. supervision) did the experimental part of the work. K.P. did the initial project plan and provided expert support and advice during the project. P.D., O.T. and M.M. provided expert support and advice during the project. All authors contributed to the writing of the manuscript.

Correspondence to Kari Peräjärvi.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Holm, P., Ihantola, S., Bogdanoff, V. et al. Novel environmental monitoring detector for discriminating fallout and airborne radioactivity. Sci Rep 13, 22550 (2023). https://doi.org/10.1038/s41598-023-48730-0

DOI: https://doi.org/10.1038/s41598-023-48730-0

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Novel environmental monitoring detector for discriminating fallout and airborne radioactivity | Scientific Reports

Cdwo4 Crystal Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.