Next Article in Journal
Acoustic Presence of Dolphins through Whistles Detection in Mediterranean Shallow Waters
Next Article in Special Issue
In-Situ Seawater Gamma Spectrometry with LaBr3 Detector at a Nuclear Power Plant Outlet
Previous Article in Journal
Verification and Validation of CFD Based Form Factors as a Combined CFD/EFD Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Radioactivity Monitoring at North Aegean Sea Integrating In-Situ Sensor in an Ocean Observing Platform

by
Christos Tsabaris
*,
Effrossyni G. Androulakaki
,
Dionysios Ballas
,
Stylianos Alexakis
,
Leonidas Perivoliotis
and
Athanasia Iona
Hellenic Centre for Marine Research, Institute of Oceanography, P.O. BOX 712, Gr-19013 Anavyssos, Greece
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2021, 9(1), 77; https://0-doi-org.brum.beds.ac.uk/10.3390/jmse9010077
Submission received: 11 December 2020 / Revised: 23 December 2020 / Accepted: 7 January 2021 / Published: 13 January 2021
(This article belongs to the Special Issue Application of Coastal/Ocean Sensors and Systems)

Abstract

:
The integration of the radioactivity spectrometer KATERINA II in a fixed station (buoy) of the POSEIDON network at the North Aegean Sea within the framework of MARRE Project is presented. The acquisition period lasted from 20 November 2019 till 22 February 2020. An intense increment of the activity concentration of radon progenies (up to an order of magnitude) was recorded during rainfall. More specifically, the 214Bi activity concentration varied from 0.09 to 0.53 Bq L−1 without rainfall and the 214Pb activity concentration varied from 0.14 to 0.81 Bq L−1. The 214Bi activity concentration during rainfall ranged from 0.4 to 5.4 Bq L−1 and of 214Pb from 0.3 to 5.3 Bq L−1. The minimum detectable activity of the KATERINA II detection system for measuring low level activities of 137Cs is optimized applying background subtraction and the full spectrum analysis technique.

1. Introduction

Radioactivity monitoring in the marine environment is widely used for monitoring natural and artificial contributions, which are attributed to radiological incidences occurring in nuclear and non-nuclear industries close to the coastline. An example of possible artificial radioactivity spread across Europe is the recent fire that took place in the Chernobyl exclusion zone.
The development of systems for the detection of radioactive contamination has been expanded the last 20 years by applying innovative cost-effective tools (hardware and measurement method) in order to efficiently support the responsible authorities and the scientific radioecological communities. Different in situ methodologies and tools regarding radioactivity monitoring have been tested in many marine areas of the world the last years [1,2,3,4,5,6,7,8,9,10,11], depending on the selected area, topography and contamination level. The continuous monitoring is the only recommended method in cases that the activity concentration of radionuclides is above the minimum detectable activity of the in-situ detection system.
Assuming that the background gamma-ray radiation level (which depends mostly on the water salinity) is constant at a specific area at the open sea, the fluctuations of the gross counting rate in seawater are mainly caused by rainfall due to the atmospheric wash out of the various natural radionuclides (in the absence of a radiological incidence). Radon as a natural radioactive inert gas and non-chemically reactive element is studied extensively in the terrestrial and marine environment concerning NORM and TENORM studies [12]. One of its isotopes, 222Rn, although it is not a gamma-emitter, can be detected via gamma-ray spectrometry from its progenies (214Bi and 214Pb). These two radionuclides have been extensively utilized as radiotracers in oceanography applications [13,14], such as water masses mixing, radioprotection purposes in non-nuclear industries (fertilizers, oil exploration, oil and construction industries), pockmarks [15], mud volcanoes [16], and submarine groundwater discharge identification [17,18]. Moreover, the monitoring of natural and artificial radioactivity variations during rainfall at the sea can provide important information in critical hazards (e.g., fires, floods). The radionuclides 214Bi and 214Pb have been used as radiotracers also in different scientific fields. Data regarding the utilization of radiotracers for meteorological studies (e.g., correlation rainfall events with wind speed direction) is scarce [6,19,20], especially in the marine environment, although it provides ideal conditions for such studies, eliminating factors adding complexity (e.g., wind, temperature variations). In these studies, strong correlations of rainfall events with radon progenies activity concentrations have been observed [6,20], depending on the passage and origin of the rainfall (if it comes from the terrestrial or marine environment). Therefore, monitoring data of radionuclide concentrations at the sea is requested as baseline information, so any possible environmental variations could be identified.
In this study, the correlation analysis between radioactivity data and rainfall parameters from neighboring stations was performed. Moreover, the phenomenon of surface salinity variations induced by rainfall dilution [21,22] was observed via radioactivity data using 40K as a tracer. The main goal of this work was to study the atmospheric wash out due to rainfalls, regarding natural (e.g., radon progenies) and artificial radioactivity (e.g., 137Cs), as well as to improve the minimum detectable activity (MDA) of the detection system in order to detect low concentrations of the artificial 137Cs. For this purpose, a new methodology is developed, which relies on spectra subtraction to minimize the measurement background prior to the analysis and utilization of the Full Spectrum Analysis (FSA) technique for the 137Cs concentration values estimation. To this end, monitoring data in the Athos region at the North Aegean Sea were obtained, for studying radioactivity fluctuations during rainfall events, in the period from 20 November 2019 to 22 February 2020.
In the next paragraph, the floating measuring system (buoy), together with the system upgrade are described. The methodology of data analysis, the variability of the activity concentration of natural radionuclides due to rainfall events and their interpretation are also described. Moreover, the results of the developed methodology for optimizing the MDA of 137Cs (artificial radionuclide) for the KATERINA II system [6] are given for long (24 h) and short-term (3 h) periods of acquisition.

2. Materials and Methods

2.1. The Fixed Mooring System (Athos Station)

