Introduction

Because of the tsunami caused by The Great East Japan Earthquake, TEPCO's Fukushima Daichi Nuclear Power Station (FDNPS) lost all electric power, followed by the severe accidents at Units 1, 2, and 3. The accident progression has been analyzed by many researches with emphasis on the thermal-hydraulics in the OECD/NEA projects such as Benchmark Study of the Accident at the Fukushima Daiichi Nuclear Power Station (BSAF)1,2,3,4,5,6,7,8 and Analysis of Information from Reactor Buildings and Contaminant Vessels of Fukushima Daiichi Nuclear Power Station (ARC-F). In addition, fission products (FPs) release, distribution of FPs in all over the plant, and source term to environment are also the topics treated in the projects. However, there are still many uncertainties such as a type of leakage, event timings, and chemical species of FPs due to lack of complete data to input calculation codes.

In order to reduce uncertainties concerning assessments of severe accidents, VERDON tests were performed at the CEA Cadarache Centre in which high burn-up UO2 and MOX fuels were heated in a furnace under various conditions to investigate FPs release and transportation9,10,11,12,13. The results indicated a rapid release kinetics for volatile (I, Cs, Te) and even for semi-volatile (Mo, Ba) FPs during oxidizing condition and re-volatilization of iodine during air injection on the contrary to Cs. Transportation of FPs are affected by not only atmosphere and temperature conditions but also chemical species. While many source term studies have assumed Cs chemical form as CsI and CsOH, recent studies pointed out possible formation of Cs2MoO4 5,6,7,14,15,16,17,18, CsBO2 19, Cs2CrO4 19,20, Cs2Si4O9 19, and CsAlSiO4 21. The formation of these chemical species would be related to the leakage pathways since some reports indicated the production of Cs2Si4O9 and CsAlSiO4 via chemisorption of Cs onto stainless steel and reaction with concrete22,23. Furthermore, the production of Cs2Si4O9 and CsAlSiO4 is possibly related to the Cs-bearing microparticle found in environmental samples after the FDNPS accident24,25,26.

There are many reports related to the distribution of 134,135,136,137Cs, 129,131I, 90Sr, 110mAg and 132Te in environment caused by the accident at the FDNPS27,28,29,30,31,32,33,34. The measured radionuclides concentrations were compared with those calculated by atmospheric transport and deposition model (ATDM) to estimate the atmospheric release of radionuclides during the FDNPS accident27. The detailed analysis of 134Cs/137Cs radioactivity ratio of the deposition indicated which reactor unit ultimately contaminated a specific area: while the contamination of most continental areas was caused by Units 2 and 3, the contamination of the northwest part near FDNPS was caused by Unit 1 because the 134Cs/137Cs radioactivity ratio was slightly lower than those in the other parts27,28,29. Furthermore, slight difference in the 135Cs/137Cs isotope ratio was reported for the area mainly contaminated by Unit 131. Source term information is important to improve such analysis.

Recently, radiochemical analysis of samples collected inside of FDNPS buildings has been planned to improve severe accident analysis and source term. In the present study, drain water sampled at the drain pit of the exhaust stack shared by Units 1 and 2 was analyzed to obtain the information about components and chemical species of radionuclides released by vent of Unit 1. Cesium-134,137 and 131I are the important FPs of severe accident analysis because of their relatively high volatile and radioactivity. However, 131I was decayed out due to its short half-life (8.0 d). Therefore, 129I which has long half-life (1.57 × 107 y) was measured to estimate the behavior of 131I.32,33 In addition, Mo isotopes and 99Tc, daughter nuclide of 99Mo, were analyzed because Mo was prone to be released as Cs2MoO4 when the reactor vessel was under oxidizing condition with sufficient water vapor. Furthermore, lower volatile FPs (90Sr and 152,154Eu) and actinide (Th, U, Pu, and Am) were analyzed to investigate the release of them as aerosols generated in reactor vessel and Molten Core Concrete Interaction.

Materials and analysis method

Drain water

