Journal of Astronomy and Space Sciences
The Korean Space Science Society
Technical Paper

Initial Results of Low Earth Orbit Space Radiation Dosimeter on Board the Next Generation Small Satellite-2

Uk-Won Nam1,https://orcid.org/0000-0003-2398-1019, Won-Kee Park1https://orcid.org/0000-0002-8292-2556, Sukwon Youn2https://orcid.org/0000-0002-2395-8700, Jaeyoung Kwak1,3https://orcid.org/0000-0001-7143-551X, Jongdae Sohn1https://orcid.org/0000-0002-6572-622X, Bongkon Moon1https://orcid.org/0000-0002-5106-0156, Jaejin Lee1https://orcid.org/0000-0002-3367-3346, Young-Jun Choi1https://orcid.org/0000-0001-6060-5851, Jungho Kim4https://orcid.org/0000-0001-9076-0019, Sunghwan Kim5https://orcid.org/0000-0001-7320-2488, Hongjoo Kim6https://orcid.org/0000-0001-9787-4684, Hwanbae Park6https://orcid.org/0000-0001-6087-2052, Sung-Joon Ye2https://orcid.org/0000-0001-8714-6317, Hongyoung Park7https://orcid.org/0009-0007-9332-642X, Taeseong Jang7https://orcid.org/0000-0001-7041-1874
1Space Science Division, Korea Astronomy and Space Science Institute, Daejeon 34055, Korea
2Department of Applied Bioengineering, Graduate School of Convergence Science and Technology, Seoul National University, Seoul 08826, Korea
3University of Science and Technology, Daejeon 34113, Korea
4Korea Research Institute of Standards and Science, Daejeon 34113, Korea
5Department of Radiological Science, Cheongju University, Cheongju 28503, Korea
6Department of Physics, Kyungpook National University, Daegu 41566, Korea
7Satellite Technology Research Center, Korea Advanced Institute of Science and Technology, Daejeon 34141, Korea
Corresponding Author : Tel: +82-42-865-3221, E-mail: uwnam@kasi.re.kr

© Copyright 2024 The Korean Space Science Society. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Aug 12, 2024; Revised: Aug 30, 2024; Accepted: Aug 31, 2024

Published Online: Sep 30, 2024

Abstract

As human exploration goals shift from missions in low Earth orbit (LEO) to long-duration interplanetary missions, radiation protection remains one of the key technological issues that must be resolved. The low Earth orbit space radiation dosimeter (LEO-DOS) instrument to measure radiation levels and create a global dose map in the LEO on board the the next generation small satellite-2 (NEXTSat-2) was launched successfully on May 25, 2023 using the Nuri KSLV-III in Korea. The NEXTSat-2 orbits the Earth every 100 minutes, in an orbit with an inclination of 97.8° and an altitude of about 550 km above sea level. The LEO-DOS is equipped with a particle dosimeter (PD) and a neutron spectrometer (NS), which enable the measurement of dosimetric quantities such as absorbed dose (D), dose equivalent (H) for charged particles and neutrons. To verify the observations of LEO-DOS, we conducted a radiation dose estimation study based on the initial results of LEO-DOS, measured from June 2023 to September 2023. The study considered four source categories: (i) galactic cosmic ray particles; (ii) the South Atlantic Anomaly region of the inner radiation belt (IRB); (iii) relativistic electrons and/or bremsstrahlung in the outer radiation belt (ORB); and (iv) solar energetic particle (SEP) events.

Keywords: low Earth orbit; space radiation; next generation small satellite-2 (NEXTSat-2); charged particles dose; neutrons dose; dose map

1. INTRODUCTION

Space radiation dosimetry is one of the most significant challenges in radiation protection due to the complex radiation fields encountered in low Earth orbit (LEO) and interplanetary space. As human exploration shifts from LEO missions to long-duration interplanetary missions, addressing radiation protection becomes a critical technological issue that must be resolved. In outer space, humans are continuously exposed to natural ionizing radiation, which can cause various harmful biological effects. Therefore, it is crucial to quantitatively assess exposure levels to estimate radiation hazards during manned space missions (Semkova et al. 2014; Wilson et al. 2002).

Low Earth orbit space radiation dosimeter (LEO-DOS) is a science payload of the next generation small satellite-2 (NEXTSat-2), which was launched successfully on May 25, 2023. The NEXTSat-2 (Fig. 1) orbits the Earth every 100 minutes, in an orbit with an inclination of 97.8° and an altitude of about 550 km (Jang et al. 2022). LEO-DOS has two types of instruments, a particle dosimeter (PD) and a neutron spectrometer (NS). The PD measures the effects of ionizing energy loss in matter due to penetrating solar energetic protons (SEP) and galactic cosmic rays (GCR), specifically in silicon solid-state detectors and after interactions with tissue-equivalent plastic (TEP), a synthetic analog of human tissue. The PD investigation quantifies the linear energy transfer (LET) spectrum in materials directly through measurements within the space radiation environment, particularly focusing on interactions of ions with energies above 40 MeV. The NS instrument is based on a stilbene organic scintillator (C14H12), enabling it to function as a neutron dosimeter for fast neutrons in the energy range of 1.6 to 8 MeV.

jass-41-3-195-g1
Fig. 1. Illustration of NEXTSat-2 spacecraft and LEO-DOS instrument. The LEO-DOS to view space is mounted below solar panel and the top deck of the satellite. The illustration of LEO-DOS. The LEO-DOS consist of a PD and NS instruments. NEXTSat-2, the next generation small satellite-2; LEO-DOS, low Earth orbit space radiation dosimeter; PD, particle dosimeter; NS, neutron spectrometer.
Download Original Figure

LEO-DOS objectives are

  • To create a global map of absorbed dose and dose equivalent of space radiation from charged particles and neutrons.

  • To study changes in space radiation and the impact of the space environment during the increasing phase of the solar activity cycle.

  • To analyze the space radiation risk to astronauts, considering neutron weighting factor in near-Earth space.

  • To assess the effect of space radiation on electronic components due to neutrons in near-Earth space.

