A dynamic model to estimate the dose rate of marine biota (K-BIOTADYN- M) and its application to the Fukushima accident

This paper describes a dynamic compartment model, K-BIOTA-DYN-M, to assess the activity concentration and dose rate of marine biota when the seawater activity varies with time, which is likely for the early phase after an accident. The model consists of seven compartments, phytoplankton, zooplankton, prey fish, benthic fish, crustacean, mollusk, and macro-algae. The phytoplankton compartment is assumed to be instantaneously in equilibrium with the seawater owing to the huge mass of the plankton in sea, and thus the activity of the phytoplankton is estimated using the equilibrium concentration ratio. The other compartments intake the radioactivity from both water and food, and lose the radioactivity by the biological elimination and radioactivity decay. Given the seawater activity, a set of ordinary differential equations representing the activity balance for biota is solved to obtain the time-variant activity concentration of biota, which is subsequently used to calculate the internal dose rate. The key parameters include the water intake rate, the daily feeding rate, the assimilation efficiency of radionuclides from food, the occupancy factor, and so on. The model has been applied to predict the activity concentration and dose rate of marine biota as a result the Fukushima nuclear accident on March 11, 2011. Using the seawater activities measured at three locations near the Fukushima NPPs, the time-variant activity concentration and dose rate during a few months after an accident for the seven model biota have been estimated. The preliminary results showed that the activity concentration of 137Cs in fish inhabiting the sea close to the Fukushima Daiichi NPP increased up to tenth-thousands of Bq/kg around the peak time of the seawater activity. This level is much higher than the food consumption restriction level for human protection; however, the estimated total dose rates (internal + external) for biota during the entire simulation time were all much less than the UNSCEAR's benchmark level to protect the aquatic biota from the ionizing radiation, indicating that the integrity of marine ecosystem would be maintained. In addition, the simulation results will be compared with the measured activity concentration of the biota as well as the results predicted by the equilibrium model. (authors)