Effect of Stokes drift on upper ocean mixing

Stokes drift is the main source of vertical vorticity in the ocean mixed layer. In the ways of Coriolis - Stokes forcing and Langmuir circulations, Stokes drift can substantially affect the whole mixed layer. A modified Mellor-Yamada 2. 5 level turbulence closure model is used to parameterize its effect on upper ocean mixing conventionally. Results show that comparing surface heating with wave breaking, Stokes drift plays the most important role in the entire ocean mixed layer, especially in the subsurface layer. As expected, Stokes drift elevates both the dissipation rate and the turbulence energy in the upper ocean mixing. Also, ilffluence of the surface heating, wave breaking and wind speed on Stokes drift is investigated respectively. Research shows that it is significant and important to assessing the Stokes drift into ocean mixed layer studying. The laboratory observations are supporting numerical experiments quantitatively.     

全球Stokes 漂流的时空分布特征研究

利用海浪模式WWIII(Wave Watch III)2008年的模拟结果对海面Stokes漂流、Stokes输运、Stokes深度以及全球Langmuir数的年平均分布特征和季节平均分布特征分别进行了详细的研究与分析。结果表明,海面Stokes漂流和Stokes输运均呈现高纬度偏大的特征,以南极绕极流海域最为突出。全球大部分海域Stokes漂流影响深度在20 m以内,呈现大洋东部偏大,西部偏小的分布特征。全球大部分海域的混合作用是剪切不稳定性和Langmuir湍效应并存的状态,甚至有些海域是以Langmuir湍效应为主。因此,在进行大尺度的海洋数值模拟时,应该考虑波浪导致的混合效应。     

失稳背景流下对称型海洋内波的生成演变

本文采用无海底地形但考虑海洋跃层和剪切背景流的二维非静力准不可压缩方程组的数值模式，开展失稳垂向剪切背景流下线性和非线性对称型海洋内波生成演变的数值实验，并对结果进行分析、比较和讨论。研究结果表明，线性内波强度随积分时间始终呈指数增长，并有内波的对称不稳定；而非线性内波强度则在发展期呈准线性增长，最终进入稳定期。线性增长比非线性增长要快得多，非线性效应具有维稳作用。对该线性和非线性对称型内波，在跃层附近位密度扰动均有大值中心，即其为跃层所俘获，这与实际观测相一致；流函数与位密度扰动两者均有很好配合，位密度扰动的正、负中心分别相应于流函数的上升、下沉运动，表明有从海底向上的斜对流发生，且以跃层为顶盖。对线性内波来说，随积分时间增加，其波形大体不变，其正、负振幅也大体相同，并有符号相反原地增长的两个倾斜环流圈，而在它们之间则有较强倾斜上升流。非线性内波波形随积分时间改变，倾斜环流圈数目也在增加，最终形成负环流强于正环流的结果，并导致流函数、位密度扰动水平梯度剧增，其可视为间断。     

Nonlinear effects in the process of propagation of internal waves in the presence of turbulence

In the Boussinesq approximation, we study the nonlinear effects observed in the process of propagation of internal waves in the presence of turbulence. The space damping factor of the waves is evaluated. The Stokes drift velocity and the Euler velocity of the mean current induced by waves due to the presence of nonlinearity are determined. It is shown that the principal contribution to the wave transfer is made by the horizontal velocity of the induced current. The Stokes drift is significant only near the bottom. The vertical component of the Stokes drift velocity obtained with regard for the turbulent viscosity is nonzero.     

Langmuir环流在雷诺平均型模式中的参数化研究进展