The sampling point was illustrated in Fig. 1. In the vent operation, the gas in the contaminant vessel was passed through the suppression chamber to reduce the radioactivity, and then, released into environment via stand by gas treatment system pipe and exhaust stack. The exhaust stack was shared by Units 1 and 2. It is considered that the vent of Unit 2 was failed since outflow side of the rupture disk in the vent line of Unit 2 was not contaminated, and the inside of the exhaust stack was mainly contaminated by vent of Unit 1. It is presumed that dew condensation water and rainwater in the stack were accumulated in the drain pit. Drain water was sampled in September 12, 2016. Salinity and pH of the drain water were preliminary measured using salinometer and pH test paper, and they were 0.01% and neutrality.

Figure 1
figure 1

Diagrammatic illustration of the sampling point.

Full size image

Analysis of γ-emitting nuclides

γ-Ray spectra were measured by a Ge semiconductor detector (GC1020S-CJT-F-1/S-2002C, Mirion Technologies Canberra). Energy and efficiency of the γ-ray spectrum were calibrated using a mixed standard source of 60Co, 133Ba, and 137Cs, which was prepared as a filter paper (Amersham plc, UK). Therefore, the sample for quantitative analysis of drain water was prepared by transfusing 30 Î¼l of sample solution into a filter paper to be the same sample form. The filter paper was doubly sealed in plastic bags. First, measurement of 134Cs and 137Cs was conducted on the drain water as it was. After that, the measurement of the nuclides excluding 134Cs and 137Cs, such as 60Co, 106Rh (daughter nuclide of 106Ru), 125Sb, and 152Eu, was conducted after Cs was eliminated using ammonium phosphomolybdate (AMP) because γ-rays from the other nuclides were considered to be hidden by the higher level of radioactivity of 134Cs and 137Cs in the drain water. The amount of AMP, HNO3 concentration and standing time were optimized to eliminate Cs and to remain the other elements before analysis of the drain water. Sample solutions for the optimization were prepared by diluting the standard solutions of B, Na, Mg, Al, Ca, Mn, Fe, Co, Ni, Se, Rb, Sr, Y, Zr, Mo, Ru, Rh, Ag, Sn, Sb, Te, Cs, Ba, Eu, Re (1 mg/ml solutions purchased from Fujifirm Wako Chemical CO.INC., Kanto Chemical CO.,INC., and Merck) and of Nb (0.1 mg/ml solution purchased from Accu Standard). Based on the optimization, 5 Î¼l of standard solutions of Co, Se, Sr, Ru, Rh, Sb, Cs, Eu and Re (1 mg/ml) were added to 5 ml of drain water and HNO3 concentration was adjusted to 2 M. The drain water sample solution was stood for overnight. After the addition of 0.02 g of AMP to the drain water sample solution, it was stirred for 1 h. AMP was filtered using a syringe and a membrane filter unit of which pore size was 0.45 Î¼m. The filtrate was assigned as Cs-eliminated drain water, and quantitative analysis of γ-ray emitting nuclides was conducted. Furthermore, γ-ray spectra of the Cs-eliminated drain water in a bottle (approximately 5.7 ml) were recorded for the qualitative analysis.

Analysis of 90Sr

Five μl of Cs standard solution was added to the Cs-eliminated drain water, and then, the Cs elimination procedure using AMP was repeated one more time with the same manner of the first time. The obtained filtrate was assigned as the twice Cs-eliminated drain water. 0.1 ml of the twice Cs-eliminated drain water was taken into a bottle, and HNO3 and pure water were added to prepare 2 ml of 8 M HNO3 sample solution. This was assigned as DW-Sr sample.

0.5 ml of the Sr resin (Eichrom Technology, LLC) swollen with ultrapure water was packed into a column and conditioned by passing through 5 ml of 8 M HNO3. Next, the DW-Sr sample was loaded into the column. The bottle which had contained DW-Sr sample was rinsed with 0.5 ml of 8 M HNO3 4 times and all rinse solutions were loaded into the column. After that, 2.5 ml of 8 M HNO3 was passed through the column 4 times. The extracted Sr was recovered by passing through 3 ml of 0.01 M HNO3. 0.2 ml of the recovery solution was sampled and mixed with 1.8 ml of 1 M HNO3 to measure Sr concentration with inductively coupled plasma mass spectrometer, ICP-MS, (Agilent 7700, Agilent Technologies Japan, Lts), and the recovery rate during chemical separation was evaluated. Ultima Gold LLT (PerkinElmer) was added to the remaining recovery solution, and β-ray spectrum was measured with liquid scintillation counter (PerkinRlmer Tri-Carb 3110TR) after reaching radioactive equilibrium between 90Sr and 90Y. Efficiency was corrected by a calibration curve obtained from the measurement of standard samples containing 20 Bq of 90Sr and 90Y, various concentration of HNO3 and Ultima Gold LLT.

