# Isotopic signature and nano-texture of cesium-rich micro-particles: Release of uranium and fission products from the Fukushima Daiichi Nuclear Power Plant

Junpei Imoto, Asumi Ochiai, Genki Furuki, Mizuki Suetake, Ryohei Ikehara, Kenji Horie, Mami Takehara, Shinya Yamasaki, Kenji Nanba, Toshihiko Ohnuki, et al.

► To cite this version:

Junpei Imoto, Asumi Ochiai, Genki Furuki, Mizuki Suetake, Ryohei Ikehara, et al. Isotopic signature and nano-texture of cesium-rich micro-particles: Release of uranium and fission products from the Fukushima Daiichi Nuclear Power Plant. Scientific Reports, 2017, 7, pp.5409. 10.1038/s41598-017-05910-z. in2p3-01700243

HAL Id: in2p3-01700243

https://hal.in2p3.fr/in2p3-01700243

Submitted on 27 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_1.png" /></figure>

OPEN

# Isotopic signature and nano-texture of cesium-rich micro-particles: Release of uranium and fission products from the Fukushima Daiichi Nuclear Power Plant

Received: 25 April 2017

Accepted: 5 June 2017

Published online: 14 July 2017

Junpei Imoto<fnref n="1" />, Asumi Ochiai<fnref n="1" />, Genki Furuki<fnref n="1" />, Mizuki Suetake<fnref n="1" />, Ryohei Ikehara<fnref n="1" />, Kenji Horie<fnref n="2" />,<fnref n="3" />, Mami Takehara<fnref n="2" />, Shinya Yamasaki<fnref n="4" />, Kenji Nanba<fnref n="5" />, Toshihiko Ohnuki<fnref n="6" />, Gareth T. W. Law<fnref n="7" />, Bernd Grambow<fnref n="8" />, Rodney C. Ewing<fnref n="9" /> & Satoshi Utsunomiya<fnref n="1" />

Highly radioactive cesium-rich microparticles (CsMPs) released from the Fukushima Daiichi Nuclear Power Plant (FDNPP) provide nano-scale chemical fingerprints of the 2011 tragedy. U, Cs, Ba, Rb, K, and Ca isotopic ratios were determined on three CsMPs (3.79–780 Bq) collected within ~10 km from the FDNPP to determine the CsMPs’ origin and mechanism of formation. Apart from crystalline Fe-pollucite, CsFeSi₂O₆·nH₂O, CsMPs are comprised mainly of Zn–Fe-oxide nanoparticles in a SiO₂ glass matrix (up to ~30 wt% of Cs and ~1 wt% of U mainly associated with Zn–Fe-oxide). The ²³⁵U/²³⁸U values in two CsMPs: 0.030 (±0.005) and 0.029 (±0.003), are consistent with that of enriched nuclear fuel. The values are higher than the average burnup estimated by the ORIGEN code and lower than non-irradiated fuel, suggesting non-uniform volatilization of U from melted fuels with different levels of burnup, followed by sorption onto Zn–Fe-oxides. The nano-scale texture and isotopic analyses provide a partial record of the chemical reactions that occurred in the fuel during meltdown. Also, the CsMPs were an important medium of transport for the released radionuclides in a respirable form.

Radionuclides with ~520 PBq initial total activity were released from the Fukushima Daiichi Nuclear Power Plant (FDNPP) as a result of the nuclear disaster that occurred after the Tohoku earthquake on March 11, 2011¹. The released radionuclides, including noble gases (Xe and Kr) and volatile fission products (I, Cs, Te, Sb, and Ag), contaminated the surface over ~14,000 km<fnref n="2" /> surrounding the FDNPP<fnref n="1" />⁻<fnref n="5" />, requiring the evacuation of some 100,000 residents. Radioactive Cs, ¹³⁴Cs and ¹³⁷Cs, are the most important radionuclides contributing to the high radiation in the environment near the FDNPP at present because of their relatively short half-lives, 2.06 and 30.07

<footnote n="1">¹Department of Chemistry, Kyushu University, 744 Motooka, Nishi-ku, Fukuoka-shi, Fukuoka, 819-0395, Japan.</footnote>

<footnote n="2">²National Institute of Polar Research, 10-3, Midori-cho, Tachikawa-shi, Tokyo, 190-8518, Japan.</footnote>

<footnote n="3">³Department of Polar Science, The Graduate University for Advanced Studies (SOKENDAI), Shonan Village, Hayama, Kanagawa, 240-0193, Japan.</footnote>

<footnote n="4">⁴Faculty of Pure and Applied Sciences and Center for Research in Isotopes and Environmental Dynamics, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki, 305-8577, Japan.</footnote>

<footnote n="5">⁵Department of Environment System Management, Faculty of Symbiotic Systems Science, Fukushima University, 1 Kanayagawa, Fukushima, 960-1296, Japan.</footnote>

<footnote n="6">⁶Laboratory for Advanced Nuclear Energy, Institute of Innovative Research, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo, 152-8550, Japan.</footnote>

<footnote n="7">⁷Centre for Radiochemistry Research, School of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK.</footnote>

<footnote n="8">⁸SUBATECH, IMT Atlantique, the University of Nantes CNRS/IN2P3, Nantes, 44307, France.</footnote>

<footnote n="9">⁹Department of Geological Sciences and Center for International Security and Cooperation, Stanford University, Stanford, CA, 94305-2115, USA. Junpei Imoto and Asumi Ochiai contributed equally to this work. Correspondence and requests for materials should be addressed to S.U. (email: utsunomiya.satoshi.998@m.kyushu-u.ac.jp)</footnote>
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Isotopic signature and nano-texture of cesium-rich micro-particles: Release of uranium and fission products from the Fukushima Daiichi Nuclear Power Plant

Received: 25 April 2017

Accepted: 5 June 2017

Published online: 14 July 2017

Junpei Imoto<fnref n="1" />, Asumi Ochiai<fnref n="1" />, Genki Furuki<fnref n="1" />, Mizuki Suetake<fnref n="1" />, Ryohei Ikehara<fnref n="1" />, Kenji Horie<fnref n="2" />,<fnref n="3" />, Mami Takehara<fnref n="2" />, Shinya Yamasaki<fnref n="4" />, Kenji Nanba<fnref n="5" />, Toshihiko Ohnuki<fnref n="6" />, Gareth T. W. Law<fnref n="7" />, Bernd Grambow<fnref n="8" />, Rodney C. Ewing<fnref n="9" /> & Satoshi Utsunomiya<fnref n="1" />

Highly radioactive cesium-rich microparticles (CsMPs) released from the Fukushima Daiichi Nuclear Power Plant (FDNPP) provide nano-scale chemical fingerprints of the 2011 tragedy. U, Cs, Ba, Rb, K, and Ca isotopic ratios were determined on three CsMPs (3.79–780 Bq) collected within ~10 km from the FDNPP to determine the CsMPs’ origin and mechanism of formation. Apart from crystalline Fe-pollucite, CsFeSi₂O₆·nH₂O, CsMPs are comprised mainly of Zn–Fe-oxide nanoparticles in a SiO₂ glass matrix (up to ~30 wt% of Cs and ~1 wt% of U mainly associated with Zn–Fe-oxide). The ²³⁵U/²³⁸U values in two CsMPs: 0.030 (±0.005) and 0.029 (±0.003), are consistent with that of enriched nuclear fuel. The values are higher than the average burnup estimated by the ORIGEN code and lower than non-irradiated fuel, suggesting non-uniform volatilization of U from melted fuels with different levels of burnup, followed by sorption onto Zn–Fe-oxides. The nano-scale texture and isotopic analyses provide a partial record of the chemical reactions that occurred in the fuel during meltdown. Also, the CsMPs were an important medium of transport for the released radionuclides in a respirable form.

Radionuclides with ~520 PBq initial total activity were released from the Fukushima Daiichi Nuclear Power Plant (FDNPP) as a result of the nuclear disaster that occurred after the Tohoku earthquake on March 11, 2011¹. The released radionuclides, including noble gases (Xe and Kr) and volatile fission products (I, Cs, Te, Sb, and Ag), contaminated the surface over ~14,000 km<fnref n="2" /> surrounding the FDNPP<fnref n="1" />⁻<fnref n="5" />, requiring the evacuation of some 100,000 residents. Radioactive Cs, ¹³⁴Cs and ¹³⁷Cs, are the most important radionuclides contributing to the high radiation in the environment near the FDNPP at present because of their relatively short half-lives, 2.06 and 30.07

<footnote n="1">¹Department of Chemistry, Kyushu University, 744 Motooka, Nishi-ku, Fukuoka-shi, Fukuoka, 819-0395, Japan.</footnote>

<footnote n="2">²National Institute of Polar Research, 10-3, Midori-cho, Tachikawa-shi, Tokyo, 190-8518, Japan.</footnote>

<footnote n="3">³Department of Polar Science, The Graduate University for Advanced Studies (SOKENDAI), Shonan Village, Hayama, Kanagawa, 240-0193, Japan.</footnote>

<footnote n="4">⁴Faculty of Pure and Applied Sciences and Center for Research in Isotopes and Environmental Dynamics, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki, 305-8577, Japan.</footnote>

<footnote n="5">⁵Department of Environment System Management, Faculty of Symbiotic Systems Science, Fukushima University, 1 Kanayagawa, Fukushima, 960-1296, Japan.</footnote>

<footnote n="6">⁶Laboratory for Advanced Nuclear Energy, Institute of Innovative Research, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo, 152-8550, Japan.</footnote>

<footnote n="7">⁷Centre for Radiochemistry Research, School of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK.</footnote>

<footnote n="8">⁸SUBATECH, IMT Atlantique, the University of Nantes CNRS/IN2P3, Nantes, 44307, France.</footnote>

<footnote n="9">⁹Department of Geological Sciences and Center for International Security and Cooperation, Stanford University, Stanford, CA, 94305-2115, USA. Junpei Imoto and Asumi Ochiai contributed equally to this work. Correspondence and requests for materials should be addressed to S.U. (email: utsunomiya.satoshi.998@m.kyushu-u.ac.jp)</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** A map showing the location of the samples. This map was produced based on Furuki *et al.*¹⁶ using the power point.</figcaption></figure>

years, respectively⁶. About 1–7% of the Cs inventory of three reactor cores was released³,⁷. Previous studies of the distribution and migration of radioactive Cs in the surface environment around Fukushima⁴,⁶,⁸–¹² indicated that an initially soluble form of radioactive Cs was released from the damaged reactors and spread over the Fukushima Prefecture and the surrounding area through dry and wet deposition³⁻⁵,⁸. Subsequently, a soluble form of Cs was tightly bound to the interlayers of clays, such as vermiculite, and remained within the top ~5 cm of the soil³, ⁹, ¹⁰. However, Cs in the contaminated soils is heterogeneous and concentrated locally at the micron scale as hot spots, as revealed through autoradiography⁹⁻¹¹, and this heterogeneity has not been fully described. A possible cause of the heterogeneity is the formation of Cs-rich microparticles (CsMPs), with a high Cs radioactivity per unit mass, found at a range of distances from the FDNPP¹². This is another important route of Cs migration in the environments. Different from the soluble Cs, the CsMPs are sparingly soluble in water. The ¹³⁴Cs/¹³⁷Cs radioactivity ratio of ~1 indicates that they originated from the FDNPP. CsMPs were initially considered as amorphous glass particles containing various elements derived from the reactors that were melted into the glass matrix¹²⁻¹⁴. However, recent studies have shown that the Cs concentrations are associated with discrete Zn–Fe-oxide nano-particles embedded in a pure SiO₂ glass matrix, as well as numerous nano-scale inclusions with a variety of fission products¹⁵, ¹⁶. The nano-texture within the CsMPs records the chemical reactions that took place during the meltdowns inside the reactors. The sparingly soluble CsMPs were identified as dominant Cs carriers during the initial Cs fallout in Tokyo, Japan, and they are expected to differ from soluble Cs in terms of environmental mobility and health impacts¹⁵. In addition to the reaction sequence of the CsMPs formation inside the reactors, trace amounts of U are present in CsMPs. The origin of U and the other nuclides remains uncertain, as insulating material in the reactors also contains trace amount of natural U. The occurrence of U derived from inside the reactors is potentially useful for understanding the reactions that nuclear fuels experienced during meltdown and even the status of the melted nuclear fuels in the damaged reactors¹⁷. This is crucial for developing an appropriate decommissioning strategy. The hypothesis of this study is that the isotopic ratio analysis in combination with atomic-resolution analysis reveals the source of U and the release processes associated with CsMPs from the FDNPP. The CsMPs are a very unique form of condensed matter that was created during the meltdown events inside the reactors at the FDNPP, very different from general concept of Cs release as a soluble Cs species, such as CsI and CsOH. Isotopic analyses of the individual CsMPs are reported for the first time. Isotopic analysis of the other nuclides, stable and radiogenic, also allow us to identify their sources, natural or fissionogenic.

## Results

**Shape, composition, and radioactivity of CsMPs.** Four CsMPs were found at three localities, which are hereafter labeled as OTZ3, OTZ10, KOI2, and OMR1 (Fig. 1). Detailed information on the sampling is provided in the method section. Figure 2 shows secondary electron scanning electron microscopy (SEM) images and energy-dispersive X-ray analysis (EDX) elemental maps of major elements of the CsMPs. In previous studies¹²⁻¹⁶, the shape of CsMPs has been described as spherical. In contrast, these particles appear as irregularly shaped aggregates rather than spheres. In addition, the radioactivity of the presently examined CsMPs ranges from 3.79
<figure><img src="image_3.png" /><figcaption>**Figure 1.** A map showing the location of the samples. This map was produced based on Furuki *et al.*¹⁶ using the power point.</figcaption></figure>

