Initial spread of 137Cs from the Fukushima Dai-ichi Nuclear Power Plant over the Japan continental shelf : a study using a high-resolution, global-coastal nested ocean model

© The Author(s), 2013. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Biogeosciences 10 (2013): 5439-5449, doi:10.5194/bg-10-5439-2013.The 11 March 2011 tsunami triggered by the M9 and M7.9 earthquakes off the Tōhoku coast destroyed facilities at the Fukushima Dai-ichi Nuclear Power Plant (FNPP) leading to a significant long-term flow of the radionuclide 137Cs into coastal waters. A high-resolution, global-coastal nested ocean model was first constructed to simulate the 11 March tsunami and coastal inundation. Based on the model's success in reproducing the observed tsunami and coastal inundation, model experiments were then conducted with differing grid resolution to assess the initial spread of 137Cs over the eastern shelf of Japan. The 137Cs was tracked as a conservative tracer (without radioactive decay) in the three-dimensional model flow field over the period of 26 March–31 August 2011. The results clearly show that for the same 137Cs discharge, the model-predicted spreading of 137Cs was sensitive not only to model resolution but also the FNPP seawall structure. A coarse-resolution (∼2 km) model simulation led to an overestimation of lateral diffusion and thus faster dispersion of 137Cs from the coast to the deep ocean, while advective processes played a more significant role when the model resolution at and around the FNPP was refined to ∼5 m. By resolving the pathways from the leaking source to the southern and northern discharge canals, the high-resolution model better predicted the 137Cs spreading in the inner shelf where in situ measurements were made at 30 km off the coast. The overestimation of 137Cs concentration near the coast is thought to be due to the omission of sedimentation and biogeochemical processes as well as uncertainties in the amount of 137Cs leaking from the source in the model. As a result, a biogeochemical module should be included in the model for more realistic simulations of the fate and spreading of 137Cs in the ocean.This project was supported by the US National Science Foundation RAPID grants No. 1141697 and No. 1141785 and the Japan Science and Technology Agency J-RAPID program. The development of Global-FVCOM was supported by NSF grants ARC0712903, ARC0732084, and ARC0804029. Z. Lai’s contribution was supported by the Natural Science Foundation of China project 41206005, China MOST project 2012CB956004, and Sun Yat-Sen University 985 grant 42000-3281301. C. Chen serves as chief scientist for the International Center for Marine Studies, Shanghai Ocean University, and his contribution was supported by the Program of Science and Technology Commission of Shanghai Municipality (09320503700)

Seasonal and interannual variability of the Arctic sea ice : a comparison between AO-FVCOM and observations

Author Posting. © American Geophysical Union, 2016. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 121 (2016): 8320–8350, doi:10.1002/2016JC011841.A high-resolution (up to 2 km), unstructured-grid, fully ice-sea coupled Arctic Ocean Finite-Volume Community Ocean Model (AO-FVCOM) was used to simulate the sea ice in the Arctic over the period 1978–2014. The spatial-varying horizontal model resolution was designed to better resolve both topographic and baroclinic dynamics scales over the Arctic slope and narrow straits. The model-simulated sea ice was in good agreement with available observed sea ice extent, concentration, drift velocity and thickness, not only in seasonal and interannual variability but also in spatial distribution. Compared with six other Arctic Ocean models (ECCO2, GSFC, INMOM, ORCA, NAME, and UW), the AO-FVCOM-simulated ice thickness showed a higher mean correlation coefficient of ∼0.63 and a smaller residual with observations. Model-produced ice drift speed and direction errors varied with wind speed: the speed and direction errors increased and decreased as the wind speed increased, respectively. Efforts were made to examine the influences of parameterizations of air-ice external and ice-water interfacial stresses on the model-produced bias. The ice drift direction was more sensitive to air-ice drag coefficients and turning angles than the ice drift speed. Increasing or decreasing either 10% in water-ice drag coefficient or 10° in water-ice turning angle did not show a significant influence on the ice drift velocity simulation results although the sea ice drift speed was more sensitive to these two parameters than the sea ice drift direction. Using the COARE 4.0-derived parameterization of air-water drag coefficient for wind stress did not significantly influence the ice drift velocity simulation.This work was supported by NSF grants OCE-1203393 for the UMASSD team and PLR-1203643 for R. C. Beardsley.2017-05-2

Surface circulation in Block Island Sound and adjacent coastal and shelf regions : a FVCOM-CODAR comparison

