Introduction

Gamma radiation emitted from naturally occurring radioisotopes, such as 40K and the radionuclides from the 232Th and 238U series and their decay products, which exist at trace levels in all ground formations, represents the main external source of irradiation to the human body [1,2,3]. The human beings are exposed mainly to natural sources of radioactivity. More than 80% of the radiation dose received by mankind is due to natural radiation sources [1]. Natural radioactivity is associated mainly to primordial radionuclides, including the isotopes 40K and the progeny of the 238U and 232Th decay series.

The annual worldwide per caput effective radiation dose from natural source is determined by adding the various components. The annual global per caput effective dose due to natural radiation sources is 2.4 mSv. However, the range of individual doses is wide. In any large population about 65% would be expected to have annual effective doses between 1 and 3 mSv, about 25% of the population would have annual effective doses less than 1 mSv and 10% would have annual effective doses greater than 3 mSv [1].

Numerous studies all over the world have been conducted to determine the activity of these radioisotopes in soils and rocks to estimate the external gamma-radiation dose to the public [1]. Several studies have been conducted in the countries of the Asian continent, to determine the activity concentrations of the naturally occurring radionuclides in soil and rock samples [4,5,6,7,8,9,10]. In most of these studies the external radiation dose was calculated. Internal radiation dose from ingestion, tobacco and building materials were also studied [11,12,13,14]. However, internal dose due to inhalation of fine aerosol particles in the surrounded Korean peninsula are not found in the literature.

The Gosan area in Jeju Island, South Korea serves as a natural background site for characterizing the air pollution in Korean peninsula since the atmosphere in the area is relatively very little affected by artificial airborne matters. Chemical composition analyses have usually been conducted for the atmospheric aerosols collected at the site. In this study, we have measured, using the Inductively Coupled Plasma-Dynamic Reaction Cell-Mass Spectrometer (ICP-DRC-MS), radioactivity concentrations of 40K, 232Th and 238U contained in the atmospheric PM10 aerosols which were collected at the Gosan during the year of 2014.

Experimental

Air sampling locations

Air sampling for PM10 aerosols was conducted at the Gosan (33°1′N, 126°10′E), which is located at the western edge of Jeju Island facing the Asian continent (Fig. 1) [15]. Jeju Island, Korea, is located at the boundary of the East China Sea and the Yellow Sea and is surrounded by mainland China, the Korean peninsula, and Japan. The Gosan sampling site is located on the hill of 72 m above sea level and is isolated from residential areas on the island. In order to understand physicochemical and radioactive properties of anthropogenic aerosols under Asian continental outflow, several international experiments have been conducted at the Gosan site, such as ACE-Asia (Aerosol Characterization Experiment-Asia) [16] and ABC-EAREX 2005 (Atmospheric Brown Cloud-East Asia Regional Experiment 2005) [17]. To understand the link between physical and chemical properties of aerosols transported from the Asian continent and regional climate change, we need to investigate the sources and formation mechanism of secondary aerosols. In addition, contributions of local effects on the Gosan site aerosols should be qualitatively and quantitatively evaluated [18].

Fig. 1
figure 1

Location of air sampling Gosan site in Jeju Island, Korea

Collection of atmospheric aerosols

The PM10 samples have been collected using PMS-103 (APM Engineering, PMS-103, Korea), which is an automatic system with Teflon filter (Pall Corporation, ZeflourTM, 47 mm/2.0 µm, USA) at 24 h basis with every 3 day intervals during the year of 2014. The air flow rate was kept to about 16.7 L/min, and total air flow was calculated from the flow rate and running time with MFC (mass flow controller). Quality assurance and quality control procedures were followed, including calibration of the sampler flow rates, collecting field and trip blanks, inspection and conditioning of the filter media before and after use.

Sample analysis

For obtaining accurate concentrations of the radionuclides of interest, complete digestion is required before ICP-MS. The addition of hydrofluoric acid strongly influences the recovery of the microwave acid digestion of environmental samples because it breaks down silicates and minerals better than other acid combinations [19]. However, due to the problems in glassware and torch damage of ICP-MS as well as toxicity of HF, the acid combination HNO3/HCl was used for the digestion process [19].

The PM10 aerosols were decomposed with acids using a START D microwave digestion system (Milestone, Italy). The aerosol filters were put in a perfluoroalkoxy Teflon vessel of the microwave digestion system together with 10 mL acid solution (5.55% HNO3/16.75% HCl). This vessel was heated at 180 °C for 15 min with 1000 W microwave radiation to digest the PM10 aerosols. The decomposed solution was transferred through 0.45 μm PVDF syringe filter (Whatman) into a 25 mL volumetric flask, and the volume was adjusted to 25 mL with 5 mL acid solution (3% HNO3/8% HCl) and ultrapure water [20]. The number of elements determined by ICP-DRC-MS (Perkin Elmer, ELAN DRC-II, USA) instruments was 3 species such as 39K, 232Th and 238U. The instrumental conditions of ICP-DRC-MS were as follows: 40 MHz RF frequency, 1500 W RF power, 15.0 L min−1 coolant, 0.9–1.05 L min−1 Ar flow, 1.2 L min−1 auxiliary. The DRC gas was used as the ammonia. The instrumental detection limits (IDL) of 39K, 232Th and 238U were 27.12 μg L−1, 2.11 ng L−1 and 2.89 ng L−1, respectively (Table 1).

