Numerical study of surface water circulation around Sekisei Lagoon, southwest Japan

This paper discusses the variability of surface currents around Sekisei Lagoon using a nested grid ocean circulation model. We developed a triple-nested grid system that consists of a coarse-resolution (1/60° or ∼1.85 km) model off Taiwan, an intermediate-resolution (1/300° or ∼370 m) model around the Yaeyama Islands, and a fine-resolution (1/900° or ∼123 m) model of Sekisei Lagoon. The nested grid system was forced by wind and heat flux calculated from six-hourly atmospheric reanalysis data and integrated over the period from May to July 2003. The coarse-resolution model was driven by lateral boundary conditions calculated from daily ocean reanalysis data to include realistic variation of the Kuroshio and mesoscale eddies with spatial scales of ∼500–700 km in the open ocean. The tidal forcing was included in the intermediate-resolution model by interpolating sea level data obtained from a data-assimilative tidal model. The results were then used to drive the fine-resolution model to simulate the surface water circulation around Sekisei lagoon. Model results show that (1) currents inside the lagoon are mainly driven by tide and wind; (2) there exists a strong southwestward current along the bottom slope in the southeast portion of the lagoon; the current is mainly driven by remote mesoscale eddies and at times intensified by the local wind; (3) the flow relaxation scheme is effective in reducing biases along the open boundaries. The simulated currents were used to examine the retention and dispersion of passive particles in the surface layer. Results show that the surface dispersion in the strong open ocean current region is significantly higher than that inside the lagoon.     

On the modifications of near-inertial waves at fronts: implications for energy transfer across scales

In the ocean, wind-generated kinetic energy (KE) manifests itself primarily in balanced currents and near-inertial waves. The dynamics of these flows is strongly constrained by the Earth’s rotation, causing the KE in balanced currents to follow an inverse cascade but also preventing wave-wave interactions from fluxing energy in the near-inertial band to lower frequencies and higher vertical wavenumbers. How wind-generated KE is transferred to small-scale turbulence and dissipated is thus a non-trivial problem. This article presents a review of recent theoretical calculations and numerical simulations that demonstrate how some surprising modifications to internal wave physics by the lateral density gradients present at ocean fronts allow for strong interactions between balanced currents and near-inertial waves that ultimately result in energy loss for both types of motion.     

Observation-based evaluation of surface wave effects on currents and trajectory forecasts

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.     

Observation-based evaluation of surface wave effects on currents and trajectory forecasts

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.

     

Wave effects in global ocean modeling: parametrizations vs. forcing from a wave model

One of the main challenges of the Copernicus Marine Service is the implementation of coupled ocean/waves systems that accurately estimate the momentum and energy fluxes provided by the atmosphere to the ocean. This study aims to investigate the impact of forcing the Nucleus for European Modelling of the Ocean (NEMO) ocean model with forecasts from the wave model of Météo-France (MFWAM) to improve classical air-sea flux parametrizations, these latter being mostly driven by the 10-m wind. Three wave-related processes, namely, wave-state-dependent stress, Stokes drift-related effects (Stokes-Coriolis force, Stokes drift advection on tracers and on mass), and wave-state-dependent surface turbulence, are examined at a global scale with a horizontal resolution of 0.25°. Three years of sensitivity simulations (2014–2016) show positive feedback on sea surface temperature (SST) and currents when the wave model is used. A significant reduction in SST bias is observed in the tropical Atlantic Ocean. This is mainly due to the more realistic momentum flux provided by the wave model. In mid-latitudes, the most interesting impact occurs during the summer stratification, when the wind is low and the wave model produces a reduction in the turbulence linked with wave breaking. Magnitudes of the large-scale currents in the equatorial region are also improved by 10% compared to observations. In general, it is shown that using the wave model reduces on average the momentum and energy fluxes to the ocean in tropical regions, but increases them in mid-latitudes. These differences are in the order of 10 to 20% compared with the classical parametrizations found in stand-alone ocean models.     

Adaptive volume penalization for ocean modeling

