NIRS external dose estimation system for Fukushima residents after the Fukushima Dai-ichi NPP accident.

The great east Japan earthquake and subsequent tsunamis caused Fukushima Dai-ichi Nuclear Power Plant (NPP) accident. National Institute of Radiological Sciences (NIRS) developed the external dose estimation system for Fukushima residents. The system is being used in the Fukushima health management survey. The doses can be obtained by superimposing the behavior data of the residents on the dose rate maps. For grasping the doses, 18 evacuation patterns of the residents were assumed by considering the actual evacuation information before using the survey data. The doses of the residents from the deliberate evacuation area were relatively higher than those from the area within 20 km radius. The estimated doses varied from around 1 to 6 mSv for the residents evacuated from the representative places in the deliberate evacuation area. The maximum dose in 18 evacuation patterns was estimated to be 19 mSv.

The great east Japan earthquake occurred at 14:46, Japan time on 11 March 2011 in the Tohoku District, northeastern part of Japan. After the quake, big tsunamis rushed to the coast of the Tohoku district including Fukushima prefecture. The tsunamis caused electric failure of Fukushima Dai-ichi Nuclear Power Plant (NPP). The cooling systems of nuclear reactors stopped, and the hydrogen explosions occurred. Subsequently, radioactive fission products leaked out to the environment12.

The interest in radiation doses due to the accident has been increasing among the Fukushima residents. The importance to estimate exposure doses of the residents has been recognized among the involved organizations. National Institute of Radiological Sciences (NIRS) started to develop the external dose estimation system for Fukushima residents at the end of March, 2011. On the other hand, Fukushima prefectural government and Fukushima Medical University decided to perform the health management survey for all Fukushima residents (about two million people) at the end of May, 20113. External dose was considered to be one of the necessary data for the health management, and then the NIRS external dose estimation system was adopted in the survey.

This paper describes the NIRS external dose estimation system and the calculation results for 18 evacuation patterns assumed by NIRS.

Results

We developed the NIRS external dose estimation system for the Fukushima residents to estimate the external effective doses for the first four months after the Fukushima Dai-ichi nuclear power plant accident. This system has been adopted in the Fukushima Health Management Survey, and the estimated result is sequentially provided to the individual Fukushima residents.

The estimated external doses of the adult residents in 18 evacuation patterns were shown in Table 1. The results show that the doses of the residents evacuated from the representative places in the deliberate evacuation area (typically 1–6 mSv) were relatively higher than those from the area within 20 km radius (around 1 mSv). The maximum dose in 18 evacuation patterns was estimated to 19 mSv.

Table 1

Estimated external effective doses (mSv) for adult of 18 patterns

Estimated external effective doses for adult (mSv)
From the area within 20 km			From the deliberate evacuation area
Pattern No.	Representative place	Max. dose place	Pattern No.	Representative place	Max. dose place
1	0.76	0.76	13	4.8	13
2	0.94	1.1	14	0.84	7.2
3	0.22	0.64	15	5.5	17
4	0.23	1.6	16	6.2	19
5	0.30	0.30	17	-	1.8
6	0.25	0.25	18	-	3.0
7	2.0	2.3
8	-	0.77
9	-	1.5
10	-	0.18
11	-	0.74
12	-	0.71

Discussion

The NIRS external dose estimation system for the Fukushima residents has been developed based on the simulation with the System for Prediction of Environmental Emergency Dose Information (SPEEDI)4 and monitoring data reported officially by the government. The uncertainties in this dose estimation depend mainly on not only those of the dose rate map but also the residents' memory on the behavior.

With regard to the dose rate map, since maps for the first 3 days are based on the SPEEDI simulation instead of the monitoring data, the uncertainty is expected to be higher than those for other days. We have no alternatives for first 3 days at present because of the scarce monitoring data. The calculated results of SPEEDI greatly depend on the source term emitted from the nuclear reactors, which is expressed as a function of time. In the current situation, we have used the available and reliable data of the source term reported by the government. However, when the accuracy improvement of the simulation is done with further analysis on the source term, recalculations with the dose rate maps base on the newly-simulated results would be necessary. On the other hand, the contributions of exposures in early stage for total estimated doses are not so much, because most residents near the NPP have evacuated early after the accident and the increases of the dose by radionuclides in the plume were limited for a few hours as shown in Fig. 1. Actually the dose by radionuclides deposited on the ground etc., which is relatively plateau for the four months from the accident, is more important. For example, even at the Minami-soma city office where the increase of dose was observed earlier than other areas as shown in Fig. 1, the contribution of the first 3 days to the total dose of four months was estimated to be about 5% by the NIRS system and about 10% when using the monitoring data. This fact indicated that the refinement of the simulation for the first 3 days would not so affect to the total dose.

