ALBERT

Description: The Tōhoku earthquake and tsunami of March 11, 2011, resulted in unprecedented radioactivity releases from the Fukushima Dai-ichi nuclear power plants to the Northwest Pacific Ocean. Results are presented here from an international study of radionuclide contaminants in surface and subsurface waters, as well as in zooplankton and fish, off Japan in June 2011. A major finding is detection of Fukushima-derived 134Cs and 137Cs throughout waters 30–600 km offshore, with the highest activities associated with near-shore eddies and the Kuroshio Current acting as a southern boundary for transport. Fukushima-derived Cs isotopes were also detected in zooplankton and mesopelagic fish, and unique to this study we also find 110mAg in zooplankton. Vertical profiles are used to calculate a total inventory of ∼2 PBq 137Cs in an ocean area of 150,000 km2. Our results can only be understood in the context of our drifter data and an oceanographic model that shows rapid advection of contaminants further out in the Pacific. Importantly, our data are consistent with higher estimates of the magnitude of Fukushima fallout and direct releases [Stohl et al. (2011) Atmos Chem Phys Discuss 11:28319–28394; Bailly du Bois et al. (2011) J Environ Radioact, 10.1016/j.jenvrad.2011.11.015]. We address risks to public health and marine biota by showing that though Cs isotopes are elevated 10–1,000× over prior levels in waters off Japan, radiation risks due to these radionuclides are below those generally considered harmful to marine animals and human consumers, and even below those from naturally occurring radionuclides.

Description: Radium has four naturally occurring isotopes that have proven useful in constraining water mass source, age, and mixing rates in the coastal and open ocean. In this study, we used radium isotopes to determine the fate and flux of runoff-derived cesium from the Fukushima Nuclear Power Plant (NPP). During a June 2011 cruise, the highest Cs concentrations were found along the eastern shelf of northern Japan, from Fukushima south, to the edge of the Kuroshio current, and in an eddy ∼ 130 km from the NPP site. Locations with the highest cesium also had some of the highest radium activities, suggesting much of the direct ocean discharges of Cs remained in the coastal zone 2–3 months after the accident. We used a short-lived Ra isotope (223Ra, t1/2 = 11.4 d) to derive an average water mass age (Tr) in the coastal zone of 32 days. To ground-truth the Ra age model, we conducted a direct, station-by-station comparison of water mass ages with a numerical oceanographic model and found them to be in excellent agreement (model avg. Tr = 27 days). From these independent Tr values and the inventory of Cs within the water column at the time of our cruise, we were able to calculate an offshore 134Cs flux of 3.9–4.6 × 1013 Bq d−1. Radium-228 (t1/2 = 5.75 yr) was used to derive a vertical eddy diffusivity (Kz) of 0.7 m2 d−1 (0.1 cm2 s−1); from this Kz and 134Cs inventory, we estimated a 134Cs flux across the pycnocline of 1.8 × 104 Bq d−1 for the same time period. On average, our results show that horizontal mixing loss of Cs from the coastal zone was ∼ 109 greater than vertical exchange below the surface mixed layer. Finally, a mixing/dilution model that utilized our Ra-based and oceanographic model water mass ages produced a direct ocean discharge of 134Cs from the FNPP of 11–16 PBq at the time of the peak release in early April 2011. Our results can be used to calculate discharge of other water-soluble radionuclides that were released to the ocean directly from the Fukushima NPP.