The POSEIDON system is an observatory research infrastructure that aims to establish a sustainable marine observing network in the Eastern Mediterranean and to provide quality and validated forecasts for the marine environment. The POSEIDON observing component consists of a network of moored buoys, two underwater gliders, a ferrybox system, a cabled observatory, a fleet of Argo floats, an HF radar system, and a repeated cruises-based ecosystem monitoring program. The system has been in operation for more than 20 years and was developed in accordance to the policy frameworks suggested by IOC/GOOS, EuroGOOS, MonGOOS and GEO. Since 2015, it has participated in the Copernicus Marine Environment Monitoring System (CMEMS—marine.copernicus.eu), providing the near real time in situ data and the wave forecasts for the Mediterranean Sea [23,24]. The POSEIDON mooring network consists of six fixed stations deployed in Aegean and Ionian Seas (Pylos station in SE Ionian, Athos station in North Aegean Sea, Mykonos station in Central Aegean, Saronikos station in the gulf that surrounds the Athens metropolitan area and two stations in South Aegean—Cretan Sea). The Athos station has operated since 2000, carrying recording sensors for the surface marine and atmospheric conditions. In the neighborhood of this station, the National Observatory of Athens owns and maintains two meteorological stations in Vatopedi area and Lemnos Island. The marine and terrestrial stations are shown in Figure 1.
A typical fixed mooring consists of steel wire or synthetic rope, which holds the integrated sensors, and is associated with one or more buoys to provide sufficient buoyancy to keep the mooring upright. Each fixed platform has different kinds of sensors at different depths, providing atmospheric and marine data in the water column. The data recorded by the moorings are transmitted in real-time using GPRS and/or satellite telecommunication links. The long-term recordings collected by the fixed stations can provide evidence on the ocean state and ocean variability, allow a better understanding on the several processes that take place in the marine environment, support the calibration and verification of the marine and weather forecasting systems and enhance the maritime safety, and the efficient planning of marine infrastructures. The mooring network is maintained regularly (twice a year) through dedicated cruises realized by the HCMR’s research vessel.
Although the power system of the station provides continuous energy for all sensors, an independent battery pack was installed for the period of measurements using the radioactivity sensor that was operated in stand-alone mode. Using the software of the installed radioactivity sensor, data are automatically saved according to the time lag of the system. A schematic view of the ocean observing system at Athos with the installed radioactivity sensor and the power supply unit are shown in Figure 2.

2.2. Acquisition and Data Analysis

The radioactivity sensor KATERINA II is integrated in the floating observing system of the POSEIDON network, and the crystal was located at a depth of 70 cm below the sea surface. KATERINA II is a compact and autonomous underwater gamma-ray spectrometer capable of monitoring natural and artificial radionuclides. A scintillation crystal of NaI(Tl) (7.62 × 7.62 cm) is used providing spectra of gamma-rays in a wide range of energy (50–2800 keV) based on a multi-channel analyzer with varied spectroscopy. The system resolution is ~6.5% for 137Cs (662 keV) detection. The dead time is less than 0.5% when the system operates in seawater. All parts of the system are tightly assembled inside a cylindrical housing made of acetal, allowing maximum operation depth at around 400 m. The KATERINA II system is designed to be integrated into measuring fixed ocean platforms. The system also has the specifications to operate in stand-alone mode as well as in sequential buffering acquisition mode.
The electronics of the shaping amplifier provide stable voltage output, compensating voltage drifts (which are the typical drawback of scintillation-based gamma-ray spectrometers). The intense peak of 40K at 1461 keV was always observed at the middle of the spectrum. Potential voltage drifts (especially in the air) are compensated via gain-stabilization procedure using the two present peaks at the energies of 50 and 1461 keV [25,26].
An important upgrade of the KATERINA II system with respect to previous versions is the standalone operation of the system providing time series of data without connection to a computer. Appropriate software was developed to perform measurements in the pre-defined time lag for a specific period of time, saving all data in a volatile memory [6]. The measured data were automatically saved mainly every 3 h of time lag in the aforementioned memory as raw measurements (as it was programmed before the acquisition start). Thus, after the end of the acquisition period, the gamma-ray spectra are automatically stored inside the memory and are retrieved every 2 months (using computer connection). The data pass first through a quality control by checking potential spikes and then a second quality control that identifies abnormal trends (with or without rainfalls).
Before the installation of the detection system at the ocean observing system (buoy), the KATERINA II system was calibrated (for energy, energy resolution and efficiency) in the energy range from energy threshold to 2800 keV for quantitative results as described elsewhere [6,25,27]. The first step in the analysis is the energy calibration as given from existing methodology at the laboratory [27,28]. The data analysis was performed by subtracting the background spectrum (the measured spectrum before rain fall) from the measured spectrum during rainfall. Several tests were made by comparing measured spectra acquired during dry periods, in order to specify a global background spectrum. This was possible due to the constant counting rate in dry periods and the stability of the measured spectra (since spectra subtraction is not applicable if voltage drifts take place during measurements). The spectrum analysis was performed using the software package SPECTRW [29].
The activity concentration of radon progenies (214Bi, 214Pb) was calculated in rainfall events, using experimental and theoretical techniques. The SPECTRW software [29] was utilized in the spectra analysis. The activity concentration of radon progenies was deduced taking into account the detection efficiency of the system, the detected counting rate, the acquisition time, and the emission probability of the detected gamma-rays [27,28]. The statistical uncertainty of the simulation runs for the detection efficiency of the system was kept below 4% in the energy regions of the photopeaks, in all cases. Thus, the total uncertainty budget is calculated using the typical propagation law of relative uncertainties taking into account the uncertainties of the system efficiency (4%) and the statistical counting rate uncertainty in the analyzed photopeaks (4–39%). The analysis is performed assuming homogenized enrichment of radon progenies in the seawater and rapid dilution process in the time lag window of the detection system (3 h).

2.3. Full Spectrum Analysis (FSA) Technique

The typical analysis technique is the photopeak analysis. However, there are some critical drawbacks due to the overlapping energy peaks and the introduction of systematic uncertainty related to background subtraction issues. The most important drawback is the inability of the method to identify unexpected radionuclides [30] that may contribute to the measured number of events in the selected energy window. Consequently, the Full Spectrum Analysis (FSA) technique was also implemented in the analysis of those spectra obtained during rainfall events. The FSA technique is based on the reproduction of standard spectra for all radionuclides of interest. Standard spectra are obtained either experimentally by performing calibration measurements [30,31], or theoretically by performing Monte Carlo (MC) simulations [32,33], including the geometry and detailed characteristics of the measurement. In this work, the standard spectra were derived using the MCNP-CP code [34]. The first step for the spectra analysis using the FSA technique, in order to derive the activity concentrations, was the subtraction of the background radiation (normalized spectrum during periods without rainfall events) from the spectra during rainfall to exclude the 40K and cosmic radiation contributions. Then, the 214Bi and 214Pb standard spectra were normalized to the experimental conditions (electronic setup, acquisition time) and were fitted to the experimental spectrum to derive the activity concentrations. Details of the FSA technique and the convergence between experimental and theoretical spectra are given in the literature [35]. Nevertheless, it is important to mention that in the simulations the radionuclides are considered to be homogeneously distributed in a large (with a radius of 50 cm to 200 cm depending on the detected gamma-ray energy) spherical volume surrounding the detector [35,36].

3. Results

3.1. Quality Control of γ-Ray Spectra