© The Author(s), 2016. This is the author's version of the work and is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Progress in Oceanography 143 (2016): 26-45, doi:10.1016/j.pocean.2016.02.005.CODAR-derived surface currents in Block Island Sound over the period of June 2000 through September 2008 were compared to currents computed using the Northeast Coastal Ocean Forecast System (NECOFS). The measurement uncertainty of CODAR-derived currents, estimated using statistics of a screened nine-year time series of hourly-averaged flow field, ranged from 3-7 cm/s in speed and 4°-14° in direction. The CODAR-derived and model-computed kinetic energy spectrum densities were in good agreement at subtidal frequencies, but the NECOFS-derived currents were larger by about 28% at semi-diurnal and diurnal tidal frequencies. The short-term (hourly to daily) current variability was dominated by the semidiurnal tides (predominantly the M2 tide), which on average accounted for ~87% of the total kinetic energy. The diurnal tidal and subtidal variability accounted for ~4% and ~9% of the total kinetic energy, respectively. The monthly-averaged difference between the CODAR-derived and model-computed velocities over the study area was 6 cm/s or less in speed and 28° or less in direction over the study period. An EOF analysis for the low-frequency vertically-averaged model current field showed that the water transport in the Block Island Sound region was dominated by modes 1 and 2, which accounted for 89% and 7% of the total variance, respectively. Mode 1 represented a relatively stationary spatial and temporal flow pattern with a magnitude that varied with season. Mode 2 was characterized mainly by a secondary cross-shelf flow and a relatively strong along-shelf flow. Process-oriented model experiments indicated that the relatively stationary flow pattern found in mode 1 was a result of tidal rectification and its magnitude changed with seasonal stratification. Correlation analysis between the flow and wind stress suggested that the cross-shelf water transport and its temporal variability in mode 2 were highly correlated to the surface wind forcing. The mode 2 derived onshore and offshore water transport, and was consistent with wind-driven Ekman theory. The along-shelf water transport over the outer shelf, where a large portion of the water flowed from upstream Nantucket Shoals, was not highly correlated to the surface wind stress.This work was supported by the NSF grants OCE-1332207 and OCE-1332666, MIT Sea Grant College Program through grant 2012-R/RC-127, and the NOAA NERACOOS program funds for NECOFS. Operational funding for the CODAR systems used in this study was provided by the Mid-Atlantic Regional Association Coastal Ocean Observing System. The development of the Global-FVCOM system has been supported by NSF grants OCE-1203393. C. Chen’s contribution was also supported by the International Center for Marine Studies at Shanghai Ocean University through the “Shanghai Universities First-class Disciplines Project”.2017-03-0

Impact of current-wave interaction on storm surge simulation : a case study for Hurricane Bob

Author Posting. © American Geophysical Union, 2013. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 118 (2013): 2685–2701, doi:10.1002/jgrc.20207.Hurricane Bob moved up the U.S. east coast and crossed over southern New England and the Gulf of Maine [with peak marine winds up to 54 m/s (100 mph)] on 19–20 August 1991, causing significant damage along the coast and shelf. A 3-D fully wave-current-coupled finite-volume community ocean model system was developed and applied to simulate and examine the coastal ocean responses to Hurricane Bob. Results from process study-oriented experiments showed that the impact of wave-current interaction on surge elevation varied in space and time, more significant over the shelf than inside the inner bays. While sea level change along the coast was mainly driven by the water flux controlled by barotropic dynamics and the vertically integrated highest water transports were essentially the same for cases with and without water stratification, the hurricane-induced wave-current interaction could generate strong vertical current shear in the stratified areas, leading to a strong offshore transport near the bottom and vertical turbulent mixing over the continental shelf. Stratification could also result in a significant difference of water currents around islands where the water is not vertically well mixed.This work was supported by the MIT Sea Grant College Program through grant 2012-R/RC-127 and the NOAA NERACOOS Program funds for NECOFS. The development of the FVCOM system has been supported by the NSF Ocean Sciences Division through grants OCE-0234545, OCE-0227679, OCE-0606928, and OCE- 0712903 and the NSF Office of Polar Programs-Arctic Sciences Division through grants ARC0712903, ARC0732084, ARC0804029, and ARC1203393. C.C.’s contribution was also supported by Shanghai Ocean University International Cooperation Program (A-2302-11-0003), the Program of Science and Technology Commission of Shanghai Municipality (09320503700), and the Leading Academic Discipline Project of Shanghai Municipal Education Commission (J50702).2013-11-3

A wet/dry point treatment method of FVCOM, part I: stability experiments