Table 1 Instrumental detection limit (IDL) and coefficient of variation (CV) for ICP-DRC-MS (n = 7)

Isotope activity concentration

Potassium, thorium and uranium were measured as trace elements in PM10 aerosol by ICP-DRC-MS. The radioactivity concentrations of the isotopes 40K, 232Th and 238U in the collected PM10 were calculated using the following Eq. [1] [21]:

$$A_{\text{i}} = \frac{\ln 2}{{T_{1/2} }} \times \frac{{\rho_{\text{i}} \times m_{\text{e}} }}{{M_{\text{i}} }} \times N$$
(1)

where Ai, T1/2, ρi, me, Mi, and N are the radioactivity concentration (Bq m−3), half-life time (s) of isotope i (40K: 1.28 × 109 years, 238U: 4.468 × 109 years, 232Th: 1.405 × 1010 years), isotopic ratio (natural abundance) of isotope i, the mass concentration of element e corresponding to isotope i (g m−3), atomic mass (g mol−1), and the Avogadro’s number (6.022 × 1023 mol−1), respectively. The isotopic ratios for 40K, 238U and 232Th were 0.000117, 0.99275 and 1.0, respectively, based on the International Union of Pure and Applied Chemistry (IUPAC) report 2009 [22].

Inhalation annual effective dose

The inhalation annual effective radiation dose (Ei) due to PM10 was calculated using the following Eq. [2] [1]:

$$E_{\text{i}} = A_{\text{i}} \times B \times d_{\text{i}} \left( {1 \, - F_{0} + F_{0} F_{\text{r}} } \right)$$
(2)

where Ai is the integrated activity concentration of radionuclide i associated with PM10 in outdoor air (Bq m−3), B is the breathing rate (m3 year−1), di is the committed dose per unit intake from inhalation or effective dose coefficient (Sv Bq−1), F0 is the indoor occupancy factor and Fr is the ratio of indoor to outdoor air concentration.

The annual effective radiation dose due to inhalation of PM10 was calculated for the six age groups identified by the International Commission on Radiological Protection (ICRP), namely 3 months, 1, 5, 10, 15 years and adults. For each of the age groups, the corresponding breathing rate and indoor occupancy factor [23], and effective dose coefficient for inhalation [24] were used. These values are summarized in Tables 2 and 3. The default modes of absorption for the isotopes were used. For the effective dose coefficient, the activity median aerodynamic diameter (AMAD) was assumed to be 1 μm. The ratio of indoor to outdoor air concentration (Fr) was assumed to be 0.3 [1, 25, 26].

Table 2 Breathing rate and indoor occupancy factor values used for calculating the inhalation annual effective radiation dose due to inhalation of PM10
Table 3 Effective dose coefficient for inhalation di (μSv Bq−1) used for calculating the inhalation annual effective radiation dose due to inhalation of PM10

Results and discussion

40K, 232Th and 238U radioactivity concentration in PM10

As an alternative method, we have measured, using the ICP-DRC-MS, radioactivity concentrations of 40K, 232Th and 238U contained in the atmospheric PM10 aerosols which were collected at the Gosan during the year of 2014. A total of 115 samples have been analyzed, of which 5 samples are those collected during Asian Dust days, 47 samples collected during normal weather days (no-event days), and the remaining samples collected during the days of haze (6 samples) and fog-mist (57 samples).

The mean mass concentration of PM10 was 47.31 μg/m3. During the study period, the mean concentrations of 40K, 232Th and 238U were 0.56, 1.02 and 0.53 mg/kg-dust, respectively. Furthermore, during the study period, the activity concentrations of the radionuclides 232Th and 238U were 0–16.86 (mean 4.16 ± 2.98 Bq/kg-dust) and 0–62.07 (mean 6.48 ± 7.43 Bq/kg-dust), respectively. They are lower than the world averages for activity concentrations of 232Th and 238U, which are 30 and 35 Bq/kg-dust, respectively [27]. The mean active concentrations of 40K, 232Th and 238U during normal days are 7.89 ± 10.89, 0.25 ± 0.37 and 0.30 ± 0.35 μBq/m3, respectively. Furthermore, they are lower than the world averages for the mean atmospheric activity concentrations of 232Th and 238U associated with dust, which are 0.5 and 1.0 μBq/m3, respectively [27]. It is worth noting that the world average is based on total suspended particulate. The low atmospheric activity concentration of 232Th and 238U might be attributed to the low atmospheric PM10 in Gosan site as compared to many of the world cities. The 232Th/238U activity concentration ratio of PM10 was 0.830.

Activity concentrations by atmospheric phenomenon