In the period from 20 November 2019 to 22 February 2020, gamma-ray spectra were carried out with the subsea detection system KATERINA II at the Athos station of POSEIDON network. The KATERINA II sensor provided a measurement (γ-ray spectrum) every 3 h from 20 November 2019 till 25 January 2020 and every 6 h from 12 February 2020 till 22 February 2020. Prior to the analysis, a quality control was performed to evaluate the collected data and search for possible malfunctions leading to poor quality data.
In this work, the system exhibited an excellent performance, leading to a total of 428 acquired spectra. A satisfactory observation was that the detection system exhibited excellent stability (less than five channels of voltage drifts), assisting effective spectrum analysis. Typical gamma-ray spectra, that were acquired in periods with rainfall (e.g., 12 December 2019) and without rainfall (defined as background spectra), are depicted in Figure 3. More specifically, the spectra obtained during and after rainfall on 12 December 2019 (3 h measurements) are depicted together with the background spectrum (3 h measurement and without rainfall).
The background spectrum was acquired in a period without rainfall close to the aforementioned date (12 December 2019). The data are depicted in arbitrary units of counts per second along with energy (channels) since there was no need for energy calibration exercises due to the stable output signals. The main characteristic of the measured spectra without rainfall is that the gross counting rate of the system varies in the limit of the statistical uncertainty. This value is considered as a background value of the gross gamma-ray intensity of the KATERINA II sensor operated at the specific site of the sea. Each variation that exceeded the statistical uncertainty was considered as a rainfall event, if the salinity (detected from the 40K peak) remained constant.
However, after a rainfall event, several photo-peaks attributed to radon progenies were identified in the spectrum (Figure 3). More specifically, photo-peaks of the gamma-ray emitters of radon (214Bi and 214Pb) are depicted as raw data along with energy (channels). The 214Pb radon progeny is observed at the gamma-ray energies at 241, 295 and 352 keV, while the 214Bi radon progeny is identified mainly at the energies of 609 and 1764 keV (gamma rays of high intensity). The aforementioned gamma-rays are clearly seen in the spectrum after rainfall. At 12 December 2019 (13:00–16:00), the spectra provide similar activity concentration (within uncertainties) of the two radon progenies, while at 12 December 2019 (16:00–19:00), the counting rate of 214Bi between the two spectra is slightly varied. At the same spectrum, the counting rate of 214Pb is reduced during the second period (16:00–19:00). This reduction of the 214Pb activity concentration between the two continuous periods of acquisition may be attributed to its decay, since rainfall intensity is also reduced.
In all spectra, an intense contribution at 511 keV is also observed from the detection system. The observed peak in this energy is attributed to the mechanism of the pair production, which may occur when highly energetic γ-rays travel along the seawater (close to the detector) or from the atmosphere (since the detector is placed relatively close to the seawater surface). A pair production event is observed through the accompanied appearance in the spectrum of the annihilation peak (511 keV). Such events are not present when the detector is deployed in high depths (above 10 m), due to the screening of the energetic γ-rays in the seawater column [6].

3.2. Gross Gamma-Ray Intensity Analysis

The observations of rainfall events (rain in mm normalized to the radioactivity sensor lag time) as given from the POSEIDON model and the two meteorological stations in Lemnos and Vatopedi [37], are depicted in Figure 4a,b, along with the gross counting rate recorded using the KATERINA II sensor. The data are depicted for the whole monitoring period and for the period from 20 November 2019 to 20 December 2019 (see Figure 4). In most cases, a high increase of the gross counting rate is observed during rainfall events.
The gross counting rate of the recorded events ranges between 11.6 events per second and 13.0 events per second in dry periods (without rain) and increases up to four times during and after rainfall events (the maximum value reaches 40 events per second). Earlier measurements [5] in the same area from May 2000 to December 2002 exhibited similar values with the total recorded event rate ranging from 9 to 12 events per second in dry periods.
During the observation period, 15 events with high gross counting rate were observed with the gross counting rate ranging from 17 events per second to 40 events per second detected during the period 20 November 2019 to 20 December 2019. In the monitoring period from 20 November 2019 to 22 February 2020, 79 events of rainfall took place, providing a gross counting rate larger than 13 events per second or cps (s−1).
The rainfall events that provide gross counting rates above the background value are correlated with rain signals from the terrestrial stations [37], although they are distant from the marine station. The rainfall events are also identified by the POSEIDON forecasting model [23,24]. Very few measurements exhibited high counting rate (cps~17s−1) in the buoy station, while the terrestrial neighboring meteorological stations did not record any signal (possibly a local low intensity rainfall took place in the area of the buoy without affecting the Vatopedi and Lemnos stations). An accurate correlation of the gross counting rate with the rainfall events is not applicable due to the lack of experimental rain intensity data from the same position of the KATERINA II sensor (since the rainfall intensity data are taken from POSEIDON model and from neighboring stations in the terrestrial environment).
An interesting result from the data analysis is that after 20 December 2019, the detected events per second (cps), as given from the KATERINA II sensor, are reduced during rainfalls. This could be associated with a change of the trend of the wind direction (preferably south winds were present). Radon as an inert gas is emanated from the Earth’s crust and then transported by wind and air advection to the observing station. Consequently, if the wind direction is not originated from the terrestrial part, the radon concentration is drastically reduced, since radon progenies’ concentrations in the seawater are negligible.

3.3. Radioactivity Analysis

