1. Introduction
The protection of first responders and the general public against radioactive contaminations caused by accidental or deliberate radiological events, including terrorist attacks and a possible war scenario, is critical. After a nuclear or radiological incident, the decision-makers must understand the spread of contamination and the affected area. This information is crucial to determine the most suitable response. In the immediate vicinity of a nuclear or radiological event, or in case of large-area ground contaminations, the remote detection of alpha and gamma radioactivity by using Unmanned Aircraft Systems (UAS) is an excellent solution to assess the contaminated area before the first responders can enter to protect them against contamination and irradiation.
Because the radionuclide inventory in such scenarios is typically unknown, there is no single dedicated detector-UAS configuration, and the operator must analyze or predict the contamination profile to select the most suitable detector-UAS configuration for mapping the contaminated area.
In the case of contaminations with gamma-emitting radionuclides, spectrometric gamma detectors that can both identify radionuclides and deduce a dose rate should be mounted on the UAS. Typical dose rate monitors, such as Geiger-Muller (GM) tubes, can also be used for a fast screening to evaluate the dose rate map for further analysis with spectrometric detectors. In the European project “Preparedness”
1 spectrometric gamma detectors were adapted and mounted on selected UAS based on their flight capabilities [
1,
2].
In case of a radiological scenario involving the dispersion of alpha-emitting radionuclides in the environment, the situation is significantly more complicated since alpha particles have a maximum range of only a few centimeters in the air and much shorter (around 100 µm) in biological tissues [
3]. In addition, ingestion or inhalation of alpha-emitting radionuclides leads to a higher associated health risk compared to gamma-emitting sources for the same intake. Alpha particles deposit all their energy in a relatively short path, usually resulting in double-strand deoxyribonucleic acid (DNA) breaks that are difficult to repair unlike single-strand DNA breaks caused by gamma-rays. [
4].
Therefore, in a scenario where alpha-emitting radionuclides are potentially involved, it is very important to know in advance the presence of such radionuclides and provide the first responders with appropriate protection equipment to avoid ingestion and inhalation, or generally, to avoid contact with such radionuclides. However, a detection system to remotely measure alpha particles is not available. Detection of alpha-emitting radionuclides with traditional detectors such as proportional counters, scintillator counters and passivated implanted planar silicon (PIPS) is highly time consuming and tedious due to the short range of alpha particles (e.g., about 3.7 cm for a 5.5 MeV particle) - this approach involves scanning very close to the surface of the contaminated area and would expose the personnel to other hazards and risks such as other types of radiation, fire, etc. Therefore, new detection techniques which keep both the detector and the operator out of the radiation field are needed.
Optical detection of alpha particle emitters in the environment is a technology that enables sensing a radiological threat at a safe distance, without putting personnel at risk or contaminating equipment. This standoff detection technique gives advanced warning of potential threats in contrast to conventional close-range handheld alpha detectors. These optical detectors are based on alpha radiation detection via radioluminescence of air. Alpha particle induces radioluminescence light when absorbed in air. Radioluminescence in air is generated mostly by the emission of molecular nitrogen (N
2), and to a much lesser extent, by trace amounts of nitric oxide (NO) in air, with wavelength spanning three ultraviolet (UV) bands and about 99% of emission occurring in the 280–440 nm spectral range [
5,
6,
7]. The diurnal background light is composed of approximately 95% UV-A (wavelengths from about 315 nm to 400 nm), 5% UV-B (wavelengths from about 280 nm to 315 nm), and almost no UV-C (wavelengths from about 100 nm to 280 nm). Therefore, the optical detection of alpha-emitting radionuclides is complex because the UV-A emissions are eclipsed by the much greater solar irradiance. Due to the low UV-C background emissions, the idea of acquiring only the UV-C (solar blind) light during daytime is a good alternative. However, the intensity of UV-C radioluminescence is very low (< 1%) which makes the detection quite challenging.
The alpha-induced air luminescence has been estimated to generate approximately 19 photons/MeV (at 1013 hPa, 22°C, and 43% relative humidity), and the dependence is found to be linear in the studied energy range from 0.3 MeV to 5.1 MeV [
6]. The emitted photons are mainly in the UV-A wavelengths, and only about 0.2 photons/MeV belong in the UV-C wavelength range [
7]. The radioluminescence signal is influenced by the characteristics of the gas atmosphere mainly due to the quenching effects of water vapor and oxygen. The UV-A photon yield can be increased by using a high-purity nitrogen atmosphere, whilst adding trace amounts of NO to the nitrogen atmosphere causes a sharp increase in the UV-C light yield. In the work of Kerst and Toivonen [
8], the intensity of the UV-C signal was maximized at 50 ppm NO diluted in pure nitrogen. Also, this effect has been illustrated in the work of Krasniqi et al. [
9] where activities as low as 300 Bq were measured in uranium samples using 3 ppm of NO.
The RemoteALPHA European project (2020-2023) [
10] had the overall objective to develop novel optical systems and methodologies for the remote detection and quantification of large-scale contamination with alpha emitters in the outdoor environment, allowing sound and quick countermeasures in the case of a radiological emergency. One of the objectives was focused on the development of an unmanned airborne monitoring system (UAMS) which integrates the UAS with a novel alpha-radioluminescence detection system to scan and obtain a map of the contaminated area in real-time.
In the framework of this project, an optical detection system based on the Poly (methyl methacrylate) (PMMA) Fresnel lens and interference filters for detecting the UV-C spectral region was built and mounted on a DJI Matrice 600 Pro. A 100 MBq 241Am source with an active area of 20 mm × 100 mm, specifically designed for the project, was used for testing the monitoring system in real flight. The source is housed inside an evacuable gas vessel with a quartz dome. This allows the space around the source to be filled with selected gases and the radioluminescence outside the chamber to be studied in the relevant wavelength range above 230 nm. The vessel was purged by high-purity nitrogen doped with 50 ppm NO to increase the UV-C light yield.
This paper focuses on the description of the developed system, its characterization, and the results obtained during real flights and simulated scenarios.
Section 2 details the materials and methods used to perform the UAS missions. This includes the hardware components, air-ground communications, and software architecture.
Section 3 discusses the detector calibration results obtained both in laboratory conditions and during flight campaigns conducted at the Drone Research Laboratory (DroneLab) of the Technical University of Catalonia (UPC) [
11]. Finally,
Section 4 provides the conclusions of this work.
2. Materials and Methods
This Section introduces the developed hardware and software to detect alpha-emitting radionuclides using an upgraded version of RIMA Software Architecture [
1]. The different components of the developed system are presented including the optical detector, the UAS, radiological sources, the laser altimeter, air-ground communications, and the on-board computer. This Section also explains the calibration method used to characterize the light collection efficiency and field of view of the detector in a controlled laboratory environment using the UV sources that simulate the photon yield of alpha sources. At last, the software developed to send and display the alpha-emitting radionuclides in real-time at the UAS Ground Control Station (GCS) is presented.
2.1. Hardware Architecture
This Section describes all the equipment used during the project to perform the alpha particle detection flights explaining the specific role of each instrument or sensor.
2.1.1. Detector
To the authors' knowledge, this is the first time that alpha particles have been detected from a UAS. There are currently no commercial systems available to measure alpha particles remotely from a UAS. Therefore, the initial step was to design and build an optical alpha detection system mountable on a commercial UAS. The main components of the detection system are light collection optics, photon filtering, and photon counting.
The radioluminescence light collection is done by a PMMA Fresnel lens with a diameter of 452.9 mm and a focal length of 424.5 mm chosen to ensure a high light throughput in a manageable form factor [
12]. The photon filtering is done with 3 Semrock FF01-260/16-25 interference bandpass filters centered at 260 nm with a bandwidth of 16 nm. To achieve the best UV-C band filtering, a set of two corrective fused silica lenses (Thorlabs LC4252 and LA4052) was installed to parallelize the collimated photon beam. An additional filtering performance was achieved by tilting two filters in the set by 3 degrees relative to the optical axis.
The photon counting is done with a PMT detector, that converts an incident light into an electrical signal analyzed by a digital pulse processing unit. The detection system uses the solar-blind Hamamatsu H11870-09 photomultiplier tube [
12] with transistor-transistor logic (TTL) output and the GBS Elektronik MCA527 OEM+ multichannel analyzer [
13] working in the scaler (gating) mode to translate these TTL signals into the timestamped count rate. During the flight, the count rate is polled and transmitted to the GCS to be displayed on a map for subsequent interpretation. If the count rate exceeds the threshold determined by the background rate, the source is present in the detector's field of view.
A carbon fiber-reinforced cone-shaped polymer frame has been built to hold the lens and the detection system. The frame materials were chosen to provide the necessary stability for the system while keeping it lightweight. The detector has a height of 790 mm and a diameter of 520 mm at the base.
Figure 1 shows an image of the detector and its different parts in a horizontal position during laboratory tests.
2.1.2. Characterization of the Detector
Investigation of the detector response was first performed under light-controlled conditions with UV LEDs to determine the field of view (FoV) and the best focusing configuration for the flights. In addition, the change in efficiency as a function of the angle of incidence and the power of the LED lamp was investigated. The results of these experiments are presented in
Section 3.
The LED lamp used to test and optimize the optical system had a 5 mm thick Teflon disc at the exit aperture and was calibrated with a
210Po activity standard. The LED lamp allows 12 different power modes, from number 0 to 11 (see
Figure 2), simulating alpha sources with energy 5.4 MeV and activities up to 7.3 GBq (setting Nr. 11).
The optical system has a sliding detector design which allows adjusting the focus according to the object distance. The efficiency of each focusing configuration, represented by detector position F (protrusion from the frame in centimeters) was tested with the LED source placed at the optical axis at various distances as illustrated in
Figure 3.
The field of view is defined as the diameter of the spot that can be seen by the optical detection system. It depends on the focal length of the used receiving optics and the object distance. While in the simple lens-PMT case, the FOV can be estimated with a formula, the use of corrective lenses and non-trivial filter allocation required experimental evaluation. To determine the field of view, the LED source was moved perpendicularly to the optical axis.
The FOV is calculated by fitting the experimental data to a Gaussian function. The amplitude A and the standard deviation σ are calculated. The FOV is defined by the integral area in [-2σ, +2σ] range that implies ca. 95% of the total light collected.
The angular field of view (AFOV) is another important metric for drone-assisted optical measurements. It is calculated from the FOV and the distance between the object and the lens:
To assess how the tilt of the UAS can affect the light collection efficiency, the radio-luminescence signal has been measured as a function of the angle between the Fresnel lens optical axis and UV LED axis using a rotation platform with 1° precision.
This was the method used to calculate the field of view of the detector, the appropriate focus setting, and the count rate values provided by the LED source.
Section 3 details the graphs obtained from these experiments in the laboratory.
2.1.3. Unmanned Aircraft System (UAS)
Given the dimensions of the detector (see
Figure 4), no commercial UAS was tall enough to natively integrate the detector under its belly, since the height of the detector is larger than UAS legs. This made it difficult to integrate the detector into a commercial UAS unless the detector itself acted as a landing gear. Therefore, four legs were attached to the detector so that the structure of the Fresnel lens is used as the landing gear of the unmanned aircraft.
Figure 4 shows the Fresnel lens detector structure with a scheme of the landing legs and their dimensions.
Table 1 displays a summary of the weight of each UAS payload element and the sum of all components.
Two UAS from the UPC fleet were selected as possible candidates to mount the payload described in
Table 1: DJI WIND 4 (W4) and DJI Matrice 600 Pro (M600P).
Figure 5 demonstrates pictures of both UAS. The authors have extensive experience with the M600P flight detector for radiological measurements, as it was used in previous projects such as “Preparedness” [
1,
2]. In addition, the M600P incorporates newer flight technology than the W4. For instance, the M600P is equipped with the A3 Pro flight controller, which improves control accuracy while increasing overall safety. On the other hand, the W4 has a better payload weight and flight autonomy than the M600P. However, the landing legs of W4 suffer from mechanical stress at high payloads. Another drawback of the W4 is the placement of the batteries, which would make integration of the detector into the UAS more difficult. Given the extensive operational experience, the improved avionics of the UAS, and easy payload integration, the M600P was chosen to mount the radioluminescence detector.
Since M600P was not equipped with a Real-time Kinematic (RTK) system that greatly enhances the positioning accuracy, especially in the vertical profile, the Lightware SF11/C laser altimeter [
14] was installed on the UAS. This altimeter has a long range (100 m), fast update rate (20 readings per second), and lightweight (35 g) laser designed specifically to be used on UAS with 0.1 m accuracy.
2.1.4. Onboard Computer and Air-Ground Communications
Each measurement carried out with the optical detector must be geo-tagged during the flight for real-time mapping. Taking this task, the onboard computer fuses real-time UAS telemetry and the detector count rate data. The optical detector can provide a count rate measurement every 100 milliseconds.
The Raspberry Pi 4 Model B was chosen as the onboard computer [
15]. This single-board computer is excellent for prototyping due to its low weight, large number of ports for connecting devices, and large support community. The Raspberry Pi 4 Model B has a Quad-core Cortex-A72 (ARM v8) 64-bit SoC with 1.8 GHz and 8 GB SRAM. The board also features 2.4 GHz and 5 GHz IEEE 802.11ac wireless LAN, Bluetooth 5.0, two USB 3.0 ports, two USB 2.0, Gigabit Ethernet, and an extended 40-pin General-Purpose Input/Output (GPIO) header where six UARTs can be found. Two of the six available ports are used. One of the UARTs is connected to the UAS autopilot to receive its telemetry and extract its current position and UAS attitude. The second UART port is used to connect the laser altimeter to obtain an accurate altitude.
The air-ground communication with the onboard computer is performed through Sixfab Raspberry Pi 4G/LTE Cellular Modem Kit [
16]. The kit provides all the components for establishing an Internet connection on the Raspberry Pi. This allows interaction with the onboard computer at any given time, in real time, and from anywhere in the world, as long as it is connected to the Internet. Through this connection, the geo-tagged measurements taken by the detector are sent to the GCS.
The onboard computer and 4G/LTE air-ground communications are powered using the UAS batteries. The M600P has a connector that provides an 18 V 7 A power output. The voltage must be converted to 5 V to power the onboard computer, for which the HOBBYWING regulator capable of providing 5 V and 3 A was used.
Figure 6 shows the integration of the Raspberry PI, the 4G/LTE model with its corresponding antennas and the voltage regulator. It also shows the connection of the autopilot and the laser altimeter (right panel).
2.1.5. The Hardware Integration in the UAS
The alpha radioluminescence detector with the landing legs, the MCA counter, the laser altimeter, and the onboard computer with a 4G/LTE modem were integrated with the M600P as shown in
Figure 7. The figure also displays the flat landing pad required for the takeoff and landing of the assembled UAMS due to the low clearance between the lens and the ground.
The following subsections describe the radiation and UV sources used during the test and validation flights.
Figure 7.
Radioluminescence detection system integrated with the DJI Matrice 600 Pro and ready to fly.
Figure 7.
Radioluminescence detection system integrated with the DJI Matrice 600 Pro and ready to fly.
2.2. UV Sources
To verify the operation of the system, it is necessary to check it in conditions where alpha radiation sources or similar are present. The problem with using alpha radiation sources lies in the stringent radiation protection protocols for both the airfield and the transport to the test site, which limits the flexibility of outdoor testing. Although alpha sources were used during the final flights, in many preliminary tests they were replaced by LED sources that simulate the UV emitted by the air radioluminescence.
2.2.1. UV-C LED Sources
The UV-C LED sources used in the flights are commercially available and commonly used for water purification, disinfection, and sterilization of home aquariums. The LED has a peak wavelength between 270 nm and 280 nm and a radiant flux of 8 mW. Teflon sheets of a few millimeters are used to moderate the radiant flux to the level representative of radioactive sources that could be found in an emergency scenario. In addition, Teflon sheets serve to diffuse the light increasing its angular spread.
As per
Section 3.1 results, the calculated FOV of the detector of only a few centimeters makes the detection of point sources significantly more complex than in the case of an extended source scenario. To analyze the detection probability of point sources, five LEDs were fabricated as shown in
Figure 8. Therefore, five simulated alpha sources were distributed over the flight area to test the performance of the system in detecting the point sources.
2.2.2. Americium Source
The reference alpha source used in the flights was manufactured by Eckert & Ziegler Cesio s.r.o (Czech Republic). The provided sealed foil of
241Am source had an activity of 100 MBq (activity tolerance: + 0%, - 30%).
241Am is an alpha and gamma emitter [
17], with the most probable alpha emission at 5485.56 keV (84.45%) and the most probable gamma emission of 59.54 keV (35.92%).
According to the manufacturer, the 241Am source consists of a thin AuPd layer where the 241Am matrix of 1 µm thickness has been homogeneously incorporated. The source has an overall thickness of 0.25 mm consisting of a silver backing and a very thin layer of gold (< 2 µm thickness) covering the alpha-emitting face. The lateral dimensions of the foil are:
To fulfil the radiation source handling requirements, the
241Am foil was mounted on an aluminum base plate and covered with a quartz dome optically transparent to the radioluminescence produced by alpha particles in the air (see
Figure 9 and
Figure 10). The chamber design allows it to be evacuated to approximately 1 hPa and purged with gas at pressures up to 1000 hPa (environmental pressure).
From the radiological point of view, this enclosed source is considered sealed, and the only concern is the external exposure to 241Am gamma emissions. The 12 cm-high dome, with a 5 mm-thick wall and nominal outer diameter of 20 cm, provided a barrier to avoid the dispersion of the radioactive material in case of an incident. The ambient dose equivalent rates, Ḣ*(10), measured for the 100 MBq source in contact with the dome wall, were 30 µSv h-1 (top of the dome) and 15 µSv h-1 (lateral position). An exclusion zone of 2 m around the sources was applied during the tests in the field to further reduce radiation doses to all personnel to the level of background ambient dose equivalent rate level (0.1 µSv h-1).
The LEDs and the 100 MBq
241Am source purged with N2+NO gas mixture have been used for the test flights with the Fresnel lens detector performed in the Drone Research Laboratory (DroneLab). The left side of
Figure 10 shows the way the gas was administered to the americium source, and the right side of the figure shows how the count rate increases when the source is purged.
Figure 10.
The N2+NO (50 ppm) gas administration to the 100 MBq 241Am source (left panel) and the corresponding increase in the detector’s count rate (right panel).
Figure 10.
The N2+NO (50 ppm) gas administration to the 100 MBq 241Am source (left panel) and the corresponding increase in the detector’s count rate (right panel).
2.3. Software Architecture
The presented hardware uses the RIMASpec software architecture already used in the previous gamma detection projects [
1,
2]. For the RemoteALPHA project, the software was upgraded to visualize a map of alpha-emitting radionuclides in the environment in real time.
Figure 11 illustrates the software architecture of the aerial and ground segments. The diagram displays how the software services cooperate to geolocate the detector measurements, record them in a CSV file, and send the data to the ground control station for display on a map. The communication between the ground and air segments is managed by a server using the MQTT protocol [
18]. The UAS operation and data visualization on the ground station are done via the developed web browser application.
2.3.1. Air segment Software Architecture
This subsection details the operation of the air segment and the role of each service. The air segment is a distributed system where each service has a specific functionality that must be fulfilled for the whole system to work. The air segment software operates on the Raspberry Pi 4 board using the Raspberry Pi OS (previously known as Raspbian), which is the officially supported operating system.
The different services exchange information through a data bus. RabbitMQ—a message-queueing software that uses the publish/subscribe paradigm [
19]—is used as a data bus. RabbitMQ facilitates the transfer of messages between connected processes: services can publish information to the bus, allowing other services interested in that information to subscribe to it.
The air segment operates as follows: the software service in charge of interacting with the UAS autopilot, called DJI Autopilot Service, publishes telemetry to the data bus at a frequency of 20Hz. At the same time, the service in charge of managing the laser altimeter, called Laser Altimeter Service, publishes the UAS altitude to the data bus with a frequency of 10Hz. Finally, the RemoteAlpha Detector Service publishes the detector’s count rate every 100 ms. The Real Time Data Processing Service subscribes to all this information and merges the data. In the process, it selects the best UAS telemetry sample to geoposition the detected count rate. Once the data has been joined together, this same service publishes the geo-tagged count rate to be saved in a CSV file by the CSV Builder Service and sent to the ground control station by the GCS Comm Service using LTE (4G) network.
2.3.2. Ground Segment Software Architecture
To visualize the count rate on a map in real time, a web application was developed. This web application has a backend and a frontend as shown in
Figure 13. The backend has been developed using the Apache Tomcat Server (
https://tomcat.apache.org/). Apache Tomcat is a popular open-source web server and Servlet container for Java code. For this work, this server was hosted on the same machine as the MQTT server. The backend is subscribed to the MQTT data through which all the information from the air segment arrives. The air segment information is processed and ready to be displayed on the different web browser clients (frontend). The communication between the backend and the frontend is done through web sockets and HTTP queries.
Figure 12.
RIMASpec ground segment software architecture.
Figure 12.
RIMASpec ground segment software architecture.
Web application allows any computer with an Internet connection and a browser to access the application using a username and a password requiring no additional software installation to run the RIMASpec ground control station. The web application has a database (MongoDB) where the users are created and where all the measurements with the drone are stored. The application has three user types: the administrator, the operator and the viewer. The viewer type is intended for users who are only interested in radiological information. Under the role of operator, the user has control over the UAS. The administrator can manage user profiles on top of the previously described functionality.
Once logged in, the screen shown in
Figure 14 appears with three panels representing application options.
Figure 13.
RIMASpec UAS control tab.
Figure 13.
RIMASpec UAS control tab.
Panel 1 has the tools needed to create and edit a flight plan for the UAS. Through this panel, the user can create or edit an operation area and save it in a file. Then, the software can generate a flight plan over that area: for example, a random flight plan, or a flight plan with different parallel tracks. Some of these flight plans are shown in
Section 3. Panel 2 presents the tools to work with the map. This panel allows zooming in and out on the map or centering the map at the UAS location. Finally, Panel 3 alerts the user when it detects a UAS that wants to connect to RIMASpec.
Figure 14.
Screenshot of RIMASpec during a test flight.
Figure 14.
Screenshot of RIMASpec during a test flight.
As mentioned above, all flight and detector data are stored in the database, and RIMASpec can replay previous flights as shown in
Figure 15 which displays a screenshot of a flight performed in December 2022. The color and height of the cylinders indicate the counts measured in 0.1 seconds. The color range and height of the cylinders can be modified to optimize the visualization of the source detection. The detected counts can be also read numerically in the right-hand data box. In this data box, the UAS telemetry is also displayed: speed, pitch, roll, yaw, longitude, latitude, and altitude. Behind the telemetry box, the commands box is placed which contains the dialog for each possible command that can be sent to the UAS. The visualization software also shows the UAS position and the location on the ground where the detector is pointed (violet line). In addition, the UAS path can also be shown on the map but hidden in
Figure 15 to better see the count rate markers.
RIMASpec software has more extensive functionality beyond the application-related features shown here, such as user management tools, layer visualization, and flight plan editing.
3. Results and Discussion
Before performing the flights with the developed sources (LED and 241Am) the Fresnel lens detector had to be calibrated in terms of the counting efficiency and the field of view. This was necessary since both these parameters are crucial to designing the optimum flight plans to localize the alpha sources.
3.1. Focus Configuration and the Field of View
The measurement setup illustrated in
Figure 3 has been used to obtain the optimal focus configuration, the detector FOV, and the efficiency. The LED source was placed at a distance d ranging from 2 to five meters from the Fresnel lens. The PMT position F defines the detector protrusion from the frame (see
Figure 3) which is easily measured with a ruler and represents the focal length according to the formula:
The position
F was varied in the range from 1 cm to 10 cm during the tests. The results of the count rates measured at constant LED power simulating 7.3 GBq are shown in
Figure 16 with the measured most effective detector position marked with a red circle.
Figure 15.
Count rate measured at different source-lens distances and different focus configurations F. Uncertainties in the figure correspond to standard deviations given with a 68% confidence interval (𝑘 = 1) and are attributed to the counting statistics of the detection system.
Figure 15.
Count rate measured at different source-lens distances and different focus configurations F. Uncertainties in the figure correspond to standard deviations given with a 68% confidence interval (𝑘 = 1) and are attributed to the counting statistics of the detection system.
Figure 16 shows the influence of the detector position
F on the detected counts for a point source located at the optical axis of the lens. At a distance of 2 m, the detector measures approximately 6000 s
-1 at the optimal position between 6 and 7 cm; at a distance of 3 m, the count rate drops to 2500 s
-1, and
F is between 3 and 5 cm. At larger distances, the optical detector position
F is approximately 3 cm. The maximum count rate at a 5 m distance is around 800 s
-1. The dimensions of the laboratory did not allow testing longer distances.
Figure 17 shows the count rates measured for every lateral position
h (perpendicular to the optical axis shown in Fig. 3) of the LED source used to calculate the FOV and AFOV of the system at different distances to the source and focus configurations F. These parameters were estimated from the Gaussian fits of corresponding lateral profiles (see
Table 2). From
Table 2 it is apparent that FOV is relatively small: for a distance of 5 m and a focus configuration at 3 cm, the FOV is 10.7 centimeters.
Figure 16.
Count rates measured at different source-lens distances and lateral positions from the optical axis.
Figure 16.
Count rates measured at different source-lens distances and lateral positions from the optical axis.
For the scenario where there are point alpha sources, the flight plan tracks should be very close together, the speed of the UAS should be very low, and the integration time should be short to maximize the probability of detecting them. Another concern with this type of scenario is that the UAS is not always completely perpendicular to the ground due to variations in wind and UAS speed resulting in changes in the roll and pitch of the whole flying system. These variations in the roll and pitch of the UAS reduce the number of counts acquired compared to a vertical position that is completely aligned with the ground. On the other hand, this variation in roll and pitch increases the probability of detecting the source because the area covered by the detector is larger than when the UAS is in constant vertical alignment.
In the case of an extended source, the detection efficiency can be obtained by laboratory calibration. In this case, the uncertainties of the detector pointing vector are not significant for the source detection and the activity calculation. In such scenarios, the detection efficiency
ε, given in units of s
-1 MBq
-1, and the FOV can be used to calculate the detection efficiency per unit area ε
s given in units of s
-1 MBq
-1 cm
2:
Table 3 and
Table 4 list the detection efficiencies and corresponding source activities calculated according to the
Table 2 data for an extended source. With relative uncertainties of the LED source and FOV of 10%, the estimated uncertainty for ε
s is about 22% [
21].
Figure 18 juxtaposes the efficiency per unit area for
F= 3 cm and
F = MAX focus configurations at different source-lens distances. The reduction in the detection efficiency per unit area can be clearly seen for the
F = MAX case: while the count rate is maximized due to the tighter focus, the FOV is reduced, resulting in a lower count rate per unit area compared to the focus configuration maximizing FOV (
F = 3 cm).
Figure 17.
Comparison of the detection efficiency per unit area for different focus configurations F.
Figure 17.
Comparison of the detection efficiency per unit area for different focus configurations F.
Finally, the decision thresholds have been calculated to decide whether there is a significant count rate increase that indicates the presence of an alpha source. Following the ISO11929-4, the decision threshold
*, with a 95% confidence level, is calculated with the following equation:
where
ũ(0) is the standard deviation of background measurements.
During the flights with high solar irradiance, typical of the Spanish region, the mean and standard deviation of the UV-C background measured with the Fresnel system is around 26 ± 28 s-1 for an integration time of 0.1 s. Therefore, the decision threshold to identify a possible alpha source is about 46 s-1. This net count rate represents a surface activity decision threshold of about 4.6 MBq cm-2 for integration times of 0.1 s. The 100 MBq 241Am source with a 20 cm2 active area has an activity per unit of area of 5 MBq cm-2, which is too close to the decision threshold. Therefore, in order to detect the americium source, the radioluminescence intensity had to be amplified. This has been done by purging the quartz dome enclosure with N2 and 50 ppm NO gas mixture.
3.2. Optimum Flight Plans for Localization and Mapping of Alpha Emitters
Following the choice of the optimum height of 5 m elaborated in the previous section, another crucial flight plan parameter, the flight speed, is considered. The flight speed is primarily determined by the field of view of the detector. With the FOV of about 10 cm, the measurement time of 100 ms necessitates flying at a speed of 1 m s-1.
The flight planner software module of the RIMASpec ground control station can generate flight plans for a UAS that cover an area with different patterns such as back and forth (
Figure 19a), random, and random with genetic algorithm (
Figure 19b).
Figure 18.
Examples of the scan patterns that can be generated by the RIMASpec GCS.
Figure 18.
Examples of the scan patterns that can be generated by the RIMASpec GCS.
For the test flights, a back-and-forth scan pattern was chosen to cover the whole area and to have a certain repeatability in the flight execution that the random algorithm does not offer. To cover the complete area, a distance between parallel lines of 10 cm needs to be used. However, the UAS does not permit placing waypoints in such close proximity, therefore the shortest distance allowed by the UAS, approximately 70 cm, was selected. Given the manufacturer-specified horizontal positioning error of 1.5 m, it will be necessary to perform several iterations of the parallel line scan pattern with a small offset to ensure that the UAS passes through areas it omitted in the previous iteration. The upper panels in
Figure 20 illustrate the scheme of two flight plans while the bottom panels show real flights conducted according to this idea.
Figure 19.
One iteration of the UAS flight plan (left panel) followed by the second iteration of the flight plan with offset (right panel).
Figure 19.
One iteration of the UAS flight plan (left panel) followed by the second iteration of the flight plan with offset (right panel).
Each cylinder in
Figure 20 corresponds to a measurement made by the detector. The markers are placed where the measurement was taken, meaning the projection is made towards the ground, taking into account the roll and pitch of the UAS.
Figure 20 demonstrates how the measurements are accumulated at the ends of the flight plan tracks where the UAS decelerates, stops, and changes its course to continue the path. It can also be observed that the decelerating before the turn causes the UAS to perform measurements outside the flight plan due to the pitch change of the UAS. In the different lines of the flight plan, due to the small variation of roll and pitch and the wind effects, the measurements are not completely aligned with the flight plan as mentioned before. These variations in roll and pitch, together with the GPS errors, make it very difficult to ensure that by following or repeating the flight plan the whole area will be scanned.
This effect can be seen even more clearly when comparing a real flight with a flight performed with the manufacturer's simulator.
Figure 21 demonstrates the real and simulated flights side by side: in a simulation, the measurements are perfectly aligned.
Figure 20.
Comparison of a real flight (left panel) with a simulated flight (right panel).
Figure 20.
Comparison of a real flight (left panel) with a simulated flight (right panel).
Considering the mentioned challenges and uncertainties associated with the complete area coverage, the flights in the DroneLab were performed using the following parameters:
Speed: 1 m s-1
Altitude above ground level: 5 m
Distance between parallel lines: 70 cm
Detector integration time: 0.1 s.
Scan pattern type: back and forth
If the sources are not detected in the first iteration of the flight plan, a second iteration with an offset will be performed.
3.3. Defining and Preparing the Flight Area
The flights were carried out at the UPC's Drone Research Laboratory (DroneLab), located in north-eastern Spain. The DroneLab is a research laboratory in the field of UAS that provides the infrastructure to carry out test flights, validation, and training for any kind of professional or research application based on UAS. The infrastructure has an 80 × 45 × 15 m3 area, which is secure and separate from the general airspace, covered by a protective net that guarantees safety. The DroneLab allows flying in a controlled environment without the need for authorization from the authorities, despite being in an urban environment and close to the Barcelona city airport.
Figure 21.
Typical configuration used for test flights with LED sources at the UPC’s DroneLab.
Figure 21.
Typical configuration used for test flights with LED sources at the UPC’s DroneLab.
Figure 22 shows the flight setup with the LEDs placement and the take-off/landing area that has been used in the different test flights during the project. The landing pad made of artificial grass protects the lens from any potential damage (e.g.,
dust particles) during the landing and takeoff.
3.4. Flight Test with LEDs and 241Am Source
As mentioned before, due to GPS uncertainties and the wind effects (roll and pitch variations), the same flight plan repeated many times does not cause UAS to fly exactly over the same place and, therefore, the measurements are scattered between consecutive runs. Without this repeatability, it is very difficult to obtain the uncertainty of the source location of the point sources.
For instance, during experiments with five LEDs placed in random positions, not all of them were detected during the same iteration. During the first flight, three of the five LEDs were detected. Using the same flight plan in the following flights, either no LEDs or all of them were detected. Concluding the flight results, the probability of detecting a point source constitutes approximately 60%. This probability can be boosted by increasing the field of view of the detection system. However, increasing the FOV also has some other consequences such as the reduction of the detector efficiency and background light leakage. Extended sources with larger lateral dimensions (e.g., on the order of GPS uncertainty), on the other hand, would not be susceptible to difficulties encountered in the case of point sources.
Figure 23 shows a flight example with 5 LEDs surrounding the americium source in the middle. In a single flight plan execution, the measurement system was able to detect the americium source and 4 out of 5 LEDs. In the figure, the colormap is set such as the yellow color (background) is set to 20 s-1 and the red color shows the count rate of 50 s-1 or more, with intermediate count rate values having a color between yellow and red.
Figure 22.
Flight performed with 5 LEDs and an americium source.
Figure 22.
Flight performed with 5 LEDs and an americium source.
4. Conclusions
The main conclusion from the results presented in this manuscript is that alpha particles can be detected remotely from an unmanned aircraft system (UAS) by using optical detectors. To the extent of the authors’ knowledge, this manifests the first time that alpha particles have been remotely detected from a UAS. The described prototype system, with certain operational improvements, can be used to assess contaminated areas without the need to put first responders at risk.
Another important conclusion of this work is that the operational experience gathered with the novel radioluminescence detection system can be readily applied to the extended alpha source detection (spread over the area), while the point source localization needs further hardware improvements. All test flights carried out during this study were performed using point sources – LEDs or the americium source – due to the difficulty of reproducing an extended source. Even though point sources are difficult to detect given the optical constraints, the prototype system was able to detect them with a reasonable probability of about 60%. Extended sources are far less dependent on alignment stability and are therefore much easier to detect.
Based on the results obtained in this study and in the framework of the project RemoteALPHA, it can be observed that:
The Fresnel-lens detection system used to measure the alpha-induced radioluminescence was relatively large. This is because maximizing radioluminescence throughput while minimizing background signal requires a compromise between the diameter of the receiving optics (which determines the geometric efficiency of the system) and the focal length (which determines the FOV, that subsequently affects both the number of detected photons and filtering [
12]
The detector's size makes the drone integration difficult which means it can only be integrated into medium- and large-sized UASs, with the sensor itself acting as a landing gear. The lens, located at the aircraft's underbelly deep below its center of gravity, resists UAS movement. Flight tests conducted in windy conditions revealed significant instabilities in UAS flight. It is not advisable to exchange the landing gear supplied by the manufacturers for incorporating the detector. Therefore, in future iterations, the compact detector design must be implemented.
The major limiting factor of the developed detection system towards the localization of point sources is its small field of view. Since the FOV of the detection system (up to 10 cm) is smaller than the uncertainty of the UAS GPS (about a meter), the detection of point sources becomes more difficult. To mitigate this issue, two clear operational improvements are therefore proposed: i) Equipping the UAS with RTK to achieve centimeter-level navigation accuracy. Results show that when the system identifies a point source during different flights, the ground station positions it at nearby locations, but not at the same location. This may be due to the GPS position uncertainties themselves. That is, with the current system we can affirm that there is alpha activity over that area, but not the exact location and how many point sources there are. With an RTK-ready UAS, the enhanced navigation would allow repeating the same flight plan: this would make passes and identify source positions more consistently and would allow estimation of the uncertainties related to the source localization; ii) Increase the FOV of the detector. This would require enlarging the PMT surface, reducing the focal length of the objective, or using mirrors with shorter focal lengths. However, this option would require a tradeoff with background light leakage [
12]. An improved FOV would open the possibility of increasing the speed of the UAS and the distance between flight tracks, thus increasing the size of the area scanned with a single flight.
Author Contributions
Conceptualization, P.R., M.L., A.V., V.D. and F.K.; methodology, P.R., A.V, T.G.,D.S, J.P., M.L., V.D. and F.K.; software, P.R., T.G., D.S, J.P. and M.L.; validation, all authors.; formal analysis, P.R, A.V., D.R and M.L.; investigation, all authors.; writing—original draft preparation, P.R. and A.V.; writing—review and editing, all authors.; All authors have read and agreed to the published version of the manuscript.
Funding
The project 19ENV02 RemoteALPHA has received funding from the EMPIR programme co-financed by the Participating States and from the European Union’s Horizon 2020 research and innovation programme.
Conflicts of Interest
The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.
References
- Royo, P.; Pastor, E.; Macias, M.; Cuadrado, R.; Barrado, C.; Vargas, A. An Unmanned Aircraft System to Detect a Radiological Point Source Using RIMA Software Architecture. Remote Sens. 2018, 10, 1712. [Google Scholar] [CrossRef]
- Vargas, A.; Costa, D.; Macia Lopez, M.; Royo, P.; Pastor, E.; Luchkov, M.; Neumaier, S.; Stöhlker, U. Luff, R. Comparison of airborne radiation detectors carried by rotary-wing unmanned aerial systems. Radiat. Meas. 2021, 145, 1–12. [Google Scholar] [CrossRef]
- Podgorsak, E.B. Radiation Physics for Medical Physicists, Springer, 2006. https://cds.cern.ch/record/840186/files/3540250417_TOC.pdf.
- Burchfield, L.A. Radiation Safety: Protection and Management for Homeland Security and Emergency Response; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
- Baschenko, S.M. Remote optical detection of alpha particle sources. J. Radiol. Prot. 2004, 24, 75–82. [Google Scholar] [CrossRef] [PubMed]
- Sand, J.; Ihantola, S.; Peräjärvi, K.; Toivonen, H.; Toivonen, J. Radioluminiscence yield of alpha particles in air. New J. Phys. 2014, 16. [Google Scholar] [CrossRef]
- Sand, J. Alpha Radiation Detection via Radioluminescence of Air. PhD thesis (Tampere University of Technology. Publication; Vol. 1449). Tampere University of Technology (2016) ISBN 978-952-15-3882-7 (printed); ISBN 978-952-15-3889-6 (PDF); ISSN 1459-2045. http://www.tut.fi/tutcris.
- Kerst, T.; Toivonen, J. Intense radioluminescence of NO/N2-mixture in solar blind spectral region. Opt. Express 2018, 26, 33764–33771. [Google Scholar] [CrossRef] [PubMed]
- Krasniqi, F.S.; Kerst, T.; Leino, M.; Eisheh, J.-T.; Toivonen, H.; Röttger, A.; Toivonen, J. Standoff UV-C imaging of alpha particle emitters. Nucl. Inst. Methods Phys. Res. A 2021, 987. [Google Scholar] [CrossRef]
- Remote and Real-Time Optical Detection of Alpha-Emitting Radionuclides in the Environment (RemoteALPHA) Project Official Web-Site. Available online: https://remotealpha.drmr.nipne.ro/ (accessed on 30 August 2023).
- Parc Mediterrani de la Tecnologia, DRONELAB Section. Available online: https://pmt.es/ca/dronlab (accessed on 30 August 2023).
- Luchkov, M.; Dangendorf, V.; Giesen, U.; Langner, F.; Olaru, C.; Zadehrafi, M.; Klose, A.; Kalmankoski, K.; Sand, J.; Ihantola, S.; et al. Novel optical technologies for emergency preparedness and response: mapping contaminations with alpha-emitting radionuclides. Nucl. Inst. Meth. A 2023, 1047, 167895. [Google Scholar] [CrossRef]
- Hamamatsu Photonics Web-Site; Photon Counting Head PMT Modules Section. Available online: https://www.hamamatsu.com/us/en/product/optical-sensors/pmt/pmt-module/photon-counting-head/H11870-09.html (accessed on 30 August 2023).
- GBS Elektronik GmbH Web Site, Multi-Channel Analyzer Section. Available online: https://www.gbs-elektronik.de/en/nuclear-measurements/multi-channel-analyzer.php (accessed on 30 August 2023).
- LightWare Lidar LLC Web-Site. Available online: https://lightwarelidar.com/products/sf11-c-100-m (accessed on 30 August 2023).
- Raspberry Pi Official Website. Available online: https://www.raspberrypi.com/products/raspberry-pi-4-model-b/ (accessed on 30 August 2023).
- Sixfab Official Website. Available online: https://sixfab.com/product/raspberry-pi-4g-lte-modem-kit/ (accessed on 30 August 2023).
- Table of Radionuclides (Vol. 5 – A = 22 to 244). M.-M. Bé et al. 2010, BUREAU INTERNATIONAL DES POIDS ET MESURES. ISBN-13 978-92-822-2234-8.
- MQTT Official Website. Available online: https://mqtt.org/ (accessed on 31 August 2023).
- Available online: https://www.rabbitmq.com/ (accessed on 31 August 2023).
- European Accreditation. EA-4/02 M: 2013 Evaluation of the Uncertainty of Measurement in Calibration. 2013, p. 75. Available online: https://www.twirpx.com/file/1996254/.
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).