Analysis of 99Tc

Aliquot of the twice Cs-eliminated drain water was adjusted to approximately 5 ml of 1 M HNO3 solution (DW-Tc). 0.5 ml of the TEVA resin (Eichrom Technology, LLC) swollen with ultrapure water was packed into a column and conditioned with 1 M HNO3. The DW-Tc was loaded into the TEVA resin column. The bottle which had contained DW-Tc sample was rinsed with ultrapure water 3 times and all rinse solutions were loaded into the column. After that, 2.5 ml of ultrapure water was passed through the column 4 times. The extracted Tc was recovered by passing through 1.5 ml of 8 M HNO3. The recovery solution was appropriately diluted to measure Re with ICP-MS. The remaining recovery solution was also appropriately diluted and mixed with Ultima Gold LLT to measure β-ray of 99Tc with liquid scintillation counter (PerkinElmer Tri-Carb 3110 TR).

Analysis of Mo isotopes

Molybdenum isotopes (mass number 92, 94, 95, 96, 97, 98 and 100) in a Mo standard solution (0.5 ng/mL, this solution was prepared by diluting a standard solution of 1 mg/ml of Mo, Fujifilm Wako Chemical Co. INC.) was measured using ICP-MS (Agilent 7700) and Sector field ICP-MS (Element 2, Thermo Fisher Scientific Co.) to study the effect of mass discrimination on the isotope ratios obtained by ICP-MS measurement. Reference values of the isotope ratio to 98Mo, R true, were divided by those of measured, R meas, and plotted against the mass difference of the isotope of interest from 98Mo, Δm.

Previously, separation method for Mo using a TEVA resin was developed to purify Mo from rubble samples35. Because the dissolved solution of the rubble sample contained large amount of Fe3+ and Cr6+, which compete with Mo extraction on the TEVA resin, ascorbic acid was added to the sample solution to decrease their extraction by reducing to Fe2+ and Cr3+. In addition, HF was added to the sample solution to avoid precipitation of Nb. However, influence of Fe, Cr, and Nb in the drain water is negligible due to lower concentration of them. Therefore, sample solution of the drain water to load the TEVA resin was adjusted to 4 M HCl using ultrapure grade HCl solution (31%, Kanto Chemical Co.,INC.), and the solution was assigned as DW-Mo sample. 0.5 ml of the TEVA resin swollen with ultrapure water was packed into a column, followed by washing with 5 ml of 1 M HNO3, 3.4 ml of ultrapure water, and 2 ml of 4 M HCl. Then, 5 ml of the DW-Mo sample was loaded into the column. The bottle which had contained DW-Mo sample was rinsed with 0.5 ml of 4 M HCl 3 times and all rinse solutions were loaded into the column. Next, the column was rinsed with 3.4 ml of 4 M HCl and 0.5 ml of ultrapure water. Finally, extracted Mo was recovered by passing through 2.5 ml of 1 M HNO3. As an operation blank solution, 4 M HCl solution was prepared and the same separation operation was carried out. Recovered Mo isotopes were measured with sector field ICP-MS. Net count rates of Mo isotopes for drain water were obtained by subtracting the measured count rates of Mo isotopes for the operation blank solution from those for the drain water. Whereas, net count rates of Mo isotopes for Mo standard solution, 0.5 ng/ml, were obtained by subtracting the measured count rates for 1 M HNO3 solution from those for the Mo standard solution.

Analysis of 129I and 129IO3