The development of various volume penalization techniques for use in modeling topographical features in the ocean is the focus of this paper. Due to the complicated geometry inherent in ocean boundaries, the stair-step representation used in the majority of current global ocean circulation models causes accuracy and numerical stability problems. Brinkman penalization is the basis for the methods developed here and is a numerical technique used to enforce no-slip boundary conditions through the addition of a term to the governing equations. The second aspect to this proposed approach is that all governing equations are solved on a nonuniform, adaptive grid through the use of the adaptive wavelet collocation method. This method solves the governing equations on temporally and spatially varying meshes, which allows higher effective resolution to be obtained with less computational cost. When penalization methods are coupled with the adaptive wavelet collocation method, the flow near the boundary can be well-resolved. It is especially useful for simulations of boundary currents and tsunamis, where flow near the boundary is important. This paper will give a thorough analysis of these methods applied to the shallow water equations, as well as some preliminary work applying these methods to volume penalization for bathymetry representation for use in either the nonhydrostatic or hydrostatic primitive equations.     

Lateral stirring of large-scale tracer fields by altimetry

Ocean surface fronts and filaments have a strong impact on the global ocean circulation and biogeochemistry. Surface Lagrangian advection with time-evolving altimetric geostrophic velocities can be used to simulate the submesoscale front and filament structures in large-scale tracer fields. We study this technique in the Southern Ocean region south of Tasmania, a domain marked by strong meso- to submesoscale features such as the fronts of the Antarctic Circumpolar Current (ACC). Starting with large-scale surface tracer fields that we stir with altimetric velocities, we determine ‘advected’ fields which compare well with high-resolution in situ or satellite tracer data. We find that fine scales are best represented in a statistical sense after an optimal advection time of ～2 weeks, with enhanced signatures of the ACC fronts and better spectral energy. The technique works best in moderate to high EKE regions where lateral advection dominates. This technique may be used to infer the distribution of unresolved small scales in any physical or biogeochemical surface tracer that is dominated by lateral advection. Submesoscale dynamics also impact the subsurface of the ocean, and the Lagrangian advection at depth shows promising results. Finally, we show that climatological tracer fields computed from the advected large-scale fields display improved fine-scale mean features, such as the ACC fronts, which can be useful in the context of ocean modelling.     

Why ripples become sweiis,and swells in shallow water decay rapidly

Abstract

In an ocean with a horizontal bottom where no wind is blowing it is shown that the spin (angular momentum) of the ocean is conserved. Thus, when energy is dissipated, at least one of three things will happen: i) Wave spectra may move towards lower frequencies. ii) The directional distribution may be changed towards long-crested waves. iii) Shear currents may be generated. By neglecting ii) and iii), the frequency shift of a spectrum is calculated due to molecular dissipation. When all energy transforming phenomena as e.g. wave breaking and turbulence generation are taken into account, the conservation of spin seems to be able to explain the frequency shift of wave spectra. In shallow water it is shown that there is energy transfer from the waves to shear currents.     

海洋表面流的高频雷达遥感

解释了利用高频雷达探测海洋表面流的基本原理，简要介绍了新研制的海态监测分析雷达OSMAR的组成及其工作流程，着重分析了利用OSMAR探测北部湾海流的实验结果．     

Turbulence avoidance and the wind-driven transport of plankton in the surface Ekman layer

Observations of turbulence avoidance in zooplankton are compared to estimates of the wind-driven turbulence in the upper ocean. Plankton that avoid wind-driven turbulence by moving deeper are no longer transported by the wind-driven Ekman currents near the surface because they are no longer near the surface. Here, a threshold level of turbulence that triggers an avoidance response is estimated, and is used to infer the wind speed and water column stratification conditions that would lead to zooplankton leaving the Ekman layer. Turbulence avoidance is argued to lead to near-shore retention in wind-driven upwelling systems, and to a reduction of the delivery of zooplankton to Georges Bank from the deeper waters of the Gulf of Maine.     

Time-optimal path planning in dynamic flows using level set equations: theory and schemes

