Predicting the long-term (137)Cs distribution in Fukushima after the Fukushima Dai-ichi nuclear power plant accident: a parameter sensitivity analysis.

Radioactive materials deposited on the land surface of Fukushima Prefecture from the Fukushima Dai-ichi Nuclear Power Plant explosion is a crucial issue for a number of reasons, including external and internal radiation exposure and impacts on agricultural environments and aquatic biota. Predicting the future distribution of radioactive materials and their fates is therefore indispensable for evaluation and comparison of the effectiveness of remediation options regarding human health and the environment. Cesium-137, the main radionuclide to be focused on, is well known to adsorb to clay-rich soils; therefore its primary transportation mechanism is in the form of soil erosion on the land surface and transport of sediment-sorbed contaminants in the water system. In this study, we applied the Soil and Cesium Transport model, which we have developed, to predict a long-term cesium distribution in the Fukushima area, based on the Universal Soil Loss Equation and simple sediment discharge formulas. The model consists of calculation schemes of soil erosion, transportation and deposition, as well as cesium transport and its future distribution. Since not all the actual data on parameters is available, a number of sensitivity analyses were conducted here to find the range of the output results due to the uncertainties of parameters. The preliminary calculation indicated that a large amount of total soil loss remained in slope, and the residual sediment was transported to rivers, deposited in rivers and lakes, or transported farther downstream to the river mouths. Most of the sediment deposited in rivers and lakes consists of sand. On the other hand, most of the silt and clay portions transported to river were transported downstream to the river mouths. The rate of sediment deposition in the Abukuma River basin was three times as high as those of the other 13 river basins. This may be due to the larger catchment area and more moderate channel slope of the Abukuma River basin than those of the other rivers. Annual sediment outflows from the Abukuma River and the total from the other 13 river basins were calculated as 3.2 × 10(4)-3.1 × 10(5) and 3.4 × 10(4)-2.1 × 10(5)ty(-1), respectively. The values vary between calculation cases because of the critical shear stress, the rainfall factor, and other differences. On the other hand, contributions of those parameters were relatively small for (137)Cs concentration within transported soil. This indicates that the total amount of (137)Cs outflow into the ocean would mainly be controlled by the amount of soil erosion and transport and the total amount of (137)Cs concentration remaining within the basin. Outflows of (137)Cs from the Abukuma River and the total from the other 13 river basins during the first year after the accident were calculated to be 2.3 × 10(11)-3.7 × 10(12) and 4.6 × 10(11)-6.5 × 10(12)Bqy(-1), respectively. The former results were compared with the field investigation results, and the order of magnitude was matched between the two, but the value of the investigation result was beyond the upper limit of model prediction.

[1]  G. Richard,et al.  Estimation of runoff critical shear stress for soil erosion from soil shear strength , 2004 .

[2]  G. Hammer,et al.  Impact of soil erosion on production in cropping systems .I. Development and validation of a simulation model , 1992 .

[3]  Masaki Kanamori,et al.  Export of 134 Cs and 137 Cs in the Fukushima river systems at heavy rains by Typhoon Roke in September 2011 , 2013 .

[4]  D. Walling,et al.  Interpreting particle size effects in the adsorption of 137Cs and unsupported 210Pb by mineral soils and sediments , 1996 .

[5]  Calvin Wyatt Rose,et al.  Application and sensitivity analysis of process-based erosion model GUEST , 1996 .

[6]  C. Adachi,et al.  CORRIGENDUM: Promising operational stability of high-efficiency organic light-emitting diodes based on thermally activated delayed fluorescence , 2014, Scientific Reports.

[7]  B. W. Simpson,et al.  A simulation study of erosion in the Emerald Irrigation Area , 1999 .

[8]  Y. Onda,et al.  Depth distribution of 137Cs, 134Cs, and 131I in soil profile after Fukushima Dai-ichi Nuclear Power Plant Accident. , 2012, Journal of environmental radioactivity.

[9]  John R. Williams,et al.  LARGE AREA HYDROLOGIC MODELING AND ASSESSMENT PART I: MODEL DEVELOPMENT 1 , 1998 .

[10]  S. Takeuchi,et al.  Development of a model to predict a radionuclide distribution based on soil migration after Fukushima Dai-ichi Nuclear Power Plant accident , 2013 .

[11]  R. Young,et al.  AGNPS: A nonpoint-source pollution model for evaluating agricultural watersheds , 1989 .

[12]  L. F. Huggins,et al.  ANSWERS: A Model for Watershed Planning , 1980 .

[13]  R. Jackson,et al.  Partitioning of strontium-90 among aqueous and mineral species in a contaminated aquifer. , 1983, Environmental science & technology.

[14]  G. R. Foster,et al.  WEPP: A new generation of erosion prediction technology , 1991 .

[15]  Yosuke Yamashiki,et al.  Initial flux of sediment-associated radiocesium to the ocean from the largest river impacted by Fukushima Daiichi Nuclear Power Plant , 2014, Scientific Reports.

[16]  R. Morgan Soil Erosion and Conservation , 1988 .

[17]  W. G. Knisel,et al.  CREAMS: a field scale model for Chemicals, Runoff, and Erosion from Agricultural Management Systems [USA] , 1980 .

[18]  Tadatoshi Yamamoto,et al.  Effect of soil characteristics on soil particle size dependence of distribution coefficient. , 1988 .

[19]  Anthony J. Jakeman,et al.  A review of erosion and sediment transport models , 2003, Environ. Model. Softw..

[20]  H. Yamazawa,et al.  Preliminary Estimation of Release Amounts of 131I and 137Cs Accidentally Discharged from the Fukushima Daiichi Nuclear Power Plant into the Atmosphere , 2011 .

[21]  Yasuyuki Shimizu,et al.  MUTUAL EFFECTS OF BED AND BANK DEFORMATION IN CHANNEL PLANE FORMATION , 2003 .

[22]  S. Churakov,et al.  Up-Scaling of Molecular Diffusion Coefficients in Clays: A Two-Step Approach , 2011 .

[23]  H. Beck Environmental gamma radiation from deposited fission products, 1960-1964. , 1966, Health physics.

[24]  W. H. Wischmeier,et al.  Predicting rainfall erosion losses : a guide to conservation planning , 1978 .

[25]  Y. Onishi,et al.  Mathematical Modeling of Radioactive Contaminants in the Fukushima Environment , 2015 .