Numerical modeling was used to provide a new estimate of the amount of 137Cs released directly into the ocean from the Fukushima Daiichi nuclear power plant (NPP) after the accident in March 2011 and to gain insights into the physical processes that led to its dispersion in the marine environment during the months following the accident. An inverse method was used to determine the time‐dependent137Cs input responsible for the concentrations observed at the NPP's two liquid discharge outlets. The method was then validated through comparisons of the simulated concentrations with concentrations measured in seawater at different points in the neighborhood of the plant. An underestimation was noticed for stations located 30 km offshore. The resulting bias in the release inventory was estimated. Finally, the maximum 137Cs activity released directly to the ocean was estimated to lie between 5.1 and 5.5 PBq (Peta Becquerel = 1015 Bq) but uncertainties remain on the amount of radionuclides released during the first few days after the accident. This estimate was compared to previous ones and differences were analyzed further. The temporal and spatial variations of the 137Cs concentration present in the coastal waters were shown to be strongly related to the wind intensity and direction. During the first month after the accident, winds blowing toward the south confined the radionuclides directly released into the ocean to a narrow coastal band. Afterwards, frequent northward wind events increased the dispersion over the whole continental shelf, leading to strongly reduced concentrations.

Abstract. We develop and implement a new method to take into account the impact of waves into the 3-D circulation model SYMPHONIE (Marsaleix et al., 2008(Marsaleix et al., , 2009a following the simplified equations of , which use glm2z-RANS theory (Ardhuin et al., 2008c). These adiabatic equations are completed by additional parameterizations of wave breaking, bottom friction and wave-enhanced vertical mixing, making the forcing valid from the surf zone through to the open ocean. The wave forcing is performed by wave generation and propagation models WAVEWATCH III ® (Tolman, 2008(Tolman, , 2009Ardhuin et al., 2010) and SWAN (Booij et al., 1999). The model is tested and compared with other models for a plane beach test case, previously tested by Haas and Warner (2009) Finally, a realistic case is simulated with energetic waves travelling over a coast of the Gulf of Lion (in the northwest of the Mediterranean Sea) for which currents are available at different depths as well as an accurate bathymetric database of the 0-10 m depth range. A grid nesting approach is used to account for the different forcings acting at different spatial scales. The simulation coupling the effects of waves and currents is successful to reproduce the powerful northward littoral drift in the 0-15 m depth zone. More precisely, two distinct cases are identified: When waves have a normal angle of incidence with the coast, they are responsible for complex circulation cells and rip currents in the surf zone, and when they travel obliquely, they generate a northward littoral drift. These features are more complicated than in the test cases, due to the complex bathymetry and the consideration of wind and non-stationary processes. Wave impacts in the inner shelf are less visible since wind and regional circulation seem to be the predominant forcings. Besides, a discrepancy between model and observations is noted at that scale, possibly linked to an underestimation of the wind stress.This three-dimensional method allows a good representation of vertical current profiles and permits the calculation of the shear stress associated with waves and currents. Future work will focus on the combination with a sediment transport model.

We develop and implement a new method to take into account the impact of waves into the 3-D circulation model SYMPHONIE (Marsaleix et al., 2008, 2009a), following the simplified equations of Bennis et al. (2011) which use glm2z-RANS theory (Ardhuin et al., 2008b). These adiabatic equations are completed by additional parameterizations of wave breaking, bottom friction and wave-enhanced vertical mixing, making the forcing valid from the surf zone through to the open ocean. The wave forcing is performed by wave generation and propagation models WAVEWATCH III® (Tolman, 2008, 2009; Ardhuin et al., 2010) and SWAN (Booij et al., 1999). The model is tested and compared with other models for a plane beach test case, previously tested by Haas and Warner (2009) and Uchiyama et al. (2010). A comparison is also made with the laboratory measurements of Haller et al. (2002) of a barred beach with channels. Results fit with previous simulations performed by other models and with available observational data. <br><br> Finally, a realistic case of energetic waves travelling over a coast of the Gulf of Lion (in the northwest of the Mediterranean Sea) for which currents are available at different depths as well as an accurate bathymetric database of the 0–10 m depth range, is then simulated. A grid nesting approach is used to account for the different forcings acting at different spatial scales. The simulation coupling the effects of waves and currents is successful to reproduce the powerful northward littoral drift in the 0–15 m depth zone. More precisely, two distinct cases are identified: when waves have a normal angle of incidence with the coast, they are responsible for complex circulation cells and rip currents in the surf zone, and when they travel obliquely, they generate a northward littoral drift. These features are more complicated than in the test cases, due to the complex bathymetry and the consideration of wind and non-stationary processes. Wave impacts in the inner shelf are less visible since wind and regional circulation seem to be the predominant forcings. Besides, a discrepancy between model and observations is noted at that scale, possibly linked to an underestimation of the wind stress. Lastly, this three-dimensional method allows a good representation of vertical current profiles and permits to calculate the shear stress associated with wave and current. Future work will focus on the combination with a sediment transport model

A phase-resolved wave model is derived from an ocean circulation model for the purpose of studying wave-current effects in nearshore zones. One challenge is to adapt the circulation model to the specificities of wave physics. This mainly concerns the consideration of nonhydrostatic effects and the parametrization of wave breaking. The non-hydrostatic pressure is calculated using the artificial compressibility method (ACM). The ACM-induced errors on wave dispersion properties are examined in detail in the context of the linear theory using idealized test cases. The possible compromise between the precision achieved on nonhydrostatic physics and the adjustable CPU cost of the ACM method is looked at in detail.The modification of the wave characteristics by the bathymetric slope and the breaking of waves are then examined from a linear slope beach laboratory experiment. Finally the model is evaluated on the issue of rip currents and their feedback on the wave field using a laboratory experiment of a beach with a bar intersected by channels.