Langmuir环流影响着海洋上层的能量输入,对海洋上混合层的形成和加深起着重要作用,对于海洋上混合层具有重要意义。近年来许多学者采用大涡模拟(LES)方法对Langmuir环流进行机制研究,并通过在雷诺平均模型中参数化Langmuir环流效应,将Langmuir环流过程引入到三维海洋环流或海洋耦合模式中,提出了一系列混合参数化方案。本文回顾了Langmuir环流在雷诺平均模式参数化中的研究进展,主要可分为以下几种方案:一种方法是用Langmuir数在KPP垂直混合参数化方案中引入湍流特征速度增强因子,并不断发展Langmuir数的定义;一种是在Mellor-Yamada2.5湍流闭合模型中增加斯托克斯漂流剪切效应项,此外还有通过修改模式中混合长方程来加入Langmuir效应等。通过在雷诺平均模式中应用的结果来看,现有的参数化方案在一定程度上改善了混合层深度和SST的模拟,肯定了Langmuir环流在加深混合等方面的作用,但仍存在一些问题需要在今后的研究中进一步改进。     

An experimental study of the spatial structures of the instabilities of longshore currents on plane beaches

A laboratory experiment on alongshore currents is conducted for two plane beaches with slopes 1:40 and 1:100 to investigate the instability of alongshore currents.The dye release experiment is also performed synchronously in surf zone.Complicated and strongly unstable motions of alongshore currents are observed in the experiment.To examine the spatial and temporal variations of the shear instabilities of longshore clearly,dye batches are released in surf zone.The deformation of the dye patch is observed efficiently and effectively with charge coupled device(CCD) system.Some essential characteristics of the shear instability are validated from the results of image analyses of the temporal variation of the dye patch.The influences of alongshore currents,Stokes drift,large-scale vorticity and the shear instabilities on the transport of dye are analyzed using the collected images.The spatial structure of the instabilities of longshore currents is studied by analyzing collected images of the dye patch.And the phase velocity of the meandering movements is obtained through measuring the movement distances of the oscillations of dye patch in alongshore direction with time.The results suggest that the propagation speed of the shear instability is approximately 50%-75% of maximum of mean alongshore currents for irregular and regular waves.The calculated propagation speed using a linear instability analysis theory is compared with the experimental results.The comparison shows agreements between them.     

Flow under standing waves: Part 1. Shear stress distribution,energy flux and steady streaming

The conditions for energy flux, momentum flux and the resulting streaming velocity are analysed for standing waves formed in front of a fully reflecting wall. The exchange of energy between the outer wave motion and the near bed oscillatory boundary layer is considered, determining the horizontal energy flux inside and outside the boundary layer. The momentum balance, the mean shear stress and the resulting time averaged streaming velocities are determined. For a laminar bed boundary layer the analysis of the wave drift gives results similar to the original work of Longuet–Higgins from 1953. The work is extended to turbulent bed boundary layers by application of a numerical model. The similarities and differences between laminar and turbulent flow conditions are discussed, and quantitative results for the magnitude of the mean shear stress and drift velocity are presented. Full two-dimensional simulations of standing waves have also been made by application of a general purpose Navier–Stokes solver. The results agree well with those obtained by the boundary layer analysis. Wave reflection from a plane sloping wall is also investigated by using the same numerical model and by physical laboratory experiments. The phase shift of the reflected wave train is compared with theoretical and empirical models.     

高频表面波对定常Ekman流解的影响

基于Jenkins(1989)建立的包含Stokes漂流、风输入和波耗散影响的修正Ekman模型,采用Paskyabi等(2012)使用的推广的Donelan等(1987)中的谱和波耗散函数,并利用Paskyabi等(2012)中修正方法给出的包含高频波的风输入函数,在粘性不依赖于水深及粘性随深度线性变化的条件下,研究了包含高频毛细重力波的随机表面波对Stokes漂流和Song(2009)导出的波浪修正定常Ekman流解的影响。结果表明高频表面波使Stokes漂流在海表面剪切加强,对定常Ekamn流解的影响通常不能忽略,但对Ekman流场的角度偏转影响很小。最后,将考虑高频表面波尾谱影响所估算的定常Ekman流解与已有观测结果以及经典Ekman解进行了比对分析。     

台风条件下朗缪尔环流对上层海洋混合的影响研究进展

