Stokes drift estimation for deep water waves based on short-term variation of wave conditions

The paper provides a simple analytical method which can be used to give estimates of the Stokes drift based on short-term variation of wave conditions. This is achieved by providing bivariate distributions of wave height and surface Stokes drift as well as wave height and volume Stokes transport for individual random waves within a sea state. The paper presents and discusses statistical aspects of these Stokes drift parameters, as well as examples of results corresponding to typical field conditions.     

Wave power statistics for individual waves

The paper provides a bivariate distribution of wave power and wave height, as well as a bivariate distribution of wave power and wave period; both bivariate distributions are for individual waves within a sea state. This is relevant for e.g. making assessments of wave power devices and their potential for converting energy from waves. The results can be applied to compare systematically the wave power potential for individual waves in a given sea state at different locations.     

Statistical properties of successive wave heights and successive wave periods

Two successive wave heights are modeled by a Gaussian copula, which is referred to as the Nataf model. Results with two initial distributions for the transformation are presented, the Næss model [Næss A. On the distribution of crest to trough wave heights. Ocean Engineering (1985);12(3):221–34] and a two-parameter Weibull distribution, where the latter is in best agreement with data. The results are compared with existing models. The Nataf model has also been used for modeling three successive wave heights.Results show that the Nataf transformation of three successive wave heights can be approximated by a first order autoregressive model. This means that the distribution of the wave height given the previous wave height is independent of the wave heights prior to the previous wave height. Thus, the joint distribution of three successive wave heights can be obtained by combining conditional bivariate distributions. The simulation of successive wave heights can be done directly without simulating the time series of the complete surface elevation.Successive wave periods with corresponding wave heights exceeding a certain threshold have also been studied. Results show that the distribution for successive wave periods when the corresponding wave heights exceed the root-mean-square value of the wave heights, can be approximated by a multivariate Gaussian distribution.The theoretical distributions are compared with observed wave data obtained from field measurements in the central North Sea and in the Japan Sea, with laboratory data and numerical simulations.     

最大熵原理应用于海浪波高分布的研究

利用最大熵原理从理论上推导出波高的最大熵分布,在此基础上研究了状态参量对波高分布和波高熵的影响。影响最大熵分布的因子是平均波高和状态参量,不同海况对应的状态参量是不同的。利用波高实测资料,得出3种不同海况下的最大熵分布,通过比较发现最大熵分布很好地符合实测数据。把最大熵分布与目前广泛应用的瑞利分布作了比较,结果表明,最大熵分布有2个优点：没有对波高作出任何限制性假定和能够描述不同海况下的波高分布。     

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.     

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.     

波高与风速联合概率分布研究

基于copula函数论述了两变量的联合概率分布方法.此方法的主要优点是边缘分布可由不同的分布函数构成,变量间可具相关性.以粤东汕尾海域极值波高与相应风速为研究实例,经分析获得以下结果:(1)优选的极值波高和风速可分别由P-Ⅲ型和GEV 分布表示:(2)拟合优度检验指标表明二者的最优连接函数为Archimedean co...     

平直斜坡上破波带内质量输移流对污染物输移影响的数值研究

破波带内外都有质量输移流存在,其对破波带内污染物输移有怎样的影响,需要进一步深入研究。本文基于实验以及考虑质量输移流的对流扩散数学模型研究了平直斜坡上破波带内质量输移流对污染物输移影响。数学模型包括波浪模型、近岸流模型以及对流扩散模型。首先建立了破波带内污染物输移数学模型,其中波浪场基于波能守恒方程来计算,波导流场基于Longuet-Higgins提出的辐射应力模拟,污染物对流扩散方程中考虑了质量输移流的影响,并利用算例验证该数学模型。其次简要介绍了平直斜坡上破波带内污染物输移实验,并分析了污染物输移特性。连续投放污染物会形成污染带,本文分析了两种波况下不同时刻污染带与岸线夹角的变化,以及污染物在垂直岸线和沿岸线方向的输移速度,结果表明对两种波况来说在初始10-40s污染团向岸线方向输移速度分别约为0.05m/s、0.017m/s,之后速度分别减小为0.001m/s、0.011m/s。数值模拟结果与实验结果比较表明：考虑质量输移流的模拟结果与实验更为吻合。因而,通过实验以及数模研究表明破波带内质量输移流对破波带内污染物在垂直岸线方向的输移有重要影响,而对沿岸方向的输移则影响较小。     

基于Clayton Copula函数的二维Gumbel模型及其在海洋平台设计中的应用