It was reported that I was extracted on the solid extraction disk, Anion-SR (Empore disk, 3 M Co.), from 3 M NaOH solution, and I and IO3 were extracted on the Anion-SR from 0.01 M HCl containing 0.1 M NaHSO3 because IO3 was reduced to I by NaHSO3 in diluted HCl solution36,37. In the present study, 0.5 ml of drain water was adjusted to 5 ml of 2.5 M NaOH solution to analyze I and to 5 ml of 0.01 M HCl solution containing 0.01 M K2S2O5 to analyze the sum of I and IO3 . Instead of NaHSO3, K2S2O5 was used because NaHSO3 is unstable in a solid state, and formula number of Na2SO3, 126, is close to mass number of stable iodine, 127. As spike samples, a known amount (approximately 50 pg, "a" in the Eqs. (1) and (2)) of 127I or 127IO3 was added to the drain water in the same solution conditions. In addition, 2.5 M NaOH solution and 0.01 M HCl solution containing 0.01 M K2S2O5 were also prepared as blank samples. The I or the sum of I and IO3 were purified by Anion-SR using the reported method36,37. Briefly, Anion-SR was sequentially conditioned with 10 ml of acetone, 10 ml of methanol, 10 ml of pure water, 10 ml of 1 M HNO3 and 30 ml of pure water. Next, the sample solution was loaded to the Anion-SR, and then, it was rinsed with 5 ml of pure water 6 times. Finally, extracted I was recovered 9.5 ml of 1 M HNO3. A known amount of Cs was added into the purified I or (I + IO3 ) solutions as an internal standard and a known amount of NaClO was added to uniform the chemical species to IO3 and to stabilize iodine in 1 M HNO3 solution. The concentrations of 127I ("x" in Eq. (1)) and 129I ("y" in Eq. (2)) were measured by using ICP-MS (Agilent 7700). The amounts of 127I and 129I were calculated by using Eqs. (1) and (2), respectively.

$${\text{x}} = \left\{ {\left( {{\text{A}}/{\text{B}}} \right) \times {\text{a}}} \right\}/\left\{ {\left( {{\text{A}}^{\prime } /{\text{B}}^{\prime } } \right) - \left( {{\text{A}}/{\text{B}}} \right)} \right\}$$

(1)

$${\text{y}} = {\text{a}}/\left\{ {\left( {{\text{A}}^{\prime } /{\text{B}}^{\prime } } \right) - \left( {{\text{A}}/{\text{B}}} \right)} \right\}$$

(2)

where A and B were the measured amounts of 127I and 129I in a sample solution, respectively, A' and B' were the measured amounts of 127I and 129I in the spike sample solutions, respectively, and a was spiked amount of 127I. The amounts of x and y were obtained for chemical species of I and the sum of I and IO3 , and the amount of IO3 was calculated by subtracting the amount of I from that of the sum of I and IO3 .

Measurement of stable isotopes, Th and U in drain water

The drain water was diluted 200 times with 1 M HNO3 and the concentrations of 59Co, 77,82Se, 88Sr, 99Tc, 101Ru, 103Rh, 121Sb, 133Cs, 153Eu, 185Re, 232Th and 238U in the drain water were determined using ICP-MS (Agilent 7700) by calibration curve method.

Measurement of Pu and Am

One (1) ml of the drain water was adjusted to 2 ml of 1 M HNO3 solution containing 0.1 M NaNO2 and assigned as DW-PuAm sample. Pulutonium-242- or 243Am-spiked drain water sample was prepared as the same manner and assigned as Pu-DW-PuAm and Am-DW-PuAm samples, respectively. 0.5 ml of TRU resin (Eichrom Technology, LLC) swollen with ultrapure water was packed into a column and conditioned with 1 M HNO3. Then, the sample solution was loaded into the column. The bottle which had contained the sample was rinsed with 1 M HNO3 3 times and all rinse solutions were loaded to the column. The column was rinsed with 3.4 ml of 1 M HNO3. Extracted Am was recovered with 3.4 ml of 4 M HCl, and next, extracted Pu was recovered with 3.4 ml of 4 M HCl containing 0.02 M TiCl3. 0.3 ml of Sm solution (1 mg/ml) and 1 ml of HF were added to the recovery solutions of Am and Pu to precipitate AmF3 and PuF3 with SmF3. The micro precipitate was collected on PTFE membrane filter and dried with a little amount of ethanol. The filter was fixed with double-stick tape on a stainless-steel dish to subject to the measurement by α-ray spectrometer (α-Ray Detector BU-020-450-AS0, α-Ray Module SOLOIST,SEIKO EG&G ORTEC). The recovery rate of the micro precipitate was normalized with the recovery of 147Sm, which is a natural alpha-emitting radioisotope (half-life: 1.06 × 1011 y, alpha energy: 2.2476 MeV).