We develop an accurate partial differential equation-based methodology that predicts the time-optimal paths of autonomous vehicles navigating in any continuous, strong, and dynamic ocean currents, obviating the need for heuristics. The goal is to predict a sequence of steering directions so that vehicles can best utilize or avoid currents to minimize their travel time. Inspired by the level set method, we derive and demonstrate that a modified level set equation governs the time-optimal path in any continuous flow. We show that our algorithm is computationally efficient and apply it to a number of experiments. First, we validate our approach through a simple benchmark application in a Rankine vortex flow for which an analytical solution is available. Next, we apply our methodology to more complex, simulated flow fields such as unsteady double-gyre flows driven by wind stress and flows behind a circular island. These examples show that time-optimal paths for multiple vehicles can be planned even in the presence of complex flows in domains with obstacles. Finally, we present and support through illustrations several remarks that describe specific features of our methodology.     

An assessment of the impact of surface currents on wave modeling in the Southern Ocean

This paper presents an assessment of the impact of the ocean circulation on modeled wave fields in the Southern Ocean, where a systematic positive bias of the modeled wave height against altimetry data has been reported. The inclusion of ocean currents in the wave model considerably reduces the positive bias of the simulated wave height for high southern latitudes. The decrease of wave energy in the presence of currents is almost exclusively related to the reduction of the relative wind, caused by an overall co-flowing current field associated with the Antarctic Circumpolar Current. Improvements of the model results are also found for the peak period and the mean period against a long-term moored buoy. At the mooring location, the effect of currents is greater for larger and longer waves, suggesting remotely generated swells are more influenced by the currents than local waves. However, an additional qualitative analysis using high-resolution currents in a finer grid nested to the global coarser grid shows that typical resolution of global hydrodynamic reanalysis is not sufficient to resolve mesoscale eddies, and as a consequence, the simulation of mesoscale wave patterns can be compromised. The results are also discussed in terms of the accuracy of forcing fields.     

Assessing the fidelity of surface currents from a coastal ocean model and HF radar using drifting buoys in the Middle Atlantic Bight

The rapid expansion of urbanization along the world’s coastal areas requires a more comprehensive and accurate understanding of the coastal ocean. Over the past several decades, numerical ocean circulation models have tried to provide such insight, based on our developing understanding of physical ocean processes. The systematic establishment of coastal ocean observation systems adopting cutting-edge technology, such as high frequency (HF) radar, satellite sensing, and gliders, has put such ocean model predictions to the test, by providing comprehensive observational datasets for the validation of numerical model forecasts. The New York Harbor Observing and Prediction System (NYHOPS) is a comprehensive system for understanding coastal ocean processes on the continental shelf waters of New York and New Jersey. To increase confidence in the system’s ocean circulation predictions in that area, a detailed validation exercise was carried out using HF radar and Lagrangian drifter-derived surface currents from three drifters obtained between March and October 2010. During that period, the root mean square (RMS) differences of both the east–west and north–south currents between NYHOPS and HF radar were approximately 15 cm s^?1. Harmonic analysis of NYHOPS and HF radar surface currents shows similar tidal ellipse parameters for the dominant M₂ tide, with a mean difference of 2.4 cm s^?1 in the semi-major axis and 1.4 cm s^?1 in the semi-minor axis and 3° in orientation and 10° in phase. Surface currents derived independently from drifters along their trajectories showed that NYHOPS and HF radar yielded similarly accurate results. RMS errors when compared to currents derived along the trajectory of the three drifters were approximately 10 cm s^?1. Overall, the analysis suggests that NYHOPS and HF radar had similar skill in estimating the currents over the continental shelf waters of the Middle Atlantic Bight during this time period. An ensemble-based set of particle tracking simulations using one drifter which was tracked for 11 days showed that the ensemble mean separation generally increases with time in a linear fashion. The separation distance is not dominated by high frequency or short spatial scale wavelengths suggesting that both the NYHOPS and HF radar currents are representing tidal and inertial time scales correctly and resolving some of the smaller scale eddies. The growing ensemble mean separation distance is dominated by errors in the mean flow causing the drifters to slowly diverge from their observed positions. The separation distance for both HF radar and NYHOPS stays below 30 km after 5 days, and the two technologies have similar tracking skill at the 95 % level. For comparison, the ensemble mean distance of a drifter from its initial release location (persistence assumption) is estimated to be greater than 70 km in 5 days.     

Synergy of satellite remote sensing and numerical modeling for monitoring of suspended particulate matter