根据Gumbel分布和Clayton Copula函数构造出二维Gumbel Clayton Copula分布,根据渤海海域某观测站测得的1970—1993年年最大波高及风速,具体介绍该二维Gumbel模型在海洋工程设计中的使用方法。通过导管架平台基底剪力计算表明,使用该二维Gumbel模型所得的50a一遇剪力值降低了37%,对于边际油田,可以降低荷载设计标准,从而减少海洋工程的投资费用。     

随机波浪下斯托克司漂移速度是否引起水面漂移物离散的研究

漂浮于自由水面的污染物的的迁移、扩散会受到天然随机海浪的影响。之前的研究（以Herterich和Hasselmann（1982）为代表）普遍认为,随机波浪作用下的斯托克司漂移速度会引起水面污染物的离散,这个离散甚至有可能跟风和海流引起的离散同一量级。本研究就随机波浪作用下的斯托克司漂移速度是否会引起水面漂移物的离散进行理论和试验探讨。从理论推导可知,随机波浪下的质量输移速度是个定常分量,因此它不会随时间变化而引起水面漂移物的离散。随后我们在实验室水槽中进行了漂移物在随机波浪（P-M谱）作用下的漂移过程的测量。试验结果也印证了随机波浪作用下的斯托克司漂移速度不会引起水面漂移物离散的结论。     

一种改进的基于逆Omega-K算法的海面场景SAR原始数据仿真方法

合成孔径雷达(SAR)海面场景原始数据仿真是研究海洋动力参数(表面波浪、风矢量和洋流)的有效工具。目前海面场景原始数据仿真方法已经基于逆Omega-K算法实现了海洋运动参数的空间变化。但是目前仅仅讨论了正侧视情况下的海面场景仿真,应用范围有限,同时没有考虑Stokes漂流以及Bragg相速度的影响,而这两者都是存在于真实海面的。通常情况下为了反演得到海面流场的二维速度矢量,雷达需要从两个不同的方位方向观察海面的同一个区域,因此这就需要考虑大斜视的雷达波束,同时Stokes漂流和Bragg相速度是SAR海表面流场观测不容忽视的两种运动。本文在不改变原有正侧视逆Omega-K算法的情况下,通过增加重新计算零方位时刻的斜视波束中心位置坐标,并据此确定SAR原始数据在多普勒域的位置来将其扩展到大斜侧视逆Omega-K算法,并通过时域Stokes漂流公式到频域内离散化Stokes漂流公式的推导来加入Stokes漂流,以及根据Bragg散射机制加入了Bragg相速度。仿真结果表明,经过聚焦成像后的SAR图像很好的体现了真实海面波浪场的形状,同时能够很好地反演出设定的雷达径向流场速度,且流速精度误差控制在6%以内。最后也证明了Bragg相速度以及Stokes漂流对于海面流场的影响不可忽视。     

Long-Term Joint Distribution of One-Tenth Large (or Significant) Wave Height with Average Period

- In this paper, by using the wave data from a few oceanographic observation stations in the coastal zone of the Yellow Sea, the East China Sea and the South China Sea, the long-term joint distribution of the one-tenth large (or significant) wave height with average period is studied. The statistical data demonstrate that the long- term distribution of the one- tenth wave height or average period fits the log-normal distribution, thus the joint distribution also fits the two-dimensional log-normal distribution. Then the conditional probability distribution of the average period is derived, and the range as well as the mode of the average wave period corresponding to a certain return period of wave height can be calculated easily.     

大亚湾风暴潮波浪传播数学模型

本文建立了大亚湾风暴潮涌浪传播数学模型,以任意多边形离散计算海区,每一多边形构成单元,以波浪运动学和动力学守恒方程模拟单元内能量传递,以风暴潮过程模拟边界入射波高过程,用风暴潮涌浪传播基本方程和波能缓坡方程结合模拟湾口巨浪向湾内的传播过程。通过分析大亚湾不同种类岸线反射系数的概率分布,并结合实测波高对模型进行率定,最终确定模型参数。将大亚湾特征点计算波高与统计推荐波高比较进行模型验证,结果显示SE向波高与H_(13%)推荐波高对应较好,可以用于大亚湾海区的波浪预报。计算当大亚湾口分别出现10年一遇、50年一遇及100年一遇的波高时,在E、ESE以及ES向入射波浪条件下大亚湾海域极值波高的分布。分别对风暴潮涌浪在不同类型岸线的爬高以及风暴潮涌浪传播对岸线的作用力进行计算。     