回顾了近10年来台风条件下朗缪尔环流影响上层海洋混合的研究进展，朗缪尔致湍流对海洋上混合层的形成和加深的重要作用已形成了基本共识，但对于朗缪尔致湍流对海洋上混合层的混合作用机制和程度仍然存在诸多不确定性。观测表明台风条件下台风眼附近的混合层平均湍流动能受到了较强的抑制，可能与台风不同位置朗缪尔致湍流的特征变异有关；台风条件下，现有的朗缪尔致湍流参数化方案在上层混合过程模拟中还有显著误差。在今后研究中，通过改进斯托克斯漂流剖面的计算方法，优化表征台风条件下海面状况的朗缪尔致湍流参数化计算方案，是进一步揭示台风条件下朗缪尔环流对海洋上层混合的影响机理的必要途径。     

Bottom friction caused by boundary layer streaming beneath random waves for laminar and smooth turbulent flow

The effect of boundary layer streaming on the sea bed shear stresses, beneath random waves, is investigated for laminar flow as well as smooth turbulent flow. It is demonstrated how bottom friction formulas for regular waves can be used to obtain the bed shear stresses resulting from steady streaming under random waves. As a result, friction factors for steady streaming under random waves are provided, and the effect of streaming versus the effect of linear waves is discussed. For laminar flow the effect of second order Stokes waves is also included. Examples are included to illustrate the applicability of the present practical method, and results are obtained using data typical for field conditions.     

大洋细结构的一种可能的机制Ⅲ.对混合过程动能耗散的估计

继第部分之后研究了惯性内波和近惯性内波由f～的作用所致的剪切不稳定引起的破碎机制。物理上,该机制很象存在由风应力所致薄表面涡旋漂流层时表面波的破碎与饱和过程。惯性内波和近惯性内波的破碎产物与小尺度湍流一起形成了混合块,它与Gregg等人(1986)的持久混合观测结果一致。依据Thorpe(1973)实验的结果作者提出了一个估计湍流动能耗散率和消衰时间的方法。结果表明,在剪切不稳定中近惯性内波在湍动耗散中起了关键作用,而惯性内波引起非常弱的湍动耗散。使用内波能量谱的标准总能量密度估计出的近惯性内波的耗散率和消衰时间与PATCHEX测量结果非常一致。文中还讨论了几个与此破碎机制有关的问题。     

Longitudinal dispersion in wave-current-vegetation flow

The flow, turbulence, and longitudinal dispersion in wave-current flow through submerged vegetation are experimentally examined. Laboratory experiments are carried out by superimposing progressive waves on a steady flow through simulated submerged vegetation. The resultant wave-current-vegetation interaction shows strong interface shear with increase in the velocity due to the wave-induced drift. The increase in turbulence in the region of vegetation is found to be about twice higher than in the no-wave case due to the additional mixing by wave motions. Solute experiments are conducted to quantify the wave-current-vegetation longitudinal dispersion coefficient (WCVLDC) by the routing method and by defining length and velocity scales for the wave-current-vegetation flow. An empirical expression for the WCVLDC is proposed. Although the increase in vertical diffusivity is observed as compared with bare-bed channels, the shear effect is stronger, which increases the value of the WCVLDC. The study can be a guideline to understand the combined hydrodynamics of waves, current, and vegetation and quantify the longitudinal dispersion therein. Published in Morskoi Gidrofizicheskii Zhurnal, No. 1, pp. 50–67, January–February, 2009.     

表面波浪对海洋上层混合的作用

A new three-dimensional numerical model is derived through a wave average on the primitive N-S equations, in which both the＂Coriolis-Stokes forcing＂ and the＂Stokes-Vortex force＂ are considered. Three ideal experiments are run using the new model applied to the Princeton ocean model （POM）. Numerical results show that surface waves play an important role on the mixing of the upper ocean. The mixed layer is enhanced when wave effect is considered in conjunction with small Langmuir numbers. Both surface wave breaking and Stokes production can strengthen the turbulent mixing near the surface. However, the influence of wave breaking is limited to a thin layer, but Stokes drift can affect the whole mixed layer. Furthermore, the vertical mixing coefficients clearly rise in the mixed layer, and the upper ocean mixed layer is deepened especially in the Antarctic Circumpolar Current when the model is applied to global simulations. It indicates that the surface gravity waves are indispensable in enhancing the mixing in the upper ocean, and should be accounted for in ocean general circulation models.     

The energy cycle of wind waves on the sea surface

The problems of wind-induced waves on the sea surface are considered. To this end, the empirical fetch laws that determine variations in the basic periods and heights of waves in relation to their fetch are used. The relation between the fetch and the physical time is found, as are the dependences of the basic characteristics of waves on the time of wind forcing. It is found that about 5% of wind energy dissipated in the near-water air layer contributes to the growth of wave heights, i.e. wave energy, although this quantity depends on the age of waves and the exponent in the fetch laws. With consideration for estimates of the probability distribution functions for the wind over the world ocean [11], it is found that the rate of wind-energy dissipation in the near-water air layer is on the order of 1 W/m². The calculations of wind waves [19] for the world ocean for 2007 have made it possible to assess the mean characteristics of the cycle of wave development and their seasonal variations. An analysis of these calculations [19] shows that about 20% of wind energy is transferred to the water surface. The remaining amount (80%) of wind energy is spent on the generation of turbulence in the near-water air layer. About 2%, i.e., one tenth of the energy transferred to water, is spent on turbulence generation due to the instability of the vertical velocity profile of the Stokes drift current and on energy dissipation in the surf zones. Of the remaining 18%, 5% is spent directly on wave growth and 13% is spent on the generation of turbulence during wave breaking and on a small-scale spectral region. These annually and globally mean estimates have a seasonal cycle with an amplitude on the order of 20% in absolute values but with a smaller amplitude in relative values. According to [19] and to the results of this study, the annually mean height of waves is estimated as 2.7 m and their age is estimated as 1.17.     

Self-localization of planetary wave structures in the ionosphere upon interaction with nonuniform geomagnetic field and zonal wind

The generation and further linear and nonlinear dynamics of planetary magnetized Rossby waves (MRWs) in the rotating dissipative ionosphere are studied in the presence of a zonal wind (shear flow). MRWs are caused by interaction with the spatially nonuniform geomagnetic field and are ionospheric manifestations of ordinary tropospheric Rossby waves. A simplified self-consistent set of model equations describing MRW-shear flow interaction is derived on the basis of complete equations of ionospheric magnetohydrodynamics. Based on an analysis of an exact analytical solution to the derived dynamic equations, an effective linear mechanism of MRW amplification in the interaction with nonuniform zonal wind is ascertained. It is shown that operators of linear problems are non-self-adjoint in the case of shear flows, and the corresponding eigenfunctions are nonorthogonal; therefore, the canonically modal approach is of little use when studying such flows; a so-called nonmodal mathematical analysis is required. It is ascertained that MRWs effectively get shear flow energy during the linear stage of evolution and significantly increase (by several orders of magnitude) their energy and amplitude. The necessary and sufficient condition of shear flow instability in an ionospheric medium is derived. Nonlinear self-localization begins with the development of shear instability and an increase in the amplitude, and the process ends with the self-organization of strongly localized isolated large-scale nonlinear vortex structures. Thus, a new degree of freedom and a way for perturbation evolution to occur appear in medium with shear flow. The nonlinear systems can be a pure monopole vortex, a vortex streets, or vortex chains depending of the shape of the sheared flow velocity profile. The accumulation of such vortices in the ionospheric medium can produce a strongly turbulent state.     

Analysis of some key parametrizations in a beach profile morphodynamical model

A numerical process-based model to forecast beach profile morphodynamics has been developed. In the present paper, an analysis of various modelling approaches and key parametrizations involved in the estimation of the wave-driven current and the suspended sediment concentration is carried out.Several resolution techniques for the 1DV horizontal (i.e., in the x-direction perpendicular to coastline) momentum equation governing the Mean Horizontal Velocity (MHV) are analysed. In the first kind of techniques, the mean horizontal velocity is computed from the momentum equation, whereas the Mean Water Level (MWL) is computed using a parametrization of the depth-averaged momentum equation. Two boundary or integral conditions are thus needed. In the second kind, both mean horizontal velocity and mean water level gradient in the x-direction are the unknowns of the momentum equation, thus, three boundary or integral conditions are needed. Various additional conditions are discussed. We show that using a technique of the first kind is equivalent to imposing the difference between the surface and the bottom shear stresses in the 1D vertical equation. Both techniques lead to results that are in good agreement with the Delta Flume experimental data, provided the Stokes drift flow discharge is imposed as an additional condition. The influence of the breaking roller model and of the turbulent viscosity parametrization are also analysed.Suspended sediment transport by the mean current and wave-induced bedload transport are taken into account in the sediment flux. Three turbulent diffusivity parametrizations are compared for suspended sediment concentration estimations. A linear profile for the turbulent diffusivity taking into account the wave bottom shear stress and the surface wave breaking turbulence production is shown to give the best results. Using experimental data, we put forward the poor estimation of the bottom sediment concentration given by the three implemented parametrizations. We thus propose a new parametrization relying on a Shields parameter based on the breaking roller induced surface shear stress. Using this new parametrization, the bottom profile used in the tests keeps its two bars which disappear otherwise. However, the morphodynamical model still overestimates the bars offshore motion, a bias already observed in other models.     

Vertical momentum transfer by internal waves when eddy viscosity and diffusion are taken into account

Free internal waves are considered in a Boussinesq approximation in the situation when horizontal eddy viscosity and diffusion in a vertically inhomogeneous flow are taken into account. The dispersion relation and wave damping factor are found in a linear approximation. The Stokes drift velocity is determined in the second order of smallness based on the wave amplitude. It has been indicated that the Stokes drift velocity, transverse with respect to the wave propagation direction, differs from zero if the flow-rate transverse component depends on the vertical coordinate. Vertical momentum fluxes differ from zero and can be comparable with or exceed the corresponding turbulent fluxes if eddy viscosity and diffusion are taken into account.     

The influence of Stokes drift on oil spills: Sanchi oil spill case

Spilled oil floats and travels across the water’s surface under the influence of wind, currents, and wave action. Wave-induced Stokes drift is an important physical process that can affect surface water particles but that is currently absent from oil spill analyses. In this study, two methods are applied to determine the velocity of Stokes drift, the first calculates velocity from the wind-related formula based upon a one-dimensional frequency spectrum, while the second determines velocity directly from the wave model that was based on a two-dimensional spectrum. The experimental results of numerous models indicated that: (1) oil simulations that include the influence of Stokes drift are more accurate than that those do not; (2) for medium and long-term simulations longer than two days or more, Stokes drift is a significant factor that should not be ignored, and its magnitude can reach about 2% of the wind speed; (3) the velocity of Stokes drift is related to the wind but is not linear. Therefore, Stokes drift cannot simply be replaced or substituted by simply increasing the wind drift factor, which can cause errors in oil spill projections; (4) the Stokes drift velocity obtained from the two-dimensional wave spectrum makes the oil spill simulation more accurate.     

Comparison of turbulence schemes for prediction of wave-induced near-bed sediment suspension above a plane bed

Based on a wave bottom boundary layer model and a sediment advection-diffusion model, seven turbulence schemes are compared regarding their performances in prediction of near-bed sediment suspension beneath waves above a plane bed. These turbulence algorithms include six empirical eddy viscosity schemes and one standard two-equation k-ε model. In particular, different combinations of typical empirical formulas for the eddy viscosity profile and for the wave friction factor are examined. Numerical results are compared with four laboratory data sets, consisting of one wave boundary layer hydrodynamics experiment and three sediment suspension experiments under linear waves and the Stokes second-order waves. It is shown that predictions of near-bed sediment suspension are very sensitive to the choices of the empirical formulas in turbulence schemes. Simple empirical turbulence schemes are possible to perform equally well as the two-equation k-ε model. Among the empirical schemes, the turbulence scheme, combining the exponential formula for eddy viscosity and Swart formula for wave friction factor, is the most accurate. It maintains the simplicity and yields identically good predictions as the k-ε model does in terms of the wave-averaged sediment concentration.