years, respectively⁶. About 1–7% of the Cs inventory of three reactor cores was released³,⁷. Previous studies of the distribution and migration of radioactive Cs in the surface environment around Fukushima⁴,⁶,⁸–¹² indicated that an initially soluble form of radioactive Cs was released from the damaged reactors and spread over the Fukushima Prefecture and the surrounding area through dry and wet deposition³⁻⁵,⁸. Subsequently, a soluble form of Cs was tightly bound to the interlayers of clays, such as vermiculite, and remained within the top ~5 cm of the soil³, ⁹, ¹⁰. However, Cs in the contaminated soils is heterogeneous and concentrated locally at the micron scale as hot spots, as revealed through autoradiography⁹⁻¹¹, and this heterogeneity has not been fully described. A possible cause of the heterogeneity is the formation of Cs-rich microparticles (CsMPs), with a high Cs radioactivity per unit mass, found at a range of distances from the FDNPP¹². This is another important route of Cs migration in the environments. Different from the soluble Cs, the CsMPs are sparingly soluble in water. The ¹³⁴Cs/¹³⁷Cs radioactivity ratio of ~1 indicates that they originated from the FDNPP. CsMPs were initially considered as amorphous glass particles containing various elements derived from the reactors that were melted into the glass matrix¹²⁻¹⁴. However, recent studies have shown that the Cs concentrations are associated with discrete Zn–Fe-oxide nano-particles embedded in a pure SiO₂ glass matrix, as well as numerous nano-scale inclusions with a variety of fission products¹⁵, ¹⁶. The nano-texture within the CsMPs records the chemical reactions that took place during the meltdowns inside the reactors. The sparingly soluble CsMPs were identified as dominant Cs carriers during the initial Cs fallout in Tokyo, Japan, and they are expected to differ from soluble Cs in terms of environmental mobility and health impacts¹⁵. In addition to the reaction sequence of the CsMPs formation inside the reactors, trace amounts of U are present in CsMPs. The origin of U and the other nuclides remains uncertain, as insulating material in the reactors also contains trace amount of natural U. The occurrence of U derived from inside the reactors is potentially useful for understanding the reactions that nuclear fuels experienced during meltdown and even the status of the melted nuclear fuels in the damaged reactors¹⁷. This is crucial for developing an appropriate decommissioning strategy. The hypothesis of this study is that the isotopic ratio analysis in combination with atomic-resolution analysis reveals the source of U and the release processes associated with CsMPs from the FDNPP. The CsMPs are a very unique form of condensed matter that was created during the meltdown events inside the reactors at the FDNPP, very different from general concept of Cs release as a soluble Cs species, such as CsI and CsOH. Isotopic analyses of the individual CsMPs are reported for the first time. Isotopic analysis of the other nuclides, stable and radiogenic, also allow us to identify their sources, natural or fissionogenic.

## Results

**Shape, composition, and radioactivity of CsMPs.** Four CsMPs were found at three localities, which are hereafter labeled as OTZ3, OTZ10, KOI2, and OMR1 (Fig. 1). Detailed information on the sampling is provided in the method section. Figure 2 shows secondary electron scanning electron microscopy (SEM) images and energy-dispersive X-ray analysis (EDX) elemental maps of major elements of the CsMPs. In previous studies¹²⁻¹⁶, the shape of CsMPs has been described as spherical. In contrast, these particles appear as irregularly shaped aggregates rather than spheres. In addition, the radioactivity of the presently examined CsMPs ranges from 3.79
<figure><img src="image_3.png" /><figcaption>**Figure 2.** SEM images of CsMPs associated with the elemental maps.</figcaption></figure>

to 780 Bq (Table 1). These values are approximately thirty three times higher than that of spherical CsMPs$^{12–16}$. The $^{134}$Cs/$^{137}$Cs radioactivity ratios of the OTZ3, KOI2, and OMR1 used for the SIMS analysis are 1.06–1.08 (average 1.07), which approximately corresponds to ~26 GWd/tU according to OrigenArp calculations$^{18}$. Although the $^{134}$Cs/$^{137}$Cs isotopic ratio values are close to those of the reactor Units #2 and #3 in the FDNPP$^{19}$, the specific reactor source could not be determined based only on the isotopic ratios, because burnups of the irradiated fuels in each reactor are heterogeneous depending on the positions of the fuel assemblies.

**Isotopic ratio based on secondary ion mass spectrometry (SIMS) and γ spectrometry.** The results of the isotopic analysis of the three CsMPs, OTZ3, KOI2, and OMR1, are summarized in Table 2. In general, the $^{235}$U/$^{238}$U isotopic ratio in nature is 0.00729, whereas the ratio in non-irradiated nuclear fuel is typically >~0.03 due to enrichment of $^{235}$U$^{20}$. As shown in Table 2, U isotopes were successfully determined to be 0.029584 and 0.029341 in OTZ3 and KOI2, respectively, whereas U was not detected in OMR1. The U isotopic ratio for the standard specimen NIST SRM610 is comparable to the value in a previous report$^{21}$, and the deviation of the
<figure><img src="image_4.png" /><figcaption>**Figure 2.** SEM images of CsMPs associated with the elemental maps.</figcaption></figure>

to 780 Bq (Table 1). These values are approximately thirty three times higher than that of spherical CsMPs$^{12–16}$. The $^{134}$Cs/$^{137}$Cs radioactivity ratios of the OTZ3, KOI2, and OMR1 used for the SIMS analysis are 1.06–1.08 (average 1.07), which approximately corresponds to ~26 GWd/tU according to OrigenArp calculations$^{18}$. Although the $^{134}$Cs/$^{137}$Cs isotopic ratio values are close to those of the reactor Units #2 and #3 in the FDNPP$^{19}$, the specific reactor source could not be determined based only on the isotopic ratios, because burnups of the irradiated fuels in each reactor are heterogeneous depending on the positions of the fuel assemblies.

**Isotopic ratio based on secondary ion mass spectrometry (SIMS) and γ spectrometry.** The results of the isotopic analysis of the three CsMPs, OTZ3, KOI2, and OMR1, are summarized in Table 2. In general, the $^{235}$U/$^{238}$U isotopic ratio in nature is 0.00729, whereas the ratio in non-irradiated nuclear fuel is typically >~0.03 due to enrichment of $^{235}$U$^{20}$. As shown in Table 2, U isotopes were successfully determined to be 0.029584 and 0.029341 in OTZ3 and KOI2, respectively, whereas U was not detected in OMR1. The U isotopic ratio for the standard specimen NIST SRM610 is comparable to the value in a previous report$^{21}$, and the deviation of the
<table><caption>**Table 1.** Radioactivity of four CsMPs described in the present study. The radioactivity was decay-corrected to March 12, 2011, 15:36 JST. The uncertainties in the radioactivity analyses are given in the parenthesis. The particle size was measured using the SEM image.</caption><thead><tr><th>Sample</th><th>Particle size (μm)</th><th>Radioactivity of <sup>134</sup>Cs (Bq)</th><th>Radioactivity of <sup>137</sup>Cs (Bq)</th><th>Ratio of <sup>134</sup>Cs/<sup>137</sup>Cs radioactivity</th></tr></thead><tbody><tr><td>OTZ3</td><td>&lt;17.3</td><td>401(±3.7)</td><td>379(±1.5)</td><td>1.06</td></tr><tr><td>OTZ10</td><td>2.2</td><td>1.85(±0.080)</td><td>1.94(±0.031)</td><td>0.955</td></tr><tr><td>OMR1</td><td>&lt;5.2</td><td>21.5(±0.50)</td><td>20.3(±0.19)</td><td>1.06</td></tr><tr><td>KOI2</td><td>&lt;4.4</td><td>18.2(±0.16)</td><td>16.8(±0.05)</td><td>1.08</td></tr></tbody></table>

three spot analyses with a standard deviation of ca. 0.00002 is 0.000001, indicating the high precision of the SIMS analysis.

Stable Cs is represented by monoisotopic ¹³³Cs (100% natural abundance), although various radioisotopes including ¹³⁴Cs, ¹³⁵Cs and ¹³⁷Cs are produced in a nuclear reactor, and the ratios of the isotopes depend on burnup. ¹³⁰Ba (0.1058 at.%), ¹³²Ba (0.1012 at.%), ¹³⁴Ba (2.417 at.%), ¹³⁵Ba (6.592 at.%), ¹³⁶Ba (7.853 at.%), ¹³⁷Ba (11.232 at.%), and ¹³⁸Ba (71.699 at.%) are stable isotopes in nature (their abundance in at.%)²². Most Ba radioisotopes (> ¹²²Ba) typically decay through isobars to the most stable isotope in various modes such as β⁻, β⁺, and electron capture. Cesium and Ba isotopes are difficult to distinguish from one another in the present SIMS analysis. However, based on gamma spectrometry, the isotopic ratio of ¹³⁴Cs/¹³⁷Cs(γ) after decay-correction to the time of the SIMS analysis were calculated to be 0.0131, 0.0135, and 0.0132 for OTZ3, KOI2, and OMR1, respectively. The ¹³⁴Cs/¹³⁷Cs(γ) isotopic ratio after decay-correction to 15:36 March 12, 2011 was ~0.073. The ¹³⁶Ba/¹³⁸Ba isotopic ratio was determined to be 0.1523, 0.1638, and 0.4293 for OTZ3, KOI2, and OMR1, respectively. The isotopic ratio of ¹³⁴Cs + Ba<sub>radiogenic</sub>/¹³⁸Ba determined by SIMS ranges from 35 to 138. The isotopic ratios ¹³⁵Cs/¹³⁷Cs + Ba<sub>radiogenic</sub> and ¹³⁷Cs + Ba<sub>radiogenic</sub>/¹³³Cs were determined to be 0.39–0.40 and 0.92–0.99, respectively. The Rb isotopic ratio, ⁸⁷Rb/⁸⁵Rb, ranges 2.2–2.4. Strontium-90 is not detected in the CsMPs, and the isotopic ratios of Ca and K reveal that these elements are natural in origin.

**Nanoscale structure and composition of CsMPs.** Two FIB thin sections were successfully prepared from the OTZ3 CsMP labeled OTZ3-1 and OTZ3-2, which were previously analyzed by SIMS. Another FIB specimen was prepared from OTZ10 CsMP that was found in the same soil sample.

A high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) image of OTZ3-1 CsMP is shown in Fig. 3a, in which the white arrow indicates the position of the particle. There are two distinct zones with different contrast (Fig. 3b). The elemental maps of the boundary between the two zones reveal that the Cs concentration is much higher in the area with bright contrast, although the other elements are distributed homogeneously (Fig. 3c, Table S1). The area with dark contrast (edx1) contains a lower Cs-content as compared with the zone with bright contrast (edx2) (Fig. 3d). The magnified HAADF-STEM image of the low Cs zone (dark contrast) is composed of nanoparticles that consist of Zn and Fe, as revealed by the elemental maps (Fig. 3e). Cesium is associated with these nanoparticles. The distribution of silicon could not be resolved due to the thickness of the specimen (Fig. 3e). Based on the HRTEM image, the nanoparticles are most likely identified as franklinite structures (ZnFe₂O₄, Fd3m, Z = 8)²³, which has also been identified in other CsMPs¹⁵,¹⁶ (Fig. 3f).

The high Cs zone (represented by edx2 in Fig. 3b) appears to contain relatively large nanoparticles, ~50 nm in size, with a high Cs-content (Fig. S1a) as compared with the franklinite nanoparticles found in the low Cs-zone. The large Cs-rich particle was subsequently identified as an Fe-rich pollucite structure (CsFeSi₂O₆·nH₂O, Ia3d, Z = 16)²⁴ based on the SAED pattern (Fig. S1b). Pollucite is a zeolite with the ideal formula (Cs,Na)(Al,Si)₃O₆·nH₂O²⁵. The Fe-pollucite is unstable under electron beam irradiation and becomes amorphous (Fig. S1c), possibly due to the loss of water. In addition, the HAADF-STEM image of the boundary between the high- and low-Cs zones displays different levels of contrast after electron beam irradiation (Fig. S1d); that is, the bright contrast of Fe-pollucite particles disappears, indicating the structural degradation of Fe-pollucite, consistent with the loss of peaks in the SAED pattern. Electron energy-loss spectroscopy (EELS) of the high-Cs zone reveals that Cs is monovalent as indicated by the Cs M-edge (Fig. S1e). Barium is closely associated with Cs-bearing phases rather than forming a separate Ba-rich phase, which is expected since almost all Ba isotopes in the CsMP are decay daughters of radioactive Cs. The iron L-edge confirms that Fe is present in the form of oxides.

The same major elements, Si, Fe, Zn, and Cs, are present in the other FIB specimen from the same aggregate, OTZ3-2, as well (Fig. 4a). The area indicated by the white arrow is magnified, and the texture is rather homogeneous, as shown by the HAADF image contrast (Fig. 4b), which is different from the spherical texture observed in the high-Cs zone (edx2 in Fig. 3b) in OTZ3-1. However, the Cs-content (Fig. 4c) is similar to that of the high-Cs zone in OTZ3-1. Also observed was a single crystal as large as ~0.5 µm, evidenced by continuous lattice fringes in the HRTEM image (Fig. 4d). Based on the diffraction patterns obtained from the two major zone axes, this large Cs-phase is an Fe-pollucite (Fig. 4e).

In order to determine the presence of U, a FIB-prepared cross-section of the OTZ10 CsMP was additionally investigated by TEM (Fig. 5a). Typically, the SAED exhibits broad diffraction maxima that correspond to diffuse scattering from amorphous domains (Fig. 5a inset). The elemental map of the entire cross-section displays a homogeneous distribution of the major constituents at the bulk scale (Fig. 5b). The Cs-concentration is 7–10 wt% (Fig. 5c, Table S1). Despite the apparently homogeneous distribution and the diffuse scattering halo in the SAED, the magnified HAADF-STEM image shows phase separation as evidenced by differences in contrast (Fig. S2a).
<table><caption>**Table 1.** Radioactivity of four CsMPs described in the present study. The radioactivity was decay-corrected to March 12, 2011, 15:36 JST. The uncertainties in the radioactivity analyses are given in the parenthesis. The particle size was measured using the SEM image.</caption><thead><tr><th>Sample</th><th>Particle size (μm)</th><th>Radioactivity of <sup>134</sup>Cs (Bq)</th><th>Radioactivity of <sup>137</sup>Cs (Bq)</th><th>Ratio of <sup>134</sup>Cs/<sup>137</sup>Cs radioactivity</th></tr></thead><tbody><tr><td>OTZ3</td><td>&lt;17.3</td><td>401(±3.7)</td><td>379(±1.5)</td><td>1.06</td></tr><tr><td>OTZ10</td><td>2.2</td><td>1.85(±0.080)</td><td>1.94(±0.031)</td><td>0.955</td></tr><tr><td>OMR1</td><td>&lt;5.2</td><td>21.5(±0.50)</td><td>20.3(±0.19)</td><td>1.06</td></tr><tr><td>KOI2</td><td>&lt;4.4</td><td>18.2(±0.16)</td><td>16.8(±0.05)</td><td>1.08</td></tr></tbody></table>

three spot analyses with a standard deviation of ca. 0.00002 is 0.000001, indicating the high precision of the SIMS analysis.

Stable Cs is represented by monoisotopic ¹³³Cs (100% natural abundance), although various radioisotopes including ¹³⁴Cs, ¹³⁵Cs and ¹³⁷Cs are produced in a nuclear reactor, and the ratios of the isotopes depend on burnup. ¹³⁰Ba (0.1058 at.%), ¹³²Ba (0.1012 at.%), ¹³⁴Ba (2.417 at.%), ¹³⁵Ba (6.592 at.%), ¹³⁶Ba (7.853 at.%), ¹³⁷Ba (11.232 at.%), and ¹³⁸Ba (71.699 at.%) are stable isotopes in nature (their abundance in at.%)²². Most Ba radioisotopes (> ¹²²Ba) typically decay through isobars to the most stable isotope in various modes such as β⁻, β⁺, and electron capture. Cesium and Ba isotopes are difficult to distinguish from one another in the present SIMS analysis. However, based on gamma spectrometry, the isotopic ratio of ¹³⁴Cs/¹³⁷Cs(γ) after decay-correction to the time of the SIMS analysis were calculated to be 0.0131, 0.0135, and 0.0132 for OTZ3, KOI2, and OMR1, respectively. The ¹³⁴Cs/¹³⁷Cs(γ) isotopic ratio after decay-correction to 15:36 March 12, 2011 was ~0.073. The ¹³⁶Ba/¹³⁸Ba isotopic ratio was determined to be 0.1523, 0.1638, and 0.4293 for OTZ3, KOI2, and OMR1, respectively. The isotopic ratio of ¹³⁴Cs + Ba<sub>radiogenic</sub>/¹³⁸Ba determined by SIMS ranges from 35 to 138. The isotopic ratios ¹³⁵Cs/¹³⁷Cs + Ba<sub>radiogenic</sub> and ¹³⁷Cs + Ba<sub>radiogenic</sub>/¹³³Cs were determined to be 0.39–0.40 and 0.92–0.99, respectively. The Rb isotopic ratio, ⁸⁷Rb/⁸⁵Rb, ranges 2.2–2.4. Strontium-90 is not detected in the CsMPs, and the isotopic ratios of Ca and K reveal that these elements are natural in origin.

**Nanoscale structure and composition of CsMPs.** Two FIB thin sections were successfully prepared from the OTZ3 CsMP labeled OTZ3-1 and OTZ3-2, which were previously analyzed by SIMS. Another FIB specimen was prepared from OTZ10 CsMP that was found in the same soil sample.

A high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) image of OTZ3-1 CsMP is shown in Fig. 3a, in which the white arrow indicates the position of the particle. There are two distinct zones with different contrast (Fig. 3b). The elemental maps of the boundary between the two zones reveal that the Cs concentration is much higher in the area with bright contrast, although the other elements are distributed homogeneously (Fig. 3c, Table S1). The area with dark contrast (edx1) contains a lower Cs-content as compared with the zone with bright contrast (edx2) (Fig. 3d). The magnified HAADF-STEM image of the low Cs zone (dark contrast) is composed of nanoparticles that consist of Zn and Fe, as revealed by the elemental maps (Fig. 3e). Cesium is associated with these nanoparticles. The distribution of silicon could not be resolved due to the thickness of the specimen (Fig. 3e). Based on the HRTEM image, the nanoparticles are most likely identified as franklinite structures (ZnFe₂O₄, Fd3m, Z = 8)²³, which has also been identified in other CsMPs¹⁵,¹⁶ (Fig. 3f).

The high Cs zone (represented by edx2 in Fig. 3b) appears to contain relatively large nanoparticles, ~50 nm in size, with a high Cs-content (Fig. S1a) as compared with the franklinite nanoparticles found in the low Cs-zone. The large Cs-rich particle was subsequently identified as an Fe-rich pollucite structure (CsFeSi₂O₆·nH₂O, Ia3d, Z = 16)²⁴ based on the SAED pattern (Fig. S1b). Pollucite is a zeolite with the ideal formula (Cs,Na)(Al,Si)₃O₆·nH₂O²⁵. The Fe-pollucite is unstable under electron beam irradiation and becomes amorphous (Fig. S1c), possibly due to the loss of water. In addition, the HAADF-STEM image of the boundary between the high- and low-Cs zones displays different levels of contrast after electron beam irradiation (Fig. S1d); that is, the bright contrast of Fe-pollucite particles disappears, indicating the structural degradation of Fe-pollucite, consistent with the loss of peaks in the SAED pattern. Electron energy-loss spectroscopy (EELS) of the high-Cs zone reveals that Cs is monovalent as indicated by the Cs M-edge (Fig. S1e). Barium is closely associated with Cs-bearing phases rather than forming a separate Ba-rich phase, which is expected since almost all Ba isotopes in the CsMP are decay daughters of radioactive Cs. The iron L-edge confirms that Fe is present in the form of oxides.

The same major elements, Si, Fe, Zn, and Cs, are present in the other FIB specimen from the same aggregate, OTZ3-2, as well (Fig. 4a). The area indicated by the white arrow is magnified, and the texture is rather homogeneous, as shown by the HAADF image contrast (Fig. 4b), which is different from the spherical texture observed in the high-Cs zone (edx2 in Fig. 3b) in OTZ3-1. However, the Cs-content (Fig. 4c) is similar to that of the high-Cs zone in OTZ3-1. Also observed was a single crystal as large as ~0.5 µm, evidenced by continuous lattice fringes in the HRTEM image (Fig. 4d). Based on the diffraction patterns obtained from the two major zone axes, this large Cs-phase is an Fe-pollucite (Fig. 4e).

In order to determine the presence of U, a FIB-prepared cross-section of the OTZ10 CsMP was additionally investigated by TEM (Fig. 5a). Typically, the SAED exhibits broad diffraction maxima that correspond to diffuse scattering from amorphous domains (Fig. 5a inset). The elemental map of the entire cross-section displays a homogeneous distribution of the major constituents at the bulk scale (Fig. 5b). The Cs-concentration is 7–10 wt% (Fig. 5c, Table S1). Despite the apparently homogeneous distribution and the diffuse scattering halo in the SAED, the magnified HAADF-STEM image shows phase separation as evidenced by differences in contrast (Fig. S2a).
<table><caption>**Table 2.** Summary of isotope analysis for three CsMPs using SIMS and γ spectrometry. A NIST standard, SRM610, which contains depleted U, was also analyzed as a reference. Ten scans were conducted for each analytical spot on the CsMPs, meaning that the sequence of ten analyses represents a depth profile of the variation in the isotopic ratios; thus, SIMS analyses provide the isotope signatures inside the CsMPs. The average values are given in the table with the standard deviation calculated for the ten analyses and the standard deviations are given in the parenthesis as 1σ. The isotopic ratios recalculated from the results of γ spectrometry are annotated as (γ). Conversion of radioactivity to isotopic ratio was conducted by decay-correcting to the time of SIMS analysis. n/d stands for not detected.</caption><thead><tr><th>Isotopic ratio</th><th>OTZ3 (1σ)</th><th>KOI2 (1σ)</th><th>OMR1 (1σ)</th><th>Natural abundance ratio</th><th>SRM610 (1σ)</th></tr></thead><tbody><tr><td><sup>235</sup>U/<sup>238</sup>U</td><td>0.029584 (0.004967)</td><td>0.029341 (0.003032)</td><td>n/d</td><td>0.00729</td><td>0.002388 (0.000001)</td></tr><tr><td><sup>134</sup>Cs/<sup>137</sup>Cs (γ)</td><td>0.02488 (0.00219)</td><td>0.03184 (0.00107)</td><td>0.03167 (0.00037)</td><td>—</td><td>—</td></tr><tr><td><sup>134</sup>Ba<sub>radiogenic</sub>/<sup>137</sup>Ba<sub>radiogenic</sub> (γ)</td><td>0.0131</td><td>0.0135</td><td>0.0132</td><td>—</td><td>—</td></tr><tr><td><sup>134</sup>Cs/<sup>134</sup>Ba<sub>radiogenic</sub> (γ)</td><td>0.515</td><td>0.517</td><td>0.520</td><td>—</td><td>—</td></tr><tr><td><sup>137</sup>Cs/<sup>137</sup>Ba<sub>radiogenic</sub> (γ)</td><td>0.190</td><td>0.190</td><td>0.191</td><td>—</td><td>—</td></tr><tr><td><sup>137</sup>Cs/<sup>137</sup>Ba<sub>radiogenic</sub> (γ)</td><td>7.47</td><td>7.46</td><td>7.47</td><td>—</td><td>—</td></tr><tr><td><sup>134</sup>(Cs + Ba<sub>radiogenic</sub>)/<sup>137</sup>(Cs + Ba<sub>radiogenic</sub>) (γ)</td><td>0.0724</td><td>0.0742</td><td>0.0727</td><td>—</td><td>—</td></tr><tr><td><sup>135</sup>Cs/<sup>137</sup>(Cs + Ba<sub>radiogenic</sub>)</td><td>0.3929 (0.0089)</td><td>0.4007 (0.0050)</td><td>0.4011 (0.0020)</td><td>—</td><td>—</td></tr><tr><td><sup>135</sup>Cs/<sup>133</sup>Cs</td><td>0.3901 (0.0121)</td><td>0.3685 (0.0062)</td><td>0.3777 (0.0025)</td><td>—</td><td>—</td></tr><tr><td><sup>137</sup>(Cs + Ba<sub>radiogenic</sub>)/<sup>133</sup>Cs</td><td>0.9941 (0.0543)</td><td>0.9199 (0.0246)</td><td>0.9417 (0.0102)</td><td>—</td><td>—</td></tr><tr><td><sup>134</sup>(Cs + Ba<sub>radiogenic</sub>)/<sup>138</sup>Ba</td><td>36.6 (3.7)</td><td>35.3 (5.9)</td><td>137.7 (21.7)</td><td>0.03371</td><td>0.03371 (0.00009)</td></tr><tr><td><sup>136</sup>Ba/<sup>138</sup>Ba</td><td>0.1523 (0.0124)</td><td>0.1638 (0.0191)</td><td>0.4293 (0.0697)</td><td>0.1095</td><td>0.1097 (0.0002)</td></tr><tr><td><sup>87</sup>Rb/<sup>85</sup>Rb</td><td>2.407 (0.014)</td><td>2.284 (0.017)</td><td>2.236 (0.120)</td><td>0.3856</td><td>0.3857 (0.00009)</td></tr><tr><td><sup>90</sup>Sr</td><td>n/d</td><td>n/d</td><td>n/d</td><td>—</td><td>—</td></tr><tr><td><sup>41</sup>K/<sup>39</sup>K</td><td>0.07217 (0.00001)</td><td>—</td><td>—</td><td>0.07216</td><td>0.07214 (0.00003)</td></tr><tr><td><sup>40</sup>K/<sup>39</sup>K</td><td>0.000127 (0.000041)</td><td>—</td><td>—</td><td>0.000125</td><td>—</td></tr><tr><td><sup>43</sup>Ca/<sup>44</sup>Ca</td><td>0.06465 (0.00115)</td><td>—</td><td>—</td><td>0.06472</td><td>0.06469 (0.000005)</td></tr></tbody></table>

The elemental map and the EDX analysis indicate that the phase in bright contrast contains Fe, Zn, Cs, and Sn; whereas, the area of dark contrast corresponds to SiO$_2$-domains (Fig. S2b and c). The HRTEM image also reveals the presence of many nanocrystals (Fig. 5d). Based on the d-spacing and the FFT image, these nanocrystals are franklinite$^{23}$. Another magnified HAADF-STEM image with an elemental map reveals that the U distribution is closely associated with Zn-Fe-oxide nanoparticles (Fig. 5e). The U concentration is 0.8 wt% (including Si) or 1.1 wt% (excluding Si), assuming U is in the form of U$_3$O$_8$ for the quantitative EDX analysis (edx2, Fig. 5f).

Discussion

The U isotopic ratios of the ten scans on the CsMPs vary with a rather large deviation, ~0.0050 for OTZ3 and ~0.0030 for KOI2, most likely because the U concentration in the CsMPs is low (ppm level in the bulk). In addition, the analyses performed are not sufficiently precise to provide the absolute concentrations as high as that of the analysis on SRM610, because the irregular surfaces of the CsMPs aggregate were analyzed without polishing. The effect of isotopic fractionation should be negligibly small during the formation of these CsMPs. Even with a large deviation in the isotopic ratios, $^{235}$U/$^{238}$U in OTZ3 and KOI2 are significantly high, 0.029584 and 0.029341, as compared with the U isotopic ratio in nature, 0.00729. As expected, the U isotopic ratios are close to that of a typical nuclear fuel enriched in $^{235}$U.

The presence of $^{138}$Ba indicates either a fissionogenic or natural origin. Natural $^{136}$Ba/$^{138}$Ba, 0.1095, is an order of magnitude greater than $^{136}$Ba/$^{138}$Ba of fissionogenic origin, ~0.0126, estimated for Units #1–3 by calculation using the ORIGEN code$^{19}$, which largely varies depending on the burnup: 0.004 for 10 MWd/kgU, 0.011 for 20 MWd/kgU, and 0.014 for 30 MWd/kgU. $^{136}$Ba/$^{138}$Ba is 0.1523, 0.1638, and 0.4293 for OTZ3, KOI2, and OMR1, respectively. These values are much higher than those of fissionogenic Ba. The $^{136}$Ba/$^{138}$Ba ratio may indicate the contribution of a low-yield fission product; the amount of $^{136}$Cs (half-life of 13.16 days) can be remarkable in some cases depending on the amount of volatilized Ba. Although the amount of fissionogenic Ba isotopes is greater than Cs isotopes, the Cs can volatilize to greater extent than Ba during the meltdown, meaning that CsMP can selectively incorporate Cs as compared with Ba. Because $^{136}$Ba is a decay product of $^{136}$Cs (half-life of ~13 days), $^{136}$Ba was the most likely incorporated in the form of Cs species, leading to the higher value of $^{136}$Ba/$^{138}$Ba. On the other hand, $^{134}$(Cs + Ba)/$^{138}$Ba determined by SIMS varies from 35 to 138, which is much higher than $^{134}$Ba/$^{138}$Ba in nature (0.03371) and that of fissionogenic origin (0.0373–0.0566). Based on results from calculations using the ORIGEN code$^{19}$, the initial amount of fissionogenic $^{134}$Cs should be about the same as $^{134}$Ba. Thus, the high $^{134}$(Cs + Ba)/$^{138}$Ba is probably due to the higher volatility of Cs$^{26}$, as well as the greater amount of Cs incorporated into the CsMPs as compared with Ba. Thus, Ba isotopes of natural and fissionogenic origin that were initially incorporated in the CsMPs are negligibly small. Rather, $^{134}$Ba and $^{137}$Ba are decay daughters derived from Cs isotopes. $^{134}$Cs/$^{137}$Cs (γ) isotopic ratio at 15:36 March 12, 2011 was ~0.073, and based on the calculation of decayed amount of $^{134}$Cs and $^{137}$Cs corresponding to that of $^{134}$Ba and $^{137}$Ba in the CsMPs, the two isotopic
<table><caption>**Table 2.** Summary of isotope analysis for three CsMPs using SIMS and γ spectrometry. A NIST standard, SRM610, which contains depleted U, was also analyzed as a reference. Ten scans were conducted for each analytical spot on the CsMPs, meaning that the sequence of ten analyses represents a depth profile of the variation in the isotopic ratios; thus, SIMS analyses provide the isotope signatures inside the CsMPs. The average values are given in the table with the standard deviation calculated for the ten analyses and the standard deviations are given in the parenthesis as 1σ. The isotopic ratios recalculated from the results of γ spectrometry are annotated as (γ). Conversion of radioactivity to isotopic ratio was conducted by decay-correcting to the time of SIMS analysis. n/d stands for not detected.</caption><thead><tr><th>Isotopic ratio</th><th>OTZ3 (1σ)</th><th>KOI2 (1σ)</th><th>OMR1 (1σ)</th><th>Natural abundance ratio</th><th>SRM610 (1σ)</th></tr></thead><tbody><tr><td><sup>235</sup>U/<sup>238</sup>U</td><td>0.029584 (0.004967)</td><td>0.029341 (0.003032)</td><td>n/d</td><td>0.00729</td><td>0.002388 (0.000001)</td></tr><tr><td><sup>134</sup>Cs/<sup>137</sup>Cs (γ)</td><td>0.02488 (0.00219)</td><td>0.03184 (0.00107)</td><td>0.03167 (0.00037)</td><td>—</td><td>—</td></tr><tr><td><sup>134</sup>Ba<sub>radiogenic</sub>/<sup>137</sup>Ba<sub>radiogenic</sub> (γ)</td><td>0.0131</td><td>0.0135</td><td>0.0132</td><td>—</td><td>—</td></tr><tr><td><sup>134</sup>Cs/<sup>134</sup>Ba<sub>radiogenic</sub> (γ)</td><td>0.515</td><td>0.517</td><td>0.520</td><td>—</td><td>—</td></tr><tr><td><sup>137</sup>Cs/<sup>137</sup>Ba<sub>radiogenic</sub> (γ)</td><td>0.190</td><td>0.190</td><td>0.191</td><td>—</td><td>—</td></tr><tr><td><sup>137</sup>Cs/<sup>137</sup>Ba<sub>radiogenic</sub> (γ)</td><td>7.47</td><td>7.46</td><td>7.47</td><td>—</td><td>—</td></tr><tr><td><sup>134</sup>(Cs + Ba<sub>radiogenic</sub>)/<sup>137</sup>(Cs + Ba<sub>radiogenic</sub>) (γ)</td><td>0.0724</td><td>0.0742</td><td>0.0727</td><td>—</td><td>—</td></tr><tr><td><sup>135</sup>Cs/<sup>137</sup>(Cs + Ba<sub>radiogenic</sub>)</td><td>0.3929 (0.0089)</td><td>0.4007 (0.0050)</td><td>0.4011 (0.0020)</td><td>—</td><td>—</td></tr><tr><td><sup>135</sup>Cs/<sup>133</sup>Cs</td><td>0.3901 (0.0121)</td><td>0.3685 (0.0062)</td><td>0.3777 (0.0025)</td><td>—</td><td>—</td></tr><tr><td><sup>137</sup>(Cs + Ba<sub>radiogenic</sub>)/<sup>133</sup>Cs</td><td>0.9941 (0.0543)</td><td>0.9199 (0.0246)</td><td>0.9417 (0.0102)</td><td>—</td><td>—</td></tr><tr><td><sup>134</sup>(Cs + Ba<sub>radiogenic</sub>)/<sup>138</sup>Ba</td><td>36.6 (3.7)</td><td>35.3 (5.9)</td><td>137.7 (21.7)</td><td>0.03371</td><td>0.03371 (0.00009)</td></tr><tr><td><sup>136</sup>Ba/<sup>138</sup>Ba</td><td>0.1523 (0.0124)</td><td>0.1638 (0.0191)</td><td>0.4293 (0.0697)</td><td>0.1095</td><td>0.1097 (0.0002)</td></tr><tr><td><sup>87</sup>Rb/<sup>85</sup>Rb</td><td>2.407 (0.014)</td><td>2.284 (0.017)</td><td>2.236 (0.120)</td><td>0.3856</td><td>0.3857 (0.00009)</td></tr><tr><td><sup>90</sup>Sr</td><td>n/d</td><td>n/d</td><td>n/d</td><td>—</td><td>—</td></tr><tr><td><sup>41</sup>K/<sup>39</sup>K</td><td>0.07217 (0.00001)</td><td>—</td><td>—</td><td>0.07216</td><td>0.07214 (0.00003)</td></tr><tr><td><sup>40</sup>K/<sup>39</sup>K</td><td>0.000127 (0.000041)</td><td>—</td><td>—</td><td>0.000125</td><td>—</td></tr><tr><td><sup>43</sup>Ca/<sup>44</sup>Ca</td><td>0.06465 (0.00115)</td><td>—</td><td>—</td><td>0.06472</td><td>0.06469 (0.000005)</td></tr></tbody></table>

The elemental map and the EDX analysis indicate that the phase in bright contrast contains Fe, Zn, Cs, and Sn; whereas, the area of dark contrast corresponds to SiO$_2$-domains (Fig. S2b and c). The HRTEM image also reveals the presence of many nanocrystals (Fig. 5d). Based on the d-spacing and the FFT image, these nanocrystals are franklinite$^{23}$. Another magnified HAADF-STEM image with an elemental map reveals that the U distribution is closely associated with Zn-Fe-oxide nanoparticles (Fig. 5e). The U concentration is 0.8 wt% (including Si) or 1.1 wt% (excluding Si), assuming U is in the form of U$_3$O$_8$ for the quantitative EDX analysis (edx2, Fig. 5f).

Discussion

The U isotopic ratios of the ten scans on the CsMPs vary with a rather large deviation, ~0.0050 for OTZ3 and ~0.0030 for KOI2, most likely because the U concentration in the CsMPs is low (ppm level in the bulk). In addition, the analyses performed are not sufficiently precise to provide the absolute concentrations as high as that of the analysis on SRM610, because the irregular surfaces of the CsMPs aggregate were analyzed without polishing. The effect of isotopic fractionation should be negligibly small during the formation of these CsMPs. Even with a large deviation in the isotopic ratios, $^{235}$U/$^{238}$U in OTZ3 and KOI2 are significantly high, 0.029584 and 0.029341, as compared with the U isotopic ratio in nature, 0.00729. As expected, the U isotopic ratios are close to that of a typical nuclear fuel enriched in $^{235}$U.

The presence of $^{138}$Ba indicates either a fissionogenic or natural origin. Natural $^{136}$Ba/$^{138}$Ba, 0.1095, is an order of magnitude greater than $^{136}$Ba/$^{138}$Ba of fissionogenic origin, ~0.0126, estimated for Units #1–3 by calculation using the ORIGEN code$^{19}$, which largely varies depending on the burnup: 0.004 for 10 MWd/kgU, 0.011 for 20 MWd/kgU, and 0.014 for 30 MWd/kgU. $^{136}$Ba/$^{138}$Ba is 0.1523, 0.1638, and 0.4293 for OTZ3, KOI2, and OMR1, respectively. These values are much higher than those of fissionogenic Ba. The $^{136}$Ba/$^{138}$Ba ratio may indicate the contribution of a low-yield fission product; the amount of $^{136}$Cs (half-life of 13.16 days) can be remarkable in some cases depending on the amount of volatilized Ba. Although the amount of fissionogenic Ba isotopes is greater than Cs isotopes, the Cs can volatilize to greater extent than Ba during the meltdown, meaning that CsMP can selectively incorporate Cs as compared with Ba. Because $^{136}$Ba is a decay product of $^{136}$Cs (half-life of ~13 days), $^{136}$Ba was the most likely incorporated in the form of Cs species, leading to the higher value of $^{136}$Ba/$^{138}$Ba. On the other hand, $^{134}$(Cs + Ba)/$^{138}$Ba determined by SIMS varies from 35 to 138, which is much higher than $^{134}$Ba/$^{138}$Ba in nature (0.03371) and that of fissionogenic origin (0.0373–0.0566). Based on results from calculations using the ORIGEN code$^{19}$, the initial amount of fissionogenic $^{134}$Cs should be about the same as $^{134}$Ba. Thus, the high $^{134}$(Cs + Ba)/$^{138}$Ba is probably due to the higher volatility of Cs$^{26}$, as well as the greater amount of Cs incorporated into the CsMPs as compared with Ba. Thus, Ba isotopes of natural and fissionogenic origin that were initially incorporated in the CsMPs are negligibly small. Rather, $^{134}$Ba and $^{137}$Ba are decay daughters derived from Cs isotopes. $^{134}$Cs/$^{137}$Cs (γ) isotopic ratio at 15:36 March 12, 2011 was ~0.073, and based on the calculation of decayed amount of $^{134}$Cs and $^{137}$Cs corresponding to that of $^{134}$Ba and $^{137}$Ba in the CsMPs, the two isotopic
<figure><img src="image_4.png" /><figcaption>**Figure 3.** (a) A HAADF-STEM image of FIB-TEM specimen of OTZ3-1 CsMP. (b) Magnified HAADF-STEM image of the area indicated by the arrow in (a). (c) HAADF-STEM image (labeled as ADF) and the elemental map of the area displayed in (b). (d) STEM-EDX spectrums of the area indicated by the yellow square in (b). (e) Magnified HAADF-STEM image of the dark contrast zone represented as edx1 in (b) associated with elemental maps. (f) HRTEM image of the Zn-Fe-oxide nanoparticle. The lattice fringe of 2.9 Å corresponds to (220).</figcaption></figure>

ratios, ¹³⁴Cs/¹³⁴Ba<sub>radiogenic</sub>(γ) and ¹³⁷Cs/¹³⁷Ba<sub>radiogenic</sub>(γ) at the time of the SIMS analysis, were ~0.19 and ~7.47, respectively (Table 2).

Since the amount of volatilized Ba is negligible and the Ba in the CsMPs is mostly radiogenic, ¹³⁵Cs/¹³³Cs determined by SIMS represents the initial isotopic ratio, because ¹³³Cs is a stable isotope and ¹³⁵Cs has a long half-life, 2.3 million years. The previous ORIGEN calculation reported ¹³⁵Cs/¹³³Cs to be 0.388, 0.344, and 0.353 for Units # 1, 2, and 3, respectively, which are close to the present SIMS results, 0.39, 0.37, and 0.38 for OTZ3, KOI2, and OMR1, respectively. The isotopic ratios ¹³⁵Cs/¹³⁷(Cs + Ba<sub>radiogenic</sub>) and ¹³⁷(Cs + Ba<sub>radiogenic</sub>)/¹³³Cs calculated by SIMS are 0.39–0.40 and 0.92–0.99, respectively, which are comparable to the initial ¹³⁵Cs/¹³⁷Cs and ¹³⁷Cs/¹³³Cs values calculated by ORIGEN, 0.40 and 0.98, 0.34 and 1.01, 0.35 and 1.01, for Units #1, 2, and 3, respectively. Moreover, ¹³⁵Cs/¹³⁷(Cs + Ba<sub>radiogenic</sub>) is also comparable to ¹³⁵Cs/¹³⁷Cs of the bulk soil samples previously analyzed
<figure><img src="image_5.png" /><figcaption>**Figure 3.** (a) A HAADF-STEM image of FIB-TEM specimen of OTZ3-1 CsMP. (b) Magnified HAADF-STEM image of the area indicated by the arrow in (a). (c) HAADF-STEM image (labeled as ADF) and the elemental map of the area displayed in (b). (d) STEM-EDX spectrums of the area indicated by the yellow square in (b). (e) Magnified HAADF-STEM image of the dark contrast zone represented as edx1 in (b) associated with elemental maps. (f) HRTEM image of the Zn-Fe-oxide nanoparticle. The lattice fringe of 2.9 Å corresponds to (220).</figcaption></figure>

ratios, ¹³⁴Cs/¹³⁴Ba<sub>radiogenic</sub>(γ) and ¹³⁷Cs/¹³⁷Ba<sub>radiogenic</sub>(γ) at the time of the SIMS analysis, were ~0.19 and ~7.47, respectively (Table 2).

Since the amount of volatilized Ba is negligible and the Ba in the CsMPs is mostly radiogenic, ¹³⁵Cs/¹³³Cs determined by SIMS represents the initial isotopic ratio, because ¹³³Cs is a stable isotope and ¹³⁵Cs has a long half-life, 2.3 million years. The previous ORIGEN calculation reported ¹³⁵Cs/¹³³Cs to be 0.388, 0.344, and 0.353 for Units # 1, 2, and 3, respectively, which are close to the present SIMS results, 0.39, 0.37, and 0.38 for OTZ3, KOI2, and OMR1, respectively. The isotopic ratios ¹³⁵Cs/¹³⁷(Cs + Ba<sub>radiogenic</sub>) and ¹³⁷(Cs + Ba<sub>radiogenic</sub>)/¹³³Cs calculated by SIMS are 0.39–0.40 and 0.92–0.99, respectively, which are comparable to the initial ¹³⁵Cs/¹³⁷Cs and ¹³⁷Cs/¹³³Cs values calculated by ORIGEN, 0.40 and 0.98, 0.34 and 1.01, 0.35 and 1.01, for Units #1, 2, and 3, respectively. Moreover, ¹³⁵Cs/¹³⁷(Cs + Ba<sub>radiogenic</sub>) is also comparable to ¹³⁵Cs/¹³⁷Cs of the bulk soil samples previously analyzed
<figure><img src="image_5.png" /><figcaption>**Figure 4.** (a) HAADF-STEM image of the FIB-TEM specimen OTZ3-2 associated with the elemental maps. (b) A magnified HAADF-STEM image of the high Cs phase indicated by the arrow in (a). (c) A representative EDX spectrum of the high Cs phase indicated by the arrow in (a). (d) HRTEM image of the high Cs phase indicated by the arrow in (a). (e) The SAED pattern of the high Cs phase from two different major zone axes.</figcaption></figure>

using thermal ionization mass spectrometry (TIMS), 0.36–0.38²⁷,²⁸. Because of the high ¹³⁷Cs/¹³⁷Ba<sub>radiogenic</sub>, the majority of 137-mass derives from ¹³⁷Cs, and the efficiency of secondary ions formation for ¹³⁷Ba should be minimal.

The values of ¹³⁴(Cs + Ba<sub>radiogenic</sub>)/¹³⁷(Cs + Ba<sub>radiogenic</sub>) (γ) determined through gamma spectrometry (0.0724–0.0742), which are the same as those of the initial ¹³⁴Cs/¹³⁷Cs, are higher than the values of ¹³⁴(Cs + Ba)/¹³⁷(Cs + Ba) for CsMPs that were determined by SIMS (0.0249–0.0318). This is attributed to the difference in ionization efficiency between Cs and Ba. Introducing the secondary ion production efficiency coefficient, ¹f=¹Ba ionization efficiency/¹Cs ionization efficiency, where i stands for the mass, the ratio by SIMS can be expressed as

$$ \begin{aligned} & (^{134}\text{Cs} + ^{134}\text{Ba}_{\text{radiogenic}} \times ^{134}\text{f})/(^{137}\text{Cs} + ^{137}\text{Ba}_{\text{radiogenic}} \times ^{137}\text{f}) \\ & = 0.029 \text{ (average for three CsMPs)} \end{aligned} \qquad (1) $$

Because the difference between ¹³⁴f and ¹³⁷f is negligible, the ratio of the ionization efficiency is simply expressed as f, and equation (1) can be rewritten as

$$ \begin{aligned} & (^{134}\text{Cs} + ^{134}\text{Ba}_{\text{radiogenic}} \times f)/(^{137}\text{Cs} + ^{137}\text{Ba}_{\text{radiogenic}} \times f) \\ & = (^{134}\text{Cs}/^{137}\text{Cs})_{\text{initial}}(0.16 + 0.84 \times f)/(0.88 + 0.12 \times f) \\ & = 0.029 \end{aligned} \qquad (2) $$
<figure><img src="image_6.png" /><figcaption>**Figure 4.** (a) HAADF-STEM image of the FIB-TEM specimen OTZ3-2 associated with the elemental maps. (b) A magnified HAADF-STEM image of the high Cs phase indicated by the arrow in (a). (c) A representative EDX spectrum of the high Cs phase indicated by the arrow in (a). (d) HRTEM image of the high Cs phase indicated by the arrow in (a). (e) The SAED pattern of the high Cs phase from two different major zone axes.</figcaption></figure>

using thermal ionization mass spectrometry (TIMS), 0.36–0.38²⁷,²⁸. Because of the high ¹³⁷Cs/¹³⁷Ba<sub>radiogenic</sub>, the majority of 137-mass derives from ¹³⁷Cs, and the efficiency of secondary ions formation for ¹³⁷Ba should be minimal.

The values of ¹³⁴(Cs + Ba<sub>radiogenic</sub>)/¹³⁷(Cs + Ba<sub>radiogenic</sub>) (γ) determined through gamma spectrometry (0.0724–0.0742), which are the same as those of the initial ¹³⁴Cs/¹³⁷Cs, are higher than the values of ¹³⁴(Cs + Ba)/¹³⁷(Cs + Ba) for CsMPs that were determined by SIMS (0.0249–0.0318). This is attributed to the difference in ionization efficiency between Cs and Ba. Introducing the secondary ion production efficiency coefficient, ¹f = ¹Ba ionization efficiency/¹Cs ionization efficiency, where i stands for the mass, the ratio by SIMS can be expressed as

$$ \begin{aligned} & (^{134}\text{Cs} + ^{134}\text{Ba}_{\text{radiogenic}} \times ^{134}\text{f})/(^{137}\text{Cs} + ^{137}\text{Ba}_{\text{radiogenic}} \times ^{137}\text{f}) \\ & = 0.029 \quad (\text{average for three CsMPs}) \end{aligned} \qquad (1) $$

Because the difference between ¹³⁴f and ¹³⁷f is negligible, the ratio of the ionization efficiency is simply expressed as f, and equation (1) can be rewritten as

$$ \begin{aligned} & (^{134}\text{Cs} + ^{134}\text{Ba}_{\text{radiogenic}} \times f)/(^{137}\text{Cs} + ^{137}\text{Ba}_{\text{radiogenic}} \times f) \\ & = (^{134}\text{Cs}/^{137}\text{Cs})_{\text{initial}}(0.16 + 0.84 \times f)/(0.88 + 0.12 \times f) \\ & = 0.029 \end{aligned} \qquad (2) $$
<figure><img src="image_6.png" /><figcaption>**Figure 5.** (a) HAADF-STEM image of the FIB-TEM specimen of the OTZ10 CsMP with the SAED as inset. (b) Elemental maps of the same area as the HAADF-STEM image in (a). (c) The EDX spectrum of the area indicated by the yellow square labeled as edx1 in (a). (d) A HRTEM image of the Zn-Fe-oxides nanoparticles accompanied with the magnified lattice image and the FFT image of the selected area indicated by the white square. The lattice fringe of 2.2 Å corresponds to (400). (e) HAADF-STEM image of the enlarged figure of the thin edge in (a) associated with elemental maps. (f) STEM-EDX point analysis on a Zn-Fe-oxide nanoparticle indicated by a red cross labeled as edx2 in (e).</figcaption></figure>

Solving equation (2): $f=0.24$.

This ionization efficiency factor is reasonably close to that determined for rhyolitic glass (SRM610); 0.288. Note, the ionization efficiency of the standard specimen is determined on the polished flat surface on the pure glass, while that for the CsMP was determined on the spherical shaped particle that is a mixture of SiO₂ glass and the Zn-Fe-oxides nanoparticles. Thus, the matrix composition is similar but not identical.

The Rb isotopic ratio, ⁸⁷Rb/⁸⁵Rb, is 2.2–2.4, which is obviously higher than the natural ⁸⁷Rb/⁸⁵Rb (0.3856), indicating that the Rb isotopes are fissionogenic. Compared with the ⁸⁷Rb/⁸⁵Rb values calculated by ORIGEN¹⁹ (~1.7), ⁸⁷Rb/⁸⁵Rb obtained by SIMS is relatively high. Nevertheless, our OrigenArp calculation gives ratio values between 2.5 and 2.7, which are closer to the ⁸⁷Rb/⁸⁵Rb values by SIMS. The difference is not ascribed to the mixture of natural Rb isotopes because the mixing with the natural isotopes would lower the ⁸⁷Rb/⁸⁵Rb isotopic ratio, but rather due to the local volatilization of Rb.

The nanoscale texture of OTZ10 and the low-Cs zone in OTZ3-1 reveal that franklinite nanoparticles are embedded in a glassy matrix of SiO₂. These textures and the Cs association with franklinite were also reported in previous studies¹⁵, ¹⁶, indicating multiple steps in the formation sequence: radioactive Cs is released to form nanoparticles and/or is present as droplets in a mist during meltdown; subsequently, numerous Zn-Fe-oxide
<figure><img src="image_7.png" /><figcaption>**Figure 5.** (a) HAADF-STEM image of the FIB-TEM specimen of the OTZ10 CsMP with the SAED as inset. (b) Elemental maps of the same area as the HAADF-STEM image in (a). (c) The EDX spectrum of the area indicated by the yellow square labeled as edx1 in (a). (d) A HRTEM image of the Zn-Fe-oxides nanoparticles accompanied with the magnified lattice image and the FFT image of the selected area indicated by the white square. The lattice fringe of 2.2 Å corresponds to (400). (e) HAADF-STEM image of the enlarged figure of the thin edge in (a) associated with elemental maps. (f) STEM-EDX point analysis on a Zn-Fe-oxide nanoparticle indicated by a red cross labeled as edx2 in (e).</figcaption></figure>

Solving equation (2): $f=0.24$.

This ionization efficiency factor is reasonably close to that determined for rhyolitic glass (SRM610); 0.288. Note, the ionization efficiency of the standard specimen is determined on the polished flat surface on the pure glass, while that for the CsMP was determined on the spherical shaped particle that is a mixture of SiO₂ glass and the Zn-Fe-oxides nanoparticles. Thus, the matrix composition is similar but not identical.

The Rb isotopic ratio, ⁸⁷Rb/⁸⁵Rb, is 2.2–2.4, which is obviously higher than the natural ⁸⁷Rb/⁸⁵Rb (0.3856), indicating that the Rb isotopes are fissionogenic. Compared with the ⁸⁷Rb/⁸⁵Rb values calculated by ORIGEN¹⁹ (~1.7), ⁸⁷Rb/⁸⁵Rb obtained by SIMS is relatively high. Nevertheless, our OrigenArp calculation gives ratio values between 2.5 and 2.7, which are closer to the ⁸⁷Rb/⁸⁵Rb values by SIMS. The difference is not ascribed to the mixture of natural Rb isotopes because the mixing with the natural isotopes would lower the ⁸⁷Rb/⁸⁵Rb isotopic ratio, but rather due to the local volatilization of Rb.

The nanoscale texture of OTZ10 and the low-Cs zone in OTZ3-1 reveal that franklinite nanoparticles are embedded in a glassy matrix of SiO₂. These textures and the Cs association with franklinite were also reported in previous studies¹⁵, ¹⁶, indicating multiple steps in the formation sequence: radioactive Cs is released to form nanoparticles and/or is present as droplets in a mist during meltdown; subsequently, numerous Zn-Fe-oxide
(franklinite) nanoparticles formed during the failure of the reactor pressure vessel (RPV), and Cs in wet form is adsorbed onto the Zn-Fe-oxide nanoparticles; then, the melted fuels break through the RPV, hit the concrete pedestal, and generate SiO gas at >2,000 K. This process is known as the molten core-concrete interaction (MCCI)²⁹, which immediately leads to the condensation of SiO₂ over the Zn-Fe-oxide nanoparticles and the incorporation of fission product nanoparticles. As shown in a previous study¹⁶, a low concentration of U (~1 wt%) is observed associated with franklinite nanoparticles that formed before the CsMP formation. Moreover, U adsorption onto the franklinite nanoparticles possibly occurs through volatilization of the oxidized form of the UO₂ fuel. In fact, UO₂ pellets can be oxidized depending on the H₂O:H₂ ratio³⁰. In general, U is considered a nonvolatile element; however, the partially oxidized form of UO₂+x can volatilize by ~10% at >~1,900 K at a heating rate of ~30 K/min; whereas, the unoxidized form UO₂ does not become volatile even at 2,700 K³¹, ³². Uranium association with microscale aerosols during a severe accident was also recognized in some experimental studies³³. The chemical form of U inside CsMPs has been determined in the OTZ10 CsMP, as well as in the other CsMP separated from the same soil sample¹⁶. The U isotopic ratio in another CsMP, OTZ3, found in the same soil sample, was determined. ²³⁵U/²³⁸U = ~0.030, as analyzed from the two CsMPs, indicating that the U incorporated in the CsMPs is not from naturally occurring trace U that has a ²³⁵U:²³⁸U ratio of 0.00729:1, which was initially considered as the source of U because trace amounts of U are included in the insulating materials surrounding reactor pressure vessel. Rather, nuclear fuel is the most likely origin of U in the CsMPs. The original isotope composition of non-irradiated fuel before operation at the FDNPP was 0.0389¹⁹. Based on calculations using the ORIGEN code, the U isotopic ratios for the irradiated fuel in Units # 1, 2 and 3 were 0.0172, 0.0193, and 0.0192¹⁹, respectively. Thus, the isotopic ratios are between the isotopic ratio calculated from burnup and that of the non-irradiated fuel. The isotopic ratios calculated by the ORIGEN code are based on the average burnup for all fuel assemblies. The intermediate values of ²³⁵U/²³⁸U are due to the fact that the burnup and the temperature are not homogeneous within the reactor core. A fuel assembly will experience different temperatures depending on its position in the core. Even within a single pellet, there is a thermal gradient; the center of the pellet can be heated as high as ~1,973 K, whereas the temperature of the rim is only ~673 K³⁴. Moreover, the burnup within a single pellet is not uniform with low burnup at the center and high burnup at the edge³⁴. The burnup of the fuel assembly is not uniform either, as the fuel rods with different burnup are assembled in a certain matter. Thus, the fuel rods with relatively lower burnup might have been heated at the temperature at which a small amount of U was volatilized, resulting in the high ²³⁵U/²³⁸U isotopic ratios for both CsMPs. The volatilized U is associated with Zn-Fe-oxide, most likely by simple adsorption, as evidenced by the U in the OTZ10 CsMP.

Dispersion of U released from the FDNPP has been reported in previous studies³⁵⁻³⁷, including ²³⁶U/²³⁸U values of ~10⁻⁹ in paddy-field water and ocean water as far as ~30 km from the FDNPP³⁵. ²³⁶U/²³⁸U values of ~10⁻⁷ were also reported in black-colored dust samples collected within a ~30-km distance³⁶. This is evidence that a small amount of U, about 150 g³⁶, was discharged from the fuels in the reactors. Nevertheless, another study measuring ²³⁵U/²³⁸U failed to identify the FDNPP as the origin due to the dilution by natural U isotopes³⁷. However, all those analyses were performed on bulk soil samples and did not include data from single CsMPs. Thus, the speciation of U was never determined. The present study convincingly identifies the U speciation as associated with Zn-Fe-oxide nanoparticles at ~1 wt%, and the nuclear fuel origin of U in CsMPs is confirmed through isotopic ratio analyses. Yet, quantitative analysis of the amount of U in the form of CsMPs in the total U that was dispersed to environment remains undetermined.

Two different occurrences of radioactive Cs have been identified: Fe-pollucite and Cs associated with franklinite nanoparticles embedded in the glassy SiO₂ matrix. The latter has been already described above and in previous studies¹⁵, ¹⁶ and is not discussed further here. On the other hand, Fe-pollucite, CsFeSi₂O₆·nH₂O, is identified for the first time in CsMPs. None of the previous studies reported pollucite during meltdown; however, a recent experimental study reported formation of pollucite and CsFeSiO₄ during the CsOH chemisorption onto stainless steel containing 5% Si³⁸. Both materials have a zeolite structure that can be synthesized by annealing a mixture of gels³⁹ and appropriate elements, such as Cs⁴⁰. Their formation by chemisorption likely occurred at the high-Cs zone during the formation of CsMPs. However, the reaction required for CsMPs formation are apparently different from the simple interaction of CsOH with stainless steel at ~1,273 K³⁸. In fact, the Cs-Cr phases that result from the reaction with stainless steel are not present in CsMPs. Rather, it is probable that CsOH is involved in the chemical reactions with Si and Fe oxide during the MCCI at >2,000 K. The clear boundary between the high- and low-Cs zones (Fig. 3b) indicates that Zn-Fe-oxide nanoparticles formed aggregates prior to the formation of Fe-pollucite, which formed at the time of MCCI.

In summary, the present study successfully determined the origin of U, Cs, Ba, Rb, K, and Ca in the CsMPs, as well as their chemical form based on the isotopic ratios and their nanoscale structure. In particular, the data obtained are critically important for delineating the reactions involving nuclear fuels, which occurred inside the reactor during the formation of the CsMPs. The detailed information on the chemical state and source provide an understanding of the source term at the FDNPP during the meltdown event. This is important information for severe accident (SA) analysis codes, such as MELCOR⁴¹, ⁴² and MAAP⁴³, since some of the phases identified in the CsMPs, such as Zn-Fe-oxide and pollucite, are not considered in these severe accident codes⁴⁴.

## Methods

**Sampling.** We collected the samples at three locations (Fig. 1). As for the sample labels, the capital letters represent the locality and the following number indicate the number in our CsMP database in Kyushu University. The OTZ3 and OTZ10 were separated from the same soil collected from the top ~1 cm of the paddy soil at Ottozawa located ~4 km west from the FDNPP in Okuma Town, Futaba County, Fukushima, on March 16, 2012. The soil was mainly composed of clay minerals, quartz, and feldspar⁹. Because it is still prohibited to enter the area due to the high radiation dose, the locality remained untouched and was not disturbed by decontamination or restoration. The radiation dose ~1 m above the ground was measured at 84 µSv/h.
(franklinite) nanoparticles formed during the failure of the reactor pressure vessel (RPV), and Cs in wet form is adsorbed onto the Zn-Fe-oxide nanoparticles; then, the melted fuels break through the RPV, hit the concrete pedestal, and generate SiO gas at >2,000 K. This process is known as the molten core-concrete interaction (MCCI)²⁹, which immediately leads to the condensation of SiO₂ over the Zn-Fe-oxide nanoparticles and the incorporation of fission product nanoparticles. As shown in a previous study¹⁶, a low concentration of U (~1 wt%) is observed associated with franklinite nanoparticles that formed before the CsMP formation. Moreover, U adsorption onto the franklinite nanoparticles possibly occurs through volatilization of the oxidized form of the UO₂ fuel. In fact, UO₂ pellets can be oxidized depending on the H₂O:H₂ ratio³⁰. In general, U is considered a nonvolatile element; however, the partially oxidized form of UO₂+x can volatilize by ~10% at >~1,900 K at a heating rate of ~30 K/min; whereas, the unoxidized form UO₂ does not become volatile even at 2,700 K³¹, ³². Uranium association with microscale aerosols during a severe accident was also recognized in some experimental studies³³. The chemical form of U inside CsMPs has been determined in the OTZ10 CsMP, as well as in the other CsMP separated from the same soil sample¹⁶. The U isotopic ratio in another CsMP, OTZ3, found in the same soil sample, was determined. ²³⁵U/²³⁸U = ~0.030, as analyzed from the two CsMPs, indicating that the U incorporated in the CsMPs is not from naturally occurring trace U that has a ²³⁵U:²³⁸U ratio of 0.00729:1, which was initially considered as the source of U because trace amounts of U are included in the insulating materials surrounding reactor pressure vessel. Rather, nuclear fuel is the most likely origin of U in the CsMPs. The original isotope composition of non-irradiated fuel before operation at the FDNPP was 0.0389¹⁹. Based on calculations using the ORIGEN code, the U isotopic ratios for the irradiated fuel in Units # 1, 2 and 3 were 0.0172, 0.0193, and 0.0192¹⁹, respectively. Thus, the isotopic ratios are between the isotopic ratio calculated from burnup and that of the non-irradiated fuel. The isotopic ratios calculated by the ORIGEN code are based on the average burnup for all fuel assemblies. The intermediate values of ²³⁵U/²³⁸U are due to the fact that the burnup and the temperature are not homogeneous within the reactor core. A fuel assembly will experience different temperatures depending on its position in the core. Even within a single pellet, there is a thermal gradient; the center of the pellet can be heated as high as ~1,973 K, whereas the temperature of the rim is only ~673 K³⁴. Moreover, the burnup within a single pellet is not uniform with low burnup at the center and high burnup at the edge³⁴. The burnup of the fuel assembly is not uniform either, as the fuel rods with different burnup are assembled in a certain matter. Thus, the fuel rods with relatively lower burnup might have been heated at the temperature at which a small amount of U was volatilized, resulting in the high ²³⁵U/²³⁸U isotopic ratios for both CsMPs. The volatilized U is associated with Zn-Fe-oxide, most likely by simple adsorption, as evidenced by the U in the OTZ10 CsMP.

Dispersion of U released from the FDNPP has been reported in previous studies³⁵⁻³⁷, including ²³⁶U/²³⁸U values of ~10⁻⁹ in paddy-field water and ocean water as far as ~30 km from the FDNPP³⁵. ²³⁶U/²³⁸U values of ~10⁻⁷ were also reported in black-colored dust samples collected within a ~30-km distance³⁶. This is evidence that a small amount of U, about 150 g³⁶, was discharged from the fuels in the reactors. Nevertheless, another study measuring ²³⁵U/²³⁸U failed to identify the FDNPP as the origin due to the dilution by natural U isotopes³⁷. However, all those analyses were performed on bulk soil samples and did not include data from single CsMPs. Thus, the specification of U was never determined. The present study convincingly identifies the U speciation as associated with Zn-Fe-oxide nanoparticles at ~1 wt%, and the nuclear fuel origin of U in CsMPs is confirmed through isotopic ratio analyses. Yet, quantitative analysis of the amount of U in the form of CsMPs in the total U that was dispersed to environment remains undetermined.

Two different occurrences of radioactive Cs have been identified: Fe-pollucite and Cs associated with franklinite nanoparticles embedded in the glassy SiO₂ matrix. The latter has been already described above and in previous studies¹⁵, ¹⁶ and is not discussed further here. On the other hand, Fe-pollucite, CsFeSi₂O₆·nH₂O, is identified for the first time in CsMPs. None of the previous studies reported pollucite during meltdown; however, a recent experimental study reported formation of pollucite and CsFeSiO₄ during the CsOH chemisorption onto stainless steel containing 5% Si³⁸. Both materials have a zeolite structure that can be synthesized by annealing a mixture of gels³⁹ and appropriate elements, such as Cs⁴⁰. Their formation by chemisorption likely occurred at the high-Cs zone during the formation of CsMPs. However, the reaction required for CsMPs formation are apparently different from the simple interaction of CsOH with stainless steel at ~1,273 K³⁸. In fact, the Cs-Cr phases that result from the reaction with stainless steel are not present in CsMPs. Rather, it is probable that CsOH is involved in the chemical reactions with Si and Fe oxide during the MCCI at >2,000 K. The clear boundary between the high- and low-Cs zones (Fig. 3b) indicates that Zn-Fe-oxide nanoparticles formed aggregates prior to the formation of Fe-pollucite, which formed at the time of MCCI.

In summary, the present study successfully determined the origin of U, Cs, Ba, Rb, K, and Ca in the CsMPs, as well as their chemical form based on the isotopic ratios and their nanoscale structure. In particular, the data obtained are critically important for delineating the reactions involving nuclear fuels, which occurred inside the reactor during the formation of the CsMPs. The detailed information on the chemical state and source provide an understanding of the source term at the FDNPP during the meltdown event. This is important information for severe accident (SA) analysis codes, such as MELCOR⁴¹, ⁴² and MAAP⁴³, since some of the phases identified in the CsMPs, such as Zn-Fe-oxide and pollucite, are not considered in these severe accident codes⁴⁴.

## Methods

**Sampling.** We collected the samples at three locations (Fig. 1). As for the sample labels, the capital letters represent the locality and the following number indicate the number in our CsMP database in Kyushu University. The OTZ3 and OTZ10 were separated from the same soil collected from the top ~1 cm of the paddy soil at Ottozawa located ~4 km west from the FDNPP in Okuma Town, Futaba County, Fukushima, on March 16, 2012. The soil was mainly composed of clay minerals, quartz, and feldspar⁹. Because it is still prohibited to enter the area due to the high radiation dose, the locality remained untouched and was not disturbed by decontamination or restoration. The radiation dose ~1 m above the ground was measured at 84 µSv/h.
The second sample (KOI2) was composed of gravels collected under the drainpipe of the assembly house in Koirino located 2.9 km southwest from the FDNPP during the same sampling campaign. The radiation dose underneath the drainpipe was extremely high compared to the surroundings; the dose at the sampling point was as high as 630 µSv/h. The gravel samples were carefully collected from the ground surface using a hand shovel and placed in plastic bags. The soil was also mainly composed of clay minerals, quartz, and feldspar. The third sample (OMR1) was collected under the drainpipe of the warehouse in Omaru located in Namie Town, Futaba County, Fukushima, ~10.5 km northwest from the FDNPP during the sampling campaign on December 20, 2012. The radiation dose ~1 m above the ground exceeded 30 µSv/h.

**Separation of CsMPs.** Prior to the procedure, both samples were sieved through a 114-µm mesh. The powder samples were dispersed on grid paper and then covered with a plastic sheet. Next, an imaging plate (IP, Fuji film, BAS-SR 2025) was placed on the samples for 5–15 min. Further, using an IP reader, autoradiograph images with pixel sizes of 100 µm were recorded. After identifying the positions of intense radioactive spots, droplets of pure water were added to these positions and then drawn using a pipette to produce suspensions with small amounts of soil particles by dilution with pure water. This procedure was repeated until the suspension did not contain a significant amount of soil particles. Subsequently, the positions containing hot spots were selected using pieces of double-stick carbon tape that were cut as small as possible with a blade. The pieces were checked by autoradiograph imaging to obtain CsMPs with maximum efficiency using scanning electron microscopy (SEM) observations. Prior to SEM analysis, the pieces were placed on an aluminum plate and coated with carbon using a carbon coater (SANYU SC-701C). The CsMPs were found and observed using two SEMs (Shimadzu, SS550 and Hitachi, SU6600) using acceleration voltages of 15–25 kV both equipped with an energy dispersive X-ray spectrometer (EDX, EDAX Genesis).

**Preparation of TEM specimens.** A focused ion beam (FIB) instrument (FEI, Quanta 3D FEG 200i Dual Beam) was utilized to prepare a thin foil of individual CsMPs with diameters of a few µm. Gallium was used as an ion source, and W deposition was used to minimize the damage from the ion bombardment. The current and accelerating voltage of the ion beam were adjusted from 100 pA to 30 nA and 5–30 kV, respectively, depending on the progress of the thinning and sample properties such as hardness and size. Each thinned piece was attached to the semilunar-shaped Cu grid for FIB and further thinned by an ion beam operating at 5 kV.

**TEM analysis.** HRTEM with EDX and high-angle annular dark-field scanning transmission electron microscopy were performed using a JEOL JEM-ARM200F and JEM-ARM200CF with an acceleration voltage of 200 kV. The JEOL Analysis Station software was used to control the STEM-EDX mapping. To minimize the effects of sample drift, a drift-correction mode was used during acquisition of the elemental map. The STEM probe size was ~0.13 nm, generating an ~140-pA current when 40 µm of the condenser lens aperture was inserted. The collection angle of the HAADF detector was ~97–256 mrad. A Gatan imaging filter system was used to conduct electron energy-loss spectroscopy in STEM mode. The convergent and collection semi-angles were ~30 mrad and ~17 mrad, respectively.

**Gamma spectrometry.** The ¹³⁴Cs and ¹³⁷Cs radioactivities of the CsMPs were determined using gamma spectrometry. The radioactivity of an additional microparticle with a size of ~400 µm obtained from surface soil in Fukushima was precisely determined at the radioisotope center in Tsukuba University, Japan, and utilized as a standard point specimen for ¹³⁴Cs and ¹³⁷Cs. The radioactivity of the point source standard was 23.9 Bq for ¹³⁴Cs and 94.6 Bq for ¹³⁷Cs as of September 29, 2015. Measurement of the radioactivity was performed on the CsMPs and the point source standard using a germanium semi-conductor detector, GMX23 (SEIKO E&G), GMX40 (SEIKO E&G) and GX6020 (Canberra) at the center for radioisotopes in Kyushu University, Japan. The acquisition times were 21,238 s for OMR1 and 84,289 s for OTZ10 using GMX23, 4,004 s for OTZ3 using GMX40, and 77, 880.8 s for KOI2 using GX6020.

**Secondary ion mass spectrometry.** Isotopic ratio analysis was performed using secondary ion mass spectrometry (SIMS, SHRIMP-II, Australian Scientific Instruments) at the National Institute of Polar Research, Tokyo. The SIMS specimens were put on Al plates or a Cu grid and fixed on 1-inch slide glass by Cu tape. They were coated with Au at a 13.5-nm thickness prior to the analysis. An O₂⁻ primary ion beam of 0.2–0.4 nA was used to sputter the specimen surface with a beam diameter of 5.0–7.0 µm. The typical mass resolution is about 4,500 (M/ΔM at 1% of peak height). The National Institute of Standards and Technology (NIST) SRM610 was used as a standard specimen, in which 461.5 mg/kg of depleted U was doped in a silicate glass matrix.

Because quantitative analysis using SIMS requires suitable standards, which have chemical compositions similar to those of the target minerals, the absolute concentrations were not obtained, but only the isotopic ratios were determined. Ten scans were completed for each analytical spot on the CsMPs, meaning that the sequence of ten analyses represents a depth profile of the variation in the isotopic ratios; thus, SIMS analyses provide the isotopic signatures inside the CsMPs. The average values are given in the Table 2 with the standard deviation calculated for the ten analyses.

**References**

1. Steinhauser, G., Brandl, A. & Johnson, T. E. Comparison of the Chernobyl and Fukushima nuclear accidents: A review of the environmental impacts. *Sci. Total Environ.* **470–471**, 800–817 (2014).

2. Steinhauser, G. Fukushima's forgotten radionuclides: A review of the understudied radioactive emissions. *Environ. Sci. Technol.* **48**, 4649–4663 (2014).

3. Koo, Y. H., Yang, Y. S. & Song, K. W. Radioactivity release from the Fukushima accident and its consequences: A review. *Prog. Nucl. Energy* **74**, 61–70 (2014).
The second sample (KOI2) was composed of gravels collected under the drainpipe of the assembly house in Koirino located 2.9 km southwest from the FDNPP during the same sampling campaign. The radiation dose underneath the drainpipe was extremely high compared to the surroundings; the dose at the sampling point was as high as 630 µSv/h. The gravel samples were carefully collected from the ground surface using a hand shovel and placed in plastic bags. The soil was also mainly composed of clay minerals, quartz, and feldspar. The third sample (OMR1) was collected under the drainpipe of the warehouse in Omaru located in Namie Town, Futaba County, Fukushima, ~10.5 km northwest from the FDNPP during the sampling campaign on December 20, 2012. The radiation dose ~1 m above the ground exceeded 30 µSv/h.

**Separation of CsMPs.** Prior to the procedure, both samples were sieved through a 114-µm mesh. The powder samples were dispersed on grid paper and then covered with a plastic sheet. Next, an imaging plate (IP, Fuji film, BAS-SR 2025) was placed on the samples for 5–15 min. Further, using an IP reader, autoradiograph images with pixel sizes of 100 µm were recorded. After identifying the positions of intense radioactive spots, droplets of pure water were added to these positions and then drawn using a pipette to produce suspensions with small amounts of soil particles by dilution with pure water. This procedure was repeated until the suspension did not contain a significant amount of soil particles. Subsequently, the positions containing hot spots were selected using pieces of double-stick carbon tape that were cut as small as possible with a blade. The pieces were checked by autoradiograph imaging to obtain CsMPs with maximum efficiency using scanning electron microscopy (SEM) observations. Prior to SEM analysis, the pieces were placed on an aluminum plate and coated with carbon using a carbon coater (SANYU SC-701C). The CsMPs were found and observed using two SEMs (Shimadzu, SS550 and Hitachi, SU6600) using acceleration voltages of 15–25 kV both equipped with an energy dispersive X-ray spectrometer (EDX, EDAX Genesis).

**Preparation of TEM specimens.** A focused ion beam (FIB) instrument (FEI, Quanta 3D FEG 200i Dual Beam) was utilized to prepare a thin foil of individual CsMPs with diameters of a few µm. Gallium was used as an ion source, and W deposition was used to minimize the damage from the ion bombardment. The current and accelerating voltage of the ion beam were adjusted from 100 pA to 30 nA and 5–30 kV, respectively, depending on the progress of the thinning and sample properties such as hardness and size. Each thinned piece was attached to the semilunar-shaped Cu grid for FIB and further thinned by an ion beam operating at 5 kV.

**TEM analysis.** HRTEM with EDX and high-angle annular dark-field scanning transmission electron microscopy were performed using a JEOL JEM-ARM200F and JEM-ARM200CF with an acceleration voltage of 200 kV. The JEOL Analysis Station software was used to control the STEM-EDX mapping. To minimize the effects of sample drift, a drift-correction mode was used during acquisition of the elemental map. The STEM probe size was ~0.13 nm, generating an ~140-pA current when 40 µm of the condenser lens aperture was inserted. The collection angle of the HAADF detector was ~97–256 mrad. A Gatan imaging filter system was used to conduct electron energy-loss spectroscopy in STEM mode. The convergent and collection semi-angles were ~30 mrad and ~17 mrad, respectively.

**Gamma spectrometry.** The ¹³⁴Cs and ¹³⁷Cs radioactivities of the CsMPs were determined using gamma spectrometry. The radioactivity of an additional microparticle with a size of ~400 µm obtained from surface soil in Fukushima was precisely determined at the radioisotope center in Tsukuba University, Japan, and utilized as a standard point specimen for ¹³⁴Cs and ¹³⁷Cs. The radioactivity of the point source standard was 23.9 Bq for ¹³⁴Cs and 94.6 Bq for ¹³⁷Cs as of September 29, 2015. Measurement of the radioactivity was performed on the CsMPs and the point source standard using a germanium semi-conductor detector, GMX23 (SEIKO E&G), GMX40 (SEIKO E&G) and GX6020 (Canberra) at the center for radioisotopes in Kyushu University, Japan. The acquisition times were 21,238 s for OMR1 and 84,289 s for OTZ10 using GMX23, 4,004 s for OTZ3 using GMX40, and 77, 880.8 s for KOI2 using GX6020.

**Secondary ion mass spectrometry.** Isotopic ratio analysis was performed using secondary ion mass spectrometry (SIMS, SHRIMP-II, Australian Scientific Instruments) at the National Institute of Polar Research, Tokyo. The SIMS specimens were put on Al plates or a Cu grid and fixed on 1-inch slide glass by Cu tape. They were coated with Au at a 13.5-nm thickness prior to the analysis. An O₂⁻ primary ion beam of 0.2–0.4 nA was used to sputter the specimen surface with a beam diameter of 5.0–7.0 µm. The typical mass resolution is about 4,500 (M/ΔM at 1% of peak height). The National Institute of Standards and Technology (NIST) SRM610 was used as a standard specimen, in which 461.5 mg/kg of depleted U was doped in a silicate glass matrix.

Because quantitative analysis using SIMS requires suitable standards, which have chemical compositions similar to those of the target minerals, the absolute concentrations were not obtained, but only the isotopic ratios were determined. Ten scans were completed for each analytical spot on the CsMPs, meaning that the sequence of ten analyses represents a depth profile of the variation in the isotopic ratios; thus, SIMS analyses provide the isotopic signatures inside the CsMPs. The average values are given in the Table 2 with the standard deviation calculated for the ten analyses.

**References**

1. Steinhauser, G., Brandl, A. & Johnson, T. E. Comparison of the Chernobyl and Fukushima nuclear accidents: A review of the environmental impacts. *Sci. Total Environ.* **470–471**, 800–817 (2014).

2. Steinhauser, G. Fukushima's forgotten radionuclides: A review of the understudied radioactive emissions. *Environ. Sci. Technol.* **48**, 4649–4663 (2014).

3. Koo, Y. H., Yang, Y. S. & Song, K. W. Radioactivity release from the Fukushima accident and its consequences: A review. *Prog. Nucl. Energy* **74**, 61–70 (2014).
4. Kinoshita, N. et al. Assessment of individual radionuclide distributions from the Fukushima nuclear accident covering central-east Japan. *Proc. Natl. Acad. Sci. USA* **108**, 19526–19529 (2011).

5. Morino, Y., Ohara, T. & Nishizawa, M. Atmospheric behavior, deposition, and budget of radioactive materials from the Fukushima Daiichi nuclear power plant in March 2011. *Geophys. Res. Lett.* **38**, L00G11 (2011).

6. Buesseler, K. et al. Fukushima Daiichi-derived radionuclides in the ocean: transport, fate, and impacts. *Annu. Rev. Mar. Sci.* **9**, 173–203 (2017).

7. Grambow, B. & Mostafavi, M. State of Fukushima nuclear fuel debris tracked by Cs137 in cooling water. *Environ. Sci.: Processes Impacts* **16**, 2472–2476 (2014).

8. Kaneyasu, N. et al. Sulfate aerosol as a potential transport medium of radiocesium from the Fukushima nuclear accident. *Environ. Sci. Technol.* **46**, 5720–5726 (2012).

9. Kaneko, M. et al. Radioactive Cs in the severely contaminated soils near the Fukushima Daiichi nuclear power plant. *Front Energy Res. 3*, 1–10 (2015).

10. Tanaka, K. et al. Vertical profiles of Iodine-131 and Cesium-137 in soils in Fukushima Prefecture related to the Fukushima Daiichi Nuclear Power station accident. *Geochem. J.* **46**, 73–76 (2012).

11. Yamasaki, S. et al. Radioactive Cs in the estuary sediments near Fukushima Daiichi Nuclear Power Plant. *Sci. Total Environ.* **551–552**, 155–162 (2016).

12. Adachi, K., Kajino, M., Zaizen, Y. & Igarashi, Y. Emission of spherical cesium-bearing particles from an early stage of the Fukushima nuclear accident. *Sci. Rep.* **3**, 2554/1–5 (2013).

13. Abe, Y. et al. Detection of uranium and chemical state analysis of individual radioactive microparticles emitted from the Fukushima nuclear accident using multiple synchrotron radiation X-ray analyses. *Anal. Chem.* **86**, 8521–8525 (2014).

14. Yamaguchi, N. et al. Internal structure of cesium-bearing radioactive microparticles released from Fukushima nuclear power plant. *Sci. Rep.* **6**, 20548 (2016).

15. Imoto, J. et al. Caesium fallout in Tokyo on 15th March, 2011 is dominated by highly radioactive, caesium-rich microparticles. *Sci. Rep.* **7**, 42118, doi:10.1038/srep42118 (2017) In press.

16. Furuki, G. et al. Caesium-rich micro-particles: A window into the meltdown events at the Fukushima Daiichi Nuclear Power Plant. *Sci. Rep.* **7**, 42731, doi:10.1038/srep42731 (2017).

17. Burns, P. C., Ewing, R. C. & Navrotsky, A. Nuclear fuel in a reactor accident. *Science* **335**, 1184–1188 (2012).

18. Bowman, S. M. & Gauld, I. C. OrigenArp Primer: How to perform isotopic depletion and decay calculations with SCALE/ORIGEN, ORNL/TM-2010/43 (2010).

19. Nishihara, K., Iwamoto, H. & Suyama, K. Estimation of Fuel Compositions in Fukushima-Daiichi Nuclear Power Plant. JAEA-Data/Code 2012–018 (2012).

20. Choppin, G. R., Liljenzin, J.-O. & Rydberg, J. *Radiochemistry and Nuclear Chemistry*, 3rd Edition. Butterworth-Heinemann, Woburn, MA, 709pp (2002).

21. Zimmer, M. M., Kinman, W. S., Kara, A. H. & Steiner, R. E. Evaluation of the homogeneity of the uranium isotope composition of NIST SRM 610/611 by MC-ICP-MS, MC-TIMS, and SIMS. *Minerals* **4**, 541–552 (2014).

22. Eugster, O., Tera, F. & Wasserburg, G. J. Isotopic analyses of barium in meteorites and in terrestrial samples. *J. Geophys. Res.* **74**, 3897–3908 (1969).

23. Pavese, A., Levy, D. & Hoser, A. Cation distribution in synthetic zinc ferrite (Zn<sub>0.97</sub>Fe<sub>2.02</sub>O<sub>4</sub>) from *in situ* high-temperature neutron powder diffraction. *Am. Mineral.* **85**, 1497–1502 (2000).

24. Kopp, O. C., Harris, L. A., Clark, G. W. & Yakel, H. L. A hydrothermally synthesized iron analog of pollucite–its structure and significance. *Am. Mineral.* **48**, 100–109 (1963).

25. Beger, R. M. The crystal structure and chemical composition of pollucite. *Zeitschrift für Kristallographie* **129**, 280–302 (1969).

26. Pontillon, Y. & Ducros, G. Behaviour of fission products under severe PWR accident conditions The VERCORS experimental programme—Part 2: Release and transport of fission gases and volatile fission products. *Nucl. Eng. Design* **240**, 1853–1866 (2010).

27. Shibahara, Y. et al. Analysis of cesium isotope compositions in environmental samples by thermal ionization mass spectrometry - 1. A preliminary study for source analysis of radioactive contamination in Fukushima prefecture. *J. Nucl. Sci. Technol.* **51**, 575–579 (2014).

28. Snow, M. S., Snyder, D. C. & Delmore, J. E. Fukushima Daiichi reactor source term attribution using cesium isotope ratios from contaminated environmental samples. *Rapid Commun. Mass Spectrom.* **30**, 523–532 (2016).

29. *Nuclear Safety in Light Water Reactors. Severe Accident Phenomenology* (ed. Sehgal, B. R.) Academic Press, Amsterdam, The Netherlands, 714pp (2012). And references therein.

30. Konings, R. J. M., Wiss, T. & Beneš, O. Predicting material release during a nuclear reactor accident. *Nature Mater.* **14**, 247–252 (2015).

31. Hiernaut, J.-P., Wiss, T., Papaioannou, D., Konings, R. J. M. & Rondinella, V. V. Volatile fission product behaviour during thermal annealing of irradiated UO₂ fuel oxidised up to U₃O₈. *J. Nucl. Mater.* **372**, 215–225 (2008).

32. Pontillon, Y. & Ducros, G. Behaviour of fission products under severe PWR accident conditions. The VERCORS experimental programme—Part 3: Release of low-volatile fission products and actinides. *Nucl. Eng. Design* **240**, 1867–1881 (2010).

33. Kissane, M. P. On the nature of aerosols produced during a severe accident of a water-cooled nuclear reactor. *Nucl. Eng. Design* **238**, 2792–2800 (2008).

34. Ewing, R. C. Long-term storage of spent nuclear fuel. *Nature Mater.* **14**, 252–257 (2015).

35. Sakaguchi, A. et al. Isotopic determination of U, Pu and Cs in environmental waters following the Fukushima Daiichi nuclear power plant accident. *Geochem. J.* **46**, 355–360 (2012).

36. Sakaguchi, A., Steier, P., Takahashi, Y. & Yamamoto, M. Isotopic compositions of <sup>236</sup>U and Pu isotopes in “Black Substances” collected from roadsides in Fukushima prefecture: fallout from the Fukushima Dai-ichi Nuclear Power Plant Accident. *Environ. Sci. Technol.* **48**, 3691–3697 (2014).

37. Zheng, J., Tagami, K. & Uchida, S. Rapid analysis of U isotopes in vegetables using ICP-MS: application to the emergency U monitoring after the nuclear accident at TEPCO's Fukushima Dai-ichi power station. *J. Radioanal. Nucl. Chem.* **292**, 171–175 (2012).

38. Di Lemma, F. G., Nakajima, K., Yamashita, S. & Osaka, M. Surface analyses of cesium hydroxide chemisorbed onto type 304 stainless steel. *Nucl. Eng. Design* **305**, 411–420 (2016).

39. Henry, P. F. & Weller, M. T. CsFeSiO₄: a maximum iron content zeotype. *Chem. Comm.* **24**, 2723–2724 (1998).

40. Bell, A. M. T. & Henderson, C. M. B. Rietveld refinement of the structures of dry-synthesized MFe<sup>III</sup>Si<sub>2</sub>O<sub>6</sub> leucites (M = K, Rb, Cs) by synchrotron X-ray powder diffraction. *Acta Crystallogr. C* **50**(10), 1531–1536 (1994).

41. Sevón, T. A MELCOR model of Fukushima Daiichi Unit 1 accident. *Annals Nucl. Energy* **85**, 1–11 (2015).

42. Fernandez-Moguel, L. & Birchley, J. Analysis of the accident in the Fukushima Daiichi nuclear power station Unit 3 with MELCOR_2.1. *Annals Nucl. Energy* **83**, 193–215 (2015).

43. Yamanaka, Y., Mizokami, S., Watanabe, M. & Honda, T. Update of the first TEPCO MAAP accident analysis of units 1, 2, and 3 at Fukushima Daiichi Nuclear Power Station. *Nucl. Technol.* **186**, 263–279 (2014).

44. Nuclear Regulation Agency, Benchmark study of the accident at the Fukushima Daiichi Nuclear Power Plant (BSAF Project) Phase I summary report. NEA/CSNI/R(2015)18.
4. Kinoshita, N. et al. Assessment of individual radionuclide distributions from the Fukushima nuclear accident covering central-east Japan. *Proc. Natl. Acad. Sci. USA* **108**, 19526–19529 (2011).

5. Morino, Y., Ohara, T. & Nishizawa, M. Atmospheric behavior, deposition, and budget of radioactive materials from the Fukushima Daiichi nuclear power plant in March 2011. *Geophys. Res. Lett.* **38**, L00G11 (2011).

6. Buesseler, K. et al. Fukushima Daiichi-derived radionuclides in the ocean: transport, fate, and impacts. *Annu. Rev. Mar. Sci.* **9**, 173–203 (2017).

7. Grambow, B. & Mostafavi, M. State of Fukushima nuclear fuel debris tracked by Cs137 in cooling water. *Environ. Sci.: Processes Impacts* **16**, 2472–2476 (2014).

8. Kaneyasu, N. et al. Sulfate aerosol as a potential transport medium of radiocesium from the Fukushima nuclear accident. *Environ. Sci. Technol.* **46**, 5720–5726 (2012).

9. Kaneko, M. et al. Radioactive Cs in the severely contaminated soils near the Fukushima Daiichi nuclear power plant. *Front Energy Res. 3*, 1–10 (2015).

10. Tanaka, K. et al. Vertical profiles of Iodine-131 and Cesium-137 in soils in Fukushima Prefecture related to the Fukushima Daiichi Nuclear Power station accident. *Geochem. J.* **46**, 73–76 (2012).

11. Yamasaki, S. et al. Radioactive Cs in the estuary sediments near Fukushima Daiichi Nuclear Power Plant. *Sci. Total Environ.* **551–552**, 155–162 (2016).

12. Adachi, K., Kajino, M., Zaizen, Y. & Igarashi, Y. Emission of spherical cesium-bearing particles from an early stage of the Fukushima nuclear accident. *Sci. Rep.* **3**, 2554/1–5 (2013).

13. Abe, Y. et al. Detection of uranium and chemical state analysis of individual radioactive microparticles emitted from the Fukushima nuclear accident using multiple synchrotron radiation X-ray analyses. *Anal. Chem.* **86**, 8521–8525 (2014).

14. Yamaguchi, N. et al. Internal structure of cesium-bearing radioactive microparticles released from Fukushima nuclear power plant. *Sci. Rep.* **6**, 20548 (2016).

15. Imoto, J. et al. Caesium fallout in Tokyo on 15th March, 2011 is dominated by highly radioactive, caesium-rich microparticles. *Sci. Rep.* **7**, 42118, doi:10.1038/srep42118 (2017) In press.

16. Furuki, G. et al. Caesium-rich micro-particles: A window into the meltdown events at the Fukushima Daiichi Nuclear Power Plant. *Sci. Rep.* **7**, 42731, doi:10.1038/srep42731 (2017).

17. Burns, P. C., Ewing, R. C. & Navrotsky, A. Nuclear fuel in a reactor accident. *Science* **335**, 1184–1188 (2012).

18. Bowman, S. M. & Gauld, I. C. OrigenArp Primer: How to perform isotopic depletion and decay calculations with SCALE/ORIGEN, ORNL/TM-2010/43 (2010).

19. Nishihara, K., Iwamoto, H. & Suyama, K. Estimation of Fuel Compositions in Fukushima-Daiichi Nuclear Power Plant. JAEA-Data/Code 2012–018 (2012).

20. Choppin, G. R., Liljenzin, J.-O. & Rydberg, J. *Radiochemistry and Nuclear Chemistry*, 3rd Edition. Butterworth-Heinemann, Woburn, MA, 709pp (2002).

21. Zimmer, M. M., Kinman, W. S., Kara, A. H. & Steiner, R. E. Evaluation of the homogeneity of the uranium isotope composition of NIST SRM 610/611 by MC-ICP-MS, MC-TIMS, and SIMS. *Minerals* **4**, 541–552 (2014).

22. Eugster, O., Tera, F. & Wasserburg, G. J. Isotopic analyses of barium in meteorites and in terrestrial samples. *J. Geophys. Res.* **74**, 3897–3908 (1969).

23. Pavese, A., Levy, D. & Hoser, A. Cation distribution in synthetic zinc ferrite (Zn<sub>0.97</sub>Fe<sub>2.02</sub>O<sub>4</sub>) from *in situ* high-temperature neutron powder diffraction. *Am. Mineral.* **85**, 1497–1502 (2000).

24. Kopp, O. C., Harris, L. A., Clark, G. W. & Yakel, H. L. A hydrothermally synthesized iron analog of pollucite–its structure and significance. *Am. Mineral.* **48**, 100–109 (1963).

25. Beger, R. M. The crystal structure and chemical composition of pollucite. *Zeitschrift für Kristallographie* **129**, 280–302 (1969).

26. Pontillon, Y. & Ducros, G. Behaviour of fission products under severe PWR accident conditions The VERCORS experimental programme—Part 2: Release and transport of fission gases and volatile fission products. *Nucl. Eng. Design* **240**, 1853–1866 (2010).

27. Shibahara, Y. et al. Analysis of cesium isotope compositions in environmental samples by thermal ionization mass spectrometry - 1. A preliminary study for source analysis of radioactive contamination in Fukushima prefecture. *J. Nucl. Sci. Technol.* **51**, 575–579 (2014).

28. Snow, M. S., Snyder, D. C. & Delmore, J. E. Fukushima Daiichi reactor source term attribution using cesium isotope ratios from contaminated environmental samples. *Rapid Commun. Mass Spectrom.* **30**, 523–532 (2016).

29. *Nuclear Safety in Light Water Reactors. Severe Accident Phenomenology* (ed. Sehgal, B. R.) Academic Press, Amsterdam, The Netherlands, 714pp (2012). And references therein.

30. Konings, R. J. M., Wiss, T. & Beneš, O. Predicting material release during a nuclear reactor accident. *Nature Mater.* **14**, 247–252 (2015).

31. Hiernaut, J.-P., Wiss, T., Papaioannou, D., Konings, R. J. M. & Rondinella, V. V. Volatile fission product behaviour during thermal annealing of irradiated UO₂ fuel oxidised up to U₃O₈. *J. Nucl. Mater.* **372**, 215–225 (2008).

32. Pontillon, Y. & Ducros, G. Behaviour of fission products under severe PWR accident conditions. The VERCORS experimental programme—Part 3: Release of low-volatile fission products and actinides. *Nucl. Eng. Design* **240**, 1867–1881 (2010).

33. Kissane, M. P. On the nature of aerosols produced during a severe accident of a water-cooled nuclear reactor. *Nucl. Eng. Design* **238**, 2792–2800 (2008).

34. Ewing, R. C. Long-term storage of spent nuclear fuel. *Nature Mater.* **14**, 252–257 (2015).

35. Sakaguchi, A. et al. Isotopic determination of U, Pu and Cs in environmental waters following the Fukushima Daiichi nuclear power plant accident. *Geochem. J.* **46**, 355–360 (2012).

36. Sakaguchi, A., Steier, P., Takahashi, Y. & Yamamoto, M. Isotopic compositions of <sup>236</sup>U and Pu isotopes in “Black Substances” collected from roadsides in Fukushima prefecture: fallout from the Fukushima Dai-ichi Nuclear Power Plant Accident. *Environ. Sci. Technol.* **48**, 3691–3697 (2014).

37. Zheng, J., Tagami, K. & Uchida, S. Rapid analysis of U isotopes in vegetables using ICP-MS: application to the emergency U monitoring after the nuclear accident at TEPCO’s Fukushima Dai-ichi power station. *J. Radioanal. Nucl. Chem.* **292**, 171–175 (2012).

38. Di Lemma, F. G., Nakajima, K., Yamashita, S. & Osaka, M. Surface analyses of cesium hydroxide chemisorbed onto type 304 stainless steel. *Nucl. Eng. Design* **305**, 411–420 (2016).

39. Henry, P. F. & Weller, M. T. CsFeSiO₄: a maximum iron content zeotype. *Chem. Comm.* **24**, 2723–2724 (1998).

40. Bell, A. M. T. & Henderson, C. M. B. Rietveld refinement of the structures of dry-synthesized MFe<sup>III</sup>Si<sub>2</sub>O<sub>6</sub> leucites (M = K, Rb, Cs) by synchrotron X-ray powder diffraction. *Acta Crystallogr. C* **50**(10), 1531–1536 (1994).

41. Sevón, T. A MELCOR model of Fukushima Daiichi Unit 1 accident. *Annals Nucl. Energy* **85**, 1–11 (2015).

42. Fernandez-Moguel, L. & Birchley, J. Analysis of the accident in the Fukushima Daiichi nuclear power station Unit 3 with MELCOR_2.1. *Annals Nucl. Energy* **83**, 193–215 (2015).

43. Yamanaka, Y., Mizokami, S., Watanabe, M. & Honda, T. Update of the first TEPCO MAAP accident analysis of units 1, 2, and 3 at Fukushima Daiichi Nuclear Power Station. *Nucl. Technol.* **186**, 263–279 (2014).

44. Nuclear Regulation Agency, Benchmark study of the accident at the Fukushima Daiichi Nuclear Power Plant (BSAF Project) Phase I summary report. NEA/CSNI/R(2015)18.
**Acknowledgements**

This study is partially supported by JST Initiatives for Atomic Energy Basic and Generic Strategic Research and by a Grant-in-Aid for Scientific Research (KAKENHI) from the Japan Society for the Promotion of Science (16K12585, 16H04634, No. JP26257402). S.U. is also supported by ESPEC Foundation for Global Environment Research and Technology (Charitable Trust) (ESPEC Prize for the Encouragement of Environmental Studies). The authors are grateful to Dr. Watanabe for her assistance on SEM analyses at the Center of Advanced Instrumental Analysis, Kyushu University. The findings and conclusions of the authors of this paper do not necessarily state or reflect those of the JST.

**Author Contributions**

S.U. conceived the idea, designed all experiments, and wrote the manuscript. J.I. performed measurements, data analysis, and wrote the first draft. A.O. conducted TEM analysis. G.F., K.H., and M.T. performed SIMS analysis. M.S. and R.I. performed SEM analysis. T.O. and K.N. provided navigation during field research in Fukushima. S.Y. performed gamma spectroscopy at Tsukuba University. G.T.W.L., B.G., and R.C.E. participated in the discussion and interpretation of the results and helped to write the paper.

**Additional Information**

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-05910-z

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017
## Acknowledgements

This study is partially supported by JST Initiatives for Atomic Energy Basic and Generic Strategic Research and by a Grant-in-Aid for Scientific Research (KAKENHI) from the Japan Society for the Promotion of Science (16K12585, 16H04634, No. JP26257402). S.U. is also supported by ESPEC Foundation for Global Environment Research and Technology (Charitable Trust) (ESPEC Prize for the Encouragement of Environmental Studies). The authors are grateful to Dr. Watanabe for her assistance on SEM analyses at the Center of Advanced Instrumental Analysis, Kyushu University. The findings and conclusions of the authors of this paper do not necessarily state or reflect those of the JST.

## Author Contributions

S.U. conceived the idea, designed all experiments, and wrote the manuscript. J.I. performed measurements, data analysis, and wrote the first draft. A.O. conducted TEM analysis. G.F., K.H., and M.T. performed SIMS analysis. M.S. and R.I. performed SEM analysis. T.O. and K.N. provided navigation during field research in Fukushima. S.Y. performed gamma spectroscopy at Tsukuba University. G.T.W.L., B.G., and R.C.E. participated in the discussion and interpretation of the results and helped to write the paper.

## Additional Information

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-05910-z

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017