Figure 1

Monitoring data reported by MEXT and Fukushima prefectural government567.

In the NIRS system, the conversion coefficients from ambient dose equivalents to effective doses for isotropic irradiation were adopted. ICRP does not provide the conversion coefficients for the special condition of the irradiation. Saito et al.8 calculated the conversion coefficients of effective dose from environmental gamma ray sources: cloud and ground shines. The coefficients of isotropic irradiation (ISO), lateral irradiation from the right side to the left side of the body (RLAT) and lateral irradiation from the left side to the right side of the body (LLAT) are close to that assumed gamma ray source on the ground according to ICRP 749 and Saito el al., and the coefficient of ISO is the most conservative. Table 2 shows the conversion coefficients from ambient dose equivalent to effective dose for adult for various radionuclides considered in this study calculated from ICRP 74 (adopted in this study) and Saito el al. This result indicated that conversion coefficient of 0.6 used in the NIRS system is reasonable.

Table 2

Conversion coefficients from ambient dose equivalent to effective dose for adult for various radionuclides

Nuclide	Xe-133	Te-129m	Te-132	I-131	I-132	Cs-134	Cs-137
Saito et al.8 (Ground)	0.41	0.61	0.51	0.57	0.62	0.61	0.61
Saito et al.8 (Air)	0.34	0.57	0.46	0.52	0.58	0.57	0.57
Present	0.44	0.58	0.48	0.53	0.59	0.58	0.57

Methods

Flow of the dose calculation

In the NIRS system, the external effective dose between March 12 and July 11, 2011 is estimated by superimposing the individual behavior data of each day on the daily dose rate map of such day. The data flow in the external dose estimation system is shown in Fig. 2. The flow is composed of three parts: 1) behavior data of residents, 2) construction of daily dose rate maps, and 3) calculation of effective dose. Each part is described below. This system utilizes two kinds of the fundamental data: daily dose rate map and behavior data of the resident. The background dose due to natural radiation sources was not included in the dose rate maps, because of the radiation protection purposes. The body size as a function of age is also considered because the external effective dose depends on body size even in the same radiation field. The dose rate maps have the data of effective dose for adults, so body size correction is considered for children less than or equal to 15 years old by using age coefficients.

Figure 2

Data flow in the NIRS external dose estimation system.

Data of behaviors of the residents

Fukushima prefectural government has been conducting the health management survey for the Fukushima residents to support their health conditions affected by the accidents in cooperation with Fukushima Medical University3. The survey consists of two kinds of investigations, which are called the basic survey and the detailed survey, respectively. In the basic survey, questionnaires regarding behavior of a resident for four months after the accident were sent by postal mail to the all residents in Fukushima prefecture, about 2 million in the population. Based on such data, the behavior data of the resident for each day can be obtained as follows;

Place to stay, i,

Time of stay in the building at i, t,

Time of stay outside the building at i, t

Type of the building at i, k,

Time of the move from i to i + 1, t, (if any)

Construction of daily dose rate maps

The dose rate maps consisting of a lot of divisions were made based on the simulated and measured data in a time series between 12 March and 11 July in 2011. The size of each division in the map is corresponding to approximately 2 km × 2 km, which is based on the size divided the second mesh (latitude distance of 7.5 min. by longitude distance of 5 min) defined by Geospatial Information Authority of Japan.

The daily dose rate maps were constructed based on two kinds of data. One is the hourly effective dose rate maps simulated by SPEEDI with the source term calculated by the MELCOR code10 by the Nuclear and Industrial Safety Agency (NISA), which was used between March 12 and 14, 201111. These data are an alternative to the monitoring data, because the number of measurement point was not sufficient to construct the dose rate maps in the period. The hourly effective dose rate maps outputted from the SPEEDI simulations, were averaged over a day, because the dose rate maps used in our system were daily maps. Also, its area was divided to those of 1 km × 1 km. Therefore, the dose rate maps generated by SPEEDI were reconstructed by dividing in proportion to the area size of our system (approximately 2 km × 2 km) with commercially-available mapping software. As a result, the effective dose (for adult) rate maps with d(h,m), where h and m are the date and the division number, respectively, were obtained between 12 and 14 March, 2011. Unfortunately, their area was limited in 98 km × 98 km, which is painted in green in Fig. 3, depending on the SPEEDI ability. At present, the doses are not added outside of the area.

Figure 3

Areas of dose rate maps used in the NIRS external dose estimation system.

Basic maps were created by using ArcGIS of ESRI.

Another data is monitoring data released by the Ministry of Education, Culture, Sports, Science and Technology (MEXT), which was used between March 15 and July 11, 2011. The monitoring data scattered in the map were converted to the spatially-continuous data by using Natural Neighbor method and then the daily dose rate in each division of approximately 2 km × 2 km was obtained by averaging the values in such division. Since the monitoring data maps by MEXT could not cover the area of a part of seaside and the boundary between Fukushima and Niigata prefectures (pink-colored areas in Fig. 3), the values of their neighbor on the right or left were used alternatively.

Monitoring data on 15 March, 2011 were not sufficient to construct the dose rate map. On the other hand, this SPEEDI results were significantly different from the monitoring data. This was expected to derive from the highest time variation of the dose rate in Fig. 1. Therefore, the dose rate map on 16 March was alternatively used for that on 15 March. Table 3 shows the comparison of effective doses on 15 and 16 March, 2011 based on the monitoring data at the points shown in Fig. 1. The differences were between −11 and +182 μSv (effective dose), so the alternative does not lead to a significant underestimation. In fact, it leads to an overestimation at most points.

Table 3

Comparison of effective dose (μSv) on 15 and 16 March, 2011 based on the monitoring data

	Effective dose (μSv)
Measurement point	15 March, 2011 (A)	16 March, 2011 (B)	Difference (B-A)
Fukushima city (Momijiyama)	70	230	160
Fukusima city (Oyama-chou)	100	265	165
Iitate village	172	354	182
Tamura city	3	5	2
Minamisouma city	43	58	14
Kawauchi village	31	23	−7
Iwaki city	65	55	−11

In the NIRS system, the background dose from natural radionuclides, etc. is not included in the estimated dose. Therefore, the background dose was estimated from the reported data before the accident. Fukushima prefectural government has reported the averaged air kerma rates per month at 23 points in the prefecture (1 in Kouno town, 4 in Naraha town, 5 in Tomioka town, 5 in Okuma town, 4 in Futaba town and 4 in Namie town) as 33–54 μGy/h12. The other measured data showed 0.04 μGy/h in Fukushima city, 0.04–0.06 μGy/h in Koriyama city, 0.04–0.05 μGy/h in Shirakawa city, 0.04–0.05 μGy/h in Aizuwakamatsu city, 0.02–0.04 μGy/h in Minami-aizu town, 0.05 μGy/h in Minami-soma town and 0.05–0.06 μGy/h in Iwaki city5. As the background dose rate, 0.04 μGy/h as the median value among them was used for the background subtractions. Actually, the background effective dose rate was set to 0.03 μSv/h by applying the conversion coefficient from air kerma to effective dose.

The monitoring data are measured in terms of ambient dose equivalent H*(10), which is an operational guantity. On the other hand effective dose E is protection quantity and used for risk assessment of human health1314. Then the conversion factors from ambient dose equivalent to effective dose were used for calculations in the system. According to the previous study13, the main radionuclides emitted from the Fukushima Dai-ichi NPP were Xe-133, Te-132 and I-132 at the early stage of the accident, and Cs-134 and Cs-137 afterwards. Adding I-131 and Te-129m to these radionuclides, the conversion coefficients of effective dose per ambient dose equivalent for each radionuclide were calculated in the following steps;

Calculate conversion factors of for isotropic irradiation (ISO) and for the energy of the i-th gamma ray emitted from a radionuclide, and , according to the ICRP publication 74 with the log-log interpolation method. Here, φ is the fluence of gamma ray.

Calculate conversion factors of for isotropic irradiation (ISO) and for each nuclide, and , by the following equations: where R is the emission ratio of the i-th gamma ray emitted from a radionuclide.

Divide by

The calculated for each nuclide is shown in Table 2. In our estimation, the value of 0.6 was adopted as for adult, rounding the value of 0.59 which is the maximum among the radionuclides. Actually, since there are differences of tissue weighting factors of effective dose between ICRP Publication 60 13 and 103 14, estimated effective doses is not the same exactly depending on the publication. However, according to the reference16, these values are similar for isotropic irradiation. Hence the correction between them was not considered in the system.

Finally, the effective dose (for adult) rate maps with d (h,m), can be obtained in a time series from 15 March to 11 July, 2011 in all parts of Fukushima prefecture and a part of neighboring four prefectures (Miyagi, Yamagata, Tochigi, and Ibaraki) painted in pink and blue in Fig. 3.

Calculation of effective dose

The external effective doses were calculated for three different situations: staying indoors/outdoors and moving from one place to another.

When one stays indoors, the dose reduction should be considered because buildings have shielding effects against radiation exposures depending on their materials and thickness of the walls. In the TECDOC 22515 by the International Atomic Energy Agency (IAEA), representative reduction factors for cloud, r, and ground sources, r, are shown. From 12 to 14 March, the radio nuclides in the plume have contributed to the dose rates in the environment. On the contrary, the radionuclides on the ground were major sources of the exposures because of the rain or snow falls on 15 March in Fukushima. In the external dose calculation system, the reduction factors for cloud source in TECDOC 225 were adopted for the calculation from 12 to 14 March. The values used in this system were 0.9, 0.6, and 0.2 for one or two story wood frame house, one or two story concrete house or building, and three or more story concrete building, respectively. For the period of 15 March and after, the factors for deposited radioactivity were applied in the calculations. The values were 0.4, 0.2 and 0.1 for one or two story wood frame house, one or two story concrete house or building, and three or more story concrete building, respectively.

Dose rates during a move may change depending on the location. The dose exposed during the move is simply calculated as the product of mean effective dose rate at the start and end points multiplying the hours of the move.

Finally, when the mesh number including the place, j is m(j), and dose reduction factor correspond to the type of the building, k is r, dose for stay on the date, h, E and dose for the move on the date, h, E express as the following equations, respectively. The effective dose from March 12 to July 11, 2011, E, can be obtained with the following equation.

The conversion coefficients from ambient dose equivalent to effective dose E/H*(10) are higher as the age is lower, because the shielding effect increases with the volumes of body stem18. When the maximum conversion coefficient was used for a conservative estimation, conversion coefficient, Cage (age coefficient) was expressed as the following linear function of age, y, between 1 and 15 years old (Fig. 4). The effective dose for children can be calculated by multiplying E by Cage.

Figure 4

Ratios of conversion coefficients for each age group to adult.

Eighteen patterns of the evacuations

For roughly grasping and estimating the doses, 18 evacuation patterns of the residents from the area within 20 km radius of the NPP, and the deliberate evacuation area2 were assumed by considering the actual evacuation information before using the survey data (Fig. 5).

Figure 5

Assumed 18 evacuation patterns of the residents.

Basic maps were created by using ArcGIS of ESRI.

Since the traffic conditions were worse compared to those in ordinary situation, the moving speeds of evacuations were assumed to be 10 km/h for the moving inside the prefecture, and 30 km/h for those from Fukushima to other prefectures. Until the end of March, the ordinary traffic conditions were considered. For the time of activities outside the houses or buildings, it is assumed that the residents stayed for 4 hours a day before the evacuations, and 2 hours a day after the evacuations.

The two places where the residents started to evacuate were assumed; the representative place such as town office, and/or the places where the doses were highest in the town until when they evacuated.

Author Contributions

M.A. conducted this project. All authors performed the conceptual design of this system, and K.A. and S.Y. have developed the current system in the Fukushima health management survey. K.A. and S.Y. wrote the text and made tables and figures. All authors reviewed the manuscript.