The results of the activity concentration of radon progenies (3h spectra) together with the meteorological data of neighboring stations are given in Table 1. Only values above 0.4 Bq L−1, for which the statistical uncertainty is tolerable, are presented using the commercial software package (SPECTRW). The FSA method was applied in cases where the activity concentration was above 1 Bq L−1, since the statistics below this value in the subtracted spectra were insufficient.
In most cases, the activity concentration of 214Bi was in good agreement (within uncertainties) with the activity concentration of 214Pb. The activity concentration of 214Bi on rainfall events varies from 0.4 to 5.4 Bq L−1, while for 214Pb from 0.3 to 5.3 Bq L−1. The activity concentration of radon progenies (214Bi and 214Pb) is higher during East winds at North Aegean Sea, in which rainfall is transferred from the terrestrial part of North Greece and not from the marine areas. The same behavior was observed also in the past at the same station [5].
The activity concentrations of radon progenies using the SPECTRW software package and those deduced using the FSA method, were found in good agreement (within uncertainties). The reproduced theoretical spectrum (using the FSA technique) together with the experimental spectrum (after the appropriate background subtraction) of a rainfall event (between 21 to 22 November 2019, summation of nine acquired spectra normalized to 3 h), are shown in Figure 5. The experimental spectrum is fitted using standard spectra of 0.3 Bq L−1 activity concentrations for 214Bi and 214Pb. The agreement between the reproduced with the experimental spectrum is satisfactory along with the detection energy taking into account the statistical fluctuations. This agreement implies that the theoretical model describes satisfactorily the measurements as given by the KATERINA II system taking into account all the mechanisms of interaction of gamma rays with water and detector materials. The observed fluctuations may be attributed to small differences between the experimental electronic setup (energy calibration, energy resolution calibration) and the corresponding values introduced in the theoretical calculations. The analysis of other radionuclides, such as 40K and 137Cs, is also performed. As concerns the 40K results, the data exhibit stable activity concentration values during dry periods representing the salinity values of the seawater.
However, the activity concentration of 40K is drastically decreased during and after strong precipitation due to the rapid mixing processes of the rainwater with the seawater. A typical example is given in Figure 6, where the influence of precipitation is clearly observed from the KATERINA II detection system by the decrease of the measured 40K activity concentration (almost 30% compared with the dry periods) during the precipitation event on 22 November 2019. The same behavior was observed during the rainfall event on 12 December 2019. This precipitation had an amount of 9.2 mm, and the rainwater was enriched with radon progenies (214Bi) concentration of (5.2 ± 0.2) Bq L−1. The calculated salinity values using the 40K activity concentrations ranged from 30 to 39 psu, while similar values ranging from 33 to 37 psu were also obtained from the buoy CTD measurements (in 1m depth) in the same station. This was expected as the 40K concentration is very well correlated with surface salinity values in the seawater [13]. The salinity data (both the values and the observed variations) agree with salinity values obtained in the past in the same station [38]. The decrease of the 40K activity concentration could be attributed to the dilution of seawater through mixing with different water masses either due to the large amount of rainwater or a passage of a seawater current mass (e.g., Black Sea Water), which has been observed in the past [38].
As concerns the 137Cs contribution at 662 keV energy, the photopeak was not clearly seen, despite the fact that the proposed method of spectrum analysis provided lower minimum detectable activity (MDA) (for 137Cs) compared to the work that has been previously done [5]. The lack of 137Cs contribution in the obtained spectra could be attributed to the high radon concentrations of rainfall and also due to the 137Cs decay process after 34 years of the Chernobyl accident [39]. The analysis procedure for the reduction of 137Cs MDA along with the results are given in the next paragraph applying both theoretical and experimental techniques.

4. Discussion

4.1. MDA of 137Cs (24 h Acquisition)

