Evaluation of 7Be and 133Xe atmospheric radioactivity time series measured at four CTBTO radionuclide stations
Introduction
Radionuclides found in trace amounts in the atmosphere have been widely used to characterise advection and transport phenomena. Although the 7Be and 133Xe radionuclides (half-lives of 53.1 d and 5.2 d, respectively) have different production mechanisms as well as different physical and chemical characteristics, they have been frequently applied in atmospheric studies. The 7Be is a cosmogenic radionuclide produced in the lower stratosphere and the upper troposphere during nuclear spallation reactions of cosmic-ray particles with oxygen and nitrogen nuclei (Yoshimori, 2005). Its half-life is long enough to be advected for large distances and be detected in ground-level air with activity concentrations depending, among others, on the latitude and the time season (Lee et al., 2015). Even though beryllium is in the atmosphere in a particulate form (attached to aerosols), hence subject to dry and wet deposition, due to its suitable physical characteristics and region of its production, it has been used to study tropopause folding and to distinguish air masses of different origin, including transport of stratospheric air to troposphere (Liu et al., 2016). It may be also coupled with detections of 210Pb which, being an intermediate product of the 238U decay chain, can be linked to air which is being advected from continental regions (Lee et al., 2007). The 133Xe is a noble gas which has been produced during atmospheric and underground nuclear weapons tests, as well as by reprocessing nuclear facilities (RNF), nuclear power plants (NPP) and isotope production facilities (IPF) during their operational releases. Xenon is chemically inert and is not subject to dry or wet deposition. Due to its suitable characteristics it has been chosen as a possible signature of nuclear weapons tests. It is interesting to mention that in addition to 133Xe, three other isotopes of Xenon are used to identify nuclear explosions (131mXe, 133mXe, 135Xe). Xenon is being monitored worldwide by the International Monitoring System (IMS), which has been set up by the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO) in order to monitor over the compliance with the CTBT when it goes into force. The monitoring process is difficult because of low concentrations expected to be released by a low-yield underground nuclear test, the high dilution that characterise atmospheric transport, and also due to the presence of legitimate sources of 133Xe emissions from RNF, NPP and IPF, unevenly distributed in space, and also emitting with both continuous or batch releases, which constitute an undesired background signal.
A comparative study of activity concentrations of 7Be and 133Xe in the ground-level air as obtained by the CTBTO/IMS radionuclide network has been carried out with the aim to contribute to a better localisation of undeclared nuclear activity sites which may occur in the future. In order to characterise the site response of different IMS stations, time series of both 133Xe and 7Be activity concentrations have been analysed, along with meteorological parameters. To achieve that, spectral analysis was first performed to quantify time series harmonic content. Then, correlation analysis between radionuclides and meteorological residuals has been carried out with the aim to establish which environmental conditions characterise radionuclide occurrence at a given station. In particular, detrended crosscorrelation analysis (DCCA) is used to establish the temporal scale where this correlation is significant. The obtained results are representative of local conditions at the measuring site. Characterisation of the complete transport process is needed, including study of emission time series. For the special case of CTBT monitoring though, analysis of really low activity concentrations being involved, establishing if detected radionuclide activity concentrations are more or less correlated to local meteorological conditions could help improve station-specific categorisation schemes (Plastino et al., 2010). The adopted methodology of analysis is explained in detail in Section 2.
Section snippets
Methodology and data analysis
If the time series data have not been excessively affected by missing data, the generalised Lomb Scargle periodogram (Lomb, 1976) can be evaluated in order to extract periodic components and quantify their percentage weight with respect to the total time series. Once peaks corresponding to periodicities have been removed, residual time series can be investigated. This is done evaluating the local Hurst exponent (Ihlen, 2012), related to noise type, and comparing time series outliers
Results and discussion
Results of the comparative analysis of beryllium and xenon time series will be presented and discussed. The harmonic and residual percentages are listed in Table 1, along with the periodicity that contributes the most. Having fully isolated the harmonic components, the residuals time series can be analysed. The local Hurst exponent obtained for 7Be and 133Xe residuals time series is shown in Fig. 1.
The Hurst exponent Ht for beryllium shows small fluctuations and is almost always between the
Conclusions
A new methodology of evaluation of 7Be and 133Xe concentrations in the atmosphere based on time series analysis, applied on atmospheric activity time series (measured at four CTBTO/IMS radionuclide stations) and to meteorological parameters have shown great variability in terms of noise structures, harmonic content, cross-correlation patterns and local Hurst exponent behaviour. Due to the very different physical and chemical characteristics of the two radionuclides, along with their different
Acknowledgements
Sampling data have been kindly provided through the virtual Data Exploitation Centre (vDEC), within the framework of an agreement between the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO) and the Department of Mathematics and Physics - Roma Tre University.
All views expressed are under the sole responsibility of the Authors and not of the Provisional Technical Secretariat CTBTO. The authors would also like to thank the anonymous reviewers for their valuable comments and suggestions
References (14)
- et al.
Multifractal detrended fluctuation analysis of nonstationary time series
Phys. A
(2002) - et al.
Time series study of a 17-year record of 7Be and 210Pb fluxes in northern Taiwan using ensemble empirical mode decomposition
J. Environ. Radioact.
(2015) Production and behavior of beryllium 7 radionuclide in the upper atmosphere
Adv. Space Res.
(2005)DCCA cross-correlation coefficient: quantifying level of cross-correlation
Phys. A: Stat. Mech. Appl.
(2011)- et al.
Ten years of development of equipment for measurement of atmospheric radioactive xenon for the verification of the CTBT
Pure Appl. Geophys.
(2010) - et al.
Long-range correlation properties of coding and noncoding DNA sequences: genbank analysis
Phys. Rev. E
(1995) - et al.
Detrended fluctuation analysis: a scale-free view on neuronal oscillations
Front. Physiol.
(2012)
Cited by (14)
Be-7 and Pb-210 in fallout and aerosols in 2000–2016 in central Europe — Deposition velocity and dependence on meteorological parameters
2022, Science of the Total EnvironmentCitation Excerpt :Ultimately, the average production rate of 810 m−2 s−1 determined by Lal and Peters (1962, 1967) was adopted by the United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) as an indicator of 7Be production in the atmosphere (UNSCEAR, 2000). 7Be transport, its dependence on meteorological conditions and the parameters of phenomena determining the isotope's fate in the atmosphere have been the subject of many studies (Feely et al., 1989, Duenas et al., 2001, Lozano et al., 2012, Pham et al., 2011, Tositti et al., 2014, Kusmierczyk-Michulec et al. 2015, Neroda et al., 2016, Bianchi et al., 2018, Zalewska et al., 2021). The 7Be concentrations observed in the ground layer of the atmosphere depend, among other things, on temperature, precipitation, humidity and atmospheric circulation, which vary significantly by latitude and thus depend on the climate zone and season.
Analysis of environmental radionuclides
2020, Handbook of Radioactivity Analysis: Volume 2: Radioanalytical Applicationstvf-EMD based time series analysis of <sup>7</sup>Be sampled at the CTBTO-IMS network
2019, Physica A: Statistical Mechanics and its ApplicationsCitation Excerpt :The methodology adopted in this paper is intended to be complementary to the one described in [14], hereafter referred to as Multistep Time-series Analysis (MTsA), and aims at extending the range of analysis to strongly nonlinear and non-stationary time series. A detailed description of MTsA software and its applications can be found in [14–17]. In this paper, cosmogenic beryllium-7 time series sampled at 28 different sites at ground level are analysed making use of the Empirical Mode Decomposition (EMD).
Analysis of trends, periodicities, and correlations in the beryllium-7 time series in Northern Europe
2019, Applied Radiation and IsotopesCitation Excerpt :However, we also used other two methods with the purpose of removing the influence of periodicities on the trends of the time series: 1) the seasonal Kendall test (Hirsch et al., 1982), which applies the Mann-Kendall trend test separately for each month and then combines the results; and 2) the trend-free pre-whitening procedure (Yue et al., 2002) applied to the seasonally adjusted monthly time series, to remove the influence of the month-to-month correlations in the significance of the trends. The methodology we used in this paper is an extension of the time series analysis carried out in Plastino et al. (2010) and has been successfully applied in investigations of radionuclides’ time series (Bianchi and Plastino, 2018; Bianchi et al., 2018a, 2018b; Longo et al., 2018). Fig. 2 presents the time evolution of the 7Be activity concentrations in the four sites.
Xenon and radon time series analysis: A new methodological approach for characterising the local scale effects at CTBT radionuclide network
2018, Applied Radiation and Isotopes