Monitoring and modeling of the distribution of suspended particulate matter (SPM) is an important task, especially in coastal environments. Several SPM models have been developed for the North Sea. However, due to waves in shallow water and strong tidal currents in the southern part of the North Sea, this is still a challenging task. In general there is a lack of measurements to determine initial distributions of SPM in the bottom sediment and essential model parameters, e.g., appropriate exchange coefficients. In many satellite-borne ocean color images of the North Sea a plume is visible, which is caused by the scattering of light at SPM in the upper ocean layer. The intensity and length of the plume depends on the wave and current climate. It is well known that the SPM plume is especially obvious shortly after strong storm events. In this paper a quasi-3-D and a 3-D SPM transport model are presented. Utilizing the synergy of satellite-borne ocean color data with numerical models, the vertical exchange coefficients due to currents and waves are derived. This results in models that for the first time are able to reproduce the temporal and spatial evolution of the plume intensity. The SPM models consist of several modules to compute ocean dynamics, the vertical and horizontal exchange of SPM in the water column, and exchange processes with the seabed such as erosion, sedimentation, and resuspension. In the bottom layer, bioturbation via benthos and diffusion processes is taken into account.Responsible Editor: Jörg-Olaf Wolff     

Coupled assimilation for an intermediated coupled ENSO prediction model

The value of coupled assimilation is discussed using an intermediate coupled model in which the wind stress is the only atmospheric state which is slavery to model sea surface temperature (SST). In the coupled assimilation analysis, based on the coupled wind–ocean state covariance calculated from the coupled state ensemble, the ocean state is adjusted by assimilating wind data using the ensemble Kalman filter. As revealed by a series of assimilation experiments using simulated observations, the coupled assimilation of wind observations yields better results than the assimilation of SST observations. Specifically, the coupled assimilation of wind observations can help to improve the accuracy of the surface and subsurface currents because the correlation between the wind and ocean currents is stronger than that between SST and ocean currents in the equatorial Pacific. Thus, the coupled assimilation of wind data can decrease the initial condition errors in the surface/subsurface currents that can significantly contribute to SST forecast errors. The value of the coupled assimilation of wind observations is further demonstrated by comparing the prediction skills of three 12-year (1997–2008) hindcast experiments initialized by the ocean-only assimilation scheme that assimilates SST observations, the coupled assimilation scheme that assimilates wind observations, and a nudging scheme that nudges the observed wind stress data, respectively. The prediction skills of two assimilation schemes are significantly better than those of the nudging scheme. The prediction skills of assimilating wind observations are better than assimilating SST observations. Assimilating wind observations for the 2007/2008 La Niña event triggers better predictions, while assimilating SST observations fails to provide an early warning for that event.     

South Atlantic circulation in a world ocean model

The circulation in the South Atlantic Ocean has been simulated within a global ocean general circulation model. Preliminary analysis of the modelled ocean circulation in the region indicates a rather close agreement of the simulated upper ocean flows with conventional notions of the large-scale geostrophic currents in the region. The modelled South Atlantic Ocean witnesses the return flow and export of North Atlantic Deep Water (NADW) at its northern boundary, the inflow of a rather barotropic Antarctic Circumpolar Current (ACC) through the Drake Passage, and the inflow of warm saline Agulhas water around the Cape of Good Hope. The Agulhas leakage amounts to 8.7 Sv, within recent estimates of the mass transport shed westward at the Agulhas retroflection. Topographic steering of the ACC dominates the structure of flow in the circumpolar ocean. The Benguela Current is seen to be fed by a mixture of saline Indian Ocean water (originating from the Agulhas Current) and fresher Subantarctic surface water (originating in the ACC). The Benguela Current is seen to modify its flow and fate with depth; near the surface it flows north-westwards bifurcating most of its transport northward into the North Atlantic Ocean (for ultimate replacement of North Atlantic surface waters lost to the NADW conveyor). Deeper in the water column, more of the Benguela Current is destined to return with the Brazil Current, though northward flows are still generated where the Benguela Current extension encounters the coast of South America. At intermediate levels, these northward currents trace the flow of Antarctic Intermediate Water (AAIW) equatorward, though even more AAIW is seen to recirculate poleward in the subtropical gyre. In spite of the model’s rather coarse resolution, some subtle features of the Brazil-Malvinas Confluence are simulated rather well, including the latitude at which the two currents meet. Conceptual diagrams of the recirculation and interocean exchange of thermocline, intermediate and deep waters are constructed from an analysis of flows bound between isothermal and isobaric surfaces. This analysis shows how the return path of NADW is partitioned between a cold water route through the Drake Passage (6.5 Sv), a warm water route involving the Agulhas Current sheeding thermocline water westward (2.5 Sv), and a recirculation of intermediate water originating in the Indian Ocean (1.6 Sv).     