We conducted a radiation dose analysis based on the initial results of LEO-DOS, measured from June 2023 to September 2023. To validate the observations of LEO-DOS, we categorized the data into four source categories as described in the reference paper (Dachev et al. 2017): (i) GCR particles; (ii) the South Atlantic Anomaly region (SAA) of the inner radiation belt (IRB); (iii) relativistic electrons and/or bremsstrahlung in the outer radiation belt (ORB); and (iv) solar energetic particle (SEP) events. In this paper, we present the design, calibration, and early flight data analysis of the LEO-DOS instrument. Detailed descriptions of the instrument configuration and design are provided in Section 2. Section 3 covers the calibration of the instrument, followed by the radiation source separation process from the measurements in Section 4. Finally, Sections 5 and 6 provide this paper’s analysis, discussion, and conclusions..

2. LOW EARTH ORBIT SPACE RADIATION DOSIMETER (LEO-DOS) INSTRUMENT DESCRIPTION

LEO-DOS has a compact, low-mass package (1.7 kg) and a total dimension of 190 × 85 × 144 mm3. The PD is designed to cover a range of 0.3–700 keV/μm LET spectra of charged particles in LEO. On the other hand, the NS consists of a combination of a stilbene organic scintillator and a plastic scintillator for VETO in a compact package that covers a range of 2–8 MeV of fast neutron energy range.

The LEO-DOS PD and NS are mounted on the satellite, with potential shielding from the electronic box, satellite body, and solar panel. However, GCR, the main dose contributor to the PD, has very high energy (several GeVs) and can penetrate most satellite materials. SEP and SAA protons also have sufficient energy so the aluminum body shields only a small fraction. Additionally, the mean chord length value used in dose calculations is the same whether the field of view (FOV) solid angle is 4π steradians or 2π steradians (refer to Eq. (2)). Radiation from behind is not effectively shielded due to the satellite structure. For NS, fast neutrons have a low cross-section for aluminum, allowing most to pass through the satellite body. Consequently, the FOVs of PD and NS are considered omni-directional.

Table 1 summarizes the main specification of the LEO-DOS instrument, including physical resources, key design values, and measurement performance metrics. These properties of the as-delivered flight model represent the optimal values resulting from science and engineering trade studies completed throughout the instrument development.

Table 1. LEO-DOS instrument specification
Parameters PD NS
Performance Properties Proton > 40 MeV
Electron > 3 MeV
Proton > 35 MeV
Electron > 2.5 MeV
Time resolved LET : 0.3–700 keV/μm Energy range : 2–8 MeV
- Flux (counts/sec/cm2)
- Absorbed dose rate (D, μGy h–1)
- Dose equivalent rate (H, μSv h–1)
- Fast neutron & gamma spectrum
- Flux (counts/sec)
- Absorbed dose rate (D, μGy h–1)
- Equivalent dose rate (H, μSv h–1)
Geometric factor 25.13 cm2 sr (2 cm2) 214.90 cm2 sr (17.10 cm2)
FOV Omni-direction Omni-direction
Shield 1.8 g cm–2 1.8 g cm–2
System Dimension 190 (W) × 85 (D) × 144 (H) mm3
Mass 1.7 kg
Power 5.5 W
Interface CAN
Transmission 350 Mb/day

LEO-DOS, low Earth orbit space radiation dosimeter; PD, particle dosimeter; NS, neutron spectrometer; LET, linear energy transfer; FOV, field of view.

Download Excel Table
2.1 Particle Dosimeter (PD)

The PD of LEO-DOS employs a silicon photodiode particle sensor, which is situated between two 10 mm layers of A-150 TEP, as illustrated in Fig. 2. This configuration, with one layer above and one below the sensor, allows for the assessment of the effective dose to the human body in accordance with ICRP Publication 103 (Charles, 2008). The particle sensor for the PD is a n-type silicon photodiode sensor with an effective area of 20 × 20 mm and 650 ± 30 μm thickness, developed for the application to the cosmic ray energetic and mass for the international space station experiment for the International Space Station (ISS) in Korea. Its properties as a particle detector are investigated in reference (Hyun et al. 2015).

jass-41-3-195-g2
Fig. 2. Schematic view of a particle dosimeter of LEO-DOS. The silicon photodiode particle sensor is sandwiched between two layers of 10 mm A-150 tissue-equivalent plastic, with one layer on the upper side and one on the lower side. An A225 charge sensitive preamplifier and shaping amplifier is placed close to the silicon sensor to reduce stray noise. LEO-DOS, low Earth orbit space radiation dosimeter.
Download Original Figure

The absorbed dose for charged particles with a silicon detector is calculated by dividing the sum of the deposited energy in the silicon sensor by the mass of the silicon sensor. The calculation formula for the absorbed dose in silicon can be expressed as follows:

D S i = 1 m i = 1 N N i ε i
(1)

where Ni is the count of the i-th channel, εi is the energy of the i-th channel, and m is the mass of the silicon sensor.

DSi is converted to DH20 by converting LET in silicon (LSi) to LET in water (LH20) through the Benton equation (Benton et al. 2010) as follows:

D H 20 = i = 1 N N i L H 20 l ¯ m H 20 , Log ( L H 20 ) = 0.2902 + 1.025 Log ( L S i )
(2)

where is a mean chord length of the PD of 919 μm obtained by dividing the thickness of the silicon sensor of PD of 650 μm by cos (45°) considering the isotropic incidence of the charged particles in space.

The biological effectiveness factor, also known as the quality factor (Q(L)) can be expressed as a function of LET (L) according to ICRP publication 60 as follows (ICRP 1991):