© The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Chen, C., Qi, J., Liu, H., Beardsley, R., Lin, H., & Cowles, G. A wet/dry point treatment method of FVCOM, part I: stability experiments. Journal of Marine Science and Engineering, 10(7), (2022): 896, https://doi.org/10.3390/jmse10070896.A 3-dimensional wet/dry point treatment method was developed for the unstructured-grid Finite-Volume Community Ocean Model (FVCOM). Analytical equations were derived to examine discretized errors that occurred during the flooding/drying process by the wet/dry point treatment for the flooding/drying process. Numerical experiments were carried out for an idealized estuary, including the inter-tidal zone. The model results show that if the ratio of internal to external mode time steps (Isplit) is appropriately selected, FVCOM was capable of simulating the flooding/drying process with sufficient accuracy to ensure the mass conservation. The up-bound limit of Isplit was restricted by the bathymetric slope of the inter-tidal zone, external mode time step, horizontal/vertical resolution, and amplitude of tidal forcing at the open boundary, as well as the thickness of the viscous layer specified in the model. Criteria for time steps via these parameters were derived from these experiments, which provide a helpful guide in selectingIsplit for applying FVCOM to realistic geometric estuaries.This research was funded by the Georgia Sea Grant College Program under grant numbers NA26RG0373 and NA66RG0282, the Georgia DNR grants 024409-01 and 026450-01, the NSF Georges Bank/Northwest Atlantic GLOBEC program under grant number NSF-OCE 02-27679, and the SMAST fishery program under the NASA grant number NAG 13-02042

Reply to comment on “Current separation and upwelling over the southeast shelf of Vietnam in the South China Sea”

Author Posting. © American Geophysical Union, 2013. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 118 (2013): 1624, doi:10.1002/jgrc.20114.2013-09-3

A new high-resolution unstructured grid finite volume Arctic Ocean model (AO-FVCOM) : an application for tidal studies

Author Posting. © American Geophysical Union, 2009. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 114 (2009): C08017, doi:10.1029/2008JC004941.A spherical coordinate version of the unstructured grid 3-D FVCOM (finite volume coastal ocean model) has been applied to the Arctic Ocean to simulate tides with a horizontal resolution ranging from 1 km in the near-coastal areas to 15 km in the deep ocean. By accurately resolving the irregular coastlines and bathymetry in the Arctic Ocean coastal regions, this model reproduces the diurnal (K1 and O1) and semidiurnal (M2 and S2) tidal wave dynamics and captures the complex tidal structure along the coast, particularly in the narrow straits of the Canadian Archipelago. The simulated tidal parameters (harmonic constituents of sea surface elevation and currents) agree well with the available observational data. High-resolution meshes over the continental shelf and slope capture the detailed spatial structure of topographic trapped shelf waves, which are quite energetic along the Greenland, Siberia, and Spitsbergen continental slope and shelf break areas. Water stratification influences the vertical distribution of tidal currents but not the water transport and thus tidal elevation. The comparison with previous finite difference models suggests that horizontal resolution and geometric fitting are two prerequisites to simulate realistically the tidal energy flux in the Arctic Ocean, particularly in the Canadian Archipelago.This research was supported by the NSF Office of Polar Programs through grants OPP ARC-0712903, ARC- 0732084, and ARC-0804029 for C. Chen, G. Gao, and G. Cowles; OPP ARC-0804010 and ARC-0712848 for A. Proshutinsky; OPP ANT-0523223, ARC0712848, NOAA Cooperative Agreement NA17RJ1223 (409) and the WHOI Smith Chair for R. C. Beardsley. J. Qi was supported by the SMAST fishery program under NOAA grants NA04NMF4720332 and NA05NMF4721131. The spherical coordinate version of FVCOM was developed with initial funds from NSF grants OCE-0606928 and OCE- 0726851. Gao was also supported by the Chinese NSF Arctic Ocean grant under contract 40476007

Studies of the Canadian Arctic Archipelago water transport and its relationship to basin-local forcings : results from AO-FVCOM

Author Posting. © American Geophysical Union, 2016. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 121 (2016): 4392–4415, doi:10.1002/2016JC011634.A high-resolution (up to 2 km), unstructured-grid, fully coupled Arctic sea ice-ocean Finite-Volume Community Ocean Model (AO-FVCOM) was employed to simulate the flow and transport through the Canadian Arctic Archipelago (CAA) over the period 1978–2013. The model-simulated CAA outflow flux was in reasonable agreement with the flux estimated based on measurements across Davis Strait, Nares Strait, Lancaster Sound, and Jones Sounds. The model was capable of reproducing the observed interannual variability in Davis Strait and Lancaster Sound. The simulated CAA outflow transport was highly correlated with the along-strait and cross-strait sea surface height (SSH) difference. Compared with the wind forcing, the sea level pressure (SLP) played a dominant role in establishing the SSH difference and the correlation of the CAA outflow with the cross-strait SSH difference can be explained by a simple geostrophic balance. The change in the simulated CAA outflow transport through Davis Strait showed a negative correlation with the net flux through Fram Strait. This correlation was related to the variation of the spatial distribution and intensity of the slope current over the Beaufort Sea and Greenland shelves. The different basin-scale surface forcings can increase the model uncertainty in the CAA outflow flux up to 15%. The daily adjustment of the model elevation to the satellite-derived SSH in the North Atlantic region outside Fram Strait could produce a larger North Atlantic inflow through west Svalbard and weaken the outflow from the Arctic Ocean through east Greenland.NSF Grant Numbers: OCE-1203393, PLR-1203643; National Natural Science Foundation of China Grant Number: 41276197; Shanghai Pujiang Program Grant Number: 12PJ1404100; Shanghai Shuguang Program2016-12-2