3-D Global Induction in the Oceans and Solid Earth: Recent Progress in Modeling Magnetic and Electric Fields from Sources of Magnetospheric,Ionospheric and Oceanic Origin

Electromagnetic induction in the Earth’s interior is an important contributor to the near-Earth magnetic and electric fields. The oceans play a special role in this induction due to their relatively high conductivity which leads to large lateral variability in surface conductance. Electric currents that generate secondary fields are induced in the oceans by two different processes: (a) by time varying external magnetic fields, and (b) by the motion of the conducting ocean water through the Earth’s main magnetic field. Significant progress in accurate and detailed predictions of the electric and magnetic fields induced by these sources has been achieved during the last few years, via realistic three-dimensional (3-D) conductivity models of the oceans, crust and mantle along with realistic source models. In this review a summary is given of the results of recent 3-D modeling studies in which estimates are obtained for the magnetic and electric signals at both the ground and satellite altitudes induced by a variety of natural current sources. 3-D induction effects due to magnetospheric currents (magnetic storms), ionospheric currents (Sq, polar and equatorial electrojets), ocean tides, global ocean circulation and tsunami are considered. These modeling studies demonstrate that the 3-D induction (ocean) effect and motionally-induced signals from the oceans contribute significantly (in the range from a few to tens nanotesla) to the near-Earth magnetic field. A 3-D numerical solution based on an integral equation approach is shown to predict these induction effects with the accuracy and spatial detail required to explain observations both on the ground and at satellite altitudes. On leave from Institute of Terrestrial Magnetism, Ionosphere and Radiowave Propagation, Russian Academy of Sciences, 142190 Troitsk, Moscow region, Russia.     

The energy and energy flux of planetary waves in an ocean of variable depth

Abstract

In a recent paper, Buchwald (1972a) has shown that besides the kinetic energy and gravitational potential energy usually associated with planetary waves in an ocean of uniform depth it is useful to define also a “spin energy”, associated with the rotation.

The present paper is basically an extension of Buchwald's result to a uniformly rotating β-plane ocean of variable depth. As in the previous work, energy conservation equations are derived and the separate energies shown to be independently conserved over the total volume of the ocean. The time-averaged energies are further shown to be propagated in the direction of the group velocity and to satisfy the equipartition rule.

Unlike Buchwald, however, we need not consider the boundary conditions in order to achieve these results. Furthermore, the use of a more realistic ocean configuration admits the possibility of a multiply connected region in the present of mean currents.

Finally, there is a physical explanation for the appearance of a spin energy in a rotating system.     

Baroclinic and topographic influences on the transport in western boundary currents

Abstract

One of the central unsolved theoretical problems of the large scale ocean circulation is concerned with explaining the very large transports measured in western boundary currents such as the Gulf Stream and the Kuroshio. The only theory up to now that can explain the size of these transports is that of non-linear recirculation in which the advective terms in the momentum equations became important near the western boundary. In this paper an alternative explanation is suggested. When bottom topography and baroclinic effects are included in a wind-driven ocean model it is shown that the western boundary current can have a transport larger than that predicted from the wind stress distribution even when the nonlinear advective terms are ignored. The explanation lies in the presence of pressure torques associated with bottom topography which can contribute to the vorticity balance in the same sense as the wind stress curl.

Three numerical experiments have been carried out to explore the nature of this process using a three dimensional numerical model. The first calculation is done for a baroclinic ocean of constant depth, the second for a homogeneous ocean with an idealized continental slope topography, and the third for a baroclinic ocean with the same continental slope topography. The nature of the vorticity balance and of the circulation around closed paths is examined in each case, and it is shown that bottom pressure torques lead to enhanced transport in the western boundary current only for the baroclinic case with variable depth.