Description: The March of 2011 earthquake and tsunami that caused a loss of power at the Fukushima nuclear power plants (FNPP) resulted in emission of radioactive isotopes into the atmosphere and the ocean. In June of 2011, an international survey of various radionuclide isotopes, including 137Cs, was conducted in surface and subsurface waters off Japan. This paper presents the results of numerical simulations aimed at interpreting these observations, investigating the spread of Fukushima-derived radionuclides off the coast of Japan and into the greater Pacific Ocean, studying the dominant mechanisms governing this process, as well as estimating the total amount of radionuclides in discharged coolant waters and atmospheric airborne radionuclide fallout. The numerical simulations are based on two different ocean circulation models, one inferred from AVISO altimetry and NCEP/NCAR reanalysis wind stress, and the second generated numerically by the NCOM model. Our simulations determine that 〉95% of 137Cs remaining in the water within ~600 km of Fukushima, Japan in mid-June 2011 was due to the direct oceanic discharge. The estimated strength of the oceanic source is 16.2 ± 1.6 PBq, based on minimizing the model-data mismatch. We cannot make an accurate estimate for the atmospheric source strength since most of the fallout cesium would have moved out of the survey area by mid-June. The model explained several features of the observed 137Cs distribution. First, the absence of 137Cs at the southernmost stations is attributed to the Kuroshio Current acting as a transport barrier against the southward progression of 137Cs. Second, the largest 137Cs concentrations were associated with a semi-permanent eddy that entrained 137Cs-rich waters collecting and stirring them around the eddy perimeter. Finally, the intermediate 137Cs concentrations at the westernmost stations were attributed to younger, and therefore less Cs-rich, coolant waters that continued to leak from the reactor in June of that year.

Description: Radium has four naturally occurring isotopes that have proven useful in constraining water mass source, age, and mixing rates in the coastal and open ocean. In this study, we used radium isotopes to determine the fate and flux of runoff-derived cesium from the Fukushima Dai-ichi Nuclear Power Plant (FNPP). During a June 2011 cruise, the highest cesium (Cs) concentrations were found along the eastern shelf of northern Japan, from Fukushima south, to the edge of the Kuroshio Current, and in an eddy ~ 130 km from the FNPP site. Locations with the highest cesium also had some of the highest radium activities, suggesting much of the direct ocean discharges of Cs remained in the coastal zone 2–3 months after the accident. We used a short-lived Ra isotope (223Ra, t1/2 = 11.4 d) to derive an average water mass age (Tr) in the coastal zone of 32 days. To ground-truth the Ra age model, we conducted a direct, station-by-station comparison of water mass ages with a numerical oceanographic model and found them to be in excellent agreement (model avg. Tr = 27 days). From these independent Tr values and the inventory of Cs within the water column at the time of our cruise, we were able to calculate an offshore 134Cs flux of 3.9–4.6 × 1013 Bq d−1. Radium-228 (t1/2 = 5.75 yr) was used to derive a vertical eddy diffusivity (Kz) of 0.7 m2 d−1 (0.1 cm2 s−1); from this Kz and 134Cs inventory, we estimated a 134Cs flux across the pycnocline of 1.8 × 104 Bq d−1 for the same time period. On average, our results show that horizontal mixing loss of Cs from the coastal zone was ~ 109 greater than vertical exchange below the surface mixed layer. Finally, a mixing/dilution model that utilized our Ra-based and oceanographic model water mass ages produced a direct ocean discharge of 134Cs from the FNPP of 11–16 PBq at the time of the peak release in early April 2011. Our results can be used to calculate discharge of other water-soluble radionuclides that were released to the ocean directly from the Fukushima NPP.