Observational and model studies of the circulation in the Gulf of Tonkin, South China Sea

Author Posting. © American Geophysical Union, 2013. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 118 (2013): 6495–6510, doi:10.1002/2013JC009455.Moored current measurements were made at one mooring site in the northern Gulf of Tonkin for about 1 year during 1988–1989. Analyses were performed to examine characteristics and variability of tidal and subtidal flows. Rotary spectra showed two peaks at diurnal and semidiurnal periods, with higher diurnal energy. Complex demodulations of diurnal and semidiurnal tidal currents indicated that the tidal current magnitudes varied significantly with seasons: more energetic in the stratified summer than in the vertically well-mixed winter. The observed subtidal currents were highly correlated with the surface wind in winter but not in summer; challenging the conceptual summertime anticyclonic circulation pattern derived using wind-driven homogenous circulation theory. The computed currents from a global ocean model were in good agreement with the observed currents. Similar to the current observations, the model-computed flow patterns were consistent with the conceptual wind-driven circulation pattern in winter but opposite in summer. Process-oriented experiments suggest that the summertime cyclonic circulation in the northern Gulf of Tonkin forms as a result of the combination of stratified wind-driven circulation and tidal-rectified inflow from Qiongzhou Strait. The interaction between the southwest monsoon and buoyancy-driven flow from Hong River can significantly intensify the cyclonic circulation near the surface, but its contribution to the vertically averaged flow of the cyclonic circulation is limited.Y. Ding has been supported by the State Scholarship Fund from the China Scholarship Council. C. Chen serves as chief scientist for the International Center for Marine Studies, Shanghai Ocean University, and his contribution has been supported by the Program of Science and Technology Commission of Shanghai Municipality (09320503700). C. Chen serves as the Zi Jiang Scholar at the State Key Laboratory for Estuarine and Coastal Research (SKLEC) of East China Normal University (ECNU) and Visiting Professor at School of Marine Sciences, Sun Yat-Sen University. C. Chen would like to credit this research to these two universities. Z. Lai’s contribution is supported by NSFC project 41206005 and Sun Yat-Sen University 985 grant 42000–3281301. The development of Global-FVCOM was funded by the US National Science Foundation Office of Polar Programs through grants ARC0712903, ARC0732084, ARC0804029, and ARC1203393.2014-06-0

Observed wintertime tidal and subtidal currents over the continental shelf in the northern South China Sea

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 119 (2014): 5289–5310, doi:10.1002/2014JC009931.Synthesis analyses were performed to examine characteristics of tidal and subtidal currents at eight mooring sites deployed over the northern South China Sea (NSCS) continental shelf in the 2006–2007 and 2009–2010 winters. Rotary spectra and harmonic analysis results showed that observed tidal currents in the NSCS were dominated by baroclinic diurnal tides with phases varying both vertically and horizontally. This feature was supported by the CC-FVCOM results, which demonstrated that the diurnal tidal flow over this shelf was characterized by baroclinic Kelvin waves with vertical phase differences varying in different flow zones. The northeasterly wind-induced southwestward flow prevailed over the NSCS shelf during winter, with episodic appearances of mesoscale eddies and a bottom-intensified buoyancy-driven slope water intrusion. The moored current records captured a warm-core anticyclonic eddy, which originated from the southwestern coast of Taiwan and propagated southwestward along the slope consistent with a combination of β-plane and topographic Rossby waves. The eddy was surface-intensified with a swirl speed of >50 cm/s and a vertical scale of ∼400 m. In absence of eddies and onshore deep slope water intrusion, the observed southwestward flow was highly coherent with the northeasterly wind stress. Observations did not support the existence of the permanent wintertime South China Sea Warm Current (SCSWC). The definition of SCSWC, which was based mainly on thermal wind calculations with assumed level of no motion at the bottom, needs to be interpreted with caution since the observed circulation over the NSCS shelf in winter included both barotropic and baroclinic components.R. Li was supported by the SOA 908 Special Project Foundation of China (908-01-ST07 and 908-01-BC10), the National High Tech Project Foundation (863) of China (2008AA09A401), the Administrator Foundation of South Branch, SOA (0683). The development of FVCOM was funded by the US NSF Office of Polar Programs through grants ARC0712903, ARC0732084, ARC0804029, and ARC1203393.2015-02-1