Search

An accurate analysis of local-scale atmospheric dispersion of radioactive materials is important for safety and consequence assessments and emergency responses to accidental release from nuclear facilities. It is necessary to predict the three-dimensional distribution of the plume in consideration of turbulent effects induced by individual buildings and meteorological conditions. In this study, first, we conducted with meteorological observations by a Doppler LiDAR and simple plume release experiments by a mist-spraying system at the site of Japan Atomic Energy Agency. Then, we developed a framework for prediction system of local-scale atmospheric dispersion based on a coupling of large-eddy simulation (LES) database and on-site meteorological observation. The LES-database was also created by pre-calculating high-resolution turbulent flows in the target site at mean wind directions of class interval 10°. We provided the meteorological observed data with the LES-database in consideration of building conditions and calculated the three-dimensional distribution of the plume with a Lagrangian dispersion model. Compared to the instantaneous shots of the plume taken by a digital camera, it was shown that the mist plume transport direction was accurately simulated. It was concluded that our proposed framework for prediction system based on a coupling of LES-database and on-site meteorological observation is effective.

The third model intercomparison project (MIP) for investigating the atmospheric behavior of atmospheric caesium-137 (137Cs) emitted from Fukushima Daiichi Nuclear Power Plant (FDNPP) (3rd FDNPP-MIP), Japan, in March 2011, was conducted. A finer horizontal grid spacing (1 km) was used than in the previous FDNPP-MIP (2nd FDNPP-MIP, Sato et al., 2018; 3 km) to evaluate the models’ performance for high-concentration events measured near FDNPP. Nine of the models used in the 2nd FDNPP-MIP were also used in the 3rd FDNPP-MIP, and all models used identical source terms and meteorological fields. The performance of the models was evaluated through a comparison with observational data. Our analyses indicated that most of the observed high atmospheric 137Cs concentrations (plumes) were reasonably well simulated by the models, and the good performance of some models improved the performance of the multimodel, highlighting the advantage of using a multimodel ensemble. The analyses also confirmed that the use of a finer grid resolution resulted in the meteorological field near FDNPP being better reproduced in the 3rd FDNPP-MIP, and the performance of the models was better than that of the 2nd FDNPP MIP. The good representation of the wind field resulted in the reasonable simulation of the narrow distribution of high deposition amount to the northwest of FDNPP and the reduction of the overestimation in deposition amount over the area to the south of FDNPP compared to the 2nd FDNPP MIP. In contrast, the performance of the models in simulating plumes observed over the Nakadori area, the northern part of Gunma, and the Tokyo Metropolitan Area (TMA) was slightly worse than in the 2nd FDNPP-MIP.

The long-lived radioactive iodine ($^{129}$I) is a useful geochemical tracer in the atmospheric environment. We recently observed clear seasonal trends in air concentration and deposition of $^{129}$I in Japan. Using these data, we developed a global atmospheric $^{129}$I transport model to reveal key processes for the global atmospheric $^{129}$I cycle. The model generally reproduced the observed seasonal change in air concentration and deposition of $^{129}$I in Japan, and the global distribution of $^{129}$I concentration in rain as presented in past literature. Numerical experiments changing the intensity of anthropogenic and natural sources were conducted to quantify the impact of anthropogenic sources on the global $^{129}$I cycle. The results indicated that the atmospheric $^{129}$I from the anthropogenic sources was deposited in winter and can be accumulated mainly in the northern part of Eurasia. In contrast, the atmospheric $^{129}$I from the natural sources dominated the deposition in summer. These results suggested that the re-emission process of $^{129}$I from the Earth's surface may be important as a secondary impact of $^{129}$I in the global-scaled environment. Furthermore, although wet deposition dominated the total deposition in the Northern hemisphere, dry deposition regionally and seasonally contributed to the total deposition over arctic and northern part of Eurasia in winter, suggesting that the dry deposition may play a key role in the seasonal change of $^{129}$I deposition in the Northern hemisphere high latitudes., 長寿命放射性ヨウ素($^{129}$I)は、大気環境における放射性核種の有用な地球化学トレーサである。本研究では、$^{129}$Iの大気濃度および沈着の観測を実施し、観測データから大気濃度および沈着の明瞭な季節変動を得た。さらに、大気中の$^{129}$I循環を支配する要因を明らかにすることを目的として、得られた観測データを用いて、移流、乱流拡散、大気沈着、光化学、ガス粒子変換、核燃料再処理工場からの$^{129}$Iの排出、海洋および陸域からの$^{129}$Iの揮発の各物理・化学過程を考慮した全球ヨウ素輸送モデルを開発した。全球ヨウ素輸送モデルは、我々が観測した$^{129}$Iの大気濃度および沈着の季節変動、そして既往文献の$^{129}$Iの降水中濃度の全球分布を良好に再現した。開発した全球ヨウ素輸送モデルを用いて人為起源と自然起源の$^{129}$Iインベントリの強度を変化させる数値実験を実施し、地球全体の$^{129}$I循環に対する人為起源の$^{129}$Iの影響を評価した。その結果、冬季においては、人為起源の$^{129}$Iが主にユーラシアの北部に沈着する可能性があることが示された。一方で、夏季においては、自然起源の$^{129}$Iが北半球中高緯度の沈着に支配的であった。これらの結果は、地球表面からの$^{129}$Iの再飛散過程が全球規模での$^{129}$I循環に重要であることを示唆している。さらに、冬季のユーラシア北部や北極域においては局所的に乾性沈着が寄与しており、乾性沈着が環境中の$^{129}$Iの季節変化に重要な影響を及ぼすことが示唆された。

When radioactive materials are released into the atmosphere due to nuclear accidents, numerical simulations that can reproduce temporal and spatial distribution of radioactive materials are useful to provide the information for emergency responses and radiological dose assessment. In this study, we attempt to improve the atmospheric dispersion simulation using an advanced meteorological data assimilation method and reconstruct the spatiotemporal distribution of radioactive materials released due to the Fukushima Daiichi Nuclear Power Station (FDNPS) accident. The atmospheric dispersion simulations were carried out by the Lagrangian particle dispersion model GEARN developed by Japan Atomic Energy Agency. To obtain meteorological fields for GEARN calculation, we used the Weather Research and Forecasting model WRF with meteorological data assimilation using four-dimensional variational method (4D-Var). GEARN calculations of the surface deposition and air concentration of radionuclides were compared with measurements. In the area close to FDNPS, the spatial distribution of the deposition of Cs-137 and I-131 simulated by GEARN agreed with the measured one. The accuracy of modeled deposition in northwest and south directions from FDNPS was particularly improved. This results were mainly attributed to the better reproducibility of wind field by using the meteorological data assimilation with 4D-Var. The improvement of the accuracy of modeled deposition distribution of Cs-137 in the East Japan area was also apparent under the meteorological fields modified by 4D-Var. The information of atmospheric dispersion processes reconstructed in this study is used for updating the existing assessment of radiological dose resulting from the FDNPS accident based on atmospheric simulations by our previous studies. It can also provide useful suggestions to make emergency response plans for nuclear facilities in Japan., 原子力施設事故によって大気中に放射性物質が放出された際に、放射性物質の時空間分布を再現可能な大気拡散シミュレーションは、緊急時対応や被ばく評価に対して有効である。本研究では、データ同化手法及び最新の領域気象モデルを用いて大気拡散シミュレーションを改良し、福島第一原子力発電所事故由来の放射性物質の時空間分布を再構築する。大気拡散シミュレーションは原子力機構が開発した大気拡散モデルGEARNによって行われた。大気拡散シミュレーションのための気象場は、4次元変分法を用いた気象データ同化解析手法と領域気象モデルWRFで計算された。シミュレーションの再現性は、放射性核種の沈着量及び大気濃度の計算値と実測値の比較から評価した。福島第一原子力発電所の付近において、観測されたCs-137とI-131の沈着量の空間分布をシミュレーションは良好に再現した。特に、発電所から北西方向及び南方向で観測されたCs-137の高い沈着量分布の再現性の向上が顕著であった。東日本の広域においても、モデル過大評価であったCs-137の沈着量の再現性が向上した。これらは、4次元変分法を用いた気象データ同化解析手法によって、再現性の高い風速場が得られたことに起因すると考えられる。本研究で再構築された大気拡散過程は、福島第一原子力発電所事故の線量影響評価の見直しや、日本国内の原子力施設における緊急時対応に対して有効な情報となる。

In this study, we developed an atmospheric dispersion model based on the finite difference method (GEARN-FDM) to simulate long-range radiological dispersion. The model includes a mass conservation and monotonic advection scheme and a horizontal diffusion scheme based on the deformation of the resolved horizontal wind. Meteorological fields in calculating dispersion were simulated by the Advanced Research Weather Research and Forecasting (WRF) model. ERA-interim reanalysis was used for the WRF nudging calculation. By using tracer gas concentrations from the European tracer experiment, the GEARN-FDM performance could be tested. The results indicated high performance with factors of 2 and 5 of 37% and 72%, respectively. To investigate the impact of horizontal diffusivity on the concentration distribution, we examined an additional simulation run using GEARN-FDM with a constant value of horizontal diffusivity. In comparison with the additional run, the plume calculated in the original GEARN-FDM run was locally and temporally diluted due to large horizontal diffusivity in a mountainous region. Therefore, geographical heterogeneity of horizontal diffusivity may affect regional-scale atmospheric dispersion simulations.

To assess the radiological dose to the public resulting from the Fukushima Daiichi Nuclear Power Station (FDNPS) accident in Japan, especially for the early phase of the accident when no measured data are available for that purpose, the spatial and temporal distributions of radioactive materials in the environment need to be reconstructed through computer simulations using the atmospheric transport, dispersion, and deposition model (ATDM). For the ATDM simulation, the source term of radioactive materials discharged into the atmosphere is essential and has been estimated in many studies. In the present study, we further refined the source term estimated in our previous study and improved the ATDM simulation with an optimization method based on Bayesian inference, which used various measurements such as air concentration, surface deposition, fallout, and newly released hourly air concentrations of 137Cs derived by analyzing suspended particulate matter (SPM) collected at air pollution monitoring stations. This optimization improved not only the source term but also the wind field in meteorological calculation, which led to the reduction of discrepancies in plume passage time at monitoring points to less than 3 h between calculations and measurements, by feeding back comparison results between the dispersion calculations and measurements of radionuclides. As a result, the total amounts of 137Cs and 131I by the present study became 1.0 × 1016 and 1.2 × 1017 Bq, respectively, and decreased by 29% and 20%, respectively, in comparison with those by previous study. The ATDM simulation successfully reproduced both the air concentrations at SPM monitoring points and surface depositions by airborne monitoring. FA10 for total samples of air concentrations of 137Cs at SPM monitoring points increased from 35.9% by the previous study to 47.3%. The deposition amount on the land decreased from 3.7 × 1015 Bq by the previous study to 2.1 × 1015 Bq, which was close to the measured amount of 2.4 × 1015 Bq. We also constructed the spatiotemporal distribution of some major radionuclides in the air and on the surface (optimized dispersion database) by using the optimized release rates and ATDM simulations. The optimized dispersion database can be used for comprehensive dose assessment in tandem with behavioral patterns of evacuees from the FDNPS accident by collaborating research group in the Japanese dose assessment project. The improvements in the present study lead to the refinement of the dose estimation.