Description: The Great East Japan Earthquake and tsunami that caused a loss of power at the Fukushima nuclear power plants (FNPP) resulted in emission of radioactive isotopes into the atmosphere and the ocean. In June of 2011, an international survey measuring a variety of radionuclide isotopes, including 137Cs, was conducted in surface and subsurface waters off Japan. This paper presents the results of numerical simulations specifically aimed at interpreting these observations and investigating the spread of Fukushima-derived radionuclides off the coast of Japan and into the greater Pacific Ocean. Together, the simulations and observations allow us to study the dominant mechanisms governing this process, and to estimate the total amount of radionuclides in discharged coolant waters and atmospheric airborne radionuclide fallout. The numerical simulations are based on two different ocean circulation models, one inferred from AVISO altimetry and NCEP/NCAR reanalysis wind stress, and the second generated numerically by the NCOM model. Our simulations determine that 〉 95% of 137Cs remaining in the water within ~600 km of Fukushima, Japan in mid-June 2011 was due to the direct oceanic discharge. The estimated strength of the oceanic source is 16.2 ± 1.6 PBq, based on minimizing the model-data mismatch. We cannot make an accurate estimate for the atmospheric source strength since most of the fallout cesium had left the survey area by mid-June. The model explained several key features of the observed 137Cs distribution. First, the absence of 137Cs at the southernmost stations is attributed to the Kuroshio Current acting as a transport barrier against the southward progression of 137Cs. Second, the largest 137Cs concentrations were associated with a semi-permanent eddy that entrained 137Cs-rich waters, collecting and stirring them around the eddy perimeter. Finally, the intermediate 137Cs concentrations at the westernmost stations are attributed to younger, and therefore less Cs-rich, coolant waters that continued to leak from the reactor in June of that year.

Description: We have found a new way to express the solutions of the RSM (Reynolds Stress Model) equations that allows us to present the turbulent diffusivities for heat, salt and momentum in a way that is considerably simpler and thus easier to implement than in previous work. The RSM provides the dimensionless mixing efficiencies Gamma-alpha (alpha stands for heat, salt and momentum). However, to compute the diffusivities, one needs additional information, specifically, the dissipation Epsilon. Since a dynamic equation for the latter that includes the physical processes relevant to the ocean is still not available, one must resort to different sources of information outside the RSM to obtain a complete Mixing Scheme usable in OGCMs. As for the RSM results, we show that the Gamma-alpha s are functions of both Ri and Rq (Richardson number and density ratio representing double diffusion, DD); the Gamma-alpha are different for heat, salt and momentum; in the case of heat, the traditional value Gamma-h = 0.2 is valid only in the presence of strong shear (when DD is inoperative) while when shear subsides, NATRE data show that Gamma-h can be three times as large, a result that we reproduce. The salt Gamma-s is given in terms of Gamma-h. The momentum Gamma-m has thus far been guessed with different prescriptions while the RSM provides a well defined expression for Gamma-m(Ri,R-rho). Having tested Gamma-h, we then test the momentum Gamma-m by showing that the turbulent Prandtl number Gamma-m/Gamma-h vs. Ri reproduces the available data quite well. As for the dissipation epsilon, we use different representations, one for the mixed layer (ML), one for the thermocline and one for the ocean;s bottom. For the ML, we adopt a procedure analogous to the one successfully used in PB (planetary boundary layer) studies; for the thermocline, we employ an expression for the variable epsilon/N(exp 2) from studies of the internal gravity waves spectra which includes a latitude dependence; for the ocean bottom, we adopt the enhanced bottom diffusivity expression used by previous authors but with a state of the art internal tidal energy formulation and replace the fixed Gamma-alpha = 0.2 with the RSM result that brings into the problem the Ri, R-rho dependence of the Gamma-alpha; the unresolved bottom drag, which has thus far been either ignored or modeled with heuristic relations, is modeled using a formalism we previously developed and tested in PBL studies. We carried out several tests without an OGCM. Prandtl and flux Richardson numbers vs. Ri. The RSM model reproduces both types of data satisfactorily. DD and Mixing efficiency Gamma-h(Ri,Rq). The RSM model reproduces well the NATRE data. Bimodal epsilon-distribution. NATRE data show that epsilon (Ri 〈 1) approximately equals 10epsilon(Ri 〉 1), which our model reproduces. Heat to salt flux ratio. In the Ri much greater than 1 regime, the RSM predictions reproduce the data satisfactorily. NATRE mass diffusivity. The z-profile of the mass diffusivity reproduces well the measurements at NATRE. The local form of the mixing scheme is algebraic with one cubic equation to solve.