The MDA value of the system KATERINA II for 137Cs was calculated using 24 h measurements during rainfall and a background spectrum (absence of precipitation events) according to Equation (1):
MDA ( Bq / m 3 ) = L D ε m × I γ × T
where, LD is the detection limit (in counts), T is the acquisition time (in s), Iγ the emission probability of the γ-ray, and εm is the marine efficiency (in m3). LD is given as a function of B (net counts) from the equation LD = 2.71 + 4.65B1/2, where B is the net counts at the corresponding gamma-ray energy (foreground with rainfall minus background counts without rain). The relative uncertainty of the MDA is given applying the typical propagation law of uncertainties and is calculated according to Equation (2):
( δ MDA MDA ) 2 = ( δ L D L D ) 2 + ( δ ε m ε m ) 2
where, the relative uncertainty of LD is propagated according to Equation (3).
( δ L D L D ) = ( 2 . 325 B 2 . 7 + 4 . 65 B R + Bg B )
where, R and Bg are the measured counts at the gamma-ray energy with and without rainfall, respectively.
Two different time periods (during rainfall) of 12 to 13 December 2019 (total counting rate 19 events per second) and 11 to 12 December 2019 (total counting rate of 15 events per second) were selected to calculate the MDA values, exhibiting the maximum and mean concentrations of radon progenies respectively, along with a background measurement. The spectrum analysis was performed using the SPECTRW software [29], while the MDA calculations were performed according to existing methodologies [28]. The calculation was first performed via peak analysis in the acquired raw data and in the same data after the subtraction of the ambient background radiation, which includes contributions from cosmic radiation, radionuclides present in the air, and the contribution of 40K, which is always prominent in the seawater. The value of the analyzed total counts drops by a factor of three, (from approximately 20,200 to 6600 counts) after the background subtraction for the measurement at 12 December 2019 and 3.5 times (from 13,300 to 3900) for the measurement at 11/12/2019, respectively.
The background subtraction process leads to a reduction of the detection limit in counts (LD) [28] of 1.75 (for the measurement at 12/12/2019) and 1.85 (for the measurement at 11/12/2019), respectively. This drastic reduction of the LD value is responsible for the MDA improvement, as the LD value is directly proportional to the MDA. The statistical uncertainty of the MDA value for 137Cs in the net spectrum is approximately 5% with highly intense rainfalls and 4% with low-intense rainfall events. The MDA statistical uncertainty value is calculated from the propagation of the statistical uncertainties in the efficiency estimation (around 4%) and the analyzed photopeak counts (using Equation (2)). The statistical uncertainty of the detection limit LD (in counts) varies from 1–4% in all the calculations (using Equation (3). The final MDA results are given in Table 2.
The MDA values are drastically reduced with the background subtraction rendering this technique a useful tool, allowing low level measurements in the marine environment using the KATERINA II detection system. The MDA results are in agreement with corresponding values of similar detection systems as reported in literature by other groups [40].

4.2. Theoretical Estimation of 137Cs MDA

The MDA of 137Cs during rainfalls was also theoretically studied applying the FSA technique, using standard spectra for 214Bi, 214Pb, and 137Cs. The reproduced FSA spectrum of the experimental data on a rainfall event during 21 to 22 November 19 (summation of nine acquired spectra normalized to 3 h), is depicted in Figure 7, introducing an activity concentration of 0.3 Bq L−1 for the radon progenies (214Pb and 214Bi). The same FSA spectrum is also depicted, introducing an activity concentration of 0.02 Bq L−1 for-137Cs. The photopeak of 137Cs becomes apparent in the spectrum for this value (0.02 Bq L−1).
The application of the FSA method, in spectra acquired during rainfall after the background subtraction, provides a methodology for low 137Cs concentration levels (~0.02 Bq L−1) detection, enabling the detection system KATERINA II to quantify low activity concentrations when it operates in the marine environment in incidences of rainfall. Due to the low resolution of the system (214Bi contribution in the 137Cs peak detection area) and the high Compton background from the energetic gamma-rays of radon progenies (such as 214Bi), the estimated MDA value drastically increases when radon progenies concentration in the marine environment increases (e.g., in case of radon progenies activity concentration above 6 Bq L−1, the 137Cs MDA becomes 0.1 Bq L−1).
The obtained MDA values for 137Cs using both experimental and theoretical techniques are drastically improved, but they are above the expected activity concentration values for 137Cs (taking into account the decay of 137Cs), according to previous measurements at the same station in similar meteorological conditions [5].

5. Summary and Conclusions

A detailed study was performed applying the KATERINA II detection system to provide quantitative data of various gamma-ray emitters which are present in the seawater before, during and after rainfall. The study of the natural radioactivity variation gives the important result that rainfall events can be detected through radon progenies concentrations, as an increment of 214Pb and 214Bi activity concentrations. The detected variation of the activity concentration of radon progenies in the absence of rainfall events, fluctuates from 0.095 (17%) to 0.53 (15%) Bq L−1 for 214Bi and from 0.14 (30%) to 0.8 (16%) Bq L−1 for 214Pb. The corresponding values in the presence of rainfall events, range from 0.4 (20%) to 5.4 (4%) for 214Bi and from 0.3 (30%) to 5.3 (4%) Bq L−1 for 214Pb. Therefore, an intense increment of the activity concentration of radon progenies (up to an order of magnitude) was recorded during rainfall.
The integration of the radioactivity sensor to the observation system gave new information for the variation of the level of the γ- activity intensity at North Aegean Sea for the period of rainfalls. During the observation period, the γ- ray intensity was not correlated with the rain intensity (as given from experimental and estimated data). However, the correlation analysis was not performed using experimental data from the same point as with the radioactivity sensor location. The estimated data of the rain intensity were taken from the Poseidon model and the experimental data from the neighbouring stations are provided from the Greek network of the National Observatory of Athens [37]. A more detailed study concerning the correlation of the radionuclide activity concentration variation with environmental parameters from the atmosphere as well as from the sea (like rain intensity, salinity, surface wind etc), is a future application.
The activity concentration of 137Cs during rainfall in Athos region, is below the MDA of the detection system, and thus, is not detected in the defined time lag of the measurement. However, in marine areas that are strongly affected by nuclear accidents (e.g., Chernobyl, Fukushima), the activity concentration of 137Cs could be detected in short acquisition periods (less than 3 h). The proposed method could also be applied in the future in areas where previous nuclear accidents took place, in order to detect variations of key radiotracers at neighboring coastal zones, as well as near aquatic ecosystems and waterways that interact with the terrestrial environment (such as streams, submarine groundwater discharge, and estuaries).
It was found that the detected activity concentration of 137Cs, after intense rainfall events at North Aegean Sea due to the atmospheric wash out process, was reduced by a factor of 3 compared to measurements in 2000 [5], where a similar detection system was integrated in the same buoy station of the POSEIDON network. Furthermore, the activity concentrations of radon progenies were found in radioactive equilibrium in most cases after such rainfall events, exhibiting high concentration values (especially if the origin and direction of the respective cloud were from a terrestrial environment). The intense rainfall events enrich the seawater mass with high concentrations of radon progenies (as tracer of rainwater mass) and lower the values of 40K (as tracer of the seawater) compared with the typical values of 40K in the sea. Thus, the dilution process of rainwater in the sea surface layer should also be studied in the future, combining the aforementioned tracers (radon progenies and 40K) which are given directly quantitatively by the KATERINA II underwater detection system.

Author Contributions

Data curation, A.I.; Formal analysis, C.T. and E.G.A.; Funding acquisition, A.I.; Investigation, C.T., E.G.A., D.B., S.A. and L.P.; Methodology, C.T., E.G.A., D.B. and S.A.; Resources, A.I.; Software, E.G.A. and S.A.; Supervision, C.T.; Validation, L.P. and A.I.; Visualization, A.I.; Writing—original draft, C.T.; Writing—review & editing, C.T., E.G.A. and L.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the MARRE project through National Strategic Reference Framework (NSRF) 2014–2020 co-financed by Greece and the European Union (European Social Fund ESF).

Institutional Review Board Statement

The study was conducted according to the MARRE project and approved for submission by the Institutional Director of the Institute of Oceanography (1/9/2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data are available on request. The data presented in this study are available on request from the corresponding author. The data are not publicly available since MARRE is an on-going project.

Acknowledgments

This work was supported by the MARRE project through National Strategic Reference Framework (NSRF) 2014–2020 co-financed by Greece and the European Union (European Social Fund ESF). The authors would like to acknowledge the POSEIDON group and especially Nikolaos Zisis for the continuous support and for the integration exercises of the sensor to the monitoring network system. The authors would like to thank the divers Vassilis Stasinos and Vassilis Mpampas for the short term cruise to access the observing station at Athos. The authors would also want to acknowledge the crew of the research vessels AEGEAO for installation exercises as well as for the system maintenance during the cruise in the frame of MARRE project (February 2020). The authors would like also to acknowledge Spiridon Velanas for designing the graph of the observing system. At last, the authors would like to thank the group at National Observation of Athens for providing meteorological data of the neighboring stations for the rainfall events during the monitoring period.

Conflicts of Interest

There is no conflict of interest.

References

  1. Aakens, U.R. Radioactivity monitored from moored Oceanographic buoys. Chem. Ecol. 1995, 10, 61–69. [Google Scholar] [CrossRef]
  2. Wedekind, C.H.; Schilling, G.; Grüttmüller, M.; Becker, K. Gamma-radiation monitoring network at sea. Appl. Radiat. Isot. 1999, 50, 733–741. [Google Scholar] [CrossRef]
  3. Povinec, P.P.; Osvath, I.; Baxter, M.S. Underwater Gamma-spectrometry with HPGe and NaI(Tl) detectors. Appl. Radiat. Isot. 1996, 47, 1127–1133. [Google Scholar] [CrossRef]
  4. Klusoň, J. Environmental monitoring and in situ gamma spectrometry. Radiat. Phys. Chem. 2001, 61, 209–216. [Google Scholar] [CrossRef]
  5. Tsabaris, C.; Ballas, D. On line gamma-ray spectrometry. Appl. Radiat. Isot. 2005, 62, 82–89. [Google Scholar] [CrossRef]
  6. Pensieri, S.; Patiris, D.; Alexakis, S.; Anagnostou, M.N.; Prospathopoulos, A.; Tsabaris, C.; Bozzano, R. Integration of underwater radioactivity and acoustic sensors into an open sea near real-time multi-parametric observation system. Sensors 2018, 18, 2737. [Google Scholar] [CrossRef] [Green Version]
  7. Byun, J.-I.; Choi, S.-W.; Song, M.-H.; Chang, B.-U.; Kim, Y.-J.; Yun, J.-Y. A large buoy-based radioactivity monitoring system for gamma-ray emitters in surface seawater. Appl. Radiat. Isot. 2020, 162, 109172. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Wu, B.; Liu, D.; Zhang, Y.; Cheng, Y. Development and deployment of an autonomous sensor for the in-situ radioactivity measurement in the marine environment. Appl. Radiat. Isot. 2018, 142, 181–186. [Google Scholar] [CrossRef]
  9. Caffrey, J.A.; Higley, K.A.; Farsoni, A.T.; Smith, S.; Menn, S. Development and deployment of an underway radioactive cesium monitor off the Japanese coast near Fukushima Dai-ichi. J. Environ. Radioact. 2012, 111, 120–125. [Google Scholar] [CrossRef] [Green Version]
  10. Sartini, L.; Simeone, F.; Pani, P.; Bue, N.L.; Marinaro, G.; Grubich, A.; Lobko, A.; Etiope, G.; Capone, A.; Favali, P.; et al. GEMS: Underwater spectrometer for long-term radioactivity measurements. Nucl. Instrum. Methods A 2011, 626–627, S145–S147. [Google Scholar] [CrossRef]
  11. Thornton, B.; Ohnishi, S.; Ura, T.; Odano, N.; Fujita, T. Continuous measurement of radionuclide distribution off Fukushima using a towed sea-bed gamma ray spectrometer. Deep-Sea Res. 2013, 79, 10–19. [Google Scholar] [CrossRef]
  12. Povinec, P.P.; Eriksson, M.; Scholten, J.; Betti, M. Marine Radioactivity Analysis. In Handbook of Radioactivity Analysis, 3rd ed.; L’Annunziata, M.F., Ed.; Academic Press: Amsterdam, The Netherlands, 2012; Chapter 12; pp. 769–832. [Google Scholar]
  13. Uddin, S.; Behbehani, M.; Aba, A.; Al-Ghadban, A.N. Naturally Occurring Radioactive Material (NORM) in seawater of the northern Arabian Gulf. Mar. Pollut. Bull. 2017, 123, 365–372. [Google Scholar] [CrossRef] [PubMed]
  14. Uddin, S.; Fowler, S.W.; Behbehani, M.; Al-Ghadban, A.N.; Swarzenski, P.W.; Al-Awadhi, N. A review of radioactivity in the Gulf region. Mar. Pollut. Bull. 2020, 159, 111481. [Google Scholar] [CrossRef] [PubMed]
  15. Patiris, D.L.; Tsabaris, C.; Schmidt, M.; Karageorgis, A.P.; Prospathopoulos, A.M.; Alexakis, S.; Linke, P. Mobile underwater in situ gamma-ray spectroscopy to localize groundwater emanation from pockmarks in the Eckernförde bay, Germany. Appl. Radiat. Isot. 2018, 140, 305–313. [Google Scholar] [CrossRef]
  16. Tsabaris, C.; Patiris, D.L.; Pappa, F.K.; Alexakis, S.; Michalopoulos, P. Preliminary investigation of olimpi field, Mediterranean Sea, using in-situ and laboratory radio-tracing methods. Deep-Sea Res. Part II Top. Stud. Oceanogr. 2020, 171, 104689. [Google Scholar] [CrossRef]
  17. Tsabaris, C.; Scholten, J.; Karageorgis, A.P.; Comanducci, J.-F.; Georgopoulos, D.; Liong Wee Kwong, L.; Patiris, D.L.; Papathanassiou, E. Underwater in situ measurements of radionuclides in selected submarine groundwater springs, Mediterranean Sea. Radiat. Prot. Dosim. 2010, 142, 273–281. [Google Scholar] [CrossRef]
  18. Eleftheriou, G.; Pappa, F.Κ.; Maragos, N.; Tsabaris, C. Continuous monitoring of multiple submarine springs by means of gamma-ray spectrometry. J. Environ. Radioact. 2020, 216, 106180. [Google Scholar] [CrossRef]
  19. Baskaran, M.; Swarzenski, P.W. Seasonal variations on the residence times and partitioning of short-lived radionuclides (234Th, 7Be and 210Pb) and depositional fluxes of 7Be and 210Pb in Tampa Bay, Florida. Mar. Chem. 2007, 104, 27–42. [Google Scholar] [CrossRef]
  20. Bottardi, C.; Albéri, M.; Baldoncini, M.; Chiarelli, E.; Montuschi, M.; Giulia, K.; Raptis, C.; Serafini, A.; Strati, V.; Mantovani, F. Rain rate and radon daughters’ activity. Atmos. Environ. 2020, 238, 117728. [Google Scholar] [CrossRef]
  21. Boutin, J.; Chao, Y.; Asher, W.E.; Delcroix, T.; Drucker, R.; Drushka, K.; Kolodziejczyk, N.; Lee, T.; Reul, N.; Reverdin, G.; et al. Satellite and In Situ Salinity: Understanding Near-Surface Stratification and Subfootprint Variability. Bull. Am. Meteorol. Soc. 2016, 97, 1391–1407. [Google Scholar] [CrossRef] [Green Version]
  22. Drushka, K.; Asher, W.E.; Jessup, A.T.; Thompson, E.J.; Iyer, S.; Clark, D. Capturing Fresh Layers with the Surface Salinity Profiler. Oceanography 2019, 32, 76–85. [Google Scholar] [CrossRef]
  23. Perivoliotis, L.; Petihakis, G.; Korres, M.; Ballas, D.; Frangoulis, C.; Pagonis, P.; Ntoumas, M.; Pettas, M.; Chalkiopoulos, A.; Sotiropoulou, M.; et al. The Poseidon System, An Integrated Observing infrastructure at the Eastern Mediterranean as a contribution to the European Ocean Observing System. Operational Oceanography serving Sustainable Marine Development. In Proceedings of the Eight EuroGOOS International Conference, Bergen, Norway, 3–5 October 2017; Buch, E., Fernández, V., Eparkhina, D., Gorringe, P., Nolan, G., Eds.; EuroGOOS: Brussels, Belgium, 2018; pp. 53–61, ISBN 978-2-9601883-3-2. [Google Scholar]
  24. Petihakis, G.; Perivoliotis, L.; Korres, G.; Ballas, D.; Frangoulis, C.; Pagonis, P.; Ntoumas, M.; Pettas, M.; Chalkiopoulos, A.; Sotiropoulou, M.; et al. An integrated open-coastal biogeochemistry, ecosystem and biodiversity observatory of the eastern Mediterranean-the Cretan Sea component of the POSEIDON system. Ocean Sci. 2018, 14, 1223–1245. [Google Scholar] [CrossRef] [Green Version]
  25. Tsabaris, C.; Vlachos, D.S.; Papadopoulos, C.T.; Vlastou, R.; Kalfas, C.A. Set up and application of an Underwater γ-ray spectrometer for radioactivity Measurements. Mediterr. Mar. Sci. 2005, 6, 35–40. [Google Scholar] [CrossRef] [Green Version]
  26. Vlachos, D.S.; Tsabaris, C. Response Function Calculation of an Underwater Gamma Ray NaI(Tl) Spectrometer. Nucl. Instrum. Methods A 2005, 539, 413–420. [Google Scholar] [CrossRef]
  27. Tsabaris, C.; Bagatelas, C.; Dakladas, T.; Papadopoulos, C.T.; Vlastou, R.; Chronis, G.T. An autonomous in situ detection system for radioactivity measurements in the marine environment. Appl. Rad. Isot. 2008, 66, 1419–1426. [Google Scholar] [CrossRef]
  28. Bagatelas, C.; Tsabaris, C.; Kokkoris, M.; Papadopoulos, C.T.; Vlastou, R. Determination of marine gamma activity and study of the minimum detectable activity (MDA) in 4pi geometry based on Monte Carlo simulation. Environ. Monit. Assess. 2010, 165, 159–168. [Google Scholar] [CrossRef]
  29. Kalfas, C.A.; Axiotis, M.; Tsabaris, C. SPECTRW: A software package for nuclear and atomic spectroscopy. Nucl. Instrum. Methods A 2016, 830, 265–274. [Google Scholar] [CrossRef]
  30. Caciolli, A.; Baldoncini, M.; Bezzon, G.P.; Broggini, C.; Buso, G.P.; Callegari, I.; Colonna, T.; Fiorentini, G.; Guastaldi, E.; Mantovani, F.; et al. A new FSA approach for in situ γ ray spectroscopy. Sci. Total Environ. 2012, 414, 639–645. [Google Scholar] [CrossRef] [Green Version]
  31. Hendriks, P.H.G.M.; Limburg, J.; de Meijer, R.J. Full-spectrum analysis of natural γ-ray spectra. J. Environ. Radioact. 2001, 53, 365–380. [Google Scholar] [CrossRef]
  32. van der Graaf, E.R.; Limburg, J.; Koomans, R.L.; Tijs, M. Monte Carlo based calibration of scintillation detectors for laboratory and in situ gamma ray measurements. J. Environ. Radioact. 2011, 102, 270–282. [Google Scholar] [CrossRef]
  33. Maučec, M.; de Meijer, R.J.; Rigollet, C.; Hendriks, P.H.G.M.; Jones, D.G. Detection of radioactive particles offshore by γ-ray spectrometry Part I: Monte Carlo assessment of detection depth limits. Nucl. Instrum. Methods 2004, 525, 593–609. [Google Scholar] [CrossRef]
  34. Berlizov, A.N. MCNP-CP a Correlated Particle Radiation Source Extension of a General Purpose Monte Carlo N Particle Transport Code; Semkov, T.M., Pommé, S., Jerome, S.M., Eds.; ACS Symposium Series 945; American Chemical Society: Washington, DC, USA, 2006; pp. 183–194. [Google Scholar]
  35. Androulakaki, E.G.; Kokkoris, M.; Tsabaris, C.; Eleftheriou, G.; Patiris, D.L.; Pappa, F.K.; Vlastou, R. In situ γ-ray spectrometry in the marine environment using full spectrum analysis for natural radionuclides. Appl. Radiat. Isot. 2016, 114, 76–86. [Google Scholar] [CrossRef] [PubMed]
  36. Tsabaris, C.; Thanos, I. An underwater sensing system for monitoring radioactivity in the marine environment. Mediterr. Mar. Sci. 2004, 5, 5–12. [Google Scholar] [CrossRef] [Green Version]
  37. Lagouvardos, K.; Kotroni, V.; Bezes, A.; Koletsis, I.; Kopania, T.; Lykoudis, S.; Mazarakis, N.; Papagiannaki, K.; Vougioukas, S. The automatic weather stations NOANN network of the National Observatory of Athens: Operation and database. Geosci. Data J. 2017, 4, 4–16. [Google Scholar] [CrossRef]
  38. Velaoras, D.; Kassis, D.; Perivoliotis, L.; Pagonis, C.; Hondronasios, A.; Nittis, K. Temperature and salinity variability in the Greek Seas based on POSEIDON stations time series: Preliminary results. Mediterr. Mar. Sci. 2013, 14, 5–18. [Google Scholar] [CrossRef] [Green Version]
  39. Vlastou, R.; Ntziou, I.T.; Kokkoris, M.; Papadopoulos, C.T.; Tsabaris, C. Monte Carlo simulation of g-ray spectra from natural radionuclides recorded by a NaI detector in the marine environment. Appl. Radiat. Isot. 2006, 64, 116–123. [Google Scholar] [CrossRef] [PubMed]
  40. Zhang, Y.; Wu, B.; Liu, D.; Lv, H.; Feng, X. Research on minimum detectable activity (MDA) of underwater gamma spectrometer for radioactivity measurement in the marine environment. Appl. Radiat. Isot. 2020, 155, 108917. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The area of study.
Figure 1. The area of study.
Jmse 09 00077 g001
Figure 2. The POSEIDON fixed station at Athos (North Aegean Sea). The sensor is represented with “K” and the power supply unit with “P”.
Figure 2. The POSEIDON fixed station at Athos (North Aegean Sea). The sensor is represented with “K” and the power supply unit with “P”.
Jmse 09 00077 g002
Figure 3. The background spectrum without rainfall and two spectra with rainfall. The acquisition period was 3 h.
Figure 3. The background spectrum without rainfall and two spectra with rainfall. The acquisition period was 3 h.
Jmse 09 00077 g003
Figure 4. The counting rate (as given from the KATERINA II system) and the rain intensity from Poseidon model and Lemnos and Vatopedi stations for two periods: 20/11/2019–22/02/2020 (b) and from 20/11/2019–20/12/2019 (a).
Figure 4. The counting rate (as given from the KATERINA II system) and the rain intensity from Poseidon model and Lemnos and Vatopedi stations for two periods: 20/11/2019–22/02/2020 (b) and from 20/11/2019–20/12/2019 (a).
Jmse 09 00077 g004
Figure 5. The experimental and reproduced spectra (after subtraction of the background spectrum) of the radon progenies (214Pb, 214Bi). The activity concentration for both progenies is 0.3 Bq L−1, and the acquisition time is 3 h.
Figure 5. The experimental and reproduced spectra (after subtraction of the background spectrum) of the radon progenies (214Pb, 214Bi). The activity concentration for both progenies is 0.3 Bq L−1, and the acquisition time is 3 h.
Jmse 09 00077 g005
Figure 6. The activity concentration of 40K and 214Bi for the measuring period (the vertical lines represent the rainfalls).
Figure 6. The activity concentration of 40K and 214Bi for the measuring period (the vertical lines represent the rainfalls).
Jmse 09 00077 g006
Figure 7. The simulated spectrum (at the energy range from threshold to 1000 keV) after the rain on 21 November 2019 (3 h spectrum with activity concentration of radon progenies 0.3 Bq L−1) and the 137Cs contribution of 0.02 Bq L−1.
Figure 7. The simulated spectrum (at the energy range from threshold to 1000 keV) after the rain on 21 November 2019 (3 h spectrum with activity concentration of radon progenies 0.3 Bq L−1) and the 137Cs contribution of 0.02 Bq L−1.
Jmse 09 00077 g007
Table 1. The activity concentration of radon progenies (214Bi and 214Pb) during rainfalls analysed with the SPECTRW software and the FSA method. The time-normalized meteorological data of two neighboring stations (Vatopedi and Lemnos) are also included.
Table 1. The activity concentration of radon progenies (214Bi and 214Pb) during rainfalls analysed with the SPECTRW software and the FSA method. The time-normalized meteorological data of two neighboring stations (Vatopedi and Lemnos) are also included.
SPECTRW AnalysisFSA AnalysisVatopedi StationLemnos Station
DateTime214Bi214Pb214Bi214PbRainWind DirectionRainWind Direction
(Bq l−1)(Bq l−1)(Bq l−1)(Bq l−1)(mm) (mm)
22/11/20191:005.2 (4)5.2 (4)5.24.89.2ESE0.8NE
25/11/201913:002.6 (4)2.6 (4)2.92.222.8ESE8ENE
25/11/201916:002.5 (4)2.5 (4)2.72.328.2E/ESE15.6ENE
26/11/201913:001.4 (4)1.2 (5)1.71.40NNE/ESE0.2SW/NW
27/11/202016:000.7 (17)0.8 (19) 0-0SE
4/12/20191:000.8 (12)0.9 (18) 11.8ESE2.2NE
4/12/20194:001.4 (10)1.3 (21)1.71.45.8SE4.4NNE
4/12/20197:001.2 (9)1.2 (14)1.41.32.4ESE4.4NE
4/12/201910:000.8 (13)0.9 (19) 1.2ESE1.4NNE
4/12/201913:000.7 (12)0.8 (16) 0.2ESE0.4NNE
8/12/20194:000.5 (16)0.5 (24) 0WNW0.4WNW
9/12/20194:000.5 (18)0.8 (16) 0WNW0.2NW
10/12/201922:000.5 (16)0.5 (27) 0ESE0ΝΕ
11/12/20191:000.7 (13)0.3 (38) 0ESE0ΝΕ
11/12/20197:000.6 (29)0.4 (32) 8.4ESE4.2ΝΕ
11/12/201910:001.0 (13)1.1 (15)1.11.210SE5.2ΝΕ
11/12/201913:000.6 (15)0.3 (39) 9.4SE0.2ΝΕ
11/12/201916:000.7 (16)0.4 (35) 3.8SE0ΝΕ
11/12/201919:000.7 (14)0.7 (19) 18.6SE0ΝΕ
11/12/201922:000.7 (12)0.6 (21) 16.8ESE0.8ΝΕ
12/12/20191:001.1 (8)0.6 (30) 8ESE/SE14.2ΝΕ
12/12/20194:000.9 (10)0.9 (16) 39.2SE2.2Ε/ΝΕ
12/12/201910:000.7 (12)0.6 (28) 43.4SSE0ΝΕ/Ε
12/12/201913:005.2 (4)5.3 (4)5.24.86.2-0Ε/ΕΝΕ
12/12/201916:005.4 (4)5.2 (4)5.24.30NW0ΝΕ
12/12/201919:001 (10)1.1 (13) 0.2NW0ΝΕ
12/12/201922:000.6 (16)0.4 (30) 0.2S0NNW
13/12/20191:000.5 (18)0.7 (22) 0S0.2NNW
14/12/201916:000.7 (14)0.8 (19) 14.6ESE/SE0ESE/SE
14/12/201919:001.4 (8)1.9 (10)1.71.93NNW0.2ENE
14/12/201922:001.4 (8)1.8 (10)1.91.80NNW0.4WNW
15/12/20191:000.7 (15)0.5 (28) 0SSW0WNW
20/12/20191:001 (20)1.2 (11) 0-0-
27/1/20209:300.6 (14)0.8 (16) 0WSW0SSW
27/1/202012:300.6 (13)0.7 (17) 0WSW2.8SSE
27/1/202015:300.5 (14)0.5 (23) 0SSW0.4SW
27/1/202018:300.5 (15)0.5 (24) 0SSW0SW
28/1/20203:300.5 (16)0.6 (21) 0WNW0-
30/1/20200:300.5 (17)0.5 (23) 0WNW0WSW
30/1/20203:300.4 (23)0.4 (23) 0WNW0WNW
30/1/20206:300.4 (20)0.5 (23) 0NW0.2WNW
30/1/20209:300.5 (15)0.5 (20) 0NNW1.2NE
Table 2. MDA values for 137Cs with and without background subtraction for two different measurements during precipitation events (24 h period with and without background subtraction).
Table 2. MDA values for 137Cs with and without background subtraction for two different measurements during precipitation events (24 h period with and without background subtraction).
Measurement PeriodMDA for 137Cs (Bq L−1)
12/12/2019
(without background subtraction)
0.041
12/12/2019
(background subtraction)
0.023
11/12/2019
(without background subtraction)
0.033
11/12/2019
(background subtraction)
0.018
Background0.031
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tsabaris, C.; Androulakaki, E.G.; Ballas, D.; Alexakis, S.; Perivoliotis, L.; Iona, A. Radioactivity Monitoring at North Aegean Sea Integrating In-Situ Sensor in an Ocean Observing Platform. J. Mar. Sci. Eng. 2021, 9, 77. https://0-doi-org.brum.beds.ac.uk/10.3390/jmse9010077

AMA Style

Tsabaris C, Androulakaki EG, Ballas D, Alexakis S, Perivoliotis L, Iona A. Radioactivity Monitoring at North Aegean Sea Integrating In-Situ Sensor in an Ocean Observing Platform. Journal of Marine Science and Engineering. 2021; 9(1):77. https://0-doi-org.brum.beds.ac.uk/10.3390/jmse9010077

Chicago/Turabian Style

Tsabaris, Christos, Effrossyni G. Androulakaki, Dionysios Ballas, Stylianos Alexakis, Leonidas Perivoliotis, and Athanasia Iona. 2021. "Radioactivity Monitoring at North Aegean Sea Integrating In-Situ Sensor in an Ocean Observing Platform" Journal of Marine Science and Engineering 9, no. 1: 77. https://0-doi-org.brum.beds.ac.uk/10.3390/jmse9010077

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop