Dose Estimation of External and Internal Exposure in Japanese Macaques After the Fukushima Nuclear Power Plant Accident

Cumulative dose of external and internal exposures following the Fukushima Daiichi Nuclear Power Plant (FNPP) accident was estimated for Japanese macaques (Macaca fuscata) around FNPP. Conversion factors for Japanese macaques modeled as ellipsoids were estimated for external exposure from contaminated ground and internal exposure uniformly distributed in the body. Conversion factors for seven radionuclides, namely, tellurium-129 (129Te), 129mTe, iodine-131 (131I), 132Te, 132I, cesium-134 (134Cs) and 137Cs were calculated using the PHITS code. The estimated factors for the seven radionuclides were consistent with those in ICRP Publication 108, using an effective radius for comparison. The external, internal and total exposures for 13 macaques in Namie Town were estimated by applying the calculated factors. The estimated cumulative exposures for the periods from the accident occurred to the sampling date, ranged from 0.26 to 1.6 Gy. The average exposure was 0.64 Gy in averaged over the 11 sampled macaques except for the 2 macaques which might be born after the FNPP accident.


Introduction
The nuclear accident at the Fukushima Daiichi Nuclear Power Plant (FNPP) was triggered by an enormous earthquake and associated tsunami (the Great East Japan Earthquake) on March 11, 2011. A large quantity of radioactive nuclides were released, resulting in the severe contamination of a wide area of the southern Tohoku region to the northern Kanto region. The main depositions occurred on March 15 and March 20-21, 2011 [1,2].
Since the accident, many studies have examined the effect of radiation on animals, insects and plants [3][4][5][6][7][8]. However, cumulative exposure in animals is quite difficult to estimate because it is challenging to account for habitat and behavioral properties. Estimates of exposure could be improved using a contamination map combined with habitat information. Dose rate and its temporal changes can be estimated for target animals in an area from the contamination map.
The National Nuclear Security Administration (NNSA) performed a rapid survey of radiation and contamination on March 17-19, 2011 using the Aerial Measuring System. These data were used for the initial NNSA March 22, 2011, report and are available for independent analysis [9]. Additionally, the Ministry of Education, Culture, Sports, Science and Technology, Japan (MEXT), conducted a contamination study using a 2-km mesh from June to August 2011 [10]. This study started 3 months after the main deposition on March 15, 2011. The radionuclides with short half-lives already decayed away; however, 134 Cs and 137 Cs remained over a wide area around FNPP.
Studies of the biological effect of radiation following the FNPP accident were initiated by Tohoku University as a part of a Comprehensive Dose Evaluation Project Concerning Animals Affected by The Fukushima Daiichi Nuclear Power Plant Accident and a Nuclear Energy Science and Technology and Human Resource Development Project by MEXT using Japanese macaques (Macaca fuscata) sampled in the Fukushima area. Macaques, like humans, belong to the order, Primates. Therefore, the biological effect on macaques might provide important information that can be applied to humans. To estimate exposure, the International Commission of Radiological Protection (ICRP) summarized conversion factors for reference animals under internal exposure to uniformly distributed radionuclides and external exposure to contaminated ground with bodies simulated as ellipsoids in ICRP Publication 108 [11]. The factors are available for 75 radionuclides from ICRP Publication 108. However, the Japanese macaque is not one of the ICRP reference animals. Therefore, we estimated the appropriate conversion factors for the species using the Monte Carlo technique. NNSA [9] and MEXT data [10] for 137 Cs contamination were combined and interpolated by the Universal Kriging method using System for Automated Geoscientific Analyses-Geographic Information System (SAGA-GIS) software [12]. A variogram was generated, which describes the variance of the difference between 137 Cs concentrations at two locations fitted by a quartic function, using SAGA-GIS [10]. The reproducibility of the variogram was 95% within 70 km from FNPP. The interpolated values were stored in 1,814 × 1,854 meshes for latitudes from 36.51611 to 38.32911 and longitudes from 139.30506 to 141.15896. The interpolated map is shown in Fig. 14

Ellipsoid Model
In ICRP Publication 108, the reference animals in terrestrial conditions are modeled as various ellipsoids [11], as listed in Table 14.1. Therefore, Japanese macaques were also modeled as ellipsoids with a density of 1 g/cm 3 according to ICRP Publication 108 assumption. Japanese macaques were categorized into three groups

Conversion Factor Calculation
The conversion factors for internal exposure to uniformly distributed radionuclides in an ellipsoid and external exposure to uniformly contaminated ground from 0 to 0.5 cm were calculated using PHITS code [13]. The calculation geometry of air occupying a 10 × 10 × 10 m area and 50-cm-deep ground soil was taken into account. The ellipsoid macaque was set at the center of the ground surface. To minimize calculation time, the four side boundaries were set to a mirror condition. Gamma-rays and β-particles were separately calculated for seven radionuclides, namely, 129 Te, 129m Te, 131 I, 132 Te, 132 I, 134 Cs and 137 Cs. The γ-ray energies and emission rates for the radionuclides were obtained from the National Nuclear Data Center [14]. Beta-particle energy spectra were obtained from the literature [15,16], and the internal conversion electrons for each radionuclide were obtained from the website of the National Nuclear Data Center [14]. The energy deposition in the ellipsoids was calculated for each γ-ray and each electron source. In addition, to consider the depth profile of radionuclides, the deposition energy in the walk model for body size L was also calculated for sources uniformly distributed in six soil layers, namely, 0.5-1 cm, 1-2 cm, 2-5 cm, 5-10 cm, 10-15 cm and 15-20 cm ( Fig. 14.2). The external exposure conversion factors were derived by weighing each deposition energy for the depth profile of the radionuclide. The factor D i ex ((mGy/day)/ (Bq/m 2 )) of the ith radionuclide from depth layer d k is expressed by where subscript j indicates emitted γ-rays or β-particles, I j is the emission rate for emitted γ-rays or β-particles, m is the mass of the macaque-simulating ellipsoid, E is the deposition energy in the ellipsoid by emitted γ-rays and β-particles and S is the calculated area expressed in m 2 . In the case of soil surface contamination, similar to the ICRP Publication 108 assumption, the depth layer d 1 is only 0-0.5 cm. To estimate the factor for the depth profile of radionuclides, the fraction of radionuclides as a function of depth, f(d k ), is introduced and summed for seven layers: The internal exposure conversion factor D i in ((mGy/d)/(Bq/kg)) for each radionuclide i was calculated using each deposition energy as follows:

Depth Profile
There is a lack of depth profile data, except for radioactive Cs. The depth profiles of 137 Cs radioactivity have been measured at various locations [17][18][19][20] [20]. The conversion factor for the 137 Cs depth-profiled soil was calculated using the relaxation mass depths from 0.5 to 2.7 according to Eq. (14.2). Additionally, the typical depth profiles for our measurement, as shown in Fig. 14.3, were tested.

Conversion Factor
The calculated conversion factors for the surface layer for seven radionuclides are listed in Table 14.3. The external exposure results for the walk model and the internal exposure results were compared with those in ICRP Publication 108 against the effective radius and are shown in Fig. 14.4. The calculated factor and those of ICRP Publication 108 agree within 20%. Therefore, if we allow for 20% uncertainty, the conversion factor in ICRP Publication 108 for the effective radius can be used.   The conversion factor ratio of external exposure for the walk to the sit model is plotted as a function of E , which is defined as the sum of mean γ-ray and β-particle energies, in Fig. 14.5. The ratio is consistent with unity within statistical uncertainty, except for the lowest E of 132 Te. The ratio for 132 Te, which emits low-energy γ-rays and β-particles, is approximately 0.83. This is due to the shielding effect by an individual's own body. Therefore, except for low-energy β-particles and γ-rays, the conversion factors for both the walk model and sit model are approximately equal.

Conversion Factor for Depth-Profiled Radionuclides
The conversion factors for the seven depth-profiled radionuclides with a relaxation depth from 0.6 to 2.6 and measured depth profiles of 134 Cs and 137 Cs in Iitate Village, Fukushima Prefecture, in comparison with the 0-0.5 cm uniform factor are shown in Fig. 14.6a, b, respectively. The relaxation depths of 137 Cs radioactivity from previous studies [17][18][19][20] are shown in light red and ranged from 0.9 to 1.5 g/cm 2 . A soil self-shielding factor introduced as a ratio of the conversion factor for the depth profiles to that for the uniformly distributed factor in surface soil characterized by the relaxation depth is summarized in Table 14.4. The soil self-shielding factors for 134 Cs and 137 Cs were 0.74 and 0.71, respectively. The shielding factor for Type 1 for the slowest measured depth profiles for farmland in Fig. 14.6b was less than that for a relaxation depth of 2.6 g/cm 2 . However, the shielding factor for Type 2 for the typical depth profile for undisturbed land is consistent with those for the relaxation depth of 0.9-1.5 g/cm 2 . The variation within the range of 0.9-1.5 g/cm 2 corresponds to the variation in the conversion factor for 134 Cs or 137 Cs, which was less than 8%.

Example of External Exposure Estimation
To consider cumulative exposure to soil contamination, temporal changes in the air dose from radionuclides deposited on the ground are good references. The air dose rate from seven radionuclides from radioactivity concentrations was estimated around FNPP in 2011 [9,10]. The air dose rate changes were estimated for several locations. For example, the estimated dose rate change in Iitate Village assuming only the physical half-life [2] is shown in Fig. 14.7. Integrating over time, the cumulative air dose can be estimated. The cumulative air dose normalized by the dose at 200 years after deposition and the radionuclide contribution as a function of time are shown in Fig. 14.8a, b, respectively. The sum of the 134 Cs and 137 Cs contributions to the cumulative air dose at 1 year after deposition or later exceeds 90%. In the sampled macaques, exposure is expected to be the same as the contributions to the air dose. Accordingly, the external exposure for some samples after 1 year following the FNPP accident can be estimated by only radioactive 134 Cs and 137 Cs, neglecting the 10% contribution of external exposure to other radionuclides.  Using the 134 Cs and 137 Cs conversion factors for the walk model, the external, the internal and the total exposures were calculated for the sampled macaques at Namie Town, Fukushima Prefecture, in 2013. Only 137 Cs concentration and 134 Cs concentration, which is the same as 137 Cs concentration in March 2011, were used in this estimation. The ground concentration of 137 Cs (interpolated values in Fig. 14.1) and the 137 Cs concentration in macaque bodies measured at Tohoku University [21] were used in this calculation. To account for the effect of the depth profile, the soil selfshielding factor of the depth profile for a relaxation depth of 1.3 g/cm 2 was assumed to be 0.74 for 134 Cs and 0.71 for 137 Cs. The macaques were caught in 2013 or later; therefore, only exposure to 134 Cs and 137 Cs was taken into account. Additionally, the ratio of the accumulated concentration of 134 Cs to 137 Cs was assumed to be unity. The estimated total exposures ranged from 0.26 to 1.6 Gy, as shown in Table 14.5. Two sampled macaques showed the body weight less than 2kg. There is a possibility that these 2 macaques were born after the FNPP accident. The average exposure was 0.64 Gy in averaged over the 11 sampled macaques except for the 2 macaques.  There is a possibility that the 2 macaques were born after the FNPP accident

Conclusions
The cumulative external and internal exposures attributed to the FNPP accident were estimated for Japanese macaques around FNPP. Conversion factors for Japanese macaques, modeled as ellipsoids, for external exposure from contaminated ground and internal exposure uniformly distributed in the body were calculated. The estimated conversion factors for seven radionuclides, namely, 129 Te, 129m Te, 131 I, 132 Te, 132 I, 134 Cs and 137 Cs, agreed within 20% with those in ICRP Publication 108 using an effective radius.
The conversion factors for depth-profiled 134 Cs and 137 Cs characterized by the relaxation depth were also estimated. The conversion factors are expressed by the factor for a uniform distribution at the 0-0.5 cm surface multiplied by a soil selfshielding factor of 0.74 for 134 Cs and 0.71 for 137 Cs and with a relaxation depth of 1.3 g/cm 2 .
The external, the internal and the total exposures for 13 sampled macaques were estimated using the calculated conversion factors. The estimated cumulative exposures for the periods from the accident occurred to the sampling date, ranged from 0.26 to 1.6 Gy. The average exposure was 0.64 Gy in averaged over the 11 sampled macaques except for the 2 macaques which might be born after the FNPP accident.
Open Access This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 chapter are included in the chapter's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the chapter'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.