Results

Analysis of γ-emitting nuclides

Before the Cs elimination from the drain water using AMP, the solution condition and AMP amount were optimized to remove Cs so that the other elements remained (Supporting Table 1 and 2), 99.9% of Cs was removed with 0.01 g of AMP from 10 ml of 1 M HNO3 containing 1 Î¼g/ml of B, Na, Mg, Al, Ca, Mn, Fe, Co, Ni, Se, Rb, Sr, Y, Zr, Mo, Ru, Rh, Ag, Sn, Sb, Te, Cs, Ba, Eu and Re, and 0.1 Î¼g/ml of Nb. Most elements remained in the solution although the rates of remaining Rb and Ag were low and those of Zr, and Nb decreased with increasing AMP amount. The concentration of Mo in the solution was increased to dozens of μg/ml because Mo in the AMP was dissolved to the solution. Remaining rates of metal ions were also affected by HNO3 concentration. Increment of HNO3 concentration tends to enhance remaining rates of metal ions. Although Zr and Nb hardly remained in 0.1 M HNO3, 94% and 67% of Zr and Nb, respectively, remained in 2 M HNO3. Therefore, 2 M HNO3 condition was selected so that larger amount of Zr and Nb remain in the solution. Furthermore, the rate of remaining Sb depended on the standing time after preparation of sample solution. Although Sb existed as SbCl3 in the standard solution, it could be changed to Sb5+ in the oxidizing solution, 2 M HNO3, to remain in the filtrate. It seemed that a part of Sb3+ was captured by AMP but Sb5+ remained in the solution. Therefore, drain water was stood for overnight after adding acid and carrier to wait for the oxidization of Sb3+ to Sb5+.

Figure 2 shows the γ-ray spectra for (a) quantitative analysis of drain water, (b) quantitative analysis of Cs-eliminated drain water, and (c) qualitative analysis of Cs-eliminated drain water. Only peaks of 134Cs and 137Cs were observed in spectrum (a) except those of natural background nuclides and backscattering. The concentrations of 134Cs and 137Cs were 2.4 × 103 and 4.8 × 104 Bq/ml, respectively. When the radioactivity was decay-corrected to March 11, 2011, they were (5.1 ± 0.02) × 104 and (5.9 ± 0.007) × 104 Bq/ml, respectively, and the 134Cs/137Cs radioactivity ratio was 0.86 ± 0.01.

Figure 2
figure 2

Gamma-ray spectra of drain water and Cs eliminated drain water. Sample solution (30 μl) was transfused into filter paper for quantitative analysis, (a) and (b). Sample solution in a bottle was measured for quantitative analysis, (c).

Full size image

The concentration of 137Cs was decreased to 63 Bq/ml for the Cs-eliminated drain water, and peaks of 125Sb showed up. The concentration of 125Sb in the drain water was 14 Bq/ml. It is reported that Sb was released by the oxidation of LWR fuel cladding by steam38. While not only 125Sb but also 60Co were observed for the qualitative analysis of the Cs-eliminated drain water, the concentration of 60Co were less than its detection limit for quantitative analysis, 0.7 Bq/ml, due to a small amount of sample, 30 Î¼l.

Analysis of β-emitting nuclides

Figure 3 shows the liquid scintillation spectrum of the Sr measurement sample separated from DW-Sr sample. The spectrum consists of β-rays emitted from 90Sr and 90Y. The maximum energies were shifted to lower energy side because of chemical quenching due to HNO3. Recovery rate of Sr for chemical separation was (92 ± 8)%. The obtained 90Sr concentration was (47 ± 4) Bq/ml at the measurement time, and (55 ± 5) Bq/ml at the time point of the accident at the FDNPS.

Figure 3
figure 3

Liquid scintillation spectrum of 90Sr-90Y for the drain water after separation and reached radioactive equilibrium.

Full size image

There is no stable Tc isotope. Because Re and Tc tend to form ReO4 and TcO4 , respectively, and show similar behavior, Re was used as a carrier of Tc. The recovery rate of Re after all chemical separation processes was 99.5%, and this result supports that the chemical separation procedure is suitable for the 99Tc separation. However, the concentration of 99Tc in the drain water was less than detection limit, 0.012 Bq/ml.

Analysis of iodine species

In this study, organic and molecular iodine were not measured because they seemed to be vaporized during long period until sampling and analysis on September 12th, 2016, and on November 17-18th, 2020, respectively. Furthermore, it was reported that the largest fraction of the radioiodine present in BWR condensate was found to be I 39, and molecular I2 was hydrolyzed and disproportionated in water to produce I and IO3 40. Therefore, these species were considered dominant species and measured. Table 1 shows the concentrations of iodine species in the drain water. As the total inorganic iodine, 20 ± 3 and 21 ± 2 ng/ml of 127I and 129I, respectively, were obtained and the values were the same order of the reported 129I concentration in the accumulated water before treatment36. It was estimated that the amounts of 127I and 129I in Unit 1 core were 2.14 × 103 and 9.44 × 103 g/core, respectively41, whereas reported 127I concentrations in the rainwater collected from Fukushima prefecture were a few ng/ml level42. Therefore, 127I in the drain water is considered a mixture of Unit 1-originated 127I and rainwater-originated 127I. Although the chemical form might be changed until the analysis and measurement, most iodine was observed as I in the present time; approximately 90% of iodine was I and 10% was IO3 in both isotopes. In case of the seawater collected at the offshore Fukushima on June 3-17th, 2011, the I/IO3 ratio for 129I different from that for nature-originated 127I was observed; the 129I originated from FDNPS existed mainly as I form but nature-originated 127I mainly as IO3 form43. However, our result shows similar I/IO3 ratio for both isotopes, indicating isotope exchange reaction reached equilibrium.

Table 1 Concentrations and percentages for each iodine chemical species.

Full size table

Table 2 Comparison of Mo isotope ratios for natural and drain water, and Unit 1 originated one and calculated on using ORIGEN2 code.

Full size table

Analysis of Mo isotopes

Stable Mo isotopes of which mass number 95, 97, 98, 100, have high fission yield, 6.497, 6.045 5.701, and 6.579%, respectively44. Furthermore, it was reported that Mo and Cs showed similar distribution within the corium after heating of MOX fuel at 2620℃ with hydrogen and steam condition, suggesting the formation of Cs2MoO4 45. Therefore, Mo is an important element to consider chemical form of radiocesium.

Mass discrimination effect was evaluated to measure Mo isotope ratios precisely. Figure 4 shows the relationship between (R true/R meas) ratio and Δm for the measurement of Mo standard solution using Agilent 7700 and Element 2 with relational expression. Noticeable mass discrimination effect was observed for the measurement using Agilent 7700, suggesting that the isotope ratios should be corrected by the mass discrimination factor obtained from the slope of the relational expression. On the other hand, mass discrimination effect for Element 2 was relatively small, and their dispersion was less than 1.5%. Therefore, Mo isotopes in the drain water sample was measured by Element 2 and the obtained isotope ratio was used without correction in this study.

Figure 4
figure 4

Relationship between mass difference (Δm) and ration of true isotope ratio to measured value (R true/ R meas).

Full size image

Table 2 shows the Mo isotope ratios in the Mo standard solution and the drain water after separation using the TEVA resin. In addition, the isotope ratios of Mo generated in the Unit 1 were estimated by ORIGEN code41,44. Isotope ratios for the Mo standard solution after separation using TEVA resin agreed with those of natural isotope ratios, indicating isotope effect is negligible during the chemical separation, whereas Mo isotope ratios for drain water was different from the natural isotope ratios and those originated from Unit 1. It is considered that nature-originated Mo and Unit 1-originated Mo were mixed in the drain water. Molybdenum-92 in the drain water is considered nature-originated isotope because the fission yield of 92Mo is approximately 10 order magnitude lower than those of 95,97,98,100Mo44. Then, count rates of nature-originated Mo isotopes in the drain water were calculated from the count rate of 92Mo and natural isotope ratio. Next, count rates of Mo isotopes originated from Unit 1 were calculated by subtracting the calculated count rates of nature-originated Mo in the drain water from the measured count rates of Mo isotopes for the drain water. Molybdenum isotope ratios originated from Unit 1 in the drain water were calculated and listed in Table 2. They were roughly agreed with those calculated by ORIGEN code, indicating that the drain water contains Unit 1-originated Mo isotopes. This is the first measurement result of the stable Mo isotopes released from FDNPS. The concentrations of the total Unit 1-originated stable Mo and nature-originated Mo in the drain water were 0.13 and 0.82 Î¼g/L, respectively. Considering that it is reported that sub-μg/L quantities of Mo are expected in rainwater46, main source of nature-originated Mo in the drain water is probably rainwater.

Analysis of Th, U, Pu and Am

Thorium and U were not detected. The limit of detection of ICP-MS measurement for Th and U were 0.0006 and 0.0013 ng/mL, respectively, and this means that the concentrations of Th and U in the drain water were less than 0.12 and 0.25 ng/mL, respectively, because the drain water was diluted 200 times for the measurement. Also, Pu and Am were not detected in the drain water. Here, the recovery rates for Pu and Am obtained from the analytical results of 242Pu or 243Am-spiked drain water samples were 93% and 101%, respectively, and the limits of detection for 242Pu and 243Am were both 0.002 Bq.

Discussion

Mole ratios of determined isotopes

The measured radioactivity concentrations of target nuclides in the drain water were decay-corrected to the time point of the accident at FDNPS, March 11th, 2011, and then, the radioactivity ratios of the target nuclides to 137Cs were calculated (see Table 3).

Table 3 Radioactivity and concentrations of target isotopes in the drain water and ratios of them to 137Cs for the measured and calculated for Unit 1 core.

Full size table

TEPCO and some researchers reported radioactivity concentrations of 134Cs and 137Cs in stagnant water in the FDNPS47,48,49. The 134Cs/137Cs radioactivity ratios for stagnant water of Unit 1 decay-corrected to March 11th, 2011 were varied from 0.89 to 0.97 and the averages reported by Komori et al. and Nishihara et al. were 0.93 ± 0.03 and 0.92 ± 0.02, respectively. Furthermore, the radioactivity ratio reported by Komori et al. was 0.91 ± 0.2 for the stagnant water sampled before October, 22th 2011 to avoid mixture of stagnant water from Unit 2 due to transportation of stagnant water from Unit 1 to Unit 2. Our value, 0.86, is slightly smaller than these reported values. It is reported that neutron spectrum and flux are not uniform throughout a reactor core50,51. In consequence, variety occurs in the Cs isotopic composition as a function of location within the core: the upper side had smaller 134Cs/137Cs radioactivity ratio. In addition, Snow et al. suggested that (i) when power excursion and/or loss of coolant happened rapidly, Cs isotope ratio would be similar to total core ratio, ii) when coolant loss happened gradually over the course of minutes to hours, radiocesium component would be changed51. Compare to the variety of Cs isotope composition in reactor core, the difference in the obtained 134Cs/137Cs radioactivity ratios is small and close to the ratio in the Unit 1 core. Therefore, power excursion and/or loss of coolant is possibly happened rapidly.

The measured 134Cs/137Cs radioactivity ratio for stagnant water of Unit 1 was clearly smaller than those for the stagnant water of Units 2 and 3 (1.00 ± 0.03), and clearly smaller 134Cs/137Cs radioactivity ratio was observed at north-northwest site near FDNPS and Oshika Peninsula27,28,29,48. Although released amounts of radionuclides by vent operation of the Unit 1 was estimated smaller than those by explosion, the 134Cs/137Cs radioactivity ratio could become an index for the source term analysis.

The 90Sr/137Cs radioactivity ratio in the drain water on March 11th, 2011 is calculated to be 9.3 × 10−4. Nishihara et al. and Asai et al. reported 90Sr radioactivity in the stagnant water of Unit 1 collected on 24th and 27th of March 201149,52. The calculated 90Sr/137Cs radioactivity ratio on 11th of March from the reported values was 1.3 × 10−4, which is one order magnitude lower than that of the drain water. On the other hand, the 90Sr/137Cs radioactivity ratio for rubble collected at the 1st floor of Unit1 on October 2013 and at just outside of Unit 1 on June 25th, 2012 and July 26th and 27th, 2012 was 1.5 × 10−3–5.0 × 10−3 and 1.4 × 10−3–5.6 × 10−3, respectively53,54, and these values were close to our value. In addition, it is reported that 90Sr concentration in the accumulated water increased gradually over time and reached almost comparable to 137Cs concentration in early 201255. Since transportation process of 90Sr during accident and after accident would be different, the sampling location and time should be noted for consideration. This complicated transportation processes of Sr may cause wide variation of 90Sr/137Cs radioactivity ratio for environmental samples56,57,58,59.

In order to discuss the chemical form of radiocesium, mole ratio is important rather than radioactivity ratio. Therefore, radioactivity concentrations were converted to mole concentrations, and the mole concentrations of the target nuclides were divided by the mole concentration of 137Cs to obtain their mole ratios to 137Cs. It is considered that the radionuclides in the drain water were originated from vent operation of the Unit 1. Therefore, the amounts of radionuclides in the Unit 1 core were considered as original inventory. The mole ratios of the target nuclides to 137Cs in the Unit 1 core were calculated by using the reported values of ORIGEN calculation and listed in Table 3. Furthermore, the mole ratios to 137Cs obtained from the measurement were divided by the mole ratios to 137Cs obtained from ORIGEN calculation (hereafter referred to as the ratio of drain water to Unit 1 core ratios). These results are summarized in Table 3. The order of ratios of the drain water to Unit 1 core ratios was 129I >137Cs≈134Cs >125Sb > Mo >90Sr. It is reported that Mo release increased under oxidizing conditions with sufficient water vapor inside reactor core through the formation of volatile species60. However, the amount of Mo was quite low compared to that of Cs in the case of this study, suggesting that the formation of Cs2MoO4 had been suppressed under steam-starved conditions.

Based on the ORIGEN calculation, isotope mole ratio of 129I to the total I and 137Cs to the total Cs in Unit 1 core were calculated to be 78.69% and 40.27%, respectively, at the accident time41 From the isotope mole ratio and measured concentration of 129I and 137Cs in the drain water, the concentration of the total I and the total Cs were calculated to be 2.1 × 10−10 and 3.4 × 10−10 mol/ml, respectively, and the ratio of the total I to the total Cs was 0.6 in the drain water. Compared to the ratio of the total I to the total Cs in the Unit1 core, 0.08, the ratio in the drain water was approximately one order magnitude larger. This trend implies that major chemical form of I was molecular I rather than CsI.

Conclusion

In order to study the behavior of radionuclides released during venting of Unit 1 at FDNPS, radiochemical analysis was conducted for the drain water sampled from the drain pit of the exhaust stack shared between Units 1 and 2.

Although γ-ray spectra of the drain water showed only peaks related to 134Cs and 137Cs, 60Co and 125Sb were detected after elimination of radiocesium using ammonium phosphomolybdate. The decay-corrected 134Cs/137Cs radioactivity ratio to March 11th, 2011 was 0.86.

In this study, Unit 1-originated stable Mo isotopes were clearly detected. Their amounts were quite low compared to Cs, suggesting that the formation of Cs2MoO4 was suppressed under the accident condition.

Inorganic iodine chemical forms in samples collected at the FDNPS were analyzed, and it was found that approximately 90% of iodine existed as I and 10% as IO3 in the samples. Although it is a big subject to extrapolate the present result to accident time, the information related to the iodine chemistry was obtained for the first time.

The ratios of the mole concentration ratios of nuclides to 137Cs in the reactor core to those in the drain water was calculated: The order was 129I >137Cs≈134Cs >125Sb > Mo >90Sr. The larger ratio for 129I implies that iodine was released as molecular iodine rather than CsI.

The results and information obtained from this study could be used as source terms to analyze contamination of environment, and would be applied for the accident analysis of Unit 1 at FDNPS to discuss chemical form of released radionuclides and reactor condition of the accident time.