On the long-term joint distribution of characteristic wave height and period and its application

By analysing the scatter diagrams of characteristic the wave height H and the period T on the basis of instrumental data from various ocean wave stations, we found that the conditional expectation and standard deviation of wave period for a given wave height can be better predicted by using the equations of normal linear regression rather than by those based on the log- normal law. The latter was implied in Ochi' s bivariate log-normal model(Ochi. 1978) for the long-term joint distribution of H and T. With the expectation and standard deviation predicted by the normal linear regression equations and applying proper types of distribution, we have obtained the conditional distribution of T for given H. Then combining this conditional P(T / H) with long-term marginal distribution of the wave height P(H) we establish a new parameterized model for the long-term joint distribution P(H,T). As an example of the application of the new model we give a method for estimating wave period associated with an extreme w     

The effects of random surface waves on the steady Ekman current solutions

The response of near-surface current profiles to wind and random surface waves are studied based on the approach of Jenkins [1989. The use of a wave prediction model for driving a near surface current model. Dtsch. Hydrogr. Z. 42, 134–149] and Tang et al. [2007. Observation and modeling of surface currents on the Grand Banks: a study of the wave effects on surface currents. J. Geophys. Res. 112, C10025, doi:10.1029/2006JC004028]. Analytic steady solutions are presented for wave-modified Ekman equations resulting from Stokes drift, wind input and wave dissipation for a depth-independent constant eddy viscosity coefficient and one that varies linearly with depth. The parameters involved in the solutions can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water, and the solutions reduce to those of Lewis and Belcher [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans. 37, 313–351] when only the effects of Stokes drift are included. As illustrative examples, for a fully developed wind-generated sea with different wind speeds, wave-modified current profiles are calculated and compared with the classical Ekman theory and Lewis and Belcher's [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans 37, 313–351] modification by using the Donelan and Pierson [1987. Radar scattering and equilibrium ranges in wind-generated waves with application to scatterometry. J. Geophys. Res. 92, 4971–5029] wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. Illustrative examples for a fully developed sea and the comparisons between observations and the theoretical predictions demonstrate that the effects of the random surface waves on the classical Ekman current are important, as they change qualitatively the nature of the Ekman layer. But the effects of the wind input and wave dissipation on surface current are small, relative to the impact of the Stokes drift.     

Numerical simulation of wave damping over porous seabeds

This paper presents a study of wave damping over porous seabeds by using a two-dimensional numerical model. In this model, the flow outside of porous media is described by the Reynolds Averaged Navier–Stokes equations. The spatially averaged Navier–Stokes equations, in which the presence of porous media is considered by including additional inertia and nonlinear friction forces, is derived and implemented for the porous flow. Unlike the earlier models, the present model explicitly represents the flow resistance dependency on Reynolds number in order to cover wider ranges of porous flows. The numerical model is validated against available theories and experimental data. The comparison between the numerical results and the theoretical results indicates that the omission or linearization of the nonlinear resistance terms in porous flow models, which is the common practice in most of analytical models, can lead to significant errors in estimating wave damping rate. The present numerical model is used to simulate nonlinear wave interaction with porous seabeds and it is found that the numerical results compare well with the experimental data for different wave nonlinearity. The additional numerical tests are also conducted to study the effects of wavelength, seabed thickness and Reynolds number on wave damping.     

A joint distribution of significant wave height and characteristic surf parameter

The paper provides a joint distribution of significant wave height and characteristic surf parameter. The characteristic surf parameter is given by the ratio between the slope of a beach or a structure and the square root of the characteristic wave steepness in deep water defined in terms of the significant wave height and the spectral peak period. The characteristic surf parameter is used to characterize surf zone processes and is relevant for e.g. wave run-up on beaches and coastal structures. The paper presents statistical properties of the wave parameters as well as an example of results corresponding to typical field conditions.     

Wave-induced drift of an elliptical surface film

This study investigates experimentally the drift velocity of an elliptical surface film advected by deep-water waves. Thin polyethylene sheets were used to simulate the inextensible surface film. The drift velocities were obtained by recording and analyzing a sequence of images captured using a video camera. The results show that the drift velocity increases with the longitudinal length of the polyethylene sheet until approximately 0.8 times the wavelength. Beyond that, further increment would not result in substantial increase in the drift velocity. The effect due to the normalized transverse width is found to be significant within a particular range. At large wave steepness, the drift velocity appears to be limited by the Stokes drift. A set of best-fit empirical equations based on the sigmoidal function is introduced for oil spill trajectory prediction.