The radioactivity concentrations of those isotopes were analyzed by atmospheric phenomenon (Asian Dust, Haze, Fog-Mist and Non-Event). During Asian Dust periods (5 samples), the mean concentrations of 40K, 232Th and 238U were 1.24, 2.81 and 0.87 mg/kg-dust, respectively. Furthermore, during the Asian Dust period, the activity concentrations of the radionuclides 232Th and 238U were 7.44–16.86 (mean 11.38 ± 3.83 Bq//kg-dust) and 9.19–13.45 (mean 10.71 ± 1.72 Bq/kg-dust), respectively. Table 4 shows also that the mean atmospheric activity concentration of 232Th and 238U of Asian Dust periods associated with PM10 was 1.47 ± 0.72 and 1.37 ± 0.55 μBq/m3, exceeding the revised world reference values. And these were highly as 6.94, 8.57 and 7.05 times, respectively, compared to the non-event periods (47 samples). The 232Th/238U ratio of Asian Dust was 1.073, which was higher than those of other atmospheric phenomenon (0.707–0.882) (Table 5).

Table 4 The associated 40K, 238U and 232Th concentrations of airborne PM10 by atmospheric phenomenon at Gosan site in Jeju Island, 2014
Table 5 The activity ratio of 40K, 232Th and 238U activity concentrations (μBq/m3) by atmospheric phenomenon at Gosan site in Jeju Island, 2014

Inflow pathways of air mass

We are analyzing five-day backward trajectories of the air inflow into the Gosan site using the HYSPLIT4 (HYbrid Single Particle Lagrangian Integrated Trajectory) model of National Oceanic and Atmospheric Administration (NOAA, USA) [28]. The frequency of air inflow into the Gosan site is 34.84% from Sector I (China continent), 22.58% from Sector II (Korean peninsula), 11.61% from Sector III (Japan) and 5.16% from Sector IV (North Pacific Ocean) during normal days. As a result of analyzing the air volume, the concentration ratio of 232Th/40K was 0.033 and 238U/40K was 0.036 in Sector 1 (Table 6). The concentration ratio of 232Th/238U in this sector was measured as 0.902. For sector 2, 238U/40K was 0.041 and 232Th/40K was 0.031. And 232Th/238U in this sector was measured as 0.750. During Asian dust days, the air inflow is dominated from Sector I (Fig. 2).

Table 6 Activity concentrations and ratio of nuclides according to the sector
Fig. 2
figure 2

5-Day backward trajectories at the Gosan Site of Jeju Island during the study period

Inhalation annual effective dose

Isotopes activity concentrations in air (μBq/m3) presented in Table 4 were used in the calculation of inhalation annual effective radiation dose (Ei) to various age groups of the public. Table 7 and Fig. 3 show the results of these calculations. For instance, ETotal was 17.56 nSv/year to the infants (< 1 year), while it was 78.45 nSv/year to the adults (male), depending on the age group. Because of the variations in air breathing rate, the total annual dose due to the natural radioactivity in airborne PM10 increases for the older age groups. These values are higher than the corresponding values calculated from the data presented in UNSCEAR 2000 [27]. It is obvious that 232Th is the main contributor to the inhalation annual effective dose. 232Th was found to be responsible for the total dose. On the other hand, 40K was found to slightly contribute to the total dose.

Table 7 The average inhalation annual effective dose (nSv/year) to various age groups (default mode F) of the public in Gosan site from 40K, 238U and 232Th associated with airborne PM10
Fig. 3
figure 3

Annual effective dose to various age groups in Jeju Island from inhalation 40K, 238U, and 232Th in PM10 aerosols (F type, Male)

Conclusions

The atmospheric PM10 aerosols (115 samples) were collected at Gosan of Jeju Island, which is one of the natural background sites of Korea, during the year of 2014. This study analyzed using ICP-DRC-MS the concentrations of potassium, thorium and uranium, and evaluated the annual effective dose by breathing from the results. The mean mass concentration of PM10 was 47.308 μg/m3. The mean radioactive concentrations of 40K, 232Th and 238U were 7.89, 0.25 and 0.30 μBq/m3, respectively. The 232Th/238U activity concentration ratio of PM10 was 0.830. The 232Th/238U ratio during Asian Dust days is 1.073, which is higher than those in other atmospheric conditions. During Asian Dust days, the air inflow is dominated from Sector I (China continent). The concentration ratio of 232Th/238U was 0.902 in Sector 1. In this study, the ratio of each nuclide was compared according to the inflow route. As a result, it was confirmed that the ratio of each nuclide was slightly different according to the inflow route. It is expected that this will be a preliminary data for analyzing the sources and various weather phenomena of fine dust on the peninsula. Jeju Island with less pollution source and low population density is also one of the best places as a background area in Asia. It is judged that the results become a preliminary data on the impact of fine dust from China, which has been recently intensified, on the Korean Peninsula. Therefore, it is judged as necessary to observe such Asian Dust phenomenon in the long term through collecting a large number of the PM10 aerosol filters.