Q ( L ) = { 1 ( 0 < L 10 ) 0.32 L 2.2 ( 10 < L 100 ) 300 / L ( L > 100 )
(3)

The dose equivalent can be expressed as the product of absorbed dose and mean quality factor (Qmean) as shown in the following equation:

H H 2 O = Q m e a n D H 2 O
(4)
2.2 Neutron Spectrometer (NS)

Fig. 3 presents a schematic view of the NS configuration. A 1.5” × 1.5” stilbene (C14H12) organic scintillator was adopted as a sensor of a fast NS. To mitigate the influence of charged particles in near-Earth space, a hollow plastic sensor (EJ-200, ELJEN Technology, Sweetwater, TX, USA) coated with EJ-510 reflective paint is installed around the stilbene sensor, serving as a VETO sensor to filter out these particles. Each scintillator sensor is optically coupled with an array of SiPMs (J-series 60035, SensL, now ON Semiconductor, Cork, Ireland). The stilbene sensor is coupled with a 4 × 4 SiPM array, while the VETO sensor is coupled with a 24-SiPM array.

jass-41-3-195-g3
Fig. 3. Schematic view of a neutron spectrometer which consist of stilbene and VETO sensor. The VETO sensor, designed with a hollow structure, employed an EJ-200 plastic detector coated with EJ-510 reflective paint.
Download Original Figure

The advantage of using stilbene for a fast neutron dosimeter is that it does not require an unfolding algorithm to obtain dosimetric parameters such as the absorbed dose (D) and equivalent dose (H) (Miller et al. 2020). Additionally, since the stilbene sensor has tissue-equivalent characteristics, the impact of neutrons in LEO on the human body can be evaluated through the stilbene absorbed dose.

The absorbed dose (Dstilbene) in stilbene is calculated by dividing the sum of the recoil proton energy in the stilbene crystal by the mass of the stilbene sensor. The calculation formula for the absorbed dose in stilbene in the energy range from 2 to 8 MeV neutrons can be expressed as follows:

D s t i l b e n e = 1 m i = 1 N N i ε p , i
(5)

where Ni is the count of the i-th channel, εp,i is the recoil proton energy of the i-th channel, and m is the mass of the stilbene crystal.

Since the stilbene scintillator in NS has a quenching effect, the measured electron equivalent energy (εee) in stilbene was converted to proton energy (εp) using the light output function to obtain the recoil proton energy (Hansen et al. 2002).

ε e e = 0.693 ε p 3 [ 1 e 0.2 ε p 0.965 ]
(6)

The equivalent dose (Hstilbene) for neutrons in low Earth orbit, measured by the stilbene sensor, can be obtained by multiplying the stilbene absorbed dose (Dstilbene) by the radiation weighting factor (wR) according to the energy of the neutron, as suggested in ICRP 60 (ICRP 1991). In ICRP 60, the radiation weighting factor is 10 for neutrons in the range of 2–20 MeV. The equivalent dose in the stilbene in the energy range from 2 to 8 MeV neutrons is as follows:

H s t i l b e n e = D s t i l b e n e × 10
(7)
2.3 Signal Processing for Low Earth Orbit Space Radiation Dosimeter (LEO-DOS)

The circuit configuration of the LEO-DOS system comprises a PD, NS, and a Space Craft Interface (SCIF) electronics, as illustrated in Fig. 4. All components used in the LEO-DOS readout electronics are commercial off-the-shelf parts, screened to withstand a total ionizing dose level of 30 krad.

jass-41-3-195-g4
Fig. 4. Schematic diagram of LEO-DOS signal processing showing a PD readout, NS readout and SCIF part. PSD, pulse shaping discrimination; LEO-DOS, low Earth orbit space radiation dosimeter; PD, particle dosimeter; NS, neutron spectrometer; SCIF, Space Craft InterFace.
Download Original Figure
2.3.1 Particle Dosimeter (PD) Signal Pulse Processing

The electronics developed for the PD signal processing consist of silicon sensor preamplifier and PD-amp (amplifer), a analog pulse processor (PD-APP) as shown in Fig. 4. The silicon sensor, positioned between two 10 mm thick A150 TEP plates, is biased at 160 V. The signal from the silicon sensor is amplified using the AMPTEK A225 (https://www.amptek.com), a preamplifier and shaping amplifier with a sensitivity of 5.2 V/pC and a peaking time of 2.4 μs. To achieve a wide range of LET spectra (0.3–700 keV/μm), we utilized two gain modes: high-gain and low-gain, with a gain ratio of 30:1. In the PD-analog pulse processor part, the signals from the high-gain and low-gain amplifiers are peak-held and digitized by a 12-bit analog-to-digital converter. The conversion time is set to 40 μs to achieve a fixed dead time. After being reduced to 9-bit (512 channels) data, the digitized data are stored as histograms using the PD-MPU (MSP432P401R, https://www.ti.com), resulting in a spectrum with 512 ADC channels.

2.3.2 Neutron Spectrometer (NS) Signal Pulse Processing

A cylindrical 1.5” × 1.5” Stilbene (C14H12) scintillator was optically coupled with a 4 × 4 Silicon photomultipliers (SiPM) array. SiPMs are a solid-state alternative to PMTs for detecting faint scintillation light. Individual SiPMs consist of arrays containing thousands of avalanche photodiode pixels operating in Geiger mode. Although SiPM parts are relatively small (e.g., the SensL J series is available in 6 × 6 mm square formats), arrays of SiPMs with uniform output can cover larger areas.

The bias voltage 27.5 V for the SiPM was supplied from the C14156 (https://www.hamamatsu.com), bias power supply with built-in temperature compensation function. The outputs of the 16 individual SiPMs were passively summed and amplified through an OPA656N(https://www.ti.com) trans-impedance amplifier. The output signals were connected to a pulse shaping discrimination (PSD) circuit, enabling the separation of gamma signals and neutron signals. For the PSD function, the total charge integrator (T), short charge integrator (S), and long charge integrator (L) are implemented in the NS-Charge-Integrator circuits. To quantitatively distinguish neutron and gamma-ray interactions, we define a PSD parameter as the charge ratio:

PSD ratio = ( L S ) L
(8)

The PSD image construction is conducted in the NS-Processor circuit, with the PSD ratio plotted against the total energy, represented on the x-axis. In the PSD image shown in Fig. 5, two distinct bands emerge for scintillators that respond differently to neutrons and gammas. The upper band corresponds to the neutron response, while the lower band corresponds to the gamma response. Often, a figure of merit (FOM) is defined as the distance between the two bands in a particular energy range, divided by the sum of the full-widths at half-maximum (FWHM) of the PSD ratio bands. A high FOM is usually interpreted as a good indicator of how well a certain scintillator device can distinguish between neutrons and gammas. The FOM is then computed by dividing the difference between the centroids by the sum of their FWHM:

jass-41-3-195-g5
Fig. 5. The PSD Image of stilbene NS. The PSD ratio is plotted against the total energy which is represented on the x-axis as channel. The PSD image is stored in the internal SRAM of the NS-MPU, with a resolution of 128 × 128 pixels. The upper band corresponds to the neutron response, while the lower band corresponds to the gamma response. PSD, pulse shaping discrimination; NS, neutron spectrometer.
Download Original Figure
FOM = x n x γ FWHM n + FWHM γ
(9)

where Xn and Xγ are the peaks of the neutron and gamma-ray bands, respectively, and FWHMn and FWHMγ are the full widths at half maximum of the two distributions.

In the NS-AMP circuit, a trigger signal is generated using the level discriminator of the SiPM output signal. The trigger signal creates integration time windows of 10, 0.2, and 0.5 μs for total integration, short integration (S), and long integration (L), respectively. The output of three integrator circuits were digitized using the 12-bit ADC embedded in the MSP432P401R, which serves as the NS-MPU. The PSD image with 128 × 128 pixel was constructed and stored in the MPU’s SRAM. A good FOM of more than 4.2 was obtained with our circuit. The VETO scintillator was also optically coupled with a 24 SiPM array. In the VETO-AMP, the summed SiPM array output is shaped and amplified, and the VETO spectrum is then obtained in the NS-MPU. Meanwhile, before processing the PSD logic in the microprocessor, the presence or absence of a VETO signal is determined. If the VETO signal exceeds the set VETO threshold value, it is considered a charged particle; if it is lower, it is classified as a neutron signal. Due to the difficulty in experimentally determining the VETO threshold value for distinguishing between neutrons and charged particles, it was set to 0.5 MeV based on simulation results (see Section 3).

2.3.3 Space Craft Interface

In the LEO-DOS Space Craft InterFace (SCIF), 28 V is supplied from the satellite and converted to 12 V DC using a DC-DC converter. This is then used to generate +3.3 V, +5.5 V, and –5.5 V DC power with good regulation quality. The SCIF-MCU, which utilizes the MSP432P401Y microprocessor, includes radiation-hardened MRAM (MR25H256ACDF) for data storage, Flash ROM (MT29F2G01ABAGSDSF-ITG) for storing system control parameters, and a watchdog (TPS3828-33DBV). The PD and NS units of the LEO-DOS are interfaced via UART (RS-232) serial communication with SCIF board. LEO-DOS data packages are transmitted to the satellite on-board computer through CAN communication according to the operation sequence. To monitor the health of the LEO-DOS, it is equipped with three temperature sensors each for the PD and NS, an additional temperature sensor for the LEO-DOS body.

2.3.4 Low Earth Orbit Space Radiation Dosimeter (LEO-DOS) Data Products

The scientific data products of LEO-DOS are categorized into three groups: primary science data, secondary science data, and housekeeping data. The primary science data products are measured at 120-second intervals and include the following: low-gain spectrum (512 channels, 2 bytes), high-gain spectrum (512 channels, 2 bytes) for PD, neutron spectrum (512 channels, 2 bytes), gamma spectrum (512 channels, 2 bytes), VETO spectrum (512 channels, 2 bytes) and PSD images (78 × 128 pixels, 3 bytes) for NS. The secondary science data products of LEO-DOS, when combined with satellite attitude information to calculate the dose, are provided as the following dataset: UT date and time, altitude (km), longitude (deg), latitude (deg), L-value, PD-D/F ratio (nGy cm–2 particle–1), PD-Flux (cm2 s–1),PD-absorbed dose rate (DSi, μGy h–1), PD-absorbed dose rate in water (DH2O, μGy h–1), PD-Q value, PD-dose equivalent rate in water (HH2O, μSv h–1), PD-temperature, NS-absorbed dose rate (Dstilbene, μGy h–1), NS-equivalent dose rate (Hstilbene, μSv h–1), and NS-temperature. Housekeeping data are measured at 3-sec intervals and, when combined with satellite information, provide the following details: UT date and time, altitude (km), longitude (deg), latitude (deg), SCIF temperature, SCIF power voltage and current, count rate (cps) of PD, count rate (cps) of NS, silicon sensor bias voltage, SiPM bias voltage and current monitoring for NS and VETO.

3. INSTRUMENT CALIBRATION

For each instrument PD and NS, the relationship between the energy deposited in the detector and the digitized height of the shaped pulse in terms of the ADC value is derived considering various radioisotopes of well-known energies in the laboratory.

The LET spectrum of PD was calibrated using a 137Cs gamma source, a 90Sr beta source and a 241Am alpha source. In order to determine the dose calibration factors of PD should be obtained for both low-LET and high-LET regions. Dose calibration is conducted using the 137Cs gamma source in KORASOL (Seoul, Korea), radiation calibration company in Korea and by using the MCNP simulation. It was confirmed that the LET spectrum and the absorbed dose rate of the 137Cs gamma source obtained through the Monte Carlo simulation using MCNP version 6.2 (Werner et al. 2017) and measurement at KORASOL were very similar. The detailed calibration results were reported in the reference (Youn et al. 2023). The calibration results of the PD are presented in Table 2.

Table 2. PD calibration results
Parameters Methods Results
Energy linearity 137Cs, 90Sr, 241Am Low LET : ylow = 0.0464x (channel) – 0.403
High LET : yhigh = 1.349x (channel) – 0.2
Electronics linearity Test pulser Linearity check in full range
LET spectrum MCNP simulation LET range : 0.3–700 keV/μm
Energy range : 0.2–460 MeV
Thermal range KASI T/V experiment 5°C–35°C

PD, particle dosimeter; LET, linear energy transfer.

Download Excel Table

For the Stilbene and VETO sensors of LEO-DOS NS, the energy calibration was carried out using three gamma peaks: the 662 keV gamma peak from a 137Cs source and the 511 and 1,274.5 keV gamma peaks from a 22Na source. A good FOMstilbene value of 4.3 for the PSD of the NS was confirmed using a 252Cf radioactive source. The VETO threshold to discriminate neutrons and charged particles was determined through the Monte Carlo simulation using MCNP. The LEO neutron and GCR proton responses of the EJ-200 plastic scintillator, which constitutes the VETO sensor, were confirmed through the simulation (ANSI 2004; Sato et al. 2006). The VETO threshold value of 0.5 MeV, which can accept 91% of neutrons and reject 98% of protons in the region below 0.5 MeV, was determined. The operating temperature range of the NS was confirmed to be 5°C to 30°C using a thermal/vacuum chamber at KASI. The calibration results of the NS are presented in Table 3.

Table 3. NS calibration results
Parameters Methods Results
Energy linearity 137Cs, 60Co, 22Na y(Energy (keV)) = 6.87x (channel) – 39.357
Energy range Analysis Eee = 0.2–3.4 MeV
En = 1.6–8.5 MeV
VETO calibration 137Cs, 60Co, 22Na Eee = 0–5.3 MeV
E(Energy (keV)) = 10.662x (channel) – 163.38
Thermal range KASI T/V experiment 5°C–30°C
Figure of merit 252Cf FOMSTILBENE : 4.3

NS, neutron spectrometer; FOM, figure of merit.

Download Excel Table

4. INITIAL OBSERVATION RESULTS

4.1 Dose Measurement during Launch and Early Orbit Phase (LEOP)

LEO-DOS aboard NEXTSAT-2 was launched on May 25, 2023, at 6:24 PM. The launch and early orbit phase (LEOP) for LEO-DOS was successfully conducted from June 8 to June 15, and normal observations have been ongoing since June 22. During the LEOP period, NEXTSat-2 maintained a stable altitude and temperature, ensuring the proper operation of LEO-DOS. Throughout the day, the temperature of the LEO-DOS increased from 15°C at initial power-on to 20°C after approximately 4 hours, and then stabilized between 20°C and 30°C, depending on the attitude of NEXTSat-2.

Fig. 6 illustrates the daily variation in the dose measured by LEO-DOS on June 25, 2023, during the LEOP period. The presented data, measured at 120-second intervals, are plotted against UT. The average altitude of the satellite was 556 ± 9.29 km over 15 orbits. It shows 9 passes through the SAA regions, 30 passes through the high-latitude GCR regions of both hemispheres, 30 passes through the geomagnetic equator. The SAA pass time was approximately 9 minutes, and the GCR region pass time was around 20 minutes.

jass-41-3-195-g6
Fig. 6. Daily variation in the dose measured by LEO-DOS on June 25, 2023, during the LEOP period. The figure shows, from top to bottom, the satellite altitude, LEO-DOS temperature, PD absorbed dose, and NS absorbed dose variations throughout the day. LEO-DOS, low Earth orbit space radiation dosimeter; LEOP, launch and early orbit phase; PD, particle dosimeter.
Download Original Figure

The SAA dose rates reached the highest value of 8,445 μGy h–1 for PD and 2.8 μGy h–1. The GCR dose rates varied from 0.6 to 15 μGy h–1 for PD and from 0.002 to 0.04 μGy h–1 for NS. The daily average GCR dose rates were 4.98 μGy h–1 for PD, 0.02 μGy h–1 for NS. This variation in dose according to the orbit showed reasonable agreement with the dose variation measured on the ISS (as shown in Fig. 4 of Dachev et al. 2015), confirming that the observations are consistent and normal.

4.2 Dose Map

The geographical distribution of the dose rate for charged particles and neutrons, measured by the LEO-DOS, is shown in the two panels of Fig. 7, where horizontal axis is the longitude, vertical axis is the latitude, and the dose rate intensity is color-coded. Measurements for the three-month period from June 2023 to September 2023 are plotted. In the case of PD, the maximal dose rates are recorded in the region of the SAA from the trapped protons. In the case of NS, most observations in the SAA region were vetoed due to the influence of strongly trapped charged particles, causing the SAA region to appear empty. Occasionally, when a SEP event occurred, traces of SEP were observed along the satellite’s trajectory in the ORB region.

jass-41-3-195-g7
Fig. 7. Dose map for 3 months by LEO-DOS. (a) Absorbed dose map for PD. (b) Absorbed dose map for NS. In the case of NS, the SAA region appears empty due to the VETO effect caused by strongly trapped charged particles. LEO-DOS, low Earth orbit space radiation dosimeter; PD, particle dosimeter; NS, neutron spectrometer.
Download Original Figure

We can now create a detailed global radiation map using LEO-DOS for radiation sources such as GCR, IRB, ORB, and SEP in LEO, including dosimetric parameters such as absorbed dose and equivalent dose for charged particles and neutrons, following the source selection process described in Section 4.3.

4.3 Dose by Source Selection

The typical radiation sources in the LEO are GCR coming from beyond the solar system, trapped protons in the SAA region of the IRB, trapped electrons from the ORB, and SEP. In this section, we treat the source selection method according to the data selection procedure established by Dachev (Dachev 2009; Dachev et al. 2015).

Fig. 8 characterizes the predominant radiation sources detected by LEO-DOS, using dose rate from flux and dose to flux (D/F) dependencies. The abscissa plots the measured flux in cm–2 s–1, while the ordinate shows the dose rate in μGy h–1 and the D/F ratio in nGy cm–2 particle–1 from July to September 2023. The red experimental points along the diagonal represent dose rate values that are linearly dependent on the flux, while the black points lying almost horizontally indicate the D/F ratio. From this figure, when the dose to flux ratio (D/F) exceeds 1.12, it can be characterized as an IRB source, whereas when the D/F ratio is less than 1.12, it can be characterized as an ORB source (Dachev 2009). However, the GCR, IRB, and ORB sources in dose rate < 15 μGy h–1 region are intermixed each other, making it difficult to distinguish between them. Regions with flux < 4 cm2 s–1 exhibit values of 15 μSv h–1, and this value was used as the criterion to differentiate between GCR and ORB sources.

jass-41-3-195-g8
Fig. 8. Analysis of the main radiation sources for LEO-DOS relies on the dose rate derived from flux and the dose-to-flux (D/F) ratio dependencies observed from July to September 2023. LEO-DOS, low Earth orbit space radiation dosimeter.
Download Original Figure

Due to the difficulty in categorizing radiation sources based on the D/F and flux relationship, the L-value was used to distinguish the radiation sources. The L-value, L-shell or McILwain L-parameter is a parameter used in space physics to describe the location of a magnetic field line in the Earth’s magnetosphere (Dachev et al. 2017). It corresponds to the radial distance from the Earth’s center to the point where the magnetic field line crosses the magnetic equator, measured in Earth radii, if the magnetospheric magnetic field is a dipole field.

Fig. 9 shows the distribution of the dose rates measured with the LEO-DOS instrument against McIlwain’s L values. Fig. 9(a) represents the SAA and GCR regions, generated during the period from August 17 to August 18, 2023. The SAA is an area where the IRB comes closer to the Earth’s surface due to a displacement of the magnetic dipole axes from the Earth’s center. The IRB is located at altitudes ranging from 0.2 to 2.0 Earth radii at the geomagnetic equator. Within the SAA region, electron energies can reach up to 10 MeV, while proton energies extend up to 700 MeV. The SAA region is characterized by L-values ranging from 1 to 2.5. These findings are in good agreement with the referenced literature (Dachev et al. 2017). In this region, the radiation dose rate of the LEO-DOS is measured at 8,774 μGy h–1. The maximal particle fluxes at an L-value of approximately 1.24 are recorded as 1,706 particles cm–2 s–1. From LEO-DOS data, SAA events are characterized by dose rate > 15 μGy h–1, L value < 2.5 and D/ F > 1.0 nGy cm–2 particle–1.

jass-41-3-195-g9
Fig. 9. Distribution of the dose rates measured with the LEO-DOS instrument against McIlwain’s L values. (a) SAA and GCR region: Galactic cosmic rays (GCR) are seen everywhere and at all times. SAA (L < 2.5) region created by high energy protons when the LEO-DOS crosses the region of the SAA. (b) ORB region is generated by sporadic high-energy electrons seen mainly at 3 < L < 7. (c) SEP region are generated by SEP event of Sun. LEO-DOS, low Earth orbit space radiation dosimeter; SAA, South Atlantic Anomaly; ORB, outer radiation belt.
Download Original Figure

GCR is charged particle that originated from sources beyond the solar system and consists of 98% proton and heavy ions (baryon component), 2% electrons. GCR are seen everywhere and at all times. In regions where L > 4, the energy deposited in the silicon sensor exhibits a spectrum similar to that of GCR (Dachev et al. 2017). The GCR plot in Fig. 9(a) was generated by selecting a period from October 17 to October 30, 2023, during which ORB events were not distinguished. From this plot, we adopt an upper GCR value of D < 15 μGy h–1 as the selection criterion.

Fig. 9(b) illustrates the ORB region generated during the period from September 21 to September 24, 2023, when ORB events were characterized by very high ORB fluxes and dose rates. Sporadic high-energy electrons are seen mainly at 3 < L < 7 where the ORB is usually situated. It primarily contains electrons with energies extending up to a few MeV. Thus, they are characterized by dose rate D < 15 μGy h–1, D/F < 1.12 nGy cm2 particle–1 and flux > 4 cm–2 s–1. These criteria are essential for understanding and mitigating radiation effects in space environments. The ORB poses challenges for spacecraft and satellite design due to its elevated radiation levels.

Fig. 9(c) depicts the dose variation caused by SEP events during period from October 5 to October 9, 2023. SEPs are primarily produced by solar flares and are associated with intense solar activity. SEPs cause a clear increase in the dose in L > 4 when an event occurs. The ion energies of SEPs typically range between a few MeV and larger per nucleon. They are characterized by dose rate > 15 μGy h–1, L > 4 and D/F > 1.12 nGy cm2 particle–1. These criteria are crucial for assessing the impact of SEPs on spacecraft, satellites, and astronaut health during space missions. Understanding and mitigating the effects of SEPs is essential for safe space exploration.

5. RESULTS AND DISCUSSION

5.1 Radiation Source Analysis

Table 4 presents the statistics for radiation sources based on selected criteria from July to September 2023. These include the average dose rate (Dsi) (μGy h–1) and (DH2O) (μGy h–1), the average quality factors (Qav), and the dose equivalent rates (H) (μSv h–1) for PD, as well as the average dose rate (D) (μGy h–1) and the equivalent dose rates (H) (μSv h–1) for NS. Additionally, four SEP events during the same period are also presented.

Table 4. Statistics for the sources period from July to September, 2023
Source No. of points Charged particles Neutrons Selecting requirements
DSi (μGy h–1) DH2O (μGy h–1) HH2O (μSv h–1) Q D (μGy h–1) H (μSv h–1) Dose rate (μGy h–1) L-value D/F (nGy cm2 part.–1) Flux (cm–2 s–1)
GCR 57,273 4.72 5.56 32.24 4.38 0.03 0.30 D < 15
GCR (L > 4) 19,438 7.98 9.47 55.51 5.35 0.04 0.40 D < 15 L > 4
IRB 4,483 1,559 1,944 2,690 1.62 0.92 0.92 D > 15 L < 2.5 D/F > 1.12
ORB 2,931 22.41 21.65 73.27 3.83 0.05 0.50 D < 15 D/F < 1.12 Flux > 4
Additional Table
SEP (7/18) 216 79.01 100.56 198.00 2.35 0.09 0.90 D > 15 L > 4 D/F > 1.12
SEP (8/5) 147 47.78 60.54 150.17 3.22 0.20 2.00 D > 15 L > 4 D/F > 1.12
SEP (8/8) 219 56.92 72.39 162.30 2.48 0.23 2.30 D > 15 L > 4 D/F > 1.12
SEP (9/1) 219 41.53 53.32 142.33 3.69 0.17 1.70 D > 15 L > 4 D/F > 1.12

GCR, galactic cosmic ray; IRB, inner radiation belt; ORB, outer radiation belt.

Download Excel Table

The average absorbed dose rate (Dsi) of GCR measured with LEO-DOS PD during no quiet conditions was 4.72 μGy h–1, which is higher than the average absorbed dose of GCR measured through Liulin onboard the ISS of 2.78 μGy h–1 (Dachev et al. 2017). This difference is likely due to the inclination of the ISS, where Liulin is mounted which orbits the Earth with an inclination of 51.6°, while the NEXTSat-2, where the LEO-DOS is mounted, orbits the Earth with an inclination of 97.8°. Due to the higher inclination of the NEXTSat-2 orbit, the LEO-DOS passes through the polar region with a higher L-value, resulting in higher radiation exposure.

The average absorbed dose rates of PD for SEP and IRB sources show a more significant difference from Liulin’s measurement value than that of GCR. The reason for this difference is that the low energy protons dominant in SEP and IRB sources can be effectively shielded by the 10 mm TEP surrounding the PD, while the GCR has very high energy, most of which easily passes the 10 mm TEP. Due to this difference, a significant difference in the dose rates occurs for sources with a high proportion of low energy protons in PD with a 10 mm TEP shield and Liulin with a thin aluminum shield. In the case of the ORB region, most of the electrons below 3 MeV, which make up the majority of the ORB environment, are shielded by the 10 mm TEP of the PD, resulting in a huge difference in dose rates from Liulin’s measurement results.

Radiation protection is crucial for understanding how harmful space radiation doses are to the human body during space travel, especially in LEO. Using the tissue equivalent sensor, LEO-DOS measured the absorbed dose (D), quality factor (Q), and dose equivalent (H) of space radiation for charged particles and neutrons.

The cumulative dose over 30 days for charged particles was approximately 0.164 Sv, significantly lower than the maximum recommended safety level of 1.5 Sv for skin in LEO by the National Council on Radiation Protection (NCRP 2000). These results are consistent with those observed by Dachev et al. (2017) in the ISS radiation environment. For neutrons, the observed value was approximately 1/100th of that for charged particles, indicating a much safer level.

5.2 Solar Energetic Particle Event

Fig. 10 illustrates the variations in SEP events from August 5 to August 12, 2023. The dose measured by LEO-DOS PD, plotted alongside the solar proton fluxes recorded by GOES-16 in geosynchronous orbit, shows the solar proton fluxes in the upper panel and the dose variation at 120-second intervals in the lower panel. The strong responses that regularly appear in the daily dose rate are associated with peaks from SAA encounters, as shown in Fig. 6. Excluding the SAA region, the dose variations measured by LEO-DOS show good agreement with the >50 MeV proton energy data from GOES-16. This is because the LEO-DOS PD has a self-shielding effect for proton energies below 40 MeV, as shown in Table 1.

jass-41-3-195-g10
Fig. 10. The variations in SEP events from August 5 to August 12, 2023. The upper panel displays the solar proton fluxes measured by GOES-16, while the lower panel presents the dose measured by LEO-DOS PD at 120 s intervals. The recurrent spikes in dose variation, associated to the SAA region as illustrated in Fig. 6. LEO-DOS, low Earth orbit space radiation dosimeter; PD, particle dosimeter; SAA, South Atlantic Anomaly.
Download Original Figure
5.3 Measurement vs Simulation for Neutron Spectrometer (NS)

In the case of NS, the measured absorbed dose (D) in the GCR region was 0.03 μGy h–1. This value was higher than the simulated value of 0.013 μGy h–1. In the MCNP simulation in this study, the absorbed dose and equivalent dose of NS were calculated using the neutron spectrum without shielding at 380 km altitude and 51.6° inclination obtained through simulation by Sato et al. (2006). The lower dose values in the MCNP simulation are due to the spectrum used in the simulation corresponding to a lower altitude and inclination compared to the LEO-DOS mission. Due to the high inclination of the LEO-DOS mission, NS passes through a polar region with low vertical cutoff rigidities, which can cause the generation of more albedo neutrons due to the higher GCR intensity (Dachev et al. 2023). Considering these points, it was concluded that the measured values showed good agreement with the values predicted by the simulation.

5.4 Consideration of Biological Effects

Radiation protection is crucial for understanding how harmful space radiation doses are to the human body during space travel, especially in LEO. Using the tissue equivalent sensor, LEO-DOS measured the absorbed dose (D), quality factor (Q), and dose equivalent (H) of space radiation.

Based on three months of observations with LEO-DOS, the measured biological effectiveness Q values for different radiation sources were 4.36 ± 0.27 for GCR and 1.73 ± 0.25 for SAA. Considering the detector geometry and the satellite’s inclination of 97.8°, these values are reasonably similar to those obtained by Dachev et al. (2017), which were 3.5 for GCR and 1.3 for SAA. Meanwhile, for the PD of LEO-DOS using a 10 mm TEP, the Q value for SEP was 2.82 ± 0.65, which is similar to the value of 2.5 measured on the ISS using a 40 mm TEP (Semkova et al. 2014).

6. SUMMARY and CONCLUSIONS

Since its successful launch on May 25, 2023, LEO-DOS has been operating normally up to 2024. This paper summarizes the results of an initial analysis of the data observed up to September 2023. One of its objectives, the creation of a global dose map for charged particles and neutrons in LEO, has been successfully achieved. Radiation source selection techniques, including GCR, IRB, ORB, and SEP, have been established. The observational results were validated by comparing them with measurements taken in the ISS orbit (Dachev et al. 2015). This has enabled quantitative analysis of the LEO-DOS observation for the space radiation environment study in LEO.

From July 2023 to September 2023, four SEP events were observed. The dose variations due to SEP were consistent with changes in proton energy > 50 MeV as measured by GOES. Neutron dose was successfully measured, with an absorbed dose of 0.03 μGy h–1 and an equivalent dose of 0.3 μSv h–1 in the GCR region. The dose for neutrons was 1/100th of the dose from charged particles.

Currently, research on the biological effects of the near-Earth space environment is being conducted through long-term LEO-DOS observation results and energy deposition analysis on the sensor. In conclusion, the LEO-DOS results will help improve models that simulate biological doses during long-term human space travel and provide accurate radiation risk estimates for future interplanetary manned flights. Additionally, the technical heritage of the LEO-DOS instrument will be applied to the development of the lunar vehicle radiation dosimeter (LVRAD) for the lunar radiation environment.

ACKNOWLEDGEMENTS

This work was supported by the National Research Foundation of Korea (NRF) Grant funded by the Korea government (MSIT) (NRF-2017M1A3A4A01077173 and NRF-2017M1A3A4A01077220). We would like to thank the SaTRec team for their efforts in ensuring the stable operation of LEO-DOS following the successful launch of NEXTSat-2.

REFERENCES

1.

American National Standards Institute [ANSI], Space environment: galactic cosmic ray model, Technical Report, ISO 15390 (2004).

2.

Benton ER, Benton EV, Frank AL, Conversion between different forms of LET. Radiat. Meas. 45, 957-959 (2010).

3.

Charles MW, ICRP Publication 103: recommendations of the ICRP. Radiat. Prot. Dosim. 129, 500-507 (2008).

4.

Dachev TP, Characterization of the near Earth radiation environment by Liulin type spectrometers. Adv. Space. Res. 44, 1441-1449 (2009).

5.

Dachev TP, Bankov NG, Tomov BT, Matviichuk YN, Dimitrov PG, et al., Overview of the ISS radiation environment observed during the ESA EXPOSE-R2 mission in 2014–2016. Space Weather. 15, 1475-1489 (2017).

6.

Dachev TP, Litvak ML, Benton E, Ploc O, Tomov BT, et al., The neutron dose equivalent rate measurements by R3DR/R2 spectrometers on the international space station. Life Sci. Space Res. 39, 43-51 (2023).

7.

Dachev TP, Semkova JV, Tomov BT, Matviichuk YN, Dimitrov PG, et al., Overview of the Liulin type instruments for space radiation measurement and their scientific results. Life Sci. Space Res. 4, 92-114 (2015).

8.

Hansen W, Richter D, Determination of light output function and angle dependent correction for a stilbene crystal scintillation neutron spectrometer. Nucl. Instrum, Methods Phys. Res. A, 476, 195-199 (2002).

9.

Hyun HJ, Anderson T, Angelaszek D, Baek SJ, Copley M, et al., Performances of photodiode detectors for top and bottom counting detectors of ISS-CREAM experiment. Nucl. Instrum, Methods Phys. Res. A, 787, 134-139 (2015).

10.

International Commission on Radiological Protection [ICRP], 1990 Recommendations of the International Commission on Radiological Protection (ICRP Publication, Ottawa, ON, 1991), 1-3.

11.

Jang TS, Park HY, Lee JS, Park MY, Kim JH, et al., Development of Small Radar Satellite NEXTSat-2, in Asia-Pacific International Symposium on Aerospace Technology, Niigata, Japan, 12-14 Oct 2022.

12.

Miller CA, Di Fulvio A, Clarke SD, Pozzi SA, Dual-particle dosemeter based on organic scintillator. Radiat. Prot. Dosim. 191, 319-327 (2020).

13.

National Council on Radiation Protection and Measurements [NCRP], Radiation protection guidance for activities in low-Earth orbit, NCRP Report No. 132 (2000).

14.

Sato T, Niita K, Iwase H, Nakashima H, Yamaguchi Y, et al., Applicability of particle and heavy ion transport code PHITS to the shielding design of spacecrafts, Rad. Meas. 41, 1142-1146 (2006).

15.

Semkova J, Dachev T, Koleva R, Bankov N, Maltchev S, et al., Observation of radiation environment in the International Space Station in 2012-March 2013 by Liulin-5 particle telescope. J. Space Weather Space Clim. 4, A32 (2014).

16.

Wilson JW, Badavi FF, Kim MY, Clowdsley MS, Heinbockel JH, et al., Natural and induced environment in low Earth orbit. NASA TM-2002-211668 (2002).

17.

Werner CJ, Bull JS, Solomon CJ, Brown FB, McKinney GW, et al., MCNP version 6.2 (Los Alamos National Laboratory, Los Alamos, NM, 2017).

18.

Youn S, Nam UW, Kim S, Kim H, Park WK, et al., Calibration and simulation of a silicon dosemeter for ambient dose equivalent in low-Earth orbit space. Radiat. Prot. Dosim. 